О регулировании состава экологического сообщества с помощью

advertisement
БИОФИЗИКА. 2004. Т. 49 № 5. С. 912-919
УДК 57.02.001.57
О РЕГУЛИРОВАНИИ СОСТАВА ЭКОЛОГИЧЕСКОГО СООБЩЕСТВА С
ПОМОЩЬЮ ИЗМЕНЕНИЯ ОТНОШЕНИЯ КОНЦЕНТРАЦИЙ РЕСУРСОВ В
СРЕДЕ
П.В. Фурсова, А.П. Левич
Биологический факультет МГУ им. М.В.Ломоносова, 119992, Москва
Изучена возможность регулирования состава сообщества, которую предоставляет
вариационная модель экологического сообщества. Для R-, S- и M- диссоциантов
бактерий Pseudomonas aeruginosa получены зависимости относительных численностей микроорганизмов от начальных отношений запасов питательных веществ
в среде. Показана теоретическая возможность управления составом сообщества.
Смена доминирующего вида при изменении отношения ресурсов, как по модельным расчетам, так и в эксперименте, проиллюстрирована на примере сообщества
протококковых водорослей.
Ключевые слова: состав сообщества, регулирование состава сообщества, вариационная модель, Pseudomonas aeruginosa.
Одной из важнейших целей количественной экологии сообществ является умение
рассчитывать численность каждой из входящих в сообщество популяций организмов как
функцию доступных ресурсов среды. Практическое применение подобных расчетов заключается в возможности управлять структурой сообщества с помощью внешних факторов – изменять видовой состав и доминирование видов, целенаправленно увеличивать или
уменьшать численности выбранных групп организмов.
В настоящей работе продемонстрированы возможности регулирования состава смешанной культуры микроорганизмов методом, предсказанным вариационной моделью экологического сообщества.
ВАРИАЦИОННАЯ МОДЕЛЬ
Предложена модель сообщества одноклеточных организмов, потребляющих ресурсы, которые не могут заменить друг друга, поскольку выполняют различные функции по
отношению к росту. В модели допустимы деление и смертность клеток, но не их слияние.
В лабораторных условиях описываемая модель соответствует культивированию без добавления или изъятия ресурсов и микроорганизмов. Изучено развитие поликультуры до
2
остановки роста, вызванной исчерпанием одного из ресурсов, но не какими-либо иными
причинами.
Стационарная стадия роста сообщества описана решением вариационной задачи
w
ì
=
H
(
n
,
n
,...,
n
)
n
ln
n
ni ln ni ® extr ,
å
w
1
2
ï
i =1
ï w
ï å qik ni £ Lk , k = 1, m,
í i =1
ïw
ïå ni = n,
ï i =1
ïîni ³ 0, i = 1, w,
(1)
где ni – численности групп организмов, образующих сообщество, w – число групп в сообществе, (L1, L2,..., Lm) – количества всех потребляемых сообществом ресурсов среды в начале опыта, qik – количество k-ого ресурса, необходимое для деления клеток группы i, в
расчете на одну клетку (потребность организма группы i в ресурсе k) [1, 2].
Целевая функция H (n1 , n2 ,..., nw ) для допустимых в модели преобразований рассчитана на основе теоретико-категорной модели [3].
Аналитические результаты исследования задачи в зависимости от начального содержания ресурсов в среде позволяют указывать те, которые будут полностью потреблены, и рассчитывать доли каждого вида и общую численность сообщества на стационарной
стадии роста.
Выявление лимитирующих (полностью потребленных) ресурсов проведено на основе теоремы стратификации [1]. Согласно этой теореме все пространство ресурсных фактоm
ров
ÕL
k
m
распадается (стратифицируется) на
- непересекающихся областей (стратов),
k =1
каждая из которых соответствует одному из подмножеств множества потребляемых сообществом ресурсов. В страте S J , где J ¹ Æ — подмножество множества ресурсов
{1,2,...,m}, выполнено:
1) решение задачи (1) ni ( L1 , L2 ,..., Lm ) зависят только от тех Lk, для которых kÎJ;
w
2) на этом решении нестрогие неравенства
åq n
i= 1
k
i i
£ Lk обращаются в строгие равенства
для всех kÎJ и в строгие неравенства для всех kÏJ.
Теорема стратификации влечёт редукцию задачи (1) к задачам:
ì H ( n1 , n2 ,..., nw ) ® extr ;
ïw
ï
j
j
íå qi ni = L , j Î J ;
ï i =1
ïn ³ 0, i = 1, w,
î i
(2)
3
формулируемым для любого J Ì {1,2,...,m}.
Решение задач (2) называется формулой видовой структуры и имеет вид:
r
ì
ü
ni LJ = n exp í-å l k qik ý ,
î kÎ J
þ
( )
w
r
где n = å ni , "вектор" LJ имеет компоненты Lj с индексами из набора J, идентифициi =1
r
рующего страт, которому принадлежит вектор LJ . Множители Лагранжа l k и полная
r
численность n как функции потребляемых полностью в страте S J ресурсов LJ ищутся из
алгебраических уравнений
ìw
ì
k kü
ïå exp í-å l qi ý = 1;
ï i =1
î kÎ J
þ
í w
ïn q j exp ì- l k q k=ü
í å
i
i ý
ïî å
i =1
î kÎ J
þ
Lj ,
j Î J.
Таким образом, в рамках вариационной модели экологического сообщества для определения его состава на стационарной стадии роста сначала необходимо выяснить, к какому страту принадлежит вектор ресурсов, задающий состав исходной среды, а затем по
формуле видовой структуры рассчитать численности групп организмов, образующих сообщество.
Возможность управления составом сообщества предоставляет следующая особенность решения вариационной задачи. Установлено, что относительные численности видов
на стационарной стадии роста зависят только от отношений начальных концентраций
полностью потребляемых сообществом ресурсов среды [4]. Именно это свойство позволяет регулировать состав сообщества: изменяя отношения запасов питательных веществ в
начальной среде, получать различные относительные численности видов на стационарной
стадии роста.
4
СОСТАВ СООБЩЕСТВА НА СТАЦИОНАРНОЙ СТАДИИ РОСТА
Для демонстрации указанных возможностей, предлагается приложение вариационной модели к описанию развития сообщества R-, S-, и M-диссоциантов штамма Pseudomonas aeruginosa К-2. В период с 1999 по 2001 гг. на кафедрах микробиологии и общей
экологии Биологического факультета МГУ им. М.В.Ломоносова были проведены опыты
по выращиванию моно- и смешанных культур (все попарные сочетания диссоциантов и
поликультура из трех диссоциантов). Частота диссоциативных переходов бактерий составляет 10-2-10-4 на одно клеточное деление, что позволяет считать сообщество диссоциантов удовлетворяющим предпосылкам вариационной модели (т.е. пренебречь явлением
диссоциации в сравнении с делением). Бактерии культивировали без пополнения запаса
ресурсов на средах с различным начальным содержанием глюкозы, нитратов и фосфатов.
По данным опытов с монокультурами рассчитали потребности диссоциантов в питательных веществах (табл.1) [5].
Таблица 1. Значения потребностей в ресурсах среды (10-12мг/кл) и их погрешности для
диссоциантов P.aeruginosa [5]
Ресурс
Диссоциант
углерод
азот
фосфор
R
129 ± 22
7.0 ± 2.5
1.0 ± 0.5
S
409 ± 55
17.5 ± 3.0
4.0 ± 1.5
M
525 ± 114
31 ± 3
6.0 ± 2.5
Затем, используя полученные значения, по разработанным алгоритмам [1, 6] построили области лимитирования. Для проверки адекватности предсказанной стратификации [7] применяли метод добавок. В момент предполагаемого достижения стационарной
стадии культуру разделяли на четыре пробирки. В три из них вносили добавки – глюкозу,
нитрат или фосфат, – четвертую оставляли без изменений (контроль). Культуру оставляли
расти еще четыре часа, затем проводили измерение оптической плотности. Если добавление компонента питания в среду приводило к возобновлению деления клеток, то этот ресурс считали лимитирующим. В случае, если внесение вещества не приводило к росту
культуры, считали, что этот фактор не ограничивал развития сообщества. Согласно полученным результатам (табл.2), в 75% случаев возобновление деления клеток происходило
после внесения питательного вещества, предсказанного моделью как лимитирующего [7].
5
Таблица 2. Влияние добавок на оптическую плотность моно- и смешанных культур [7]
номер лимитирующий реопыта
сурс по модели
1
С
2
N
3
P
4
C
5
C
6
N
7
P
8
C
9
N
10
C
11
C
12
C
13
C
14
P
добавка
R
S
M
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
C
N
P
121
91
62
220
102
100
220
97
106
140
80
115
227
73
100
110
96
94
101
131
100
117
119
125
217
74
109
93
93
98
94
95
134
95
143
95
212
162
235
129
86
86
131
76
87
62
100
114
100
112
104
122
100
100
129
100
93
101
68
57
RS
RM
SM
RSM
92
104
96
142
75
75
82
59
100
96
104
92
160
100
100
89
84
78
104
116
100
160
90
90
121
79
83
87
98
81
96
128
100
150
71
71
120
100
80
103
81
73
103
88
94
120
93
87
143
102
92
109
81
95
100
68
105
90
100
97
160
90
110
97
66
89
143
86
100
119
81
72
112
91
108
Возможной причиной отсутствия влияния добавок на рост бактерий является недостаточный четырехчасовой срок экспозиции культур с добавками.
Стратификация пространства ресурсов при фиксированном уровне фосфора в начальной среде для смешанной культуры трех и двух (RS) диссоциантов представлена на
113
75
100
129
91
82
6
рис.1а и 1б (для культуры из двух диссоциантов областью лимитирования трех ресурсов
является отрезок на границе IV, V и VI областей – см. рис.1б).
а)
б)
Рис.1. Стратификация пространства ресурсных факторов для смешанной культуры трех
(а) и двух (б) диссоциантов. В области I лимитирует углерод (С), в II – азот (N), в III –
фосфор (P), в IV – углерод и азот, в V – углерод и фосфор, в VI – азот и фосфор, в VII –
все три ресурса.
С C
Вид областей лимитирования в координатах æç ; ö÷ для смеси двух диссоциантов
èN Pø
изображен на рис.2.
Рис.2. Стратификация пространства ресурсов в координатах æç ; ö÷ для смеси двух дисèN Pø
С C
социантов (RS). Горизонтальные пунктирные линии соответствуют значениям отношения
7
C/P, для которых проведены расчеты зависимости относительных обилий диссоциантов
как функций отношения C/N (см. рис.3).
В каждом из стратов, используя формулу видовой структуры, рассчитали относительные численности групп организмов, образующих сообщество (вывод соответствующих формул содержится в одной из предшествующих работ [6]). В работе полностью
представлен случай R- и S-диссоциантов в смеси. Зависимость относительных численностей видов в сообществе от отношения ресурсов в среде можно проследить по рис.3.
а)
б)
в)
г)
Рис.3. Зависимость относительных численностей (отн.числ.) R- и S-диссоциантов (r, s) от
отношения углерода к азоту (C/N) в исходной среде. а) C/P
C/P = 112.5; г) C/P = 112.9.
Î (0,110.83]; б) C/P = 111; в)
8
Каждый из рис. 3а, 3б, 3в, 3г (все возможные конфигурации состава культуры в зависимости от соотношений питательных веществ в начале культивирования) соответствует фиксированному отношению C/P. В областях I, II, III относительные численности постоянны. Численные значения задаются выражениями r = x0qR , s = x0qS , где x0 – корень
k
k
k
уравнения x qR + x qS = 1 , r и s – относительные численности R- и S-диссоциантов, индекс k
k
обозначает ресурс, а параметры qRk , qSk – соответствующие потребности (значения из
табл.1). В области IV графиками зависимостей относительных численностей от отношения C/N являются участки гипербол, задаваемых уравнениями:
s=
r=
-17.5x + 409
,
-10.5x + 280
7x - 129
C
, где x =
. В области ограничения роста азотом и фосфором (область
-10.5x + 280
N
VI) решение соответствующей редуцированной задачи (2) задается формулами
r=
4z - 17.5
-z + 7
N
, s=
, где z =
. Для построения графиков этих зависимостей на
3z - 10.5
3z - 10.5
P
плоскости (C/N; отн.числ.) (рис.3) значения z , соответствующее области VI, рассчитывается исходя из заданного уровня отношения
ношением
C
C C
, а именно, как частное : . В области V постоянные относительные чисN
P N
ленности при фиксированном отношении
s=
C
и интервала значений, принимаемых отP
4h - 409
C
определяются по формулам r =
,
P
3h - 280
-h + 129
C
, подстановкой в них конкретного значения h = .
3h - 280
P
В серии экспериментов, проведенных со смешанными культурами бактерий
P.aeruginosa среды были ограничены по одному из питательных веществ. В соответствующих областях однофакторного лимитирования предсказываемый моделью состав
культуры сравнивали с данными эксперимента [8]. Анализ результатов (табл.3) показывает, что из 36 наборов данных в 29 расчетный состав поликультур отличается от экспериментального не более, чем на 12 %. Максимальное отклонение предсказанных моделью
значений от опытных составляет 17%. Полученные различия могут быть объяснены
ошибками в значениях потребностей, применявшихся в расчетах (табл.1), о чем свидетельствует проведенный анализ чувствительности модели к вариациям параметров [9].
БИОФИЗИКА. 2004. Т. 49 № 5. С. 912-919
Таблица 3. Состав смешанных культур на стационарной стадии роста [8]
номер
1
2
3
культура
R
S
S
M
R
S
R
M
S
M
R
S
M
R
S
S
M
R
S
M
R
S
M
состав по результатам эксперимента (в процентах)
состав по модельным расчетам (в
процентах)
67
33
69
31
62
38
71
29
50
50
43
41
16
70
30
50
50
67
11
22
43
41
16
66
34
54
46
65
35
74
26
60
40
60
29
11
69
31
54
46
60
29
11
60
29
11
номер
4
5
6
культура
R
M
S
M
R
S
R
M
S
M
R
S
M
R
M
S
M
R
S
M
состав по результатам эксперимента (в процентах)
состав по модельным расчетам (в
процентах)
60
40
42
58
61
39
62
38
50
50
62
32
6
78
22
43
57
68
17
15
73
27
54
46
69
31
73
27
54
46
62
23
15
73
27
54
46
62
23
15
номер
7
8
9
культура
R
S
R
M
S
M
R
S
M
R
S
M
R
S
R
M
S
M
R
S
M
состав по результатам эксперимента (в процентах)
состав по модельным расчетам (в
процентах)
79
21
87
13
52
48
75
12
13
68
19
13
60
40
86
14
65
35
60
27
13
69
31
73
27
54
46
62
23
15
62
23
15
72
28
78
22
57
43
68
22
10
номер
10
11
12
культура
R
S
R
M
S
M
R
S
M
R
S
R
M
R
S
M
R
M
состав по результатам эксперимента (в процентах)
состав по модельным расчетам (в
процентах)
62
38
84
16
72
28
63
22
15
70
30
80
20
45
31
24
75
25
72
28
78
22
57
43
68
22
10
65
35
74
26
60
29
11
74
26
БИОФИЗИКА. 2004. Т. 49 № 5. С. 912-919
Провести опытную проверку предсказываемых моделью значений во всем диапазоне изменений отношений ресурсов не представляется возможным. Это связано, во-первых, с
узостью интервалов Lk/Ls (k, s = 1, 2, 3), задающих области многофакторного лимитирования (см. рис.2 и рис.3), и с точностью, с которой могут быть составлены исходные среды,
а во-вторых, с незначительными различиями состава сообщества в стратах в сравнении с
погрешностью определения долей диссоциантов.
Регулирование состава возможно лишь в задаваемых моделью пределах. Как видно
из рис.3, для видов, обладающих полученными значениями потребностей, нельзя изменить доминирующую группу организмов. Однако, при "подходящих" значениях потребностей, варьируя соотношение питательных веществ в среде, на стационарной стадии роста можно получать различные распределения видов по рангам (от наименее до наиболее
многочисленного). В качестве демонстрации указанного свойства приведем расчеты состава сообщества, состоящего из четырех протококковых водорослей, потребляющих азот
(N) и фосфор (P) [10]. На рис.4 представлены относительные численности видов
pi , i = 1, 4 , исходя из следующих экспериментальных значений потребностей (размерность
109 мг/кл):
q1N = 3.94; q 2N = 0.91; q3N = 1.69; q 4N = 0.41;
q1P = 0.79; q 2P = 0.28; q3P = 0.22; q 4P = 0.21.
Рис.4. Относительные численности pi , i = 1, 4 для экспериментальных потребностей как
функции от N/P [10].
11
По результатам лабораторных экспериментов смену доминирующего вида можно
проследить на примере культивирования двух видов из порядка протококковых Ankistrodesmus falcatus и Scenedesmus quadricauda [10, 11]. Результаты исследования роста сообщества на средах с различными начальными отношениями азота и фосфора (N/P = 1.3; 57)
проиллюстрированы на рис.5.
Рис.5. Зависимость относительных численностей видов от отношения азота к фосфору в
лабораторных экспериментах с двухвидовыми культурами.: 1 – S. quadricauda, 2 – A. falcatus [10].
В заключении отметим еще одно свойство вариационной задачи, имеющее отношение к управлению (или регулированию) составом сообщества: относительная численность
заданного вида принимает наибольшее значение при отношениях в среде ресурсов, равных отношениям потребностей в них заданного вида [4]. Данные экспериментов подтверждают этот теоретический факт. Так, например, полученное из опытов отношение потребностей в азоте и фосфоре для A. falcatus равно 7, а для S. quadricauda – 25. Следовательно, на средах с отношениями ресурсов в среде 1.3 и 57 должен преобладать анкистродесмум и сценедесмус, соответственно. Именно это и наблюдалось в эксперименте (рис.5)
(Более подробное описание указанного опыта, а также других экспериментов содержится
в книге, посвященной управлению структурой сообществ фитопланктона [10]). Таким образом, сформулированное свойство относительных численностей сообщества также может быть использовано при целенаправленном влиянии на структуру сообщества.
12
Работа поддержана Российским фондом фундаментальных исследований (гранты 0204-48085 и 03-04-06044).
13
СПИСОК ЛИТЕРАТУРЫ
1. Левич А.П., Алексеев В.Л., Никулин В.А. Математические аспекты вариационного
моделирования в экологии сообществ // Математическое моделирование. 1994. Т. 6. №
5. С. 55-76.
2. Levich A.P. Variational modelling theorems and algocoenoses functioning principles // Ecological Modelling. 2000. V. 131. № 2 - 3. P. 207-227.
3. Левич А.П. Теория множеств, язык теорий категорий и их применение в теоретической
биологии. Учебное пособие. М.: Изд-во Моск. ун-та, 1982.
4. Левич А.П., Алексеев В.Л., Рыбакова С.Ю. Оптимизация структуры экологических сообществ: модельный анализ // Биофизика. 1993. Т. 38. Вып. 5. С. 877-885.
5. Фурсова П.В., Милько Е.С., Опарина И.А., Максимов В.Н., Левич А.П. Определение
потребностей диссоциантов Pseudomonas aeruginosa в углероде, азоте и фосфоре по
результатам культивирования без пополнения запасов ресурсов // Микробиология.
2004 (в печати).
6. Фурсова П.В. Расчеты видовых обилий и областей лимитирования в вариационной модели экологического сообщества // Математическое моделирование. 2003б. Т.15. № 5.
С. 115-128.
7. Фурсова П.В., Милько Е.С., Опарина И.А., Левич А.П. Выявление компонентов питания, ограничивающих рост моно- и смешанных культур диссоциантов бактерий
Pseudomonas aeruginosa // Вестник Моск. ун-та. Серия биология. 2004 (в печати).
8. Фурсова П.В., Милько Е.С., Опарина И.А., Левич А.П. Подходы к управлению составом сообщества диссоциантов Pseudomonas aeruginosa: экспериментальные данные и
модельные расчеты // Биотехнология. 2004 (в печати).
9. Фурсова П.В. Потребности в глюкозе, нитратах и фосфатах и их вариации в анализе
смешанной культуры диссоциантов Pseudomonas aeruginosa // Изв. АН. Сер. биол.
2003а. № 1. С. 122 – 127.
10. Левич А.П., Максимов В.Н., Булгаков Н.Г. Экспериментальная и теоретическая экология фитопланктона: управление структурой и функциями сообществ. М.: Изд-во НИЛ,
1997. 188 с.
11. Левич А.П., Булгаков Н.Г. О возможности регулирования видовой структуры лабораторного альгоценоза // Изв. РАН. Серия биол. 1993. № 4. С. 569-578.
12
ON ECOLOGICAL COMMUNITY STRUCTURE REGULATION BY RESOURCES
CONCENTRATION RATIO IN ENVIRONMENT CHANGING
P.V. Fursova, A.P. Levich
Biological faculty, MSU, 119992, Moscow
The opportunity of community structure regulation, which is given by variational model of
ecological community, is investigated. For R-, S- and M- dissociants of bacteria Pseudomonas aeruginosa microorganisms relative abundances dependencies on the initial resources concentration ratio in environment are received. The theoretical opportunity of
community structure control is shown. The change of a dominant species under resources
concentration ratio change both on modelling accounts and in experiment is illustrated on
an example of Chlorococcales community.
Key words: community structure, community structure regulation, variational model,
Pseudomonas aeruginosa.
Download