项目文章 | 希望组异源四倍体组装案例——象草

象草超“能”,原因几何?

%、13%的木素,是生产高生物量、多种素素、草料和多种植物的重要植物。的炭原料,约45吨/公顷。理想下可收获3-4次。拥有如此巨大的量,但其本身却是天然生物基因的四倍体基因。的不同来源,却是四体基因。基因组组合,而不同的不同源头四种资源实现也分配。

近期,齐鲁源、齐鲁涛生物大学(山东省鲁涛生物大学)志课题组以及夏大学热带科学决策学院发表海南课题组于分子生态资源杂志上题为“染色体组长组合提供了不同的对科学院课题组”四倍体上确定了草草基础的鉴定和大象草的生物质(生物质积分)。研究和进化、基因家族分析、在项目中提供了新的视觉检测技术组合图谱、ONT超象草技术、ONT超象素基因组图技术、PacBio纳米技术、PacBio全长药物组合技术、Hi-D C技术以及下一个诺组件服务。

01 研究思路

02 研究结果

1. 源四倍体象草基因组的从头组件及注释

通过K-mer分析和流式细胞仪,象草的基因组排列预测为2.0 Gb和2.13 Gb。另外,通过K-mer分析的综合比例为1.3%。最终为2.22 Gb的组成序列当中,N50组合为 25 Mb8con。4N5M 组合 25 Mb。去 Gbtig 获得 N.46 组合后,con2.2.4 组合 N587 组合谱。90 Mb。10 倍的 Bioscaffold 以及 2.0.0 辅助的组合组合,通过了 N50 的基因组列表。 。 加载系数为 6.88% 的 N50146.844,BUSCO 97.8% M-9% 的评估为 C 77,19 帮助预测的一个。

图1 异源象草基因组圈图

2.异源四倍体象草的进化分析

禾本科植物和树种基因有不同的来源和类似的基因系统(图2)。三。禾本科植物的组合中线发生融合22马,植物形成组中线(图2b)。 2的编码基因’B亚的重建和染色体家族的后裔,1条染色体家族的共同支持,母体之间有2个基因组之间的关系。一个“基因组基因组整合”的源体组合(图2c和)似乎出现了3次多倍组化事件。中,一个亚组11492显着发现,1719个组,而象草B亚0031这些基因家族有2个发现有1个下来)。在进化过程中保留家族链(图)基因家族的家族基因组丰富的或能力亚基具有重要的分析能力。以及过氧化物的家族基因。

图2 狼蛛相关的尾草进化属中象草

3. 象草的生物量特性相关分析

在A’和B亚组中(图3b)的两个结果中,组和B亚组中的两个和三个重复分别出现的结果,A’B亚组的基因组序列相同。积分。

图3 象草的合成基因表达和速热图

136个基因的基因发现的木质素合成相关,个数源可能是由于对四象体形成了5亚种的加成。(图4b,通过WGCNA分析)鉴定了664个木质素合成的有加。

图4 象草和木素热合成的基因表达图

03 小结

本文利用三代长基因组基因技术以及辅助组件技术,象草生物量巨大的潜在成分组,为研究象草组的提供了组合、组合、组合、组合、多种组合化的重要和进化也为等方面提供了基础,其他开发利用的资源。

平台升级 | 希望组HiFi产能扩大50%,交付快人一步!

为满足市场需求,给客户带来更高效更经济的服务,2022年4月22日希望组再次引进PacBio公司的最新Sequel IIe测序平台,HiFi产能扩大50%,目前,希望组已经累计构建PacBio Sequel/Sequel I I/IIe达12台。测序平台产能瓶颈解除之后,希望组将持续助力更多的动植物基因组完成图、三代群体重测序研究以及医学诊断等领域的研究与应用。

Sequel IIe真容

相较于以往Sequel II机型,Sequel IIe测序仪在测序原理上并无差别,但在硬件和软件上进行了一些改进和优化,用户可在Sequel IIe系统上直接获得高品质的HiFi 测序数据,总体上能够缩短用户近70%的下游分析时间,同时能够减少近90%的数据传送时间及节约相应的存储空间。与此同时,基于Sequel IIe测序仪的Read population可视化优化后可直接看到reads在长度和准确度上的空间分布情况及subread passes数对最终reads质量的影响情况,更快一步看到样本测序初步质量。

既快又好,希望组的HiFi测序数据质量及相应基因组组装结果是这样的!

1. 高产保质的HiFi测序

HiFi数据量和长度统计

近半个月来,希望组承接的40个HiFi测序项目中,除去个别特殊处理样本,每个cell的HiFi数据产出量均达25G以上!并且统计结果显示,产出数据的N50长度均能达到15kb左右,部分样本甚至可达20kb!

2. 快速高效的HiFi基因组组装

HiFi基因组组装统计

近期希望组完成的HiFi基因组组装案例,覆盖了真菌及各种大小型动植物,且组装的各项指标均达到预期!对最终组装结果进行BUSCO评估,用以评估基因组组装、基因集和转录本的准确性和完整性,各项目的BUSCO评估结果显示Complete BUSCOs 比例均达95%以上,相关基因组组装的完整性非常高!

3. HiFi宏基因组组装——HiFi测序应用的另一块高地

目前宏基因组测序方式主要是短读长的二代测序和长读长的三代测序。二代测序获得的宏基因组组装由于测序读长的限制,分析结果往往出现基因信息不完整、种群分类分辨率偏低、只能到属水平、组装连续性小等不理想状况。新升级的Pacbio Sequel II及Sequel IIe HiFi reads,兼顾长读长和高精准度,完全可以解决以上问题,通过长且准的测序方式提升组装完整性和连续性。HiFi宏基因组,在保证准确度和组装质量的前提下,提升了组装速度,为后续深入分析研究奠定基础。目前希望组正在基于华为云平台开发自动化分析流程,将进一步缩短分析流程,加快项目周期,为科研工作者和分析人员提供更便捷精准的平台,敬请期待!

项目文章 | Nature Plants封面-苏铁基因组发布

种子植物包括裸子植物(gymnosperms)和被子植物(angiosperms),裸子植物分为四大类,即苏铁类(cycads)、银杏类(Ginkgo)、松柏类(conifers)和买麻藤类(gnetophytes)。裸子植物基因组较大,重复序列含量高,结构复杂,迄今为止,现存最原始种子植物苏铁分支尚缺少完整的基因组图谱。

2022年4月18日,由22个机构65位科学家联合在Nature Plants发表了题为“The Cycas genome and the early evolution of seed plants”的封面文章,报道了现存最原始种子植物苏铁参考基因组,填补了种子植物基因组研究的空白。苏铁基因组的发布,代表着种子植物基因组演化研究中的最后一块拼图完成,为后续比较基因组学的开展奠定了基础。希望组参与了本研究项目中攀枝花苏铁的测序、组装及初步注释服务。

苏铁基因组测序材料选取苏铁类的基部类群、也是整个苏铁类分布纬度最北的种类攀枝花苏铁(Cycas panzhihuaensis)。基于长片段测序与MGI-SEQ测序,苏铁基因组组装大小为10.5 Gb,contig N50为12Mb,结合Hi-C数据,挂载为11条染色体。其中共注释32,353个蛋白编码基因,BUSCO评估完整度为91.6%,是目前裸子植物中最高质量的大基因组图谱。

裸子植物具有4大分支,1118种,关于裸子植物内部大分支之间的系统演化关系一直有不同的学术观点。基于15个维管植物基因组3282个直系同源低拷贝核基因、90个种子植物转录组1569个直系同源基因、72种维管植物叶绿体和线粒体基因组数据的系统分析结果表明,苏铁单独(线粒体数据)、或和银杏一起(核基因、叶绿体数据)构成其它所有裸子植物的姐妹群。

基因组加倍是植物演化适应的重要驱动力,关于裸子植物共同祖先是否经历了全基因组加倍事件一直存在争议。研究者采用对重复基因同义替代分析和系统发育基因组学方法,并使用基因组内共线性区域进行比较验证,发现现存裸子植物的最近共同祖先可能经历了一次古老的全基因组复制事件(命名为ω,图一a)。伴随着种子植物起源,许多关键创新性状如种子发育、花粉、次生生长相关的基因家族均发生了创新或扩张。在种子植物的祖先节点共发现663个新获得的基因家族和368个扩张的基因家族。其中,106个新获得和55个显著扩张的基因家族与种子生理发育有关,包括调控胚胎早期发育、种子休眠和萌发、种子能量和营养代谢,种皮形成以及种子的免疫和应激反应等(图一b)。

图一、a) 基于系统发育关系推断种子植物的全基因组加倍事件;b) 种子植物的基因家族创新和扩张。

最显著扩张的种子生理相关家族是cupin蛋白家族。攀枝花苏铁编码一类新的vicilin-like贮藏蛋白 vicilin-like antimicrobial peptides(v-AMP),在基因组中呈串联基因阵列分布,多在授粉胚珠后期和受精胚珠时期表达,而后逐渐降低,暗示v-AMP基因在种子发育过程特定时期发挥重要作用。LAFL家族(LEC1ABI3LEC2FUS3)是种子发育核心调控基因,苏铁等裸子植物的FUS3LEC2基因可构成一个新的进化枝,定义为FUS3 / LEC2-like类型,与被子植物的FUS3LEC2形成姐妹分支关系。FUS3 / LEC2-like类别是裸子植物特有的。在攀枝花苏铁授粉后,其会表现出明显的表达,表明可能在裸子植物胚胎发生早期发挥特定作用 。

苏铁类起源于古生代二叠纪早期,距今已有至少2亿7千万年历史。在经历大量灭绝以后,现代苏铁多是近期几次辐射演化的后代。如今苏铁具有2科10属。研究者基于现存苏铁目339种植物的转录组数据,重建了苏铁类自身的系统发育关系。分子钟分析表明,现存苏铁的多样化同步发生于距今1100至2000万年之间,是中新世以来气候剧烈变化的结果(图二)。

图二、苏铁目系统发育树支持现存苏铁是辐射演化的结果

雌雄性别分化是一种进化性状。裸子植物已报道的1118个物种中,65%的种类都是雌雄异株, 而在被子植物中雌雄异株只占到6%。苏铁类植物均为雌雄异株(图三),由于其生长缓慢,以往只能在开花时才能判断性别,而成株树龄多在十年以上以及在适宜的环境才会开花。研究苏铁性别决定的分子机制,可以在植株进入花期前就确定性别,对于苏铁类植物的就地和迁地保护和园林培育具有重要意义。苏铁的性别控制基因一直困扰科学界。研究者通过对源于四川攀枝花苏铁国家级保护区62株雌雄苏铁群体测序,表达差异分析,和雄性Y染色体的组装,找到雌雄表达差异最大的一个基因来自雄株的Y染色体,该基因编码一个MADS-box转录因子,推测其调控雌雄苏铁的性器官发育,揭示了苏铁性别决定的遗传机制。该转录因子的同源基因也仅能在雄株基因组中检测到,说明了该性别决定机制在苏铁类植物中的保守性。

图三、a) 攀枝花苏铁;b) 攀枝花苏铁雄株和雌株的孢子体。

早期维管植物的精子都是有鞭毛,可以游动的。随着演化,鞭毛丢失。在现生种子植物中仅苏铁和银杏保留精子具鞭毛的特征。研究人员发现,苏铁和银杏均保留了大量鞭毛组装所需基因,但与苏铁相比,银杏中RSP类基因有一定的丢失(RSP2RSP3RSP 9, 和RSP11等)。此外,与鞭毛行使功能密切相关的外周致密纤维合成基因(ODFs),只在苏铁和银杏基因组中存在,在其它种子植物中则全部丢失。进一步证实了苏铁在种子植物演化中古老的地位。
水平基因转移泛指不同物种间的基因交流,在陆生植物适应性进化过程中起到推动作用。研究者在攀枝花苏铁基因组中发现一种细胞毒素蛋白基因(fitD),这种基因起源于细菌,通过水平基因转移的方式转移到真菌和苏铁中(图四a)。基于苏铁类339个物种转录组数据,研究人员发现该毒素蛋白只在苏铁属物种中存在。同时,fitD基因在种子和根部高表达,这可能也是苏铁种子和根部具有毒性的原因之一。基于基因重组技术,在大肠杆菌表达的毒素蛋白产物对小菜蛾和棉铃虫有显著致死性(图四b-f),显示出毒素蛋白具有一定的农业应用前景。

 图四、a)苏铁水平转移毒蛋白基因的演化历史。b-f) 苏铁毒蛋白基因表达,及对昆虫毒性实验。

该研究由深圳华大生命科学研究院、深圳市仙湖植物园、中国科学院昆明植物研究所、兰州大学、中国环境科学研究院、河南大学和南京林业大学等22个机构65位科学家联合完成。该论文第一作者为深圳华大生命科学研究院刘阳、王思博、李林洲、杨婷、魏桐,深圳仙湖植物园董珊珊,兰州大学武生聃等为共同第一作者,分别在基因组不同的领域贡献了自己的专业力量。深圳仙湖植物园张寿洲,深圳华大生命科学研究院刘欢,中科院昆明植物所龚洵,美国佛罗里达大学Douglas E. Soltis,比利时根特大学Yves Van de Peer为文章共同通讯作者。该项目得到深圳市城市管理和执法局科研专项、国家重点研发计划、生态环境部生物多样性调查与评估等基金支持。

项目文章 | 华中农大联合深圳基因所组装康乃馨染色体级别基因组并结合多组学数据解析花色、花型和花香等重要性状调控密码

康乃馨(Dianthus caryophyllus)是石竹科石竹属多年生植物。因其花色绚丽,花型独特,香气馥郁,深受世界各地人民的喜爱,被称作世界‘四大切花’之一,具有极高的观赏价值和经济价值。

2022年3月5日,华中农业大学园艺林学学院傅小鹏研究团队,联合中国农业科学院农业基因组研究所(简称‘基因组所’)武志强研究团队,在国际著名学术期刊Plant Biotechnology Journal上发表了题为“Integrated multi-omic data and analyses reveal the pathways underlying key ornamental traits in carnation flowers”的研究论文。希望组在本文章中提供了基因组测序及Hi-C测序服务!

该研究首次组装了康乃馨染色体级别的基因组,并结合比较基因组,转录组,代谢组等方法,系统地研究了康乃馨基因组的进化,并对康乃馨的花色,花型和花香等重要观赏性状的分子机理进行解析,鉴定了参与这些性状形成的关键基因。本研究的开展为定向改良康乃馨的观赏性状,培育康乃馨新品种奠定了基础。

研究团队利用第三代ONT、二代Illumina和HIC技术,对康乃馨栽培品种‘斯嘉丽’( D. caryophyllus, cultivar ‘Scarlet Queen’)进行了全基因组测序和组装。康乃馨(2n=30)的基因全长为636.30 Mb,重复序列占70.62%,共注释到43,925个基因,Contig N50 和Scaffold N50分别为14.67 Mb和38.55M。研究发现康乃馨与藜科在~64.07 MYA发生了分化,并经历了一次近期的WGT事件,这次加倍事件促进了康乃馨特殊香气丁香酚合成相关的基因的扩张。

图1.康乃馨基因组的组装与进化

基于本基因组,对康乃馨红边花瓣进行代谢组分析,发现天竺葵素的积累促使花瓣红边的成色,叶黄素和多个类黄酮物质的共同积累,促使花瓣基部成黄色;结合转录组分析,筛选到几个MYBs, bHLHs 和WRKY44等基因共同作用于ANS调控康乃馨复色的形成。同时,利用GC_MS以及转录组,对香石竹丁香酚在不同时期花瓣的积累过程进行解析,发现EGS基因存在外显子上的结构变异,其提前终止可能是康乃馨丁香酚丢失的原因。另外,研究发现康乃馨汇总A、C类基因的异位表达可能是影响重瓣型康乃馨形成的重要因素。

图2. 康乃馨花瓣花边颜色的形成

图3.潜在的康乃馨中丁香酚合成通路

 图4.康乃馨重瓣花形成的遗传机制

本研究为开展康乃馨重要观赏性状的调控和遗传改良提供了宝贵的数据信息和坚实的理论基础,对推动康乃馨分子育种育种具有里程碑式的意义。

华中农业大学博士张晓妮(现为基因组所武志强组博士后)为论文第一作者,华中农业大学傅小鹏教授和基因组所武志强研究员为共同通讯作者。该研究是在包满珠教授和BENDAHMANE教授全程指导和帮助下完成的。同时,该研究得到了国家自然基金、国家青年科学基金、深圳市科创委以及基因组所启动基金等的支持。

原文链接:https://onlinelibrary.wiley.com/toc/14677652/0/ja

文章来源:植物生物技术Pbj公众号

进展迅速丨高质量人类二倍体基因组组装策略来临

自T2T联盟于2021年5月在BioRixv预印版上发表了首个人类基因组完成图及相关文章之后,近日人类泛参考基因组联盟紧跟步伐在BioRixv预印版上发布了题为“Automated assembly of high-quality diploid human reference genomes”的文章。

样本来源

目前常用的人类参考基因组GRCh38,是科研团队花费20多年的努力而完成的高质量的人类基因组。该基因组极大程度上推进了生命科学相关领域的发展。但是它仍然存在着许多gaps(gap数已从原来的150,000个降低至了现今的995个)和错误。因为GRCh38是20位匿名志愿者的混合样本,所以并不代表单个个体的人类基因组。2021年,T2T联盟利用最新的长读长技术与自主研发的生信算法组装出了高质量端粒到端粒的参考基因组CHM13(首个人类基因组完成图)。由于它来源于一个带有重复基因的葡萄胎细胞系,因而几乎是纯合基因组。在组装难度上,纯合基因组比二倍体基因组组装相对较低。为了解决这些原有人类基因组序列的局限性,研究团队成立了人类泛参考基因组联盟(Human Pangenome Reference Consortium,HPRC)。该联盟旨在整合高质量、低成本、二倍体人类基因组,用以构建代表人类遗传多样性的人类泛参考基因组。

人类泛参考基因组联盟使用的是二倍体基因组,组装难度会更高。在样本选择上,选择了广泛使用的HG002的永生化淋巴细胞GM24385,用于ONT测序和Bionano光学图谱测序。二代illumina测序和PB测序数据使用的样本是由大量GM24385细胞而来的NIST标准样本8391 DNA样本。父系样本HG003和母系样本HG004分别来源于公开的GM24149和GM14143细胞系,并通过二代illumina测序对NIST标准样本8392(含HG002、HG003和HG004)的DNA进行测序。

组装策略及结果

本文是以HPRC的名义发表的首篇文章,该文章确定了当前基因组测序和automated组装方法的组合可以在人工干预最少的情况下产生最完整、最准确和最经济的二倍体基因组组装。在组装过程中,团队选择了最优的高精度长读长reads和父母本及子代数据进行单倍体分型组装(trio-based assembles)的方法。第一个高质量二倍体人类参考基因组(XY型),平均每条染色体仅有~4个Gaps(范围在0-12),大多数染色体长度和CHM13相比仅有±1%的差异。将近四分之一的蛋白编码基因在单倍型间存在同义氨基酸变化,而其中着丝粒区域变异的密集度最高。该研究为构建涵盖从单个核苷酸到大型结构重排的所有遗传变异的人类泛参考基因组奠定了基础。

测序及组装方法的测试总结

HPRC Trio pipeline v1.0组装流程图

小结

使用目前已有的方法和技术,研究者们可将二倍体HG002基因组的组装达到整体完整度的98.5%。利用这些方法组装不同的二倍体人类基因组,构建人类泛参考基因组,将会更全面地了解人类遗传多样性,提高精准医疗的准确性,并对生物基因组学有更深入的了解。

美国学者Erich D. Jarvis为论文的第一作者,中南大学王建新团队、中山医眼科肖传乐团队、中国农科院阮珏团队、昆明动物所张国捷团队以及华大基因杨琛涛团队均参与该项项目研究。

希望组一直致力于推动各个物种基因组完成图的组装及应用。近期还为动植物基因组完成组量身定制包含测序及组装在内的全套解决方案,敬请期待!

原文链接:https://doi.org/10.1101/2022.03.06.483034

基因组T2T完成图新工具丨BAC-long(>150Kb)测序技术正式推出!

2022年伊始,Nature杂志发布年度前沿技术展望,其中人类基因组T2T(端粒到端粒)完成图成为7大技术展望之首。

人类基因组T2T完成图

当Karen Miga和Adam Phillippy在2019年启动端粒到端粒(T2T)联盟时,大约有十分之一的人类基因组还没完全绘制出来,而今该数字已降至为零。该联盟2021年5月在BioRixv上发布的预印版中报道了第一个人类T2T基因组序列,有将近2亿个新的碱基对加入了广泛应用的人类基因组参考序列中,成为了一个全新的T2T-CHM13基因组版本,撰写了人类基因组计划的最终章。

从人类基因组端粒到端粒的完成图发布,到已发布的拟南芥基因组完成图,物种基因组完成图时代已经来临。

科学的进步一直伴随着技术的突破,人类基因组T2T项目的完成也不例外。近年来,长读长测序技术持续改进,陆续打破了传统的测序技术极限,成为了基因组T2T完成图组装的必备工具。其中,PacBio公司开发的HiFi测序技术,可以得到超过20Kb,准确度在99.9%以上的测序序列,成为了基因组组装领域的必备工具;ONT公司开发的Ultralong测序技术,可以将DNA测序长度扩展到N50 100Kb的极限,成为基因组T2T组装的首选技术。

但是,从人类基因组到模式生物基因组,以及逐步拓展到非模式生物基因组T2T组装,将面临更大的挑战:

1)二倍体基因组全染色体分型问题;
2)二倍体/多倍体基因组高杂合度问题;
3)超大基因组(>10Gb)带来的高重复问题;
4)超高/超长重复序列组装问题;
5)新进化产生的大片段重复(Segmental Duplications)识别问题;
6)一些复杂物种的大规模STR识别问题;
7)一些植物的高占比的着丝粒和端粒的识别问题;

这些问题的解决,也许需要采用比现有的人类基因组T2T项目更前沿的技术。

为了全面推动基因组组装进入T2T完成图时代,希望组(GrandOmics)历经近5年时间的持续研发,在逾10,000张ONT测序芯片的实战基础上,特别测试了超过1,200张超长DNA测序。在此基础上,研发出希望组自主知识产权的BAC-long(>150Kb)试剂盒,并正式面向全球提供BAC-long(>150Kb)测序服务。

测序技术里程碑——BAC-long测序

短读长组装到BAC-long组装

BAC文库构建可以插入平均150Kb的长片段DNA,是早期基因组组装的主要工具。2001年,人类基因组草图完成,国际人类基因组联盟主要采用的即是BAC-Based技术。每完成一个BAC的测序,都需要付出大量的人力、物力、时间和经费,最终整个项目整整持续了11年,花费超过了27亿美元。

2021年发布的CHM13-T2T人类基因组序列是研究者们耗时数年,利用长读长测序技术(PacBio平台的HiFi技术以及ONT平台的Ultralong(>100Kb)技术)完成的第一个零GAP人类基因组组装序列。但是,测序样本CHM13本身是加倍的单倍体,拥有几乎为零的基因组杂合度,难度远远低于普通的人类基因组二倍体,以及其他模式和非模式物种的T2T组装。

为了将真正的T2T基因组组装,从人类基因组拓展到更为广阔的领域,包括动植物基因组、农业育种、医学研究、濒危物种保护等,将T2T基因组组装树立为新的行业标准,希望组正式推出BAC-long试剂盒——直接测序即可获得近乎BAC插入片段长度,提供测序读长N50超过150Kb的BAC-long reads,从而可以让测序长度更长,更容易跨过端粒、着丝粒等复杂重复区域,在基因组组装中的表现更优越,为每个物种的基因组完成图助力。

以下为近期希望组BAC-long项目中某植物、某动物的测序reads分布图,显示reads长度更长。

想要获得基因组完成图,除了需要BAC-long reads,跨越复杂重复区域,降低基因组组装拼接难度外,还需要相匹配的完成图组装软件。希望组特地为用户匹配了NextDenovo 2.0系列组装软件,彻底解决完成图组装的后顾之忧。

敬请关注希望组升级软件NextDenovo 2.0!

PJ | 中南林科大、华中农大等联合破译狭叶油茶基因组,为解析油茶自交不亲和、油脂合成等重要性状的形成与调控提供重要参考

2022年3月20日,中南林业科技大学袁德义/张琳油茶科研创新团队,联合华中农业大学金双侠课题组等单位,在经典的国际植物学TOP期刊The Plant Journal上发表了题为“Chromosome-level genome of Camellia lanceoleosa provides a valuable resource for understanding genome evolution and self-incompatibility”的研究论文。希望组在本文章中提供了基因组测序及Hi-C测序服务!

油茶是我国重要的木本食用油料树种,在推进地方经济绿色增长、维护国家食用油安全和乡村振兴中占有十分重要的地位。茶油中不饱和脂肪酸含量达90%以上,还含有丰富的角鲨烯、维生素E、谷甾醇等物质,具有软化血管、降低血脂和血压的作用,是目前国内外最优质的植物食用油。目前,油茶主栽品种主要为多倍体的普通油茶(Camellia oleifera),由于油茶基因组大且亚基因组间的同源异源构成非常复杂,使得多倍体油茶基因组的解析非常困难,严重阻碍了油茶的分子遗传改良。狭叶油茶(Camellia lanceoleosa)是油茶组唯一的二倍体野生种,和多倍体普通油茶亲缘关系最近,破译狭叶油茶基因组不仅可以深入挖掘油茶资源中的优异性状,而且为油茶重要功能基因挖掘利用奠定了坚实基础,使得油茶育种不再是盲人摸象,从而开启了油茶分子育种育种时代。

研究团队利用三代Nanopore、二代Illumina和HIC技术,获得了狭叶油茶(2n=30)高质量染色体水平基因组。基因组大小约为3.00 Gb,杂合率高达2.2%,91.85%的序列被挂载到15条染色体上,共注释到54,172个基因。狭叶油茶重复序列占基因组的80.63%,其中转座子占比78.53%。深入分析发现,长末端反转录转座子家族长期而缓慢的扩增及在过去2百万年内的爆发式扩张,加之缺少快速有效的DNA删除机制,最终导致狭叶油茶基因组变得庞大。狭叶油茶与茶叶均为山茶科山茶属二倍体植物。比较基因组研究发现,狭叶油茶与茶叶共享最近的一次WGD事件,并在6-7百万年前发生分化,2号和11号染色体在狭叶油茶与茶叶间存在较大的倒位,这可能是油茶和茶叶在基因组结构上的一个重要变异。

图1 狭叶油茶基因组的组装

图2 狭叶油茶转座子的插入与移除

图3 狭叶油茶比较基因组分析

基于狭叶油茶基因组,结合GC_MS以及转录组分析,发现油脂合成的关键基因ACC的扩张及DGAT、GPD、SAD在种子的偏向性表达是狭叶油茶高油脂和高油酸含量的重要原因。儿茶素、茶氨酸、咖啡碱是茶叶品质的重要组成成分,基于UPLC-MS/MS技术的代谢组检测到狭叶油茶叶片也富含儿茶素和茶氨酸,而咖啡碱主要富集在种皮和根中,多组学分析结果表明,SAM-dependent N-methyltransferases与咖啡碱的积累与分布密切相关。细胞学分析显示,狭叶油茶也是后期自交不亲和植物。结合亲和性相关基因的表达、结构特征及染色体定位,解析了狭叶油茶的自交不亲和性特征。狭叶油茶的油脂含量、脂肪酸比例以及次生代谢物分布与含量都与普通油茶类似,以上性状的解析为理解普通油茶油脂合成等重要经济性状的形成与调控提供了重要参考。

图4 狭叶油茶自交不亲和特征解析

中南林业科技大学龚文芳副教授、肖诗鑫讲师及在读硕士生王林凯为论文共同第一作者,中南林业科技大学袁德义教授、张琳教授和华中农业大学金双侠教授为论文共同通讯作者,黄冈师范学院朱华国教授、胡孝明教授及中国农业科学院农业基因组研究所廖振阳博士后也参与了本项研究。本研究得到了国家重点研发计划项目(2018YFD1000603)、湖南省自然科学基金(2020JJ5968)的联合资助。

原文链接:

https://onlinelibrary.wiley.com/doi/epdf/10.1111/tpj.15739

技术突破丨希望组“10X+ONT-Q20+”长读长单细胞测序实测通过!!!

2021年ONT在科研团体大会(Nanopore Community Meeting)上发布了”Q20+”试剂产品。这是利用包含最新“Q20+”化学试剂并支持“双链”(Duplex)测序的新试剂盒LSK112和最新芯片R10.4相结合的测序技术,可以实现>99%的原始reads(单链)准确度,或约Q30的双链(Duplex)准确度。

希望组(GrandOmics)作为中国最早开展ONT测序服务的公司之一,一直走在技术应用的前沿,率先将高通量单细胞分选平台10X和ONT最新技术进行组合创新,给火热的单细胞测序领域,提供了全新的技术解决方案。

测试背景

2020年8月于Nature Communication杂志发表的题为“High throughput, error corrected Nanopore single cell transcriptome sequencing”的研究论文,将单细胞与Nanopore全长转录组进行了技术整合,使得在单细胞水平上检测差异RNA剪接和RNA编辑成为可能。利用Illumina对10X Genomics文库进行短读长测序,为每个基因和基因组区域定义相关的Cell Barcode,随后为每个细胞、基因或基因组区域定义相关的UMIs组合。进而利用上述信息指导Cell Barcode和UMI分配到已比对至基因组上的Nanopore reads,通过对相同UMI的reads进行分组,校正和拆分Nanopore测序reads,用于开展下游的isform分析和RNA编辑。考虑到ONT-Q20+测序即将成为行业新的长读长测序技术标准,我们迫不及待的进行了新的尝试。

测试方式

10X单细胞平台获得的全长cDNA分成两份,一份进行常规二代测序,另外一份进ONT Q20+R10.4测序,利用二代测序数据的UMI指导进行 ONT测序的数据拆分。

测试样本

某人类PBMC样本,采用最新GPU进行Super高精准度Basecalling。尽管大大提高了计算消耗,但预期可以得到最佳的单碱基质量。

测试结果

一、下机数据统计

下机数据质控统计表

Reads长度分布图

单张测序芯片的产量为36.7Gb,pass产量为29Gb。尽管产量仍有提升空间,但通过芯片及算法的升级,basecalling的准确度有了大幅度的提升,reads平均分可达到15.98,最终获得了20 Mb的reads进行下游分析,和其他长读长平台所产生的reads数相比,有着明显的产量优势。该芯片样本为单链测序样本,相信在”Duplex”的模式下,样本的测序准确度会有更进一步的提升。

2. 与参考基因组比对identity

左图Q20+R10.4的比对identity分布图,右图为LSK110 R9芯片的比对identity分布图

从图中可以发现,同个样本不同的测序试剂和芯片,对于最终比对identity的准确性有着显著的差异。LSK110 R9 base-calling fast模式的比对identity中位值在87.1%,Q20+R10.4的比对identity中位值达到了97.9%。整体数据的准确度有了明显的提升。

3.数据分解率

随着数据准确性的不断提高,利用二代UMI数据指导拆分三代全长转录组数据的比例也会有相应的提升空间,内部数据表明二代UMI数据指导拆分三代全长转录组数据的拆分率可达70%,远远高于Nature Communication文章里50%的拆分率。相信在“Duplex”模式下的数据表现会有更加亮眼。敬请关注希望组后续的相关报道!

当10x单细胞样本遇上Oxford Nanopore Q20试剂及R10.4芯片,有效提升UMI的数据拆分率,拥有更多有效数据进行下游分析。尽管“Q20+”目前处于试测阶段,但相信未来在基因组、转录组上的应用场景会不断增加。希望组早在2020年年底已经推出了利用纳米孔平台进行单细胞全长转录组的产品SCAN-seq。现在有了Q20试剂的加持,拓宽了单细胞全长转录组的应用场景和需求,可以为广大科研工作者提供更好、更优质、更前沿的技术!

项目文章丨PNAS封面!李家堂团队等揭示“飞蛙”滑翔的遗传机制

近日,中国科学院成都生物研究所李家堂研究团队等以云南西双版纳地区的“飞蛙”——黑蹼树蛙为研究对象,通过多维度研究手段,阐明了其攀爬和滑翔行为相关表型的遗传机制。相关研究在《美国科学院院刊》以封面论文形式发表。希望组在本研究中提供了三代测序组装、注释以及Hi-C测序、挂载服务。

动物复杂性状是动物长期适应演化的结果,是动物多样性形成的重要基础。自然界物种采取各种各样的进化策略以适应不同的栖息地,如高原、海洋、荒漠等。部分类群通过演化出了特殊的表型以适应树栖生活。树栖生活拓展了这些物种对垂直空间资源的利用,有助于它们躲避天敌,获取丰富的食物资源等。但森林环境复杂的立体结构也对动物的运动能力提出了严苛的要求。

在白垩纪末期大规模物种灭绝事件后,无尾目多个科的物种独立演化出了攀爬和滑翔的相关表型,并成功拓殖树栖生态位,黑蹼树蛙就是其中的代表性物种。黑蹼树蛙隶属于树蛙科,常年生活在热带雨林树冠层,是典型的树栖蛙类。据报道,其最高栖息高度达57米,为目前树栖蛙类停留高度的最高记录。黑蹼树蛙具有强大的滑翔能力,因此又被称为“飞蛙”,为揭开两栖动物树栖适应之谜提供了良好的动物模型。

 黑蹼树蛙,饶涛绘

近日,中国科学院成都生物研究所李家堂研究团队等以云南西双版纳地区的“飞蛙”——黑蹼树蛙为研究对象,通过多维度研究手段,阐明了与攀爬和滑翔行为相关表型的遗传机制。

研究团队选择了蹼发达的黑蹼树蛙(满蹼)和宝兴树蛙(微蹼)开展行为学实验。研究人员分别设置了1,1.5和2米的高台,让两种树蛙分别落下并记录运动轨迹。

结果发现
黑蹼树蛙在下降过程中会尽力撑开四肢,依靠满蹼维持空气动力学平衡,使其身体与水平面夹角始终小于宝兴树蛙。而宝兴树蛙则更像是在“跳楼”。这一结果提示树蛙的蹼在滑翔过程中起到重要作用,是支撑滑翔行为的重要性状。

两种树蛙形态学比较与行为学实验

两种树蛙蝌蚪发育时期的肢体形态学比较提示两者之间存在不同的发育模式。以两种树蛙高质量基因组为基础,结合蝌蚪四肢发育过程中的转录组数据,通过时序基因共表达网络分析,发现黑蹼树蛙蝌蚪肢体在发育阶段特异共表达一系列与Wnt信号通路和血管重构相关的基因。这种发育时期特异表达模式可能通过参与指和指间区域生长速率的调控对蹼足的形成起到关键作用。

两种树蛙蝌蚪肢体发育过程形态学比较与时序基因共表达网络分析

树蛙中调控角蛋白和细胞骨架形成的PPL基因受到正选择,且存在树蛙属内保守的氨基酸替换,这可能有助于树蛙攀爬相关性状——吸盘的形成。壁虎的刚毛主要由β角蛋白扩张形成,而树蛙指/趾尖吸盘的主要结构蛋白为α角蛋白,而α角蛋白同样是哺乳动物和人类毛发纤维最主要的成分之一。

“这些祖先角蛋白可能在早期四足动物需要皮肤强化的区域中表达,而随后分化为支持两栖动物和哺乳动物不同的适应性结构”。

本研究通过前期大量的野外观察,整合形态学、行为学和组学等学科研究手段,系统解析了“飞蛙”滑翔的遗传机制,为人类认识动物有趣行为提供新的视角,对动物特殊功能的仿生研究及人类并指症等相关疾病的防治有重要基础科学价值。

该研究成果以“Genomic adaptations for arboreal locomotion in Asian flying treefrogs”为题在《美国科学院院刊》以封面论文形式发表(点击左下角阅读原文)。本研究得到中国科学院B类先导科技专项、国家自然科学基金等项目的资助。