CC BY
130
УДК 575.174.015.3:582.475.2 Н.В. Орешкова, М.М. Белоконь
ГЕНЕТИЧЕСКАЯ ДИФФЕРЕНЦИАЦИЯ ПОПУЛЯЦИЙ ЛИСТВЕННИЦЫ СИБИРСКОЙ (LARIXSIBIRICA LEDEB.) ПО МИКРОСАТЕЛЛИТНЫМ ЛОКУСАМ*
Авторами статьи впервые исследована и применена на практике методика использования ядерных микросателлитных маркеров в анализе генетической изменчивости и уровня дифференциации лиственницы сибирской.
Ключевые слова: лиственница сибирская, микросателлитные локусы, генетическое разнообразие и дифференциация, сохранение генофонда.
N.V. Oreshkova, M.M. Belokon GENETIC DIFFERENTIATION IN SIBERIAN LARCH (LARIXSIBIRICA LEDEB.) POPULATIONS BY MICROSATELLITE LOCI
The technique of the nuclear microsatellite marker use in the analysis of genetic variability and differentiation level of the Siberian larch is researched and put into practice by the authors of the article for the first time.
Key words: Siberian larch, microsatellite loci, genetic diversity and differentiation, gene pool conservation.
Введение. Виды рода Larix Mill. являются важнейшими компонентами бореальных лесов северного полушария. Один из наиболее распространенных представителей этого рода - лиственница сибирская (Larix sibirica Ledeb.). Основные ее массивы приурочены к континентальным районам Южной Сибири, югозападной окраине Средне-Сибирского плоскогорья и лесотундры Западной Сибири, а также к северной части Монголии [Ирошников, 2004].
Ранее проблемы генетического разнообразия и дифференциации вида решались лишь с помощью изоферментных маркеров [Аллозимный полиморфизм..., 1989; Семериков, Матвеев, 1995; Semerikov, 1999; Ларионова, Яхнева (Орешкова), Кузьмина, 2003; Орешкова, Ларионова, 2006, 2007; Орешкова, 2010; Орешкова, Барченков, 2010; Биоразнообразие лиственниц..., 2010 и др.]. Однако для адекватной количественной оценки генетического разнообразия этих данных не хватает, поскольку из-за специфики данного класса маркеров ускользает основная часть генома, а полученные показатели генетической изменчивости существенно занижены [Сулимова, 2004].
В последние годы для исследований генетического полиморфизма хвойных начали активно использовать представителей нового поколения ДНК-маркеров - микросателлиты (варьирующие участки ДНК, состоящие из повторяющихся фрагментов длиной от 1 до 6 пар оснований), которые обладают высоким уровнем полиморфизма, большой частотой встречаются в геноме и достаточно равномерным распределением по длине генома.
Цель исследований. Используя ядерные молекулярные маркеры нового поколения, получить качественно новую информацию генетической изменчивости и уровня дифференциации лиственницы сибирской.
Объекты и методы исследований. В качестве объектов исследований взяты выборки из шести популяций лиственницы сибирской, произрастающие на территории Таймырского муниципального района Красноярского края и в Республиках Алтай и Тывы. Названия популяций и их местоположение представлены в табл. 1.
* Работа выполнена в рамках проектов РФФИ (№ 10-04-00786-а; № 11-04-00478-а; №11-04-92112-ЯФ_а) Интеграционного проекта СО РАН № 76 и ФЦП «Научные и научно-педагогические кадры инновационной России» (направление 1.1, Госконтракты 02.740.11.0281, 14.740.11.0164, 14.740.12.0826).
Таблица 1
Географическое расположение и характеристика исследованных популяций лиственницы сибирской
Популяция Численность Район расположения Географические координаты Высота над у.м., м
Республика Алтай
Черга 32 Шебалинский район Республики Алтай, окр. с. Черга 5Г29'с.ш. 85032'в.д. 600-650
Перевал 33 Шебалинский район Республики Алтай, окр. автодор. на Кукуй-ском перевале 5Г27'с.ш. 85015'в.д. 930-1200
Республика Тыва
Улар-1 30 Эрзинский район Республики Тыва, на север 32 км от с. Эрзин 50027'с.ш. 95030'в.д. 1873
Улар-2 30 Эрзинский район Республики Тыва, в 16-19 км севернее п. Морен 50031'с.ш. 95027'в.д. 2181
Таймырский район Красноярского края
Боганидское 30 Таймырский район Красноярского края, 26 км от г. Дудинка 69020'с.ш. 86052'в.д. 53
Талнах 30 Таймырский район Красноярского края, окр. п. Талнах .ш. .д. 'с. 'в. 00 S 23 00 с» SS 68 70
Для выделения ДНК использовали хвою. Выделения проводили по стандартному протоколу для растительных тканей с применением цетилтриметиламмониумбромида (CTAB].
Полное отсутствие экспериментальных данных по изменчивости ядерных микросателлитных локусов у лиственницы сибирской явилось основной трудностью при подборе праймеров. Последовательности использованных в работе праймеров были взяты из литературных данных, где они использовались для анализа других видов рода Larix. лиственницы японской (L kaempfeii Sarg.J - группа bcLK [Isoda, Watanabe, 2006], лиственниц альпийской (L lyallii Pari.J и западной (L occidental's Nutt.) группы UAKLY [Isolation, characterization..., 2000; Contrasting microsateiiite variation., 2006] и UBCLX [Development and characterization of microsateiiite, 2009]. Для проведения ПЦР использовали готовые реакционные смеси производства ООО «Лаборатория Изоген».
Характеристики микросателлитных локусов и условия ПЦР-амплификации приведены в табл. 2.
Таблица 2
Характеристика микросателлитных локусов, отобранных для анализа генетической изменчивости
лиственницы сибирской
Локус Мотив t°C отжига Число аллелей Размер фрагмента Источник литературы
bcLK056 (AG)20 Touchdown 63-53°С 21 142-196 Isoda, Watanabe, 2006
bcLK066 (TG)12 8 143-157
bcLK189 (AG)17AT(AG)6 23 158-202
bcLK224 (AG)17 8 128-148
bcLK225 (GA)20 28 160-224
bcLK232 (AG)19 8 133-149
bcLK260 (TG)14(AG)9 9 94-110
bcLK235 (TC)9(AC)2AG(AC)14 58°С 21 172-220
UBCLXtet_1-22 (TATC)9(TA)12 6 176-186 Chen et al., 2009
UAKLly6 (GT)17 23 214-262 Khasa et al., 2000, 2006
Продукты амплификации разделяли путем электрофореза в 6 %-м полиакриламидном геле с использованием Трис-БОТД-боратного электродного буфера. Гели окрашивали в растворе бромистого этидия и визуализировали с помощью системы гель-документации «Kodak». Молекулярный вес фрагментов определяли путем сопоставления со стандартными маркерами в программе Photo-Capt. Обработку полученных данных производили в программе GenAlEx 6.2 [Peakall, Smouse, 2006]. Для анализа связи между географическими и генетическими расстояниями использовали тест Мантела [Mantel, 1967].
Результаты исследований и их обсуждение. Из 10 проанализированных ядерных микросателлит-ных локусов у лиственницы сибирской все оказались полиморфными. Причем наиболее высокополиморфными оказались локусы bcLK056, bcLK189, bcLK225, bcLK235 и UAKLly6, у которых наблюдалось от 21 до 28 аллелей. У остальных локусов (bcLK066, bcLK224, bcLK232, bcLK260, UBCLXtet_1-22) выявленный полиморфизм был сравнительно ниже. У них обнаружено от 6 до 9 аллельных варианта (табл. 2). В исследовании, включившем 185 проанализированных особей, по перечисленным локусам нуль-аллели выявлены не были. Всего при анализе шести популяций лиственницы сибирской было выявлено 155 аллельных варианта, из них около 22 % оказались уникальными, встречающиеся только в какой-либо одной из популяций.
При анализе показателей генетической изменчивости исследованных популяций лиственницы (табл. 3) выявили, что популяции из Республики Тыва в среднем характеризовались наиболее высоким аллельным разнообразием и уровнем наблюдаемой и ожидаемой гетерозиготности (Na=9,6; Ho=0,488±0,046; He=0,673±0,054) по сравнению с алтайскими (Na=9,0; Ho=0,409±0,054; He=0,655±0,062) и таймырскими (Na=7,8; Ho=0,478±0,035; He=0,628±0,041).
Таблица 3
Показатели генетической изменчивости лиственницы сибирской, рассчитанные по результатам микросателлитного анализа
Популяции N Na Ne Ho He F
Черга 32 9,2 5,310 0,416±0,073 0,663±0,088 0,347
Перевал 33 8,8 4,938 0,401 ±0,083 0,646±0,091 0,342
В среднем для района 9,Q 5,124 Q,4Q9±0,054 Q,655±0,062 Q,344
Улар-1 30 9,6 4,788 0,510±0,058 0,695±0,061 0,261
Улар-2 30 9,6 5,341 0,467±0,074 0,651 ±0,093 0,239
В среднем для района 9,6 5,Q64 Q,4BB±0,046 Q,673±0,054 Q,25Q
Боганидское 30 7,7 3,777 0,493±0,066 0,618±0,068 0,209
Талнах 30 7,9 3,459 0,463±0,051 0,639±0,052 0,262
В среднем для района 7,B 3,61B Q,47B±0,035 Q,62B±0,041 Q,236
В среднем для изученных популяций 8,8 ±0,599 4,602 ±0,393 0,458 ±0,026 0,652 ±0,030 0,277 ±0,026
Примечание: Nа - среднее число аллелей на локус; Не - эффективное число аллелей на локус; Но - наблюдаемая гетерозиготность; Не - ожидаемая гетерозиготность; Р - индекс фиксации.
Установленные значения основных показателей генетического полиморфизма свидетельствуют о достаточно высоком уровне генетического разнообразия лиственницы сибирской в исследованных регионах (табл. 3) и находятся в пределах значений аналогичных показателей, рассчитанных для лиственниц японской, альпийской и западной [Contrasting microsatellite variation., 2006; Isoda, Watanabe, 2006; Development and characterization., 2009].
Сопоставление наблюдаемой и ожидаемой гетерозиготности показало, что во всех популяциях наблюдался дефицит гетерозиготных генотипов по большинству изученных микросателлитных локусов. Наиболее высокие значения индекса фиксации Райта (F) (табл. 3) были выявлены у популяций из Республики Алтай (F=0,344), а наименьшие у популяций из Таймыра (F=0,236). Это, по-видимому, можно объяснить большой фрагментированностью и низкой численностью популяций, а также самоопылением, приводящим к высокой степени инбридинга.
Исследование популяционной структуры с помощью F-статистик Райта [Wright, 1965; Guries, Ledig, 1982] показало, что каждое дерево в изученных популяциях лиственницы сибирской обнаруживает в среднем 29 %-й дефицит гетерозиготных генотипов (Fis=0,293) относительно популяции и почти 36 % дефицит гетерозигот (Fit=0,358) относительно вида (табл. 4).
Таблица 4
Значения показателей F-статистик Райта
Локус Fis Fit Fst
bcLK232 0,372 0,408 0,058
bcLK056 0,240 0,303 0,083
bcLK224 0,405 0,424 0,033
bcLK225 0,324 0,370 0,068
bcLK066 0,150 0,183 0,039
bcLK189 0,378 0,401 0,037
bcLK260 0,311 0,395 0,122
bcLK235 0,187 0,237 0,061
UBCLXtet_1-22 0,257 0,487 0,310
UAKLly6 0,304 0,371 0,097
Среднее 0,293±0,026 0,358±0,029 0,091 ±0,026
Чуть более 9 % всей наблюдаемой изменчивости приходится на межпопуляционную (Рзт=0,091). Внутри популяций сосредоточено 90,9 % всего генетического разнообразия. Максимально популяции дифференцированы по локусам иВ01Х1е_1-22, Ьс1К260, иЛШув и Ьс1К056.
На основании частот аллелей выявленных ядерных микросателлитных локусов проанализирована межпопуляционная дифференциация исследованных популяций лиственницы сибирской. В табл. 5 приведены показатели стандартных генетических расстояний М. Нея [№, 1972]. Генетическое расстояние 0 между популяциями колеблется от 0,067 до 0,469, составляя в среднем 0,254. Значения 0 слабо коррелировали с географическими расстояниями согласно тесту Мантела (г=0,552, Р=0,07).
Таблица 5
Генетические расстояния D между изученными популяциями лиственницы сибирской
Черга Перевал Улар-1 Улар-2 Боганидское Перевал
0,072 ***
0,189 0,184 *** Улар-1
0,311 0,336 0,115 *** Улар-2
0,414 0,422 0,201 0,124 *** Боганидское
0,469 0,454 0,280 0,176 0,067 Талнах
Генетическая дифференциация исследованных популяций соответствует их географической удаленности друг от друга. Наименьшие генетические отличия установлены внутри каждого района исследований: между популяциями из Таймырского района (Рэт=0,019; 0=0,067), Республик Алтай (Рэт=0,016; 0=0,072) и Тыва (Рэт=0,028; 0=0,115). Наибольшая дифференциация выявлена между популяциями из разных районов: таймырские и алтайские выборки (Рэт=0,095; 0=0,439), таймырские и тувинские (Рэт=0,043; 0=0,195), а также алтайские и тувинские (Рэт=0,058; 0=0,255). Анализ генетических дистанций с помощью многомерного шкалирования (РОД) хорошо иллюстрирует это (рис.).
0,150 0,100
СЧ
то 0,050 со
I 0,000
Ч а
§ -0,050 -0,100 -0,150
-0,300 -0,200 -0,100 0,000 0,100 0,200 0,300
Координата 1
Проекция изученных популяций лиственницы сибирской на плоскости двух координат по данным РСА-анализа (компьютерная программа вепАЕх У.6.2) матрицы генетических расстояний М. Нея
Таким образом, анализ генетических расстояний между популяциями показал отсутствие тесной взаимосвязи между географическим положением выборок и степенью их генетической дифференциации. Это объясняется тем, что лиственница сибирская произрастает в различных географических зонах с разнообразными экологическими условиями, образуя большое количество разновидностей, экотипов и климатипов. Влияние же экологических факторов перекрывает показатели генетической дифференциации. Подобные выводы были сделаны и по результатам изоферментного анализа [Орешкова, 2010; Орешкова, Барченков, 2010].
Заключение. В результате проведенных исследований с использованием десяти ядерных микроса-теллитных локусов нами впервые были получены оценки уровня генетического разнообразия и дифференциации лиственницы сибирской из разных районов ее распространения по данным маркерам. Представленные данные являются лишь начальным этапом огромной работы по изучению генетического разнообразия не только вида, но и рода в целом, выявления диагностических локусов, маркирующих выборки, выделения генетических резерватов, изучения межвидовой гибридизации лиственниц.
Благодарности. Авторы выражают огромную благодарность научному сотруднику лаборатории техногенных лесных экосистем ИЛ СО РАН, канд. биол. наук Д.Н. Орешкову за предоставленный материал из Республики Тыва и Таймырского района Красноярского края.
Литература
1. Биоразнообразие лиственниц Азиатской России / А.П. Абаимов, И.Ю. Адрианова, Е.В. Артюкова [и др.]. - Новосибирск: ГЕО, 2010. - 160 с.
2. Ирошников А.И. Лиственницы России. Биоразнообразие и селекция. - М.: ВНИИЛМ, 2004. - 182 с.
3. Ларионова А.Я., Яхнева Н.В. (Орешкова), Кузьмина Н.А. Генетическая изменчивость лиственницы сибирской в Нижнем Приангарье // Лесоведение. - 2003. - № 4. - С. 17-22.
4. Орешкова Н.В., Ларионова А.Я. Аллозимный полиморфизм лиственницы сибирской ^эпх вШов 1_е-с1еЬ.) // Лесные экосистемы Северо-Восточной Азии и их динамика: мат-лы междунар. конф. -Владивосток, 2006. - С. 220-223.
5. Орешкова Н.В., Ларионова А.Я. Генетическое разнообразие лиственницы сибирской в Ужурской лесостепи (Красноярский край) // Вестн. СВНЦ ДВО РАН. - 2007. - № 3. - С. 50-55.
6. Орешкова Н.В. Генетическая дифференциация лиственницы сибирской ^впх вШов 1_еСеЬ.) в Средней Сибири // Хвойные бореальной зоны. - 2010. - № 1-2. - С. 147-153.
7. Орешкова Н.В., Барченков А.П. Генетические особенности и морфологическая изменчивость лиственницы сибирской в Алтае-Саянской горной области // Вестн. КрасГАУ. - 2010. - № 10. - С. 59-64.
8. Семериков В.Л., Матвеев А.В. Изучение генетической изменчивости лиственницы сибирской ^впх вШов 1_СЬ.) по изоферментным локусам // Генетика. - 1995. - Т. 31. - № 8. - С. 1107-1113.
9. Сулимова Г.Е. ДНК-маркеры в генетических исследованиях: типы маркеров, их свойства и области применения // Успехи современной биологии. - 2004. - Т. 124. - № 3. - С. 260-271.
♦ Улар-2 ♦ Улар-1
♦ Боганидское ♦ Черга
♦
Перевал
♦ Талнах
10. Аллозимный полиморфизм лиственницы сибирской (Larix sibirica Ledeb.) / А.В. Шурхал, А.В. Подогас,
B.Л. Семериков [и др.] // Генетика. - 1989. - Т. 25. - № 10. - С. 1899-1901.
11. Development and characterization of microsatellite loci in western larch (Larix occidental's Nutt.) / C. Chen,
C. Liewlaksaneeyanawin, T. Funda [at el.] // Molecular Ecology Resources. - 2009. - Vol. 9. - P. 843-845.
12. A genetic linkage map for Pinus radiata based on RFLP, RAPD, and microsatellite markers / M.E. Devey, J.C. Bell, D.N. Smith [at el.] // Theor. Appl. Genet. - 1996. - Vol. 92. - № 6. - P. 673-679.
13. Guries R.P., Ledig F.T. Gene diversity and population structure in pitch pine (Pinus rigida Mill.) // Evolution. -1982.-Vol. 36. - P. 3537-402.
14. Isoda K., Watanabe A. Isolation and characterization of microsatellite loci from Larix kaempferi // Molecular Ecology. - 2006. - Vol. 6. - P. 664-666.
15. Isolation, characterization, and inheritance of microsatellite loci in alpine larch and western larch / D.P. Kha-sa, C.H. Newton, M.H. Rahman [at el.] // Genome. - 2000. - № 43 (3). - P. 439-448.
16. Contrasting microsatellite variation between subalpine and western larch, two closely related species with different distribution patterns / D.P. Khasa, J.P. Jaramillo-Correa, B. Jaquish [at el.] // Molecular Ecology. -2006. - Vol. 15. - P. 3907-3918.
17. Mantel N. The detection of disease clustering and a generalized regression approach // Cancer Research. -1967. - № 27. - P. 209-220.
18. Nei M. Genetic distance between populations // Amer. Natur. - 1972. - Vol. 106. - P. 283-291.
19. Peakall R., Smouse P.E. GenAlEx V6: Genetic analysis in Excel. Population genetic software for teaching and research // Molecular Ecology Notes. - 2006. - Vol. 6. - № 1. - P. 288-295.
20. Semerikov V.L., Semerikov L.F., Lascoux M. Intra- and interspecific allozyme variability in Eurasian Larix Mill. species // Heredity. - 1999. - Vol. 82. - P. 193-204.
21. Wright S. The interpretation of population structure by F-statistics with special regard to systems of mating // Evolution. - 1965. - Vol. 19. - P. 355-420.
УДК 630.907.1 Е.М. Рунова, В.А. Савченкова, О.И. Гринько
ОСОБЕННОСТИ СОСТОЯНИЯ ОСНОВНЫХ ЛЕСООБРАЗУЮЩИХ ДРЕВЕСНЫХ ПОРОД В ПЕРИОД СМЫКАНИЯ КРОН
В статье рассмотрена зависимость, показывающая динамику численности молодняка на вырубках разного возраста, выживаемость молодняков основных лесообразующих древесных пород в зависимости от возраста. Построены модели, позволяющие ее прогнозировать.
Ключевые слова: молодняк лесообразующих древесных пород, сомкнутость крон, тип вырубок.
E.M. Runova, V.A. Savchenkova, O.I. Grinko PECULIARITIES OF THE BASIC FOREST FORMING TREE SPECIES CONDITION IN THE CROWN CLOSURE PERIOD
The dependence which shows the young growth number dynamics on the cuttings of different age, young growth survival rate of the basic forest forming tree species, depending on age is considered in the article. The models which allow to forecast it are constructed.
Key words: young growth of the forest forming tree species, crown closure, type of cuttings.
Введение. Наиболее экономичным способом, позволяющим полнее использовать производительные силы природы, предотвратить нежелательную смену древесных пород и ускорить процесс выращивания леса на территории Приангарья, является содействие естественному возобновление леса путем сохранения в процессе рубки леса подроста и молодняка (тонкомера). Но полностью избежать трудоемкого процесса
