CC BY
52
Вестник Томского государственного университета. Биология
2010 № 3 (11)
ЦИТОЛОГИЯ И ГЕНЕТИКА
УДК 575:582.475
Г.Г. Гончаренко1, А.А. Сурков1, В.Г. Митрофанов2, Л.И. Корочкин2
1 Гомельский государственный университет им. Ф. Скорины (г. Гомель, Республика Беларусь)
2Институт биологии развития им. Н.К. Кольцова РАН (г. Москва)
УРОВЕНЬ ГЕННОГО ПОТОКА У DROSOPHILA LITTORALIS (DIPTERA: DROSOPHILIDAE) В ПРИРОДНЫХ ПОПУЛЯЦИЯХ ВОСТОЧНОЙ ЕВРАЗИИ
Исследование выполнено при финансовой поддержке Национальной академии наук Беларуси (грант № 2003865) и Министерства образования Республики Беларусь (грант № 2005562).
На основании проведенного генетического анализа 15 природных популяций Drosophila littoralis Meigеn Восточной Евразии с использованием 14 генов, кодирующих изоферменты, были установлены основные показатели генетической подразделен-ности и генного потока. Величина генного потока (Nm) для всех исследованных нами популяций Восточной Евразии составила 4,21. Это говорит о том, что популяции D. littoralis обмениваются генетическим материалом в среднем с интенсивностью 4,2 мигранта за поколение. Для отдельных групп популяций эта величина возрастает до 5,43. Полученные генетические данные указывают на связь величины генного потока с географической удаленностью популяций друг от друга.
Ключевые слова: Drosophila littoralis; изоферменты; генетическая структура популяций; генный поток.
Работа выполнялась в рамках научного сотрудничества между Институтом биологии развития Н.К. Кольцова РАН (г. Москва, Россия) и Гомельским государственным университетом им. Ф. Скорины (г. Гомель, Республика Беларусь).
Представители двукрылых видов-двойников Drosophila группы virilis успешно использовались в качестве модельной системы для изучения процессов видообразования [1—3], генетики развития [4-5], молекулярной эволюции [6-8], а также таксономии и систематики [9-10]. Основным микроэволюци-онным фактором, сглаживающим действие естественного отбора, дрейфа генов и мутационного процесса, лежащих в основе генетической дифференциации природных популяций, является генный поток. Именно генный поток путем обмена наследственным материалом между популяциями выравнивает их генетическую структуру, позволяя виду сохранять единый генофонд. Величина генного потока зависит от сложного взаимодействия различных мик-роэволюционных сил и может серьезно различаться в изолированных или непрерывных популяциях одного вида.
В последние десятилетия с появлением генетических маркеров у исследователей впервые возникла возможность точно оценивать величину генного потока. Однако исследования, направленные на его измерение у высших насекомых в Восточной Евразии непосредственно в природных популяциях, практически отсутствуют. Удобной моделью для оценки генного потока в природных популяциях насекомых является хорошо генетически изученная Drosophila littoralis Ме1§еп - восточно-евразийский представитель двойниковых видов, входящих в группу virilis.
Целью нашей работы была оценка уровня генного потока в 15 природных популяциях у вида D. littoralis в Восточной Евразии на основе использования в качестве молекулярно-генетических маркеров 14 генов, кодирующих изоферменты.
Материалы и методы исследования
D. littoralis обитает вблизи незагрязненных лесных водоемов. Её взрослые особи были отловлены в 15 природных популяциях (рис. 1), которые находятся в восточной части ареала распространения данного вида на территории Восточной Евразии [8, 9, 11].
Название популяций, места обитания и год отлова исследованных особей D. littoralis: 1) «Мукачево» - вблизи г. Мукачево, Украина, 1986; 2) «Латорица» -р. Латорица в 10 км восточнее Мукачево, Украина, 1985; 3) «Воловец» - г. Воло-вец, Украина, 1986; 4) «Черновцы» - вблизи г. Черновцы, Украина, 1985; 5) «Днепровская» - вблизи г. Речица, Белоруссия, 1985, 2003-2005; 6) «Гомель ручейная» - ручей в 10 км южнее г. Гомель, Белоруссия, 1984-1987, 2003-2005; 7) «Гомель болотная» - Галые болота восточнее г. Гомель, Белоруссия, 1981— 1984, 2003-2005; 8) «Орша» - вблизи г. Орша, Белоруссия, 1985; 9) «Латвия» -вблизи г. Цесие, Латвия, 1989; 10) «Кропотово» - вблизи п. Кропотово, Московская обл., Россия, 1980; 11) «Финляндия» - вблизи г. Оулу, Финляндия, 1980;
12) «Карелия» - Карельское побережье Белого моря, Карелия, Россия, 1985;
13) «Новосибирск» - р. Обь в 100 км южнее г. Новосибирск, Новосибирская обл., Россия, 1987-1989; 14) «Алтай» - р. Катунь в 80 км южнее г. Бийск, Алтайский край, Россия, 1989; 15) «Талас» - вблизи г. Талас, Кыргызстан, 1986-1987. Местоположение проанализированных популяций показано на рис. 1.
Изученные популяции были разделены на группы: европейская группа -популяции 1-12, сибирская - 13-14 и тянь-шаньская - 15. В свою очередь 12 чисто европейских популяций подразделяются на южно-европей-ские - 1-4 (Украина), восточно-европейские - 5-10 (Белоруссия, Россия, Латвия) и североевропейские - 11-12 (Финляндия, Карелия).
Взрослые особи вида D. littoralis исследовались методом электрофореза. Каждая особь гомогенизировалась в 25 мкл дистиллированной воды или гелевого буфера. Электрофоретическое фракционирование гомогенизированных экстрактов индивидуальных особей проводилось нами по 11 ферментам в 13-14%-ном крахмальном геле с использованием двух буферных систем: А) трис-ЭДТА-боратная, рН 8,6; В) трис-цитрат, рН 6,2. Все параметры электрофоретического фракционирования, а также методики экстракции и гисто-
химического выявления ферментов подробно приведены нами ранее [1, 12,
13]. Обозначение выявленных электрофоретических вариантов дано по общепринятой номенклатуре Пракаша с соавт. [14].
Для выявления значения параметров генетической подразделённости популяций использовали показатели Бзт [15] и Оэт [16]. Величина генного потока (Кет) определялась из соотношения Кеш = (1-Рзт)/4Рзт, где Рзт - коэффициент подразделенности популяций [15].
Рис. 1. Места обитания проанализированных популяций и распространение Б. ШотаИз на территории Восточной Евразии:
«•» - популяции, «—» - граница распространения
Результаты исследования и обсуждение
В ходе электрофоретического исследования особей Б. ¡ШогаНз из 15 природных популяций по 11 ферментным системам удалось выявить 41 различный электрофоретический вариант, находящийся, как было показано ранее [17-18], под генетическим контролем 14 локусов.
Для оценки генетической структуры были рассчитаны частоты встречаемости аллелей в каждой из 15 исследованных популяций Б. ¡ШогаЫз. Аллельные частоты по 14 проанализированным генам представлены в табл. 1.
Как видно из табл. 1, полностью мономорфными в исследованных популяциях являются четыре локуса: Бит, а-ОрЛ, с-МЛ, т-МЛ, поскольку по каждому из них найден только один аллель. Наибольшая изменчивость обнаружена по генам, кодирующим а-эстеразу-3, В-эстеразу-2, кислую фосфатазу-1 и октанолдегидрогеназу (табл. 1).
Т а б л и ц а 1
Аллельные частоты в исследованных природных популяциях Б. ШХотаНь
Локус, Популяции
аллели Мук Лт Вол Чер Дн Г-р Г-б Орш Лат Кр Фин Кар Нов Алт Тал
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
РОМ
п
0,40
0,80
1,00
1,20
101 18 10 36 62 106 156 4 8 4 3 10 4 12 36
0,029 0,000 0,000 0,027 0,032 0,009 0,006 0,000 0,000 0,250 0,000 0,000 0,000 0,000 0,027
0,971 1,000 1,000 0,918 0,952 0,973 0,969 0,750 1,000 0,750 1,000 1,000 0,750 0,750 0,973
0,000 0,000 0,000 0,055 0,016 0,009 0,019 0,250 0,000 0,000 0,000 0,000 0,250 0,167 0,000
0,000 0,000 0,000 0,000 0,000 0,009 0,006 0,000 0,000 0,000 0,000 0,000 0,000 0,083 0,000
МЕ
п
1,10
1,20
136 18 10 36 62 106 156 4 8 4 3 10 4 12 36
0,000 0,000 0,000 0,000 0,000 0,000 0,006 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
1,000 1,000 1,000 1,000 1,000 1,000 0,994 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
НК-1
п
1,00
1,40
1,80
2,00
90 18 10 36 62 106 136 4 8 4 3 10 4 12 36
0,011 0,055 0,000 0,000 0,032 0,000 0,014 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
0,989 0,945 1,000 1,000 0,919 0,982 0,979 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
0,000 0,000 0,000 0,000 0,032 0,009 0,007 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
0,000 0,000 0,000 0,000 0,016 0,009 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
НК-8
п
1,00
1,05
90 18 10 36 62 106 136 4 8 4 3 10 4 12 36
1,000 1,000 1,000 1,000 0,968 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
0,000 0,000 0,000 0,000 0,032 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
ким
п
1,00
90 20 10
1,000 1,000 1,000
36
1,000
62
1,000
106
1,000
144 4
1,000 1,000
0
0,000
4 3 10 0 0 36
1,000 1,000 1,000 0,000 0,000 1,000
а-Е8Т-3
п
0,90
0,95
1,00
1,02
1,10
1,14
132 25 10 36 68 116 149 4 0 4 4 10 0 0 36
0,000 0,000 0,000 0,000 0,000 0,008 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
0,190 0,200 0,100 0,194 0,132 0,155 0,255 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,083
0,023 0,040 0,100 0,166 0,059 0,086 0,042 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
0,257 0,320 0,400 0,306 0,324 0,302 0,315 0,250 0,000 0,250 0,250 0,700 0,000 0,000 0,734
0,333 0,200 0,300 0,223 0,324 0,379 0,248 0,500 0,000 0,250 0,750 0,200 0,000 0,000 0,128
0,197 0,240 0,100 0,111 0,162 0,069 0,140 0,250 0,000 0,500 0,000 0,100 0,000 0,000 0,055
Р-Е8Т-2
п 134 22 10 35 63 112 154 4 8 4 3 10 4 10 36
1,36 0,000 0,000 0,000 0,000 0,000 0,000 0,007 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
1,39 0,000 0,000 0,000 0,000 0,016 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
1,44 0,045 0,000 0,200 0,000 0,111 0,125 0,104 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
1,46 0,023 0,000 0,000 0,000 0,016 0,000 0,032 0,000 0,000 0,500 0,000 0,200 0,000 0,000 0,000
1,48 0,701 0,454 0,600 0,886 0,587 0,572 0,565 0,500 0,750 0,500 1,000 0,700 0,500 0,600 0,083
1,51 0,231 0,545 0,200 0,114 0,270 0,303 0,292 0,500 0,250 0,000 0,000 0,100 0,500 0,400 0,917
АСРН-1
п
1,03
1,08
1,14
1,20
135 22 10 36 62 106 152 4 9 4 3 8 4 12 36
0,007 0,000 0,000 0,000 0,000 0,038 0,040 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
0,918 0,955 1,000 0,805 0,823 0,840 0,809 1,000 0,556 1,000 1,000 0,750 0,750 0,500 1,000
0,053 0,000 0,000 0,056 0,064 0,094 0,092 0,000 0,333 0,000 0,000 0,125 0,000 0,250 0,000
0,023 0,045 0,000 0,139 0,113 0,028 0,059 0,000 0,111 0,000 0,000 0,125 0,250 0,250 0,000
АЭН
п 136 22 10 36 62 80 112 4 8 4 3 6 4 12 36
1,00 0,000 0,000 0,000 0,000 0,016 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
1,40 1,000 1,000 1,000 1,000 0,984 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000 1,000
О к о н ч а н и е т а б л. 1
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
ODH
n 28 4 4 0 30 52 12 0 0 2 1 0 0 0 0
0,95 0,035 0,000 0,000 0,000 0,000 0,019 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
1,00 0,895 1,000 1,000 0,000 0,901 0,943 0,917 0,000 0,000 0,500 1,000 0,000 0,000 0,000 0,000
1,05 0,035 0,000 0,000 0,000 0,000 0,038 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
1,10 0,035 0,000 0,000 0,000 0,066 0,000 0,083 0,000 0,000 0,500 0,000 0,000 0,000 0,000 0,000
1,20 0,000 0,000 0,000 0,000 0,033 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000
a-GPDH
n 1,00 2000002080004 12 6 1,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000 1,000 0,000 0,000 0,000 1,000 1,000 1,000
c-MDH
n 1,91 2000002080004 12 6 1,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000 1,000 0,000 0,000 0,000 1,000 1,000 1,000
m-MDH
n 1,00 2000002080004 12 6 1,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000 1,000 0,000 0,000 0,000 1,000 1,000 1,000
IDH
n 0,98 1,00 2000002080004 12 6 1,000 0,000 0,000 0,000 0,000 0,000 1,000 0,000 0,875 0,000 0,000 0,000 1,000 0,833 1,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,000 0,125 0,000 0,000 0,000 0,000 0,167 0,000
Примечание. Мук - Мукачево, Лт - Латорица, Вол - Воловец, Чер - Черновцы, Дн -Днепровская, Г-р - Гомель ручейная, Г-б - Гомель болотная, Орш - Орша, Лат - Латвия, Кр - Кропотово, Фин - Финляндия, Кар - Карелия, Нов - Новосибирск, Алт - Алтай, Тал - Талас; п - число проанализированных геномов.
На основании аллельных частот, используя Р-статистики Райта (Рзт) и О-статистики Неи (Оэт), мы предприняли попытку оценить состояние равновесия и степень подразделенности исследованных природных популяций Б. littoralis Восточной Евразии.
Показатель подразделенности Рэт для всех множественных аллелей подсчитывался как средневзвешенный по всем исследованным популяциям и варьировал в полиморфных локусах от 0,082 (а-Б81-3) до 0,148 (Р-Б81-2) (табл. 2). Среднее значение Рэт составило 0,056. Это говорит о том, что 94,4% всей изменчивости находится внутри популяций Б. littoralis и только 5,6% приходится на меж-популяционную изменчивость. Относительно близкое среднее значение, равное 0,076 (табл. 2), было установлено и по другому показателю, определяющему подразделенность, - Оэт, который, как было показано в работе Неи [16], эквивалентен параметру Рзт. Полученные значения Рзт и Оэт, приведенные в табл. 2, позволяют говорить об определенной неоднородности генетической структуры, по крайней мере, в изученной нами части ареала Б. littoralis, что может объясняться географической удаленностью и изолированностью ряда проанализированных популяций Восточной Евразии (см. рис. 1).
Генный поток вычислялся для всех исследованных групп популяций Б. Ш-toralis Восточной Евразии. Результаты сведены в табл. 3. Рассчитав значение Рэт для каждой группы популяций (табл. 3), мы определили величину генного потока (Кеш), которая оказалась для европейско-сибирско-тянь-шаньской группы популяций равной 4,21. Это говорит о том, что изученные популяции Б. littoralis обмениваются генетическим материалом в среднем с интенсивно-
стью 4,2 мигранта за поколение. При исключении популяции Тянь-Шаня величина Кеш для европейско-сибирской группы увеличивается до 4,5 мигрантов за поколение. Значение генного потока для европейской группы популяций при исключении популяций Алтая и Новосибирска увеличивается до 4,8 мигранта за поколение. Для юго-восточно-европейской - при исключении популяций Карелии и Финляндии - величина генного потока составила 5,2 мигранта за поколение и для восточно-европейских при исключении Карпатских популяций - 5,4 мигранта за поколение (табл. 3).
Т а б л и ц а 2
Значения показателей F- и G-статистик по 14 локусам в исследованных природных популяциях Б. НИотаОь
Локус GST
РОМ 0,130 0,145
МЕ 0,006 0,006
НК-1 0,027 0,035
НК-8 0,030 0,030
Рим 0,000 0,000
а-ЕЭт-3 0,082 0,107
Р-ЕБт-2 0,148 0,212
АСРН-1 0,114 0,150
АБН 0,015 0,015
ОБН 0,132 0,260
а-ОРБН 0,000 0,000
с-МБН 0,000 0,000
ш-МБН 0,000 0,000
1БН 0,105 0,105
Среднее 0,056 0,076
Т а б л и ц а 3
Показатели коэффициента генетической подразделенности и генного потока в группах популяций Б. НИотаОь
Группы популяций в« Nem
Европейско-сибирско-тянь-шаньские 0,056 4,21
Европейско-сибирские 0,053 4,47
Европейские 0,050 4,75
Юго-восточно-европейские 0,046 5,18
Восточно-европейские 0,044 5,43
Полученные данные однозначно указывают на достаточно интенсивный обмен генетическим материалом между исследованными популяциями Б. littoralis в Восточно-Евразийском регионе, несмотря на наличие между ними географических преград, таких как Карпатские, Уральские, Алтайские и Тянь-шаньские горы. Степень отличия в показателе Кеш хорошо соответствует характеру распределения и взаимосвязи популяций у Б. littoralis и напрямую связана с географической удаленностью популяций друг от друга.
Таким образом, в работе на основании проведенного генетического анализа 15 природных популяций Б. littoralis Восточной Евразии с использованием
14 генов были установлены основные показатели генетической подразделен-ности и генного потока. Величина генного потока (Nem) для всех исследованных нами популяций составила 4,21. Показано, что для отдельных групп популяций эта величина возрастает до 5,43. Полученные нами генетические данные однозначно указывают на зависимость величины генного потока от географической удаленности популяций друг от друга.
Литература
1. Гончаренко Г.Г., Митрофанов В.Г., Корочкин Л.И., Савицкий Б.П. Первый этап видооб-
разования у двух подвидов Drosophila группы virilis // ДАН СССР. 1989. Т. 304, № 2. С. 448-451.
2. Майр Э. Зоологический вид и эволюция. М.: Мир, 1968. 462 с.
3. Patterson S.T., Stone W.S. Evolution in the genus Drosophila. N.Y.: McMillan, 1952. 212 р.
4. Корочкин Л.И. Биология индивидуального развития (генетический аспект). М.: МГУ,
2002. 264 с.
5. Korochkin L.I. Genetic control and development expression of esterase isozymes in Droso-
phila of the virilis group // Developmental biology. London: Academic Press Inc., 1975. P. 99-117.
6. Nei M. Interspecific gene differences and evolutionary time estimated from electrophoretic
data on protein identity // Amer. Natur. 1971. Vol. 105. P. 385-398.
7. Spicer G.S., Bell C.D. Molecular phylogeny of the Drosophila virilis species group (Diptera:
Drosophilidae) inferred from mitochondrial 12S and 16S ribosomal RNA gene // Genes Ann. Entomol. Soc. Am. 2002. Vol. 95. P. 156-161.
8. Throckmorton L.H. The virilis species group // The genetics and biology of Drosophila / Eds.
by M. Ashburner, E. Novistky. London: Academic, 1982. Vol. 3B. P. 227-297.
9. Гончаренко Г.Г., Емельянов И.М. Электрофоретический ключ для типировки взрослых
особей двойниковых видов Drosophila группы virilis, обитающих в Палеарктике // Доклады АН СССР. 1990. Т. 313, № 2. С. 448-452.
10. Goncharenko G.G., Emelianov I.M. An electrophoretic key to adult members of the sibling species belonging to the Drosophila virilis group (Diptera, Drosophilidae) inhabiting Soviet Union and adjacent countries // Z. zool. Syst. Evolut.-forsch. 1992. Vol. 30. P. 281-286.
11. Lakovaara S., Saura A., Lankinew P., Pohjola L., Lokki P. The use of isoenzymes in tracing evolution and in classifying Drosophilidae // Zool. Scr. 1976. Vol. 5. P. 173-179.
12. Гончаренко Г.Г., Митрофанов В.Г., Катохин А.Н. Изучение биохимического полиморфизма у Drosophila imeretensis в природных популяциях Краснодарского края // Генетика. 1984. Т. ХХ, № 4. С. 620-627.
13. Сурков А.А., Гончаренко Г.Г., Митрофанов В.Г., Корочкин Л.И. Методический подход к исследованию генофондов короткоусых двукрылых Drosophila группы virilis в природных популяциях Беларуси // Известия ГГУ им. Ф. Скорины. 2003. № 5. С. 50-54.
14. Prakash S., Lewontin R.C., Hubby J.L. A molecular approach to the study of genic heterozygosity in natural populations. IV. Patterns of genic variation in central, marginal and isolated populations of Drosophilapseudoobscura // Genetics. 1969. Vol. 61. P. 841-858.
15. Wright S. The interpretation of population structure by F-statistics with special regards to systems of mating // Evolution. 1965. Vol. 19. P. 395-420.
16. NeiM. Molecular Population Genetics and Evolution. Amsterdam: Holland Press, 1975. 278 p.
17. Гончаренко Г.Г. Аллозимная диагностика видов-двойников Drosophila группы virilis // ДАН СССР. 1987. Т. 295, № 4. С. 976-980.
18. Гончаренко Г.Г., Сурков А.А., Митрофанов В.Г., Корочкин Л.И. Генетикоэволюционные и таксономические взаимоотношения у видов-двойников Drosophila группы virilis Палеарктики // Известия ГГУ им. Ф. Скорины. 2004. № 3. С. 144-157.
Поступила в редакцию 20.05.2010 г.
Grigory G. Goncharenko1, Alexander A. Surkov1,
Vladimir G. Mitrophanov2, Leonid I. Korochkin2
1Francisk’s Skorina Gomel State University, Gomel, Republic of Belarus 2N.K. Koltsov’s Institute of Developmental Biology of Russian Academy of Sciences, Moscow, Russia
THE LEVEL OF GENE FLOW IN NATURAL POPULATIONS OF DROSOPHILA LITTORALIS (DIPTERA: DROSOPHILIDAE)
OF EAST EURASIA
Fifteen natural populations of Drosophila littoralis Meigеn from East Eurasia, were investigated by starch-gel electrophoresis. A total of 41 alleles were observed at 14 loci. Interpopulations genetic diversity (FST, GST) was 5,6-7,6% of the total genetic diversity and the level of gene flow (Nem) was 4,2 migrants per generation. The genetic data obtained suggest that the level gene flow is connected with geographical distance of populations of D. littoralis from each other.
Key words: Drosophila littoralis; isoenzymes; populations structure; gene flow; inheritance.
Received May 20, 2010
