Teknik mühendislik makalesi incelemesi

DOI: 10.1186/s12859-023-05385-y

PacBio Yüksek Doğruluklu Okumalardan Mitokondriyal Genom Montajı İçin Python Boru Hattı: MitoHiFi

Biyoçeşitlilik genomik araştırmalarında, PacBio HiFi gibi yüksek doğruluklu uzun okuma teknolojileri karmaşık genomların montajında devrim yaratmıştır. Ökaryotik hücrelerde mitokondriyal genomlar nükleer genomla birlikte yüksek kapsama oranlarında dizilenmektedir.

Yazar
Uliano-Silva M et al.
Yayın
Tarih
2023-07-18
01

Makalenin Kapsamı ve Akademik İçeriği

Küresel biyoçeşitlilik genomik projelerinde yüksek kaliteli referans genomların üretilmesi temel bir gaye haline gelmiştir. Ökaryotik organizmalarda hücresel solunum ve enerji metabolizmasının anahtarı olan mitokondri genomları, maternal kalıtım ve haploid yapıları gereği evrimsel araştırmalarda sıklıkla kullanılmaktadır. Bununla birlikte, uzun okumalı yeni nesil dizileme teknolojileri ortaya çıkmadan önce kullanılan kısa okumalı dizileme veya PCR tabanlı yöntemler, mitokondriyal genomlarda yer alan geniş yinelemeli dizileri ve nükleer genoma entegre olmuş mitokondriyal kopyaları (NUMT) tam anlamıyla ayırt edememekteydi. PacBio HiFi okumalarının sunduğu yüksek doğruluk ve uzun okuma boyutu organel genom dizilemesinde yeni bir dönem başlatmış, fakat bu özel veri türünü doğrudan ve yüksek ölçekte işleyen uzmanlaşmış boru hatlarının eksikliği hissedilmiştir.

Bu sorunu çözmek amacıyla geliştirilen MitoHiFi yazılımı, PacBio HiFi okumalarını veya pre-monte edilmiş kontig dizilerini ham veri olarak kabul eden Python tabanlı açık kaynaklı bir biyoenformatik araçtır. Boru hattı, yakın akraba bir türe ait referans mitokondri genomunu kullanarak HiFi okumalarını hizalamakta, referans genom uzunluğundan belirgin şekilde uzun olan okumaları NUMT adayı olarak süzmektedir. Ardından seçilen kaliteli okumalar Hifiasm yazılımı ile monte edilmekte, BLAST+ ve circularizationCheck modülleri ile dairesellik, kapsama oranı ve uzunluk kriterlerine göre doğrulanmaktadır. Kodlayıcı gen ve RNA anetasyonları MitoFinder veya MITOS entegrasyonuyla gerçekleştirilmekte, tRNA-Phe lokusu başlangıç noktası kabul edilerek standart oryantasyona dönüştürülmektedir. Böylece süreç uçtan uca otomatize edilmiştir.

Yazılımın deneysel uygulaması ve doğrulaması kapsamında Darwin Tree of Life (DToL), Vertebrate Genomes Project (VGP) ve Aquatic Symbiosis Genomics (ASG) projeleri bünyesinde 368 Metazoa ve 6 Fungi türünü kapsayan toplam 374 mitokondriyal genom başarıyla monte edilmiştir. Daha önce kamuya açık veritabanlarında referans dizisi bulunan 60 tür üzerinde yapılan karşılaştırmalı analizlerde, MitoHiFi ile elde edilen genomların 16 türde önceki referanslara kıyasla daha uzun olduğu saptanmıştır. Detaylı nokta grafik (dotplot) ve dizi hizalama analizleri, bu uzunluk farkının önceki kısa okumalı montajlarda çöken tandem yineleme dizilerinin HiFi okumaları sayesinde eksiksiz çözümlenmesinden kaynaklandığını ortaya koymuştur.

Sonuç olarak MitoHiFi, yüksek doğruluklu uzun okuma verilerinden organel genomu montajında yüksek ölçeklenebilirlik, otomasyon ve hassasiyet sağlayan güçlü bir boru hattı olarak biyoenformatik literatürüne kazandırılmıştır. Çeşitli taksonomik gruplarda test edilen yazılım, metazoan türlerinde mükemmel başarım gösterirken bitki organel genomlarındaki karmaşık izomer yapıları ve aşırı uzun yinelemeler için ek grafik tabanlı yaklaşımlara ihtiyaç duyulduğunu da göstermektedir. Bu çalışma, organel genomiği, taksonomi ve popülasyon genetiği alanında çalışan araştırmacılar için güvenilir, tekrarlanabilir ve Docker/Singularity konteyner desteği sunan temel bir biyoenformatik kaynağı niteliğindedir.

Araştırmanın amacı

Bu çalışmanın amacı, PacBio yüksek doğruluklu (HiFi) dizi okumalarını kullanarak ökaryotik canlıların mitokondriyal genomlarını otomatik, yüksek doğruluklu ve ölçeklenebilir biçimde monte eden ve genelleştiren MitoHiFi adlı açık kaynaklı Python yazılım hattını geliştirmek ve geniş biyolojik taksonlarda geçerlemektir(M et al., 2023).

İnceleme ve modelleme yöntemi

MitoHiFi yazılımı, PacBio HiFi ham dizi okumalarını veya önyüklemeli kontigleri girdi olarak alan, yakın akraba mitokondriyal referans genom eşleşmesiyle NUMT dizilerini filtreleyen, Hifiasm ile montaj ve MitoFinder/MITOS ile gen anotasyonu gerçekleştiren modüler bir Python boru hattı yaklaşımı sunmaktadır(M et al., 2023, s. 3).

Kuramsal ve kavramsal çerçeve

Çalışmanın teorik çerçevesi, biyoenformatik algoritma tasarımı, dairesel konsensus sekanslama (PacBio HiFi) ilkeleri, organel genomik anatomisi, dizileme hataları ile heteroplazmik varyasyonların ayırt edilmesi ve dizi hizalama tabanlı filtreleme yöntemlerine dayanmaktadır.

Makalenin ele aldığı literatür

Mevcut organel montaj araçlarının ağırlıklı olarak kısa okumalı verilere göre tasarlanması nedeniyle organel genomlarındaki geniş yineleme yapılarını ve heteroplazmiyi çözmede yetersiz kalması, nükleer aktarımları (NUMT) gerçek mitokondri dizilerinden ayıramaması ve yüksek ölçekli biyoçeşitlilik projelerinde otomasyon eksikliği literatürdeki temel sorunlar olarak ele alınmaktadır.

Araştırma Tasarımı

Evren veya inceleme alanı
Darwin Tree of Life (DToL), Vertebrate Genomes Project (VGP) ve Aquatic Symbiosis Genomics (ASG) projeleri kapsamında dizilenen ökaryotik canlı türlerinin biyolojik veri havuzu.
Örneklem veya araştırma materyali
Metazoa şubesinden 368 tür ve Fungi aleminden 6 tür olmak üzere toplam 374 ökaryotik türe ait PacBio HiFi yüksek doğruluklu uzun okuma verileri.
Veri toplama süreci
Hedef canlı türlerinden elde edilen nükleer ve organel DNA materyallerinin PacBio dairesel konsensüs (HiFi) dizileme teknolojisi ile yaklaşık 25 kat kapsama derinliğinde dizilenmesi.
Veri analiz yöntemi
Hifiasm montaj mimarisi, BLAST+ hizalama algoritması, circularizationCheck dairesellik doğrulama modülü ve MitoFinder/MITOS anotasyon araçlarının Python3 diliyle entegre otomasyonu.
Mitokondriyal genom montajıPacBio HiFi okumalarıNUMT filtrelemeHeteroplazmiYinelemeli dizilerGen anotasyonu

Araştırmanın Bulguları

01

MitoHiFi yazılımı, DToL, VGP ve ASG projeleri çerçevesinde 39 Metazoa takımı ile 6 Fungi türünü kapsayan toplam 374 ökaryotik türün mitokondriyal genomunu yüksek doğrulukla başarıyla monte etmiştir(M et al., 2023, s. 7).

02

Kamuya açık INSDC veri tabanlarında daha önce yayımlanmış referans genomu bulunan 60 türün incelemesinde, MitoHiFi tarafından oluşturulan 16 genomun önceki dizilimlerden daha uzun olduğu ve bu uzunluk farkının önceden saptanamayan yinelemeli dizilerin tam çözümlenmesinden kaynaklandığı tespit edilmiştir(M et al., 2023, s. 8).

03

Yazılım hattı, dairesel genom kopyalarındaki tekrarları belirleyip dizinin başlangıç noktasını tRNA-Phe geni konumuna göre otomatik olarak döndürerek biyoinformatik standartlara uygun çıktılar üretmektedir(M et al., 2023, s. 3).

04

MitoHiFi, nükleer mitokondriyal transferlerden (NUMT) türetilmiş okumaları uzunluk eşikleri ve hizalama filtreleri sayesinde etkili biçimde eleyerek organel genomlarının yanlış dizilenmesini ve birleştirilmesini engellemektedir(M et al., 2023, s. 9).

02

Bulguların Tartışılması ve Araştırmanın Sonucu

Bulguların Yazarlar Tarafından Yorumlanışı

MitoHiFi yazılımının yüksek doğruluklu PacBio HiFi verilerini kullanarak 374 farklı ökaryotik türde eksiksiz mitokondriyal genomlar oluşturduğu ve kamu veritabanlarında mevcut 60 türün 16'sında önceden saptanamamış yinelemeli dizileri ortaya çıkardığı gösterilmiştir.

Yazarlar, daha önce dizilenmiş genomlarda eksik olan bu yinelemeli dizilerin varlığını, geleneksel kısa okumalı dizileme veya PCR tabanlı yaklaşımların tandem tekrarları birleştirememesi ve nükleer mitokondriyal transferleri (NUMT) ayırt edememesi olgusuyla açıklamaktadır.

Mevcut organel montaj yazılımlarının kısa okuma verilerinin sınırlılıklarına göre tasarlandığı vurgulanarak, MitoHiFi'ın HiFi okumalarının sunduğu uzunluk ve yüksek doğruluk avantajıyla Formenti ve arkadaşlarının işaret ettiği karmaşık yinelemeli mitokondri yapılarını çözmede belirgin bir üstünlük sağladığı literatür bağlamında ifade edilmektedir.

Sonuç olarak araç, biyoçeşitlilik genom projelerinde organel genomlarının ölçeklenebilir, otomatik ve yüksek doğrulukla monte edilmesini olanaklı kılmakta, kamu veri tabanlarında yer alan hatalı veya eksik referans mitokondri genomlarının güncellenmesine doğrudan zemin hazırlamaktadır.

Bulguların Önceki Araştırmalarla Karşılaştırılması

Önceki kısa okumalı veya PCR tabanlı organel montaj araçlarının tandem tekrarları ve NUMT yapılarını çözmedeki yetersizliklerine karşın, MitoHiFi uygulamasının HiFi uzun okuma verisini işlemedeki üstünlüğü vurgulanmış ve omurgalı organel genomlarındaki yinelemeli dizilerin doğru montajı doğrulanmıştır(M et al., 2023, s. 9).

Sonuç Bölümünün Sunduğu Çıkarımlar

MitoHiFi, PacBio HiFi verilerini kullanarak geniş bir filogenetik yelpazede yer alan ökaryotik türlerin mitokondriyal genomlarını otomatik, hızlı ve güvenilir biçimde monte eden ve gen anotasyonunu gerçekleştiren açık kaynaklı bir Python aracıdır. Büyük biyoçeşitlilik projelerinde 374 tür üzerinde başarıyla uygulanmış olup, biyoinformatik analizlerde standart bir organel montaj hattı sağlama potansiyeline sahiptir(M et al., 2023, s. 10).

Literatür Taraması İçin Sunduğu Çerçeve

Bu çalışma, biyoçeşitlilik genomiğinde organel genomlarının montajında karşılaşılan teknik ve biyolojik zorluklara yüksek doğruluklu uzun okuma (PacBio HiFi) verilerini temel alan otomatik bir çözüm sunarak literatürde işlev görür. Ökaryotik genom projelerinde mitokondriyal diziler genellikle nükleer DNA ile birlikte dizilenmekte, ancak geleneksel biyoenformatik araçlar nükleer mitokondriyal transferleri (NUMT) ve tandem tekrarları tam çözememektedir. Makale, geliştirilen MitoHiFi boru hattı sayesinde organel veri setlerinin güvenilir şekilde ayrıştırılması ve montajı için uçtan uca biyoinformatik standart getirmektir.

Çalışma, GetOrganelle ve MitoFinder gibi daha önce kısa okumalı dizileme verileri veya hedefli zenginleştirme yaklaşımları için geliştirilmiş organel montaj yazılımlarının sınırlılıklarını aşma amacı taşır. Kısa okumalara dayalı araçların tandem tekrarları birleştirme yetersizliği ve NUMT dizilerini gerçek organel genomundan ayırt edememe sorunları, Formenti ve arkadaşlarının ortaya koyduğu karmaşık yineleme yapılarıyla belirginleşmiştir. MitoHiFi, HiFi okumalarının sunduğu yüksek doğruluk (Q20 üzeri) ve uzun parça boyutlarını (15-20 kb) Hifiasm algoritmik altyapısıyla harmanlayarak önceki literatürü yöntemsel açıdan ileri taşır.

Literatür taraması yapacak araştırmacılar için bu çalışma, organel genomiğinde veri kalitesi ve montaj doğruluğu standartlarının yeniden tanımlanmasına katkı sağlar. Kamu veri tabanlarında bulunan referans mitokondri genomlarının önemli bir kısmının eksik tandem tekrarlar içerebileceğini ve heteroplazmi ile NUMT karmaşasından etkilenebileceğini gösterir. Araştırmacılar, sistematik taramalarda mevcut referans dizilerin temsil yeteneğini sörgülarken bu makalenin sunduğu karşılaştırmalı bulguları ve algoritmik filtreleme kriterlerini teorik zemin olarak kullanabilir.

Kaynakça

  1. Formenti, G., Rhie, A., Balacco, J., Haase, B., Mountcastle, J., Fedrigo O., & Jarvis, E. D. (2021). Complete vertebrate mitogenomes reveal widespread repeats and gene duplications. Genome Biology, 22, 120.
03

Tez ve Makale Araştırmalarında Kullanım

Makalenin doğrulanmış yöntemi, bulguları ve kuramsal çerçevesi; tez ve makale çalışmalarında kullanılacağı bölümle birlikte aşağıdaki üç araştırma taslağına dönüştürülmüştür.

01

Tez veya makalede kullanım · method

PacBio HiFi okumaları kullanılarak yürütülen organel genom montajlarında nükleer mitokondriyal transferler (NUMT) filtrelenebilir ve dairesel genomlar tRNA-Phe konumuna göre otomatik döndürülebilir.

Araştırmacılar biyoinformatik yöntem bölümünde MitoHiFi hattının filtreleme parametrelerini (-p, --max-read-len) ve tRNA tabanlı döndürme algoritmasını referans göstererek kendi dizileme veri işleme protokollerini gerekçelendirebilir(M et al., 2023).

Böcek Takımlarında Mitokondriyal Genom Yinelemelerinin HiFi Okumaları ile Karşılaştırmalı Analizi

Araştırma sorusu
Farklı böcek takımlarına ait türlerde mitokondriyal tandem yineleme oranları kısa okumalı referans genomlara kıyasla HiFi tabanlı montajlarda ne oranda artış göstermektedir?
Yöntem taslağı
Kamu veri tabanlarında kısa okuma referansı bulunan 30 farklı böcek türünün PacBio HiFi ham veri setleri seçilir. Örneklem büyüklüğü filogenetik temsiliyet esasına göre belirlenir. MitoHiFi hattı çalıştırılarak yeni montajlar elde edilir ve FastANI ile Dotter araçları kullanılarak eski referanslarla dizi uzunluğu ve yineleme kopya sayıları yönünden karşılaştırılır.

Makaleyle bağlantısı: Makalede 16 türde tespit edilen yeni tandem yineleme kopyalarının böcek takımlarındaki yaygınlığını ve evrimsel dinamiklerini sistematik biçimde sınama ihtiyacından türetilmiştir(M et al., 2023, s. 8).

02

Tez veya makalede kullanım · literatüre-review

Kısa okumalı dizileme yöntemleriyle kamu veri tabanlarına aktarılan geleneksel mitokondri referans dizileri tandem yineleme bölgelerini eksik yansıtabilmektedir.

Literatür taraması bölümünde kamu veritabanlarındaki referans mitokondri genomlarının doğruluk ve uzunluk sınırlılıkları tartışılırken, MitoHiFi'ın 60 türlük karşılaştırmalı analizi ampirik kanıt olarak sunulabilir(M et al., 2023, s. 8).

Deniz Omurgasızlarında NUMT Kontaminasyonunun Filtreleme Algoritmalarıyla Engellenmesi

Araştırma sorusu
MitoHiFi yazılımının maksimum okuma uzunluğu filtre parametreleri deniz mollusk dizileme verilerinde NUMT yanlış montajlarını ne oranda engellemektedir?
Yöntem taslağı
Mollusca şubesinden seçilen 15 türün PacBio HiFi verileri işlenir. Farklı -p (dizi hizalama oranı) ve --max-read-len eşik değerleri uygulanarak elde edilen contig'ler çekirdek genom dizileriyle BLAST kullanılarak karşılaştırılır. Nükleer aktarımların organel montajına sızma oranları nicel olarak hesaplanır.

Makaleyle bağlantısı: Yazarların tartışma bölümünde Mollusca gibi değişken mitokondri yapısına sahip takımlarda %50 hizalama eşiğinin önemini vurgulamasından yola çıkılmıştır(M et al., 2023, s. 8).

03

Tez veya makalede kullanım · Tartışma

Kısa okuma teknolojilerine dayalı eski mitokondri montaj araçları nükleer DNA'ya aktarılmış organel dizilerini (NUMT) ayırt etmekte yetersiz kalmaktadır.

Tartışma bölümünde elde edilen genom montajlarının kalitesi değerlendirilirken eski araçların sınırlılıkları ve HiFi okumalarının sunduğu üstünlükler bu bulguyla desteklenebilir(M et al., 2023, s. 8).

Makrofungal Türlerde Mitokondriyal Genom Boyut Değişkenliği ve Intron Yapılarının Değerlendirilmesi

Araştırma sorusu
Fungi alemindeki mitokondriyal genom boyut farklılıkları MitoHiFi ile birleştirilen yeni dizilerde hangi genetik elemanların varlığıyla açıklanmaktadır?
Yöntem taslağı
Basidiomycota ve Ascomycota şubelerinden 20 mantar türünün PacBio HiFi verileri MitoHiFi -a fungi parametresi ile monte edilir. Elde edilen genom boyutları ve gen modelleri MITOS ve MitoFinder çıktıları üzerinden karşılaştırmalı intron-ekzon analizi ile incelenir.

Makaleyle bağlantısı: Makalenin mantar türlerinde 43 kb ile 133 kb arasında değişen mitokondri genomlarının başarıyla monte edildiğini gösteren ampirik bulgusuna dayanmaktadır(M et al., 2023, s. 10).