低输入PacBio测序生成高质量个体果蝇基因组并表征突变过程

Low-input PacBio sequencing generates high-quality individual fly genomes and characterizes mutational processes.

作者信息Hangxing Jia, Shengjun Tan, Yingao Cai, Yanyan Guo, Jieyu Shen, Yaqiong Zhang, Huijing Ma, Qingzhu Zhang, Jinfeng Chen, Gexia Qiao, Jue Ruan, Yong E Zhang
PMID38969648
发布时间2024-07-05
DOI10.1038/s41467-024-49992-6

实验完整度

高

包含体外建库方法开发、两只果蝇个体基因组与Wolbachia基因组组装、SV/SNP及基因转换等多层验证

主要模型

黑腹果蝇ISO1品系单只雄性成虫 内共生细菌Wolbachia wMel 果蝇家庭 multiplexed 测序样本

重点核对

单只雄性果蝇ISO1-1和ISO1-2的gDNA输入量为100 ng PacBio Sequel II平台每个个体使用一个测序cell,产生7.7–8 GB CCS reads LILAP建库中Tn5转座体用量0.5 μl,T4 DNA连接酶缓冲液为通用缓冲液 SV和SNP分别通过组装水平分析和CCS reads支持(≥3条reads)进行验证 两个果蝇家庭multiplexed文库用于验证SV/SNP在家庭中的存在

摘要

以PacBio为代表的长读长测序通过克服重复序列等挑战正在彻底改变基因组学。然而,其对DNA的高需求(>1微克)使得小生物难以企及。我们开发了一种低输入(100纳克)、低成本且无扩增的PacBio测序文库构建方法(LILAP),该方法在同一管内利用基于Tn5的标签化和DNA环化。我们用两只黑腹果蝇个体测试LILAP,生成了近乎完整的基因组,超越了先前已有的单只果蝇基因组。通过分析这两个基因组中的变异,我们表征了突变过程:复杂转座(转座子插入伴随额外重复和/或缺失)偏好以非B DNA结构为特征的区域,转座子的基因转换在DNA和RNA水平上均发生。同时,我们为这些果蝇中的内共生细菌Wolbachia生成了两个完整组装,并同样检测到转座子转换。因此,LILAP为果蝇及其共生生物的突变研究以及其他小生物或珍贵样本的探索提供了广泛采用PacBio测序的前景。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
能否开发一种低DNA输入、低成本且无扩增的PacBio文库构建方法,以实现小生物个体及其共生菌的高质量基因组组装,并据此研究短读长难以解析的突变过程?
核心机制
复杂转座(转座子插入伴随额外重复和/或缺失)偏好非B DNA结构区域,可能由易错DNA修复(如MMBIR)介导;转座子的基因转换在DNA和RNA水平发生,贡献了22.5%的聚集SNP。
主要证据
使用100 ng gDNA和LILAP生成两只ISO1雄性果蝇的PacBio HiFi数据,组装出QV 63–65、覆盖99.6%主要染色体的基因组;分析337个SV和1756个SNP,鉴定9个复杂转座和58个TE转换事件,并通过aISO1-Anno数据验证大部分事件;同时组装两个完整Wolbachia基因组并检测到TE转换。
研究意义
LILAP有望广泛用于小生物或珍贵样本的PacBio测序,推动果蝇及其共生生物的突变研究,并为其他小生物基因组探索提供途径。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

LILAP文库构建方法开发

建立低输入、低成本、无扩增的PacBio环化文库构建流程

在同一管内完成Tn5介导的标签化、外切酶消化多余接头、DNA环化与损伤修复,使用发夹接头和T4 DNA连接酶缓冲液作为通用缓冲液

2

单只果蝇基因组测序与组装

评估LILAP能否从单只果蝇的100 ng DNA生成高质量基因组

对两只ISO1雄性成虫提取gDNA,构建LILAP文库,PacBio Sequel II HiFi测序,使用hifiasm组装并抛光,评估QV、基因组覆盖度和基因覆盖度

3

结构变异鉴定与验证

鉴定两个个体基因组中的SV并验证其真实性

使用svmu等工具鉴定SV,通过CCS reads和aISO1-Anno扩增文库验证,分析TE插入类型及复杂转座事件

4

SNP与基因转换分析

鉴定SNP并检测TE基因转换事件及其对cSNP的贡献

使用PAV鉴定SNP,通过CCS reads过滤,定义cSNP并搜索TE转换供体,分类 proximal/distant/interarm 转换事件

5

Wolbachia基因组组装与突变分析

评估LILAP能否同时组装共生菌基因组并检测其突变

对饥饿处理后的果蝇gDNA测序数据组装Wolbachia完整基因组,与wMel比较鉴定突变,分析ltrA基因转换

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
蛋白酶KTIANGENRT403
磁珠动物组织基因组DNA提取试剂盒TIANGENDP341
Tn5转座酶Mei5 BiotechnologyMF650-01
T4 DNA连接酶缓冲液NEBB0202S
外切核酸酶INEBM0568L
λ外切核酸酶NEBM0262S
T4 DNA聚合酶NEBM0203L
T4 DNA连接酶NEBM0202L
DNA聚合酶INEBM0209S
尿嘧啶DNA糖基化酶NEBM0280S
甲酰胺基嘧啶DNA糖基化酶NEBM0240S
核酸内切酶VIIINEBM0299S
核酸内切酶IVNEBM0304S
T4嘧啶二聚体糖基化酶NEBM0308S
基因组DNA清洁与浓缩-10试剂盒ZYMO RESEARCHD4010
T7外切核酸酶NEBM0263L
外切核酸酶IIINEBM0206L
RecJf核酸外切酶NEBM0264L
外切核酸酶INEBM0568L

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
LILAP文库构建
gDNA输入量100 ng;Tn5转座体0.5 μl;T4 DNA连接酶缓冲液为通用缓冲液;各步温度和时间
阅读提示:Methods中“Library preparation for the single fly”及Fig. 1c
单只果蝇基因组测序
ISO1雄性成虫1–3日龄;饥饿6 h;每个个体一个Sequel II测序cell;CCS reads中位长度约4.5 kb
阅读提示:Methods中“Genomic DNA extraction”和“Library preparation for the single fly”及Fig. 2a
基因组组装与抛光
hifiasm v0.12默认参数;抛光流程Racon、merfin、BCFtools;Sniffles2支持reads数≥3
阅读提示:Methods中“Genome assembly, polishing, and assessment”
SV鉴定与验证
svmu检测SV;aISO1-Anno验证80个共享SV和9个复杂转座;CCS reads支持≥3条
阅读提示:Methods中“SV analysis”及Supplementary Data 5
SNP与基因转换分析
PAV鉴定SNP;VAF阈值常染色体≥0.3、性染色体≥0.9;cSNP距离bins;转换事件至少2个连锁SNP
阅读提示:Methods中“SNP analysis”和“Detection of gene conversion”
Wolbachia基因组分析
饥饿处理6 h;与wMel参考基因组比较;ltrA_7位点长2301 bp
阅读提示:Methods中“Wolbachia genome and meta-genomic analyses”及Fig. 5