Біоінформатичний аналіз послідовностей основних генів вірусу імунодефіциту великої рогатої худоби (gag, env, pol)

The sequences bioinformatic analysis results of the of the bovine immunodeficiency virus main genes (gag, env, pol), and updated data on the phylogenetic relationships of the pathogen are presented in the article. The existence of at least two genetic groups of BIV, which have a common origin from a common ancestor, has been established. In addition, the complexity of conducting a full-fledged and comprehensive sequence analysis for a complete understanding of the evolution and virus phylogenetic relationships was noted. It was determined that the lack of a single methodological approach to BIV sequencing in different world scientific communities complicates the analysis and interpretation of the obtained data. The application of whole -genome sequencing technologies is proposed to solve this issue

Keywords: phylogenetic relationships, Lentivirus, Retroviridae

Виклад основного матеріалу

Вірус імунодефіциту великої рогатої худоби (Bovine immunodeficiency virus, BIV), що належить до роду Lentivirus сімейства Retroviridae, спричиняє прогресуюче та стійке інфікування великої рогатої худоби [1]. Геном вірусу складається з двох копій одноланцюгової позитивної РНК, яка може інтегруватися в клітинну ДНК. За генетичною та антигенною структурою BIV споріднений з вірусом імунодефіциту людини (ВІЛ) та вірусом імунодефіциту мавп, тому його використовують як тваринну модель для вивчення інших ретровірусів [2]. BIV уперше був виділений у 1969 році від тільної молочної корови з ознаками виснаження, персистуючого лімфоцитозу, лімфаденопатії та ураженням центральної нервової системи за відсутності інфікування вірусом лейкозу ВРХ [3]. Вірус імунодефіциту ВРХ розповсюджений в усьому світі, а серопозитивну молочну та м'ясну худобу реєструють в багатьох країнах, включаючи США, Канаду, Іран, Німеччину, Францію, Швейцарію, Австралію, Велику Британію, Японію, Італію, Корею, Мексику, Індію та Бразилію [4].

У тварин, інфікованих BIV, не було виявлено ознак важкого прояву хвороби. Але, ураховуючи подібність цього вірусу до інших патогенних лентівірусів, є припущення, що BIV також здатен спричиняти клінічний прояв в організмі великої рогатої худоби, однак це не відбувається через високу швидкість оберту виробничих тварин та пов'язаний з цим короткий термін життя тварини. Тому, цілком можливо, що BIV-інфекція може бути пов'язана зі зниженням надою молока, втратою ваги в інфікованих тварин і сприяти вторинним інфекціям в їхньому організмі [5].

Геномна мінливість є важливою особливістю РНК-вірусів, зокрема і вірусу імунодефіциту ВРХ, який навіть в межах одного вірусного ізоляту демонструє значні геномні варіації [6].

В структурі BIV виявлено три основні рамки зчитування (ORF), що кодують гени gag, env і pol [1,4, 7].

Для молекулярно-генетичних досліджень найбільш придатними та інформативними є гени env і pol [8]. Ген env, який зазнає найбільшого селективного тиску у геномі ретровірусів, використовують для досліджень молекулярної епідеміології збудника. Ген pol, який є більш консервативним, використовують для створення філогенетичних дерев і встановлення еволюційних зв'язків ізолятів [4].

У доступній нам літературі є публікації щодо філогенетичного аналізу тих чи інших генів, проведеного науковцями з метою встановлення філогенетичних зв'язків збудника, що виділений в їхніх країнах. Ці дані демонструють більш ніж 98-99 % нуклеотидної гомології зі штамом R-29, що є референтним [8-11].

Мета досліджень

Провести біоінформатичний аналіз наявних у міжнародній базі даних GenBank нуклеотидних послідовностей основних генів BIV (gag, env, pol) та отримати оновлені дані щодо філогенетичних зв'язків збудника.

Матеріали та методи

Для формування баз даних нуклеотидних послідовностей основних генів вірусу імунодефіциту ВРХ проводили відповідний інформаційний пошук у міжнародній базі даних GenBank (www.ncbi.nml.nig.gov). Отримані послідовності розподіляли до відповідних підбаз з метою подальшого проведення філогенетичного аналізу.

Множинне вирівнювання обраних послідовностей проводили за допомогою програми BioEdit v. 7.0.0, модулю ClustalW [12].

Біоінформатичний аналіз отриманих послідовностей проводили з використанням алгоритму BLAST-online інтернет-ресурсу GenBank [13]. Для проведення філогенетичного аналізу використовували програму MEGA 11 [14]. Філогенетичний аналіз секвенованих послідовностей основних генів BIV та побудову дендрограм проводили за методом найближчих сусідів (Neighbour joining, bootstrap 1000, модель попарних дистанцій). Аналіз філогенетичного дерева здійснювали шляхом візуальної оцінки його топології та попарних відстаней між компонентами вибірки.

Результати досліджень

На основі інформації, представленої у міжнародній базі даних GenBank, було сформовано базу даних повних геномів та частково секвенованих послідовностей основних генів gag, env і pol ізолятів вірусу імунодефіциту ВРХ, виявленого у різних географічних регіонах світу. Послідовності були розподілені на окремі підбази, кожна з яких містила інформацію тільки про гени gag, env і pol. Усі послідовності були отримані у форматах FASTA (.fasta) або GenBank (*.gb).

Сформована база даних містила 52 послідовності геномних РНК BIV довжиною від 170 до 8 482 пар нуклеотидів (п. н.). З них: 2 повногеномні послідовності вірусу довжиною 8 482 п. н., 8 послідовностей гена gag довжиною від 171 до 393 п. н., 16 послідовностей гена env довжиною від 170 до 1 857 п. н., 22 послідовності гена pol довжиною від 176 до 451 п. н.

Після множинного вирівнювання послідовностей кожного з генів проводили пошук консервативних регіонів і варіабельних ділянок з метою їх аналізування.

Результати множинного вирівнювання нуклеотидних послідовностей генів gag та env демонструють відсутність єдиного методологічного підходу щодо секвенування цих генів, що унеможливлює їх аналізування за допомогою стандартних біоінформатичних пакетів.

Множинне врівнювання нуклеотидних послідовностей гена pol демонструвало відносну однорідність вибірки, що дозволило нам сформувати матрицю для аналізу в програмному пакеті MEGA 11. За результатами множинного вирівнювання було встановлено, що максимальне перекриття вирівняних ділянок становило 176 п. н.

Після множинного вирівнювання отриманих послідовностей були побудовані філогенетичні дерева за алгоритмом найближчого сусіда для кожного з генів BIV.

Аналіз філодендрограми, побудованої для гена gag (рис. 1), свідчить про те, що існує щонайменше дві генетичні гілки BIV, які походять від спільного пращура та у певний час відокремились одна від одної.

Рис. 1 Філогенетичний аналіз ізолятів BIV на основі послідовностей гена gag

За результатами філогенетичного аналізу гена env також встановлено існування двох генетичних гілок BIV (рис. 2).

Рис. 2 Філогенетичний аналіз ізолятів BIV на основі послідовностей гена env

Існування двох геногруп було також підтверджено під час побудови дендрограми для гена pol (рис. 3).

Рис. 3 Філогенетичний аналіз ізолятів BIV на основі послідовностей гена pol

Таким чином, біоінформатичний аналіз нуклеотидних послідовностей генів gag, env і pol BIV демонструє існування щонайменше двох генетичних груп BIV і вказує на їхнє спільне походження від єдиного пращура. Ці дані узгоджуються з даними, отриманими іншою науковою групою дослідників, які у своїх дослідженнях використали ген pol [4].

Однак, мала вибірка доступних для біоінформатичного аналізу послідовностей, що опубліковані в міжнародній базі даних GenBank, та відсутність єдиного методологічного підходу щодо секвенування цих генів не дозволяє провести повноцінний і всебічний аналіз та оцінити потенційний ризик впливу цього вірусного агента. Крім того, еволюційна мінливість цього РНК-вірусу, його генетична й антигенна спорідненість з вірусом імунодефіциту людини робить його потенційно небезпечним у концепті Єдиного здоров'я. Що, враховуючи негативний досвід людства щодо інших РНК-вірусів (пташиний грип, COVID-19, ВІЛ та ін.), вимагає більш детального вивчення та ретельного аналізування усіх можливих ризиків, пов'язаних з цим.

Висновки

1. Біоінформатичний аналіз нуклеотидних послідовностей генів gag, env і pol BIV демонструє існування щонайменше двох генетичних груп вірусу та вказує на їхнє спільне походження від єдиного пращура.

2. Молекулярно-генетичний аналіз BIV значно ускладнений відсутністю єдиного методологічного підходу щодо секвенування цього збудника у різних лабораторіях світу.

3. Вирішенням цього питання може стати застосування технологій повногеномного секвенування, що дасть більш чітке та конкретизоване розуміння філогенетичних зв'язків збудника.

Перспективи використання отриманих результатів. Отримані результати можуть стати підґрунтям для проведення більш детального вивчення циркуляції BIV в Україні, а результати біоінформатичного аналізу взято за основу для проведення філогенетичного аналізу українських ізолятів після їх секвенування.

Список літератури