Genetic and morphological diversification in gastropods of the Baicaliidae family

Full Text

Введение Эндемичное семейство брюхоногих переднежаберных моллюсков Baicaliidae (Fischer, 1885) - самое многочисленное в видовом отношении семейство моллюсков оз. Байкал. В настоящее время известен 41 вид, принадлежащий к этому семейству из немногим более 110 эндемичных видов байкальских моллюсков (Ситникова, 2004). На данный момент семейство Baikaliidae подразделяется на восемь родов, семь из которых ранее относили к единственному роду Baicalia, кроме сильно отличающегося морфологически рода Liobaicalia (Кожов, 1936). Механизмы формирования такой разнообразной в морфологическом плане группы остаются во многом непонятными, и их исследование может оказаться важным для эволюционной теории в целом (Michel, 1994; Martens, 1997; Shcherbakov, 1999). Систематика байкалий основана на морфологических признаках и анатомии репродуктивной системы и неоднократно пересматривалась на протяжении более чем столетней истории изучения семейства (Ситникова, 2004). Этой группе посвящено несколько молекулярно-филогенетических исследований, относящихся к различным уровням их эволюционных взаимосвязей с другими моллюсками (Дарикова, Щербаков, 2009; Зубаков и др., 1997; Перетолчина и др., 2007; Hausdorf et al., 2001), однако движущие силы видовой радиации внутри этого очень быстро эволюционирующего семейства до сих пор остаются во многом неясными. Одной из важных причин возникших сложностей, очевидно, является очень высокая скорость эволюции байкалий. Действительно, первые же результаты сравнения нуклеотидных последовательностей фрагмента митохондриального гена первой субъединицы цитохром с-оксидазы (СО1) разных представителей сем. Baicaliidae заставили сделать предположение о взрывном видообразовании, приведшем к формированию «букета» видов (Зубаков и др., 1997), начавшемся не более 3,5 миллионов лет назад. Несмотря на молодость современных видов байкалий, их эволюционая дистанция до ближайших групп гастропод оказалась очень велика. На основании сравнения последовательностей митохондриальных генов 12S рРНК, 16S рРНК и фрагмента гена СО1 ряда представителей суперсемейства Rissoidea и представителей трех видов байкалий показано, что байкалии отделились от их общего с Bithyniidae и Amnicolidae предка около 140 млн лет назад (Hausdorf et al., 2001). Детальное исследование генетического полиморфизма нескольких видов рода Baicalia показало, что видообразование в этой группе происходило довольно быстро, виды этого рода достаточно молоды и размах внутривидового полиморфизма часто сравним с генетическими расстояниями между представителями хорошо морфологически и экологически различающихся видов (Перетолчина и др., 2007). Поскольку митохондриальные маркеры отражают эволюционную судьбу митохондрионального генома, которая может существенно отличаться от эволюционной истории соответствующих видов, большое распространение приобретают ядерные некодирующие (интронные) EPIC-маркеры (exon-primed intron-crossing), скорость эволюции которых может быть сравнима с таковой для митохондриальных маркеров (Feher et al., 2012). Ранее нами была предпринята попытка использования последовательности интрона ядерного гена фосфофруктокиназы для установления межродовых отношений байкалий. Отмечена достаточно высокая «разрешающая способность» интронного маркера, древо построенное по данной последовательности для представителей семи видов байкалий по топологии соответствовало древу, полученному по последовательности фрагмента гена СО1, но оба они содержали политомии (Дарикова, Щербаков, 2009). Таким образом, к настоящему времени накопился довольно большой объем информации не только по сравнению видов байкалий по морфологическим и анатомическим признакам, но и по генетическому полиморфизму различных маркеров как митохондриальных, так и ядерных. С увеличением массива разнородных сведений о байкалиях возникает возможность синтеза всей этой информации для получения максимально полной картины эволюции в рамках этого семейства. В настоящей работе на основании сравнения полиморфизма последовательностей интрона α-субъединицы АТФ-синтазы, ранее опубликованных нуклеотидных последовательностей митохондриальных генов, а также морфологических описаний соответствующих видов, сделана первая попытка получения филогенетического древа байкалий, объединяющего все признаки, основанная на едином байесовском подходе (Rieppel, 2005; Ronquist et al., 2012). Материалы и методы Сбор и фиксация образцов Образцы одиннадцати видов моллюсков сем. Baicaliidae были собраны в летний период 2009 г. с помощью водолазов или драги с глубин 5-90 м в разных районах оз. Байкал (табл. 1, рис. 1). Пробы фиксировали 80 %-м этиловым спиртом, через сутки фиксатор заменяли на 70 %-й спирт. Экстракция ДНК ДНК экстрагировали по модифицированному методу Дойла (Doyle, Doyle, 1987) из кусочка мышечной ткани ноги моллюска. Амплификация и определение нуклеотидных последовательностей Амплификацию фрагмента ядерного гена α-субъединицы АТФ-синтазы проводили с использованием универсальных праймеров, фланкирующих единственный интрон этого гена (Jarman et al., 2002). Праймеры имеют следующую структуру: ATPSa F: 5′-GAGCCMATGCAGACTGGTATTAAGGCYGT-3′ ATPSa R: 5′-TTGAANCKCTTCTGGTTGATGATGGTGTC-3′, где R - G/A; K - G/T; Y - C/T; M - C/A; D - A/G/T; H - A/C/T; N - любой нуклеотид (по номенклатуре IUPAC). Параметры при 35 циклах амплификации были следующими: преденатурация ДНК при 94 °С - 2 мин, денатурация ДНК при 94 °С - 20 с, отжиг праймеров при 54 °С - 1 мин, элонгация нуклеотидной цепи при 72 °С - 1 мин (+5 минут на последнем цикле). Секвенирование продуктов амплификации проводилось в ЗАО «Синтол» (г. Москва). Выравнивание нуклеотидных последовательностей с учетом вторичной структуры проводили с помощью программы Мafft v.6 (Katoh et al., 2009). Для установления возможных рекомбинантных событий использовали методы, представленные в программe RDP 3 (RDP, GENECONV, Bootscan, Chimera, SciScan) (Martin et al., 2010). Филогенетический анализ Весь набор данных был разбит на 5 групп: интронная последовательность, каждое положение кодона митохондриального гена СО1 отдельно и морфологические признаки. Морфологические признаки (табл. 2, 3) превращали в наборы бинарных признаков с помощью программы Factor из пакета программ Phylip-3.69 (Felsenstein, 1989) с учетом гипотез о вероятных изменениях этих признаков (Кожов, 1936; Ситникова, 1991; Sitnikova et al., 2001). Нуклеотидные последовательности фрагмента гена СО1 байкалий соответствующих видов с номерами AF445328.1, FN185835.1, FN185867.1, HM543405.1, HQ113296.1, Z92998.1, Z93000.1, Z93002.2, Z93003.2, Z93005.2 были взяты из базы данных NCBI, за исключением вида Korotnewia angigyra, для которого данная информация отсутствует. Филогенетические древа были построены байесовским методом с помощью программы MrBayes 3.2.1 (Ronquist et al., 2012). Для байесовского анализа нуклеотидных данных была выбрана модель нуклеотидных замен GTR + gamma. Для построения каждого из древ запускались 1 «холодная» и 3 «горячих» марковских цепи длиной 5 млн шагов, для построения финального древа отбиралось каждое тысячное сгенерированное древо, первые 25 % древ не учитывались. Результаты Анализ нуклеотидных последовательностей Определены нуклеотидные последовательности интрона α-субъединицы АТФ-синтазы для представителей одиннадцати видов из семи родов эндемичных байкальских моллюсков сем. Baicaliidae. Полученные последовательности экспонированы в Genbank с номерами KF201695-KF201709. Длина этого интрона варьирует от 509 до 524 н. п. у моллюсков разных видов благодаря многочисленным, но не протяженным вставкам и делециям (табл. 4). Короткие (1-2 н. п.) индели чаще всего связаны с мононуклеотидными трактами длиной около 5 н. п. В нашем случае четыре из пяти более крупных делеций (от 5 нуклеотидов) связаны с потерей одного из повторов в тандемно повторенных минисателлитных или минисателлит-подобных последовательностях. Блоки, содержащие индели, были исключены из последующего филогенетического анализа, в результате длина анализируемых последовательностей составила 482 н. п., 40 из которых были вариабельными. Количество трансверсий в данном интроне на 30 % выше количества транзиций. Филогенетический анализ Для всех молекулярных данных использовалась модель нуклеотидных замен GTR+gamma (значение gamma определялось независимо для каждого из четырех наборов нуклеотидных данных); для морфологических признаков - броуновская модель (Brown, Lemmon, 2007). Древо полученное на основании только интронной последовательности является более разрешенным, по сравнению с ранее полученным древом по фрагменту гена СО1, несмотря на две политомии (Зубаков и др., 1997). Представители части родов формируют отдельные клады, однако представители родов Godlewskia, Korotnewia и Parabaikalia объединяются в общую кладу (рис. 2). Минимальное генетическое расстояние (0,2 % нуклеотидных замен) - между представителями сестринских видов Baicalia carinata и Baicalia rugosa, максимальное (3,5 %) - между видами Maackia herderiana и Parabaikalia florii. В результате объединения данных по нуклеотидным последовательностям фрагмента гена СО1 и интрона α-субъединицы АТФ-синтазы, и морфологических признаков было получено дихотомическое филогенетическое древо представителей семи родов моллюсков сем. Baicaliidae (рис. 3). Самые короткие длины ветвей отмечаются для пар сестринских видов: Baicalia carinata и Baicalia rugosa; Parabaikalia kobeltiana и Parabaikalia florii. Виды родов Korotnewia, Parabaikalia и Pseudobaikalia кластеризуются вместе. Представители остальных родов формируют отдельные клады. Обсуждение Интрон α-субъединицы АТФ-синтазы, наряду с митохондриальными молекулярными маркерами, ранее уже использовался как филогенетический маркер в ряде исследований (Feher et al., 2012; Freya, Vermeij, 2008; Hedtke et al., 2011; Keever et al., 2009; Puritz et al., 2012; Teske et al., 2011). Показано, что этот маркер применим для исследований популяционной структуры некоторых видов беспозвоночных (морских звезд и двустворчатых моллюсков). В частности, был отмечен высокий уровень гетерозиготности в популяциях морской звезды Patiria miniata, связанный главным образом с присутствием многочисленных коротких инделей, при относительно небольшом количестве нуклеотидных замен (Keever et al., 2009). Также данный маркер успешно использовался на более высоком таксономическом уровне. При сравнении топологий древ, полученных по нуклеотидной последовательности данного интрона и древ, построенных по фрагментам митохондриальных генов (16S рРНК и СО1) для брюхоногих морских моллюсков рода Nerita незначительное несоответствие топологии полученных древ было отмечено лишь для комплекса вида N. undata, который, предположительно появился в начале плиоцена, около 4,3 млн лет назад, и является одной из самых «молодых» групп, из рассматриваемых в данной работе (Freya, Vermeij, 2008). В интроне α-субъединицы АТФ-синтазы моллюсков сем. Baicaliidae отсутствуют протяженные индели (табл. 4), небольшие вариации длины чаще всего связаны с изменением числа повторов в микро- и минисателлит-подобных последовательностях. Известно, что наличие инделей может повышать вероятность нуклеотидных замен в близлежащих регионах ДНК (на расстоянии до нескольких сотен нуклеотидных пар), преимущественно трансверсий в гетерозиготах, содержащих обе аллели, которыми представлена индель. Хотя большое значение имеет «региональный» эффект, то есть сайты более подверженные нуклеотидным заменам одновременно подвержены и инделям (McDonald, 2011). Несмотря на преобладание среди нуклеотидных замен трансверсий, можно отметить, что их распределение практически не связано с присутствием инделей, и накопление трансверсий главным образом происходит у видов морфологически более далеких родов. Филогенетическое древо построенное по последовательности интрона α-субъединицы АТФ-синтазы байкалий содержит меньше политомий по сравнению с древом, полученным ранее по последовательности митохондриального гена СО1. Интронный маркер позволяет получить статистически достоверные клады не только на межвидовом таксономическом уровне, но и на популяционном уровне, в частности для полиморфного вида Baicalia carinata (рис. 2). Ранее популяционная структура этого вида была установлена с использованием последовательности фрагмента митохондриального гена СО1, в то время как изменчивость внутреннего транскрибируемого спейсера (ITS1) не отражала популяционной структуры. При этом показано, что восточная популяция вида B. carinata дала начало виду B. rugosa, причем по маркеру ITS1, в отличие от митохондриального маркера, генетические расстояния между видами B. carinata и B. rugosa были больше, чем генетические расстояния на уровне популяций B. carinata (Перетолчина и др., 2008). Ранее на основании морфологических признаков вид Baicalia carinata также разделялся на два подвида, обитающие на разных: северо-восточном и юго-западном ареалах (Ситникова, 2004). Полученное древо по последовательности интрона α-субъединицы АТФ-синтазы главным образом соответствует морфологическим сведениям о рассматриваемых видах, за исключением положения вида Godlewskia wrzesniowski, который кластеризуется с родами Korotnewia и Parabaikalia, морфологически будучи сильно отличным от них. В результате использования всех данных известных для представителей рассматриваемых видов сем. Baicaliidae было получено дихотомическое древо (рис. 3), в отличие от крайне политомического древа, полученного ранее по фрагменту гена цитохром-с-оксидазы (Sherbakov, 1999; Щербаков, 2003). На полученном древе видов представители разных видов образуют статистически достоверные клады. Однако виды родов Korotnewia, Parabaikalia и Pseudobaikalia кластеризуются вместе, и если роды Korotnewia и Parabaikalia достаточно близки как экологически, так и морфологически, то положение вида Pseudobaikalia pulla tenuicosta оказалось достаточно неожиданным, т. к. род Pseudobaikalia имеет ряд морфологических признаков значительно отличающих его от других родов сем. Baicaliidae. Однако ранее на основании данных по интронной последовательности ядерного гена фосфофруктокиназы было показано, что другой представитель рода Pseudobaikalia (Pseudobaikalia zachwatkini) и представитель вида Korotnewia semenkewitschi имеют генетическую дистанцию, сходную с межвидовыми дистанциями рода Baikalia (Дарикова, Щербаков, 2009). На древе, построенном только по митохондриальному маркеру СО1, также очень короткой оказалась генетическая дистанция между видами Korotnewia korotnewia gracilis и Parabaikalia elata (Щербаков, 2003). Интересно, что генетически близкими оказались виды со сходными биотопами: исследованные виды родов Korotnewia, Parabaikalia и Pseudobaikalia - зарываются в песок; Teratobaikalia duthiersii и Maackia herderiana чаще встречаются на каменистых грунтах; вид Godlewskia wrzesniowski предпочитает песчаные грунты, но редко зарывается; а виды рода Baicalia обитают как на песчаных, так и каменистых грунтах (Ситникова, 2004). Несоответствия морфологических и молекулярных данных представляют большой интерес для дальнейших исследований. Возможно, такая топология полученного древа может быть связана с быстрой морфологической эволюцией родов Korotnewia, Parabaikalia и Pseudobaikalia.

About the authors

Mariya Vladimirovna Kovalenkova

Limnological Institrute, Siberian Division, Russian Academy of Sciences

Email: kovalenkovam@mail.ru
postgraduatе, рrincipal engineer. Laboratory of Molecular Systematics

Tatyana Yakovlevna Sitnikova

Limnological Institrute, Siberian Division, Russian Academy of Sciences

dr., senior researcher. Laboratory of Biology of aquatic invertebrates

Dmitriy Yuryevich Shcherbakov

Limnological Institrute, Siberian Division, Russian Academy of Sciences

Email: sherb@lin.irk.ru
dr., Head of Laboratory of Molecular Systematics

References

Supplementary files