กลับไปหน้าบทความ

อ่าน 4 นาที

การบีบอัดข้อมูลลำดับจีโนม

เปลี่ยนทางจากการเคลื่อนไหว

เทคโนโลยี การจัดลำดับจีโนมความเร็วสูงส่งผลให้ต้นทุนการจัดลำดับจีโนมลดลงอย่างมาก และทำให้ข้อมูลจีโนมเพิ่มขึ้นอย่างรวดเร็วอย่างน่าทึ่ง...

การบีบอัดข้อมูลลำดับจีโนม

เทคโนโลยี การจัดลำดับจีโนมความเร็วสูงส่งผลให้ต้นทุนการจัดลำดับจีโนมลดลงอย่างมาก และทำให้ข้อมูลจีโนมเพิ่มขึ้นอย่างรวดเร็วอย่างน่าทึ่ง เทคโนโลยีเหล่านี้ช่วยให้โครงการจัดลำดับจีโนมขนาดใหญ่ เช่นโครงการ 1000 Genomes Projectและโครงการ 1001 ( Arabidopsis thaliana ) Genomes Project สามารถดำเนินการได้ การจัดเก็บและการถ่ายโอนข้อมูลจีโนมจำนวนมหาศาลกลายเป็นปัญหาหลัก ซึ่งกระตุ้นให้เกิดการพัฒนาเครื่องมือบีบอัดข้อมูลประสิทธิภาพสูงที่ออกแบบมาโดยเฉพาะสำหรับข้อมูลจีโนม ความสนใจที่เพิ่มขึ้นอย่างมากในการพัฒนาอัลกอริทึมและเครื่องมือใหม่ๆ สำหรับการจัดเก็บและจัดการข้อมูลการจัดลำดับจีโนมซ้ำเน้นย้ำถึงความต้องการที่เพิ่มขึ้นสำหรับวิธีการที่มีประสิทธิภาพในการบีบอัดข้อมูลจีโนม

แนวคิดทั่วไป

แม้ว่าเครื่องมือบีบอัดข้อมูลมาตรฐาน (เช่น zip และ rar) จะถูกนำมาใช้ในการบีบอัดข้อมูลลำดับ (เช่นฐานข้อมูลไฟล์แบนGenBank ) แต่แนวทางนี้ถูกวิพากษ์วิจารณ์ว่าฟุ่มเฟือยเกินไป เนื่องจากลำดับจีโนมมักมีเนื้อหาที่ซ้ำกัน (เช่นลำดับไมโครแซทเทลไลต์ ) หรือลำดับจำนวนมากมีความคล้ายคลึงกันสูง (เช่น ลำดับจีโนมหลายลำดับจากสายพันธุ์เดียวกัน) นอกจากนี้ คุณสมบัติทางสถิติและทฤษฎีสารสนเทศของลำดับจีโนมอาจถูกนำมาใช้ประโยชน์ในการบีบอัดข้อมูลลำดับได้[ 1 ] [ 2 ] [ 3 ]

รูปที่ 1:ขั้นตอนหลักของเวิร์กโฟลว์สำหรับการบีบอัดข้อมูลการจัดลำดับจีโนมใหม่: (1) การประมวลผลข้อมูลการจัดลำดับดั้งเดิม (เช่น การลดชุดข้อมูลดั้งเดิมให้เหลือเพียงความแปรผันที่สัมพันธ์กับลำดับอ้างอิงที่ระบุ) (2) การเข้ารหัสข้อมูลที่ประมวลผลแล้วให้อยู่ในรูปแบบไบนารี และ (3) การถอดรหัสข้อมูลกลับเป็นรูปแบบข้อความ

ตัวแปรพื้นฐาน

เมื่อมีเทมเพลตอ้างอิงแล้ว จะต้องบันทึกเฉพาะความแตกต่าง (เช่น การแทนที่นิวคลีโอไทด์เดี่ยวและการแทรก/การลบ) เท่านั้น ซึ่งจะช่วยลดปริมาณข้อมูลที่ต้องจัดเก็บลงได้อย่างมาก แนวคิดเรื่องการบีบอัดเชิงสัมพัทธ์นั้นชัดเจน โดยเฉพาะอย่างยิ่งในโครงการจัดลำดับจีโนมใหม่ ซึ่งมีเป้าหมายเพื่อค้นหาความแปรผันในจีโนมแต่ละบุคคล การใช้แผนที่โพลีมอร์ฟิซึมของนิวคลีโอไทด์เดี่ยว ( SNP ) อ้างอิง เช่นdbSNPสามารถใช้เพื่อปรับปรุงจำนวนตัวแปรที่จะจัดเก็บให้ดียิ่งขึ้นได้[ 4 ]

พิกัดจีโนมสัมพัทธ์

แนวคิดที่มีประโยชน์อีกอย่างหนึ่งคือการจัดเก็บพิกัดจีโนมแบบสัมพัทธ์แทนพิกัดแบบสัมบูรณ์[ 4 ]ตัวอย่างเช่น การแสดงฐานตัวแปรลำดับในรูปแบบ ' Position1Base1Position2Base2… ', '123C125T130G' สามารถย่อให้เหลือ '0C2T5G' ได้ โดยที่จำนวนเต็มแสดงถึงช่วงระหว่างตัวแปร ค่าใช้จ่ายคือการคำนวณทางคณิตศาสตร์เล็กน้อยที่จำเป็นในการกู้คืนพิกัดแบบสัมบูรณ์ บวกกับการจัดเก็บปัจจัยการแก้ไข ('123' ในตัวอย่างนี้)

ข้อมูลเบื้องต้นเกี่ยวกับจีโนม

สามารถลดจำนวนลงได้อีกหากทราบตำแหน่งที่เป็นไปได้ทั้งหมดของการแทนที่ในกลุ่มลำดับจีโนมล่วงหน้า[ 4 ]ตัวอย่างเช่น หากทราบตำแหน่งทั้งหมดของ SNP ในประชากรมนุษย์แล้ว ก็ไม่จำเป็นต้องบันทึกข้อมูลพิกัดของตัวแปร (เช่น '123C125T130G' สามารถย่อเป็น 'CTG' ได้) อย่างไรก็ตาม แนวทางนี้ไม่ค่อยเหมาะสมนัก เนื่องจากข้อมูลดังกล่าวมักจะไม่สมบูรณ์หรือไม่พร้อมใช้งาน

การเข้ารหัสพิกัดจีโนม

มีการใช้รูปแบบ การเข้ารหัสเพื่อแปลงจำนวนเต็มพิกัดให้เป็นรูปแบบไบนารีเพื่อให้ได้ประโยชน์ในการบีบอัดเพิ่มเติม รูปแบบการเข้ารหัส เช่นรหัส Golombและรหัส Huffmanได้ถูกนำมาใช้ในเครื่องมือบีบอัดข้อมูลจีโนม[ 5 ] [ 6 ] [ 7 ] [ 8 ] [ 9 ] [ 10 ] แน่นอนว่า รูปแบบการเข้ารหัสจำเป็นต้องมี อั ลกอริทึมการถอดรหัสควบคู่ไปด้วย การเลือกรูปแบบการถอดรหัสอาจส่งผลต่อประสิทธิภาพของการดึงข้อมูล ลำดับ

ทางเลือกในการออกแบบอัลกอริทึม

วิธีการบีบอัดข้อมูลทางพันธุกรรมแบบสากลอาจไม่ใช่แนวทางที่ดีที่สุดเสมอไป เนื่องจากวิธีการเฉพาะอาจเหมาะสมกว่าสำหรับวัตถุประสงค์และเป้าหมายที่เฉพาะเจาะจง ดังนั้น การพิจารณาตัวเลือกการออกแบบหลายอย่างที่อาจส่งผลต่อประสิทธิภาพการบีบอัดจึงเป็นสิ่งสำคัญ

ลำดับอ้างอิง

การเลือกซีเควนซ์อ้างอิงสำหรับการบีบอัดแบบสัมพัทธ์อาจส่งผลต่อประสิทธิภาพการบีบอัด การเลือกใช้ซีเควนซ์อ้างอิงแบบฉันทามติแทนซีเควนซ์อ้างอิงที่เฉพาะเจาะจงกว่า (เช่น ซีเควนซ์อ้างอิงเคมบริดจ์ ฉบับปรับปรุง ) อาจส่งผลให้มีอัตราส่วนการบีบอัดที่สูงขึ้น เนื่องจากซีเควนซ์อ้างอิงแบบฉันทามติอาจมีอคติในข้อมูลน้อยกว่า[ 4 ]อย่างไรก็ตาม ความรู้เกี่ยวกับแหล่งที่มาของซีเควนซ์ที่กำลังถูกบีบอัดอาจถูกนำมาใช้เพื่อให้ได้ผลลัพธ์การบีบอัดที่ดียิ่งขึ้น แนวคิดในการใช้ซีเควนซ์อ้างอิงหลายซีเควนซ์ได้รับการเสนอ[ 4 ] Brandon et al. (2009) [ 4 ]ได้กล่าวถึงการใช้แม่แบบซีเควนซ์อ้างอิงเฉพาะกลุ่มชาติพันธุ์ โดยใช้การบีบอัด ข้อมูลตัวแปร ดีเอ็นเอไมโท คอน เดรียเป็นตัวอย่าง (ดูรูปที่ 2) ผู้เขียนพบ การกระจาย แฮพลอไท ป์ที่มีอคติ ใน ซีเควนซ์ ดีเอ็นเอไมโทคอน เดรียของชาวแอฟริกัน ชาวเอเชีย และชาวเอเชีย - ยุโรป เมื่อเทียบกับซีเควนซ์อ้างอิงเคมบริดจ์ ฉบับปรับปรุง ผลลัพธ์ของพวกเขาชี้ให้เห็นว่าลำดับอ้างอิงเคมบริดจ์ ที่แก้ไขแล้ว อาจไม่เหมาะสมเสมอไป เนื่องจากจำเป็นต้องจัดเก็บตัวแปรจำนวนมากขึ้นเมื่อใช้กับข้อมูลจากบุคคลที่มีเชื้อชาติห่างไกลกัน นอกจากนี้ ลำดับอ้างอิงสามารถออกแบบโดยอิงจากคุณสมบัติทางสถิติ[ 1 ] [ 4 ]หรือออกแบบ[ 11 ] [ 12 ]เพื่อปรับปรุงอัตราส่วนการบีบอัด

รูปแบบการเข้ารหัส

มีการสำรวจการประยุกต์ใช้รูปแบบการเข้ารหัสประเภทต่างๆ เพื่อเข้ารหัสเบสที่แปรผันและพิกัดจีโนม[ 4 ]รหัสคงที่ เช่นรหัส Golombและรหัส Riceเหมาะสมเมื่อการกระจายตัวของตัวแปรหรือพิกัด (แสดงเป็นจำนวนเต็ม ) ถูกกำหนดไว้อย่างชัดเจน รหัสแปรผัน เช่นรหัส Huffmanให้รูปแบบการเข้ารหัสเอนโทรปีทั่วไปมากขึ้นเมื่อการกระจายตัวของตัวแปรและ/หรือพิกัดพื้นฐานไม่ได้ถูกกำหนดไว้อย่างชัดเจน (โดยทั่วไปจะเป็นกรณีในข้อมูลลำดับจีโนม)

รายชื่อเครื่องมือบีบอัดข้อมูลการจัดลำดับจีโนมใหม่

อัตราส่วนการบีบอัดของเครื่องมือบีบอัดข้อมูลจีโนมที่มีอยู่ในปัจจุบันมีตั้งแต่ 65 เท่าถึง 1,200 เท่าสำหรับจีโนมมนุษย์[ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] [ 9 ] [ 10 ] [ 13 ]ตัวแปรหรือการแก้ไขที่ใกล้เคียงกันมากของจีโนมเดียวกันสามารถบีบอัดได้อย่างมีประสิทธิภาพมาก (ตัวอย่างเช่น มีรายงานอัตราส่วนการบีบอัด 18,133 [ 6 ]สำหรับการแก้ไขสองแบบของจีโนม A. thaliana เดียวกัน ซึ่งเหมือนกัน 99.999%) อย่างไรก็ตาม การบีบอัดดังกล่าวไม่ได้บ่งชี้ถึงอัตราส่วนการบีบอัดทั่วไปสำหรับจีโนมที่แตกต่างกัน (แต่ละบุคคล) ของสิ่งมีชีวิตเดียวกัน รูปแบบการเข้ารหัสที่พบบ่อยที่สุดในบรรดาเครื่องมือเหล่านี้คือการเข้ารหัส Huffmanซึ่งใช้สำหรับการบีบอัดข้อมูลแบบไม่สูญเสีย

เครื่องมือบีบอัดข้อมูลลำดับจีโนมที่เข้ากันได้กับรูปแบบไฟล์ลำดับจีโนมมาตรฐาน (BAM และ FASTQ)
ซอฟต์แวร์ คำอธิบาย อัตราส่วนการบีบอัด ข้อมูลที่ใช้ในการประเมิน แนวทาง/รูปแบบการเข้ารหัส ลิงก์ ใบอนุญาตการใช้งาน อ้างอิง
เพตาสวีท เครื่องมือบีบอัดแบบไม่สูญเสียข้อมูลสำหรับไฟล์ BAM และ FASTQ.gz; อ่านข้อมูลกลับแบบเรียลไทม์อย่างโปร่งใสผ่านไฟล์เสมือน BAM และ FASTQ.gz 60% ถึง 90% ลำดับจีโนมมนุษย์จากโครงการ 1000 Genomes Project https://petagene.comทางการค้า [ 14 ]
จีโนซิป โปรแกรมบีบอัดไฟล์ข้อมูลทางพันธุกรรมแบบอเนกประสงค์ – รองรับการบีบอัดไฟล์ FASTQ, SAM/BAM/CRAM, VCF/BCF, FASTA, GFF/GTF/GVF, PHYLIP, BED และไฟล์จาก 23andMe [ 15 ] [ 16 ]ลำดับจีโนมมนุษย์จากโครงการ 1000 Genomes Project เฟรมเวิร์กที่ขยายได้ของ Genozip http://genozip.comใช้เพื่อการค้าได้ แต่ใช้ได้ฟรีสำหรับการใช้งานที่ไม่ใช่เชิงพาณิชย์ [ 17 ]
การบีบอัดจีโนม (G-SQZ) เครื่องมือบีบอัดข้อมูลแบบไม่สูญเสียคุณภาพ ออกแบบมาเพื่อจัดเก็บและวิเคราะห์ข้อมูลลำดับนิวคลีโอไทด์ 65% ถึง 76% ลำดับจีโนมมนุษย์จากโครงการ 1000 Genomes Project การเข้ารหัสฮัฟฟ์แมน http://public.tgen.org/sqz-ไม่ได้ระบุ- [ 8 ]
CRAM (ส่วนหนึ่งของSAMtools ) การบีบอัดข้อมูลลำดับอ้างอิงที่มีประสิทธิภาพสูงและปรับแต่งได้ [ 18 ]หอจดหมายเหตุนิวคลีโอไทด์แห่งยุโรป ลดลมและrANShttp://www.ebi.ac.uk/ena/software/cram-toolkitอะปาเช่-2.0 [ 19 ]
ตัวบีบอัดจีโนม (GeCo) เครื่องมือที่ใช้โมเดล Markov หลายแบบผสมผสานกันเพื่อบีบอัดลำดับข้อมูลอ้างอิงและลำดับข้อมูลที่ไม่อ้างอิง ลำดับจีโนมนิวเคลียร์ของมนุษย์ การเข้ารหัสเลขคณิต http://bioinformatics.ua.pt/software/geco/หรือhttps://pratas.github.io/geco/จีพีแอลวี3 [ 13 ]
ตัวแปลงสัญญาณ GenomSys การบีบอัดแบบไม่สูญเสียข้อมูลของไฟล์ BAM และ FASTQ ในรูปแบบมาตรฐาน ISO/IEC 23092 [ 20 ] (MPEG-G) 60% ถึง 90% ลำดับจีโนมมนุษย์จากโครงการ 1000 Genomes Project การเข้ารหัสเลขคณิตไบนารีแบบปรับตามบริบท (CABAC) https://www.genomsys.comทางการค้า [ 21 ]
ฟาสตาฟการบีบอัดไฟล์ FASTA / UCSC2Bit ให้เป็นไฟล์เก็บถาวรแบบบีบอัดที่เข้าถึงได้แบบสุ่ม ชุดเครื่องมือสำหรับเชื่อมต่อไฟล์ FASTA, ดัชนี และไฟล์พจนานุกรมเสมือน ซึ่งช่วยให้สามารถผสานรวมระบบไฟล์ได้อย่างเรียบร้อย (คล้าย API) โดยไม่จำเป็นต้องคลายการบีบอัดไฟล์เก็บถาวรทั้งหมดเพื่อการเข้าถึงแบบสุ่ม/บางส่วนไฟล์ FASTAการเข้ารหัสแบบฮัฟฟ์แมนตามที่นำมาใช้โดย Zstdhttps://github.com/yhoogstrate/fastafsจีพีแอล-วี.2.0[ 22 ]
เครื่องมือบีบอัดข้อมูลลำดับจีโนมไม่รองรับรูปแบบไฟล์ลำดับจีโนมมาตรฐาน
ซอฟต์แวร์คำอธิบายอัตราส่วนการบีบอัดข้อมูลที่ใช้ในการประเมินแนวทาง/รูปแบบการเข้ารหัสลิงก์ใบอนุญาตการใช้งานอ้างอิง
ตัวบีบอัดความแตกต่างของจีโนม (GDC)เครื่องมือแบบ LZ77 สำหรับบีบอัดจีโนมหลายชุดของพืชชนิดเดียวกัน180 ถึง 250 เท่า / 70 ถึง 100 เท่าลำดับจีโนมนิวเคลียร์ของมนุษย์และยีสต์Saccharomyces cerevisiaeการเข้ารหัสฮัฟฟ์แมนhttp://sun.aei.polsl.pl/gdcจีพีแอลวี2[ 5 ]
การจัดลำดับจีโนมใหม่ (Genome Re-Sequencing หรือ GRS)เครื่องมือที่ใช้ลำดับอ้างอิงโดยไม่ขึ้นอยู่กับแผนที่ SNP อ้างอิงหรือข้อมูลความแปรผันของลำดับ159 เท่า / 18,133 เท่า / 82 เท่าลำดับจีโนมนิวเคลียร์ของมนุษย์, อาราบิโดปซิส ทาเลียนา (จีโนมเดียวกันแต่มีการแก้ไขหลายเวอร์ชัน) และโอไรซา ซาติวาการเข้ารหัสฮัฟฟ์แมนhttps://web.archive.org/web/20121209070434/http://gmdd.shgmo.org/Computational-Biology/GRS/ใช้งานได้ฟรีสำหรับการใช้งานที่ไม่ใช่เชิงพาณิชย์[ 6 ]
การเข้ารหัสลำดับจีโนมใหม่ (GReEN)เครื่องมือที่ใช้แบบจำลองการคัดลอกเชิงความน่าจะเป็นสำหรับการบีบอัดข้อมูลการจัดลำดับใหม่โดยใช้ลำดับอ้างอิงประมาณ 100 เท่าลำดับจีโนมนิวเคลียร์ของมนุษย์การเข้ารหัสเลขคณิตhttp://bioinformatics.ua.pt/software/green/-ไม่ได้ระบุ-[ 7 ]
ดีเอ็นเอซิปชุดเครื่องมือบีบอัดประมาณ 750 เท่าลำดับจีโนมนิวเคลียร์ของมนุษย์การเข้ารหัสฮัฟฟ์แมนhttp://www.ics.uci.edu/~dnazip/-ไม่ได้ระบุ-[ 9 ]
จีโนมซิปการบีบอัดข้อมูลโดยอ้างอิงจากจีโนมต้นแบบ สามารถใช้ฐานข้อมูลภายนอกเกี่ยวกับความแปรผันทางพันธุกรรมได้ (เช่น dbSNP)ประมาณ 1200 เท่าลำดับจีโนมนิวเคลียร์ของมนุษย์ (วัตสัน) และลำดับจากโครงการ 1000 จีโนมการเข้ารหัสเอนโทรปีสำหรับการประมาณค่าการแจกแจงเชิงประจักษ์https://sourceforge.net/projects/genomezip/-ไม่ได้ระบุ-[ 10 ]
ดึงข้อมูลมาจาก " https://en.wikipedia.org/w/index.php?title=Compression_of_genomic_sequencing_data&oldid=1334916733 "

สรุปเนื้อหา

ข้อมูลสำคัญจากบทความ

ข้อมูลสำคัญเกี่ยวกับ การบีบอัดข้อมูลลำดับจีโนม

เทคโนโลยี การจัดลำดับจีโนมความเร็วสูงส่งผลให้ต้นทุนการจัดลำดับจีโนมลดลงอย่างมาก และทำให้ข้อมูลจีโนมเพิ่มขึ้นอย่างรวดเร็วอย่างน่าทึ่ง...

แนวคิดทั่วไป

แม้ว่าเครื่องมือบีบอัดข้อมูลมาตรฐาน (เช่น zip และ rar) จะถูกนำมาใช้ในการบีบอัดข้อมูลลำดับ (เช่น ฐานข้อมูลไฟล์แบน GenBank ) แต่แนวทางนี้ถูกวิพากษ์วิจารณ์ว่าฟุ่มเฟือยเกินไป เนื่องจากลำดับจีโนมมักมีเนื้อหาที่ซ้ำกัน (เช่นลำดับไมโครแซทเทลไลต์ )...

ตัวแปรพื้นฐาน

เมื่อมีเทมเพลตอ้างอิงแล้ว จะต้องบันทึกเฉพาะความแตกต่าง (เช่น การแทนที่นิวคลีโอไทด์เดี่ยวและการแทรก/การลบ) เท่านั้น ซึ่งจะช่วยลดปริมาณข้อมูลที่ต้องจัดเก็บลงได้อย่างมาก แนวคิดเรื่องการบีบอัดเชิงสัมพัทธ์นั้นชัดเจน โดยเฉพาะอย่างยิ่งในโครงการจัดลำดับจีโนมใหม่...

พิกัดจีโนมสัมพัทธ์

แนวคิดที่มีประโยชน์อีกอย่างหนึ่งคือการจัดเก็บพิกัดจีโนมแบบสัมพัทธ์แทนพิกัดแบบสัมบูรณ์ [ 4 ] ตัวอย่างเช่น การแสดงฐานตัวแปรลำดับในรูปแบบ ' Position1Base1Position2Base2… ', '123C125T130G' สามารถย่อให้เหลือ '0C2T5G' ได้ โดยที่จำนวนเต็มแสดงถึงช่วงระหว่างตัวแปร...