跟踪实验进化人群中结构变异的动态
Summary
我们开发了一种具有成本效益的方法来跟踪非单核苷酸多态性等位基因动力学，该方法可以轻松适应实验进化冷冻档案。将三重PCR技术与自动平行毛细管电泳相结合，以量化实验进化过程中插入等位基因的相对频率。
Abstract
结构变异（SV）（即缺失、插入、重复和倒置）现在已知在表型变异中起重要作用，因此在疾病确定或适应新环境等过程中起重要作用。然而，单核苷酸变异比SV受到更多的关注，可能是因为它们更容易检测，并且它们的表型效应更容易预测。短读长深度测序技术的发展极大地改善了SV的检测，但从汇集测序（poolseq）数据中量化其频率在技术上仍然复杂且昂贵。
在这里，我们提出了一种相当简单且便宜的方法，该方法允许研究人员跟踪SV等位基因频率的动态。作为应用的一个例子，我们跟踪细菌实验进化群体中插入序列（IS）的频率。该方法基于围绕结构变异边界的引物三重体设计，使得通过复制野生型（WT）和衍生等位基因产生的扩增子在大小上相差至少5%，并且它们的扩增效率相似。然后通过平行毛细管电泳确定每个扩增子的数量，并将其归一化为校准曲线。该方法可以很容易地扩展到其他结构变异（缺失、重复和倒置）的频率定量，以及自然人群（包括患者体内病原体群体）的池-序列方法。
Introduction
结构变异（SV）是基因组序列的改变，通常影响50 bp或更多。所描述的 SV 的四类是大插入、大删除、反转和重复。直到最近，就其表型效应及其作为疾病遗传决定因素的作用或其对适应的贡献而言，对单核苷酸变异（SNV）的关注比对结构变异的关注更多。这可能是因为检测SNV和预测其表型效应更容易。然而，短读长深度测序技术极大地改善了SV的检测，至少在单个个体或克隆基因组中是这样1。同时，它们的表型效应得到了更好的表征，并且已经记录了许多关于它们作为人类疾病2，3或适应新环境4的遗传决定因素的例子。
缺失和插入通常是由于移动遗传元件（MGE）插入引起的，比单核苷酸多态性（SNP）更具破坏性，并导致移码突变和蛋白质结构修饰。基因内的缺失和MGE插入几乎总是导致基因失活，当插入序列（IS）包含启动子或终止序列时，插入到非编码区域会导致相邻基因的抑制或组成性表达5。虽然必需基因的敲除对细菌适应性有明显的不利影响，但在某些情况下，非必需基因的缺失是有益的。尽管有其固有的成本，但重复也可能是有利的，并参与适应，因为它们导致基因剂量的变化;根据条件6，特定蛋白质活性的增加可能是有利的。
微生物实验进化种群通常从克隆开始。这种遗传多样性的最初缺乏，加上试管的“封闭环境”特征，导致通过水平基因转移和重组获得基因的进化潜力非常有限。在这些特定条件下，对缺失、重复和基因组内 MGE 插入的适应的贡献尤为重要;细菌通常通过功能丧失突变（主要是由于缺失或MGE插入）来适应，影响在稳定，通常营养丰富，单一栽培的人工环境中无用的基因7。在运行时间最长的大 肠杆菌 进化实验中，IS150插入在50，000代后进化的人群中特别频繁，IS元件占突变的35%，这些突变在保留其祖先点突变率的种群中达到高频率8。
进化和重新测序研究将实验进化和下一代测序（NGS）技术相结合，以研究细菌如何在表型和基因组水平上适应不同的环境条件和压力，例如不同的碳和能源，抗生素和渗透应激9，10，11.这些研究通常仅在实验终点获得进化群体或克隆的基因组信息，在某些情况下，在许多中间时间点12，13，14。这些数据提供了对适应给定环境所涉及的基因和途径的见解，但很少允许研究人员随着时间的推移跟踪从头出现和扫荡等位基因的动态。
遵循这些动态的一种方法是选择有限数量的分离感兴趣的等位基因（因为它们影响的基因的功能，因为它们在独立群体中平行扫描等），并使用扩增子测序来量化等位基因比例，在同一测序运行中汇集许多时间点15。该方法已成功用于跟踪实验16 和自然17 微生物种群中小尺寸变异（SNP或1 bp插入缺失）的动力学。然而，在较大的插入缺失或MGE插入的情况下，扩增子的大小差异会引起PCR效率差异，从而扭曲了读取和等位基因比例之间的关系。在某些情况下，两个等位基因之间的大小差异优于扩增子的经典长度。在这里，我们将三重态PCR技术与自动平行毛细管电泳相结合，以根据大小判别量化插入等位基因的相对频率。这种方法允许利用未充分利用的实验时间点来确定新兴突变等位基因的动力学，并以具有成本效益的方式跟踪其固定或丢失的频率。我们应用这种方法来跟踪通过IS10插入突变的新兴 mutS–等位基因，为突变的基因型提供超突变表型。
这种方法需要两个靶等位基因，大小差异为≥5%。首先，引物三联体被设计成产生相似大小的片段，它们共享一个共同的引物。其次，优化PCR条件，并使用野生型（WT）和突变gDNA的混合物产生校准曲线。最后，通过PCR扩增样品，并通过平行定量毛细管电泳定量每个等位基因的相对频率。
Protocol
Representative Results
Discussion
在这里，我们提出了一种具有成本效益的方法，该方法可以跟踪实验进化群体中新兴适应性SV等位基因的动力学。该方法将经典PCR技术和自动平行毛细管电泳相结合，可以确定两个等位基因的相对数量。一旦建立，它允许并行定量许多样品中的等位基因比例，并且比WGS便宜得多。该方法可被视为等效于非SNP突变的扩增子测序，并可解决大SV扩增子测序的技术局限性。我们已经表明，该方法使用克服…
Divulgaciones
The authors have nothing to disclose.
Acknowledgements
这项工作得到了ERC HGTCODONUSE（ERC-2015-CoG-682819）对S.B.的支持，这项工作中使用的数据（部分）是通过蒙彼利埃进化科学研究所的GenSeq技术设施在ANR“Investissements d’avenir”计划（ANR-10-LABX-04-01）的支持下产生的。
Materials
| 96 Well Skirted PCR Plate | 4titude | 4Ti – 0740 | PCR |
| Agarose molecular biology grade | Eurogentec | EP-0010-05 | Agarose gel electrophoresis |
| Agilent DNF-474 HS NGS Fragment Kit Quick Guide for the Fragment Analyzer Systems | Agilent | PDF instruction guide | |
| Buffer TBE | Panreac appliChem | A4228,5000Pc | Agarose gel electrophoresis |
| Calibrated Disposable Inoculating Loops and Needles | LABELIANS | 8175CSR40H | Bacterial culture |
| Dneasy Blood and Tissue Kit | Qiagen | 69506 | DNA extraction |
| Electrophoresis power supply | Amilabo | ST606T | Agarose gel electrophoresis |
| Fragment Analyzer Automated CE System | Agilent | Parallel capillary electrophoresis | |
| Fragment DNA Ladder | Agilent | DNF-396, range 1-6000bp | Parallel capillary electrophoresis |
| GENTAMICIN SULFATE SALT BIOREAGENT | Sigma-Aldrich | G1264-1G | Bacterial culture |
| High Sensitivity diluent marker | Agilent | DNF-373 | Parallel capillary electrophoresis |
| High Sensitivity NGS quantitative analysis kit | Agilent | DNF-474 | Parallel capillary electrophoresis |
| Ladder quick load 1 kb plus DNA ladder | NEB | N0469S | Agarose gel electrophoresis |
| LB Broth, VegitoneNutriSelect Plus | Millipore | 28713 | Bacterial culture |
| Master Mix PCR High Fidelity Phusion Flash | Thermo Fisher Scientific | F548L | PCR |
| Primers | Eurogentec | PCR | |
| Prosize data analysis software v.4 | Agilent | V.4 | Parallel capillary electrophoresis |
| Qubit assays | Invitrogen | MAN0010876 | DNA quantification |
| Qubit dsDNA HS Assay Kit | LIFE TECHNOLOGIES SAS | Q32854 | DNA quantification |
| Thermocycler | Eppendorf | Ep gradients | PCR |
| UVbox, eBOX VX5 | Vilber Lourmat | Agarose gel electrophoresis visualisation | |
| Water for injectable preparation | Aguettant | PROAMP | PCR |
Referencias
- Mahmoud, M., et al. Structural variant calling: the long and the short of it. Genome Biology. 20 (1), 246 (2019).
- Bragg, D. C., et al. Disease onset in X-linked dystonia-parkinsonism correlates with expansion of a hexameric repeat within an SVA retrotransposon in TAF1. Proceedings of the National Academy of Sciences. 114 (51), 11020-11028 (2017).
- Stransky, N., Cerami, E., Schalm, S., Kim, J. L., Lengauer, C. The landscape of kinase fusions in cancer. Nature Communications. 5, 4846 (2014).
- Tenaillon, O., et al. The molecular diversity of adaptive convergence. Science. 335 (6067), 457-461 (2012).
- Vandecraen, J., Chandler, M., Aertsen, A., Van Houdt, R. The impact of insertion sequences on bacterial genome plasticity and adaptability. Critical Reviews in Microbiology. 43 (6), 709-730 (2017).
- Andersson, D. I., Gene Hughes, D. amplification and adaptive evolution in bacteria. Annual Review of Genetics. 43, 167-195 (2009).
- Bailey, S. F., Bataillon, T. Can the experimental evolution programme help us elucidate the genetic basis of adaptation in nature. Molecular Ecology. 25 (1), 203-218 (2016).
- Consuegra, J., et al. Insertion-sequence-mediated mutations both promote and constrain evolvability during a long-term experiment with bacteria. Nature Communications. 12 (1), 980 (2021).
- Burch, C. L., Romanchuk, A., Kelly, M., Wu, Y., Jones, C. D. Genome-wide determination of barriers to horizontal gene transfer. bioRxiv. , (2022).
- Slomka, S., et al. Experimental evolution of Bacillus subtilis reveals the evolutionary dynamics of horizontal gene transfer and suggests adaptive and neutral effects. Genética. 216 (2), 543-558 (2020).
- Choudhury, D., Saini, S. Evolution of Escherichia coli in different carbon environments for 2,000 generations. Journal of Evolutionary Biology. 32 (12), 1331-1341 (2019).
- Tenaillon, O., et al. Tempo and mode of genome evolution in a 50,000-generation experiment. Nature. 536 (7615), 165-170 (2016).
- Behringer, M. G., et al. Escherichiacoli cultures maintain stable subpopulation structure during long-term evolution. Proceedings of the National Academy of Sciences. 115 (20), 4642-4650 (2018).
- Voordeckers, K., et al. Adaptation to high ethanol reveals complex evolutionary pathways. PLoS Genetics. 11 (11), 1005635 (2015).
- Levy, S. F., et al. Quantitative evolutionary dynamics using high-resolution lineage tracking. Nature. 519 (7542), 181-186 (2015).
- Bruger, E. L., Marx, C. J. A decade of genome sequencing has revolutionized studies of experimental evolution. Current Opinion in Microbiology. 45, 149-155 (2018).
- Grubaugh, N. D., et al. An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar. Genome Biology. 20 (1), 8 (2019).
- Bedhomme, S., et al. Evolutionary changes after translational challenges imposed by horizontal gene transfer. Genome Biology and Evolution. 11 (3), 814-831 (2019).
- Tenaillon, O., Toupance, B., Le Nagard, H., Taddei, F., Godelle, B. Mutators, population size, adaptive landscape and the adaptation of asexual populations of bacteria. Genética. 152 (2), 485-493 (1999).
