CC BY
29
УДК 577.1:663.13:576.8
РАСЧЕТ МОРФОЛОГИЧЕСКОГО ПАРАМЕТРА В ПРОЦЕССЕ РОСТА ПОПУЛЯЦИИ МИКРООРГАНИЗМОВ
© А.А. Арзамасцев, Д.В. Слетков, И.В. Исаева
Arzamastsev A.A., Slctkov D.V., Isayeva I.V. Calculating the morphological parameter in the process of micro-organisms’ population growth. The article looks at the dependence of the estimated fractional dimensions (i.e. morphological parameter) of yeast culture representation on the main indices of its growth. Since it takes quite little time to calculate the fractional dimensions of a digital representation, the specified technology can be utilised to control the process of microorganisms growing.
Зарождение фрактальной геометрии связано с трудами Пуанкаре, Фату, Жюлиа, Кантора, Хаусдорфа и Мандельброта [I]. В настоящее время ее использование позволило сделать существенные шаги, связанные с компьютерным представлением форм, структур поверхностей, анимацией и др. [1, 2, 4]. Фрактальная геометрия используется для объяснения формообразования различных объектов в биологии, минералогии, физике твердого тела [3].
При работе с биологическими объектами часто требуется поставить в соответствие изображению такого объекта некоторое числовое значение. Подобная проблема имеет место, например, при управлении процессом культивирования микроорганизмов, когда по изображению популяции в поле микроскопа необходимо сделать выводы о ее морфологическом и физиологическом состояниях или при наблюдении за ростом и изменением формы у макрообъектов, таких как мхи, кустарники, кроны деревьев и т. д. Удобным математическим аппаратом для решения данной проблемы является подсчет фрактальной размерности изображений. Однако результат таких вычислений не во всех случаях совпадает с фрактальной размерностью реального объекта, поэтому в данной работе мы будем использовать название «морфологический параметр». На возможность использования подобного подхода при наблюдении за микроорганизмами указано в статьях [5, 6].
Целью данной работы явилось определение зависимости между морфологическим параметром изображений популяции дрожжей БассИаготусез се/ч?уша<? и основными показателями, наблюдаемыми в процессе их роста.
МАТЕРИАЛЫ И МЕТОДЫ
В качестве объекта исследования выбраны дрожжи БассИаготусеБ сегеУ1з1ае. По культуральным и морфологическим признакам используемые микроорганизмы относятся к семейству Енс1отусе1асеае.
В качестве питательной среды для роста микроорганизмов выбрано мелассное сусло, которое готовили следующим образом. Мелассу разводили в воде до концентрации сухих веществ 16-17 %, добавляли по
5 мл Н3Р04 и НС1 до pH = 4,8 на 10 л среды. Хорошо перемешивали, разливали в колбы объемом 5 л, стерилизовали в автоклаве 20 мин при 2 атм. Такая среда часто используется в микробиологической промышленности для производства дрожжей из отходов этанола. Засев дрожжей производили при температуре 30 °С.
В ходе эксперимента наблюдение вели за следующими показателями процесса: концентрациями клеток микроорганизмов, сухих веществ питательной среды, этанола и объемом двуокиси углерода, выделяемой культурой за некоторый промежуток времени (»10 мин). По ним рассчитывали физиологические показатели: скорости роста биомассы, потребления субстратов, выделения метаболитов и дыхания.
Концентрацию клеток оценивали с помощью подсчета количества клеток в камере Горяева (глубина 0,1 мм; площади малой и большой ячейки 1/400 и 1/25 мм2 соответственно) [7]. Число клеток рассчитывали по формуле: N= а-5-\0-b-1000, где /V- количество клеток в 1 мл; а - количество клеток, подсчитанных в 5-16 = 80 малых квадратах; 5 - коэффициент пересчета на 1 мм2; 10 - коэффициент пересчета на 1 мм1 (при глубине камеры в 0,1 мм); Ь - число, выражающее, во сколько раз разведена проба (обычно использовали разведение в 10 раз); 1000 - коэффициент перевода в см3 или в мл.
Определение содержания сухих веществ по плотности проводили ареометром при температуре 20 °С. Содержание спирта в бражке определяли методом перегонки, с измерением концентрации спирта спиртометром при температуре 20 °С. Объем выделившегося углекислого газа в результате брожения определяли методом вытеснения. Измерения производили через часовые интервалы в течение 10 часов. Скорость продукции углекислого газа вычисляли по формуле: 0 = = Vg / № ■ I), где <2 - скорость продукции углекислого газа; У% - объем выделившегося газа; Ух - объем ме-лассно-дрожжевого сусла; / - время измерения.
Технология получения и обработки изображений заключалась в следующем. Фотографирование объекта производили через оптический микроскоп (увеличение 140-150) с дополнительным набором фотографических колец один раз в час.
Полученные цветные фотографии сначала преобразовывали в монохромный формат, устраняли дефекты освещенности, контрастировали и выделяли контур. Для обработки использовали программу Ас1оЬе РЬоихяИор.
Расчет морфологического параметра размерности изображений (£)) проводили по точечному методу [8, 9] с минимальным размером покрытия (АА/ш), равным 3, и максимальным размером покрытия (ЬМах), равным 51. Алгоритм реализован в виде программы, которая позволяет производить расчет й для монохромных изображений (формат *.ВМР) [10]. Программа работает с изображениями любого разрешения. Время расчета морфологического параметра зависит от размера холста, сложности изображения, производительности компьютера и занимает обычно 5-10 с. Данный параметр коррелирует с морфологическими характеристиками биологических объектов, что дает право использовать его в дальнейшем анализе подобных систем [10].
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
На рис. 1 показано типичное изображение популяции дрожжей БассИаготусез сеге\ч.ч1ае, полученное непосредственно с микроскопа (рис. 1а) и после его обработки в программе Ас1оЬе РИоЮзИор (рис. 1Ь). Видно, что изображение, полученное с микроскопа, нуждается в коррекции освещенности, удалении дефектов съемки и контрастировании. Эти операции в той или иной степени должны быть применены ко всем снимкам, полученным в ходе экспериментов. Методика обработки изображений была адаптирована с помощью настройки различных фильтров таким образом, что это обеспечивало наименьшее искажение морфологических характеристик объекта.
Поскольку в работе использовали монохромное представление изображений, окончательная обработка каждого снимка заключалась в выделении контура
(рис. 1Ь). Из рис. 1 видно, что в процессе обработки морфологические характеристики (количество, взаимное расположение клеток, их размеры) остаются практически неизменными, что позволило использовать технологию для определения морфологического параметра в дальнейших экспериментах.
Предварительное тестирование методики определения морфологического параметра проводили с использованием микроорганизмов и техники съемки, описанной выше.
На рис. 2 показано изменение плотности популяции и ее морфологических характеристик в процессе культивирования. Видно, что в начальных стадиях популяция в основном состоит из одиночных и двойных почкующихся клеток (рис. 2а), в дальнейшем наблюдается увеличение количества клеток и образование больших клеточных блоков (рис. 2Ь, 2с).
а) Ь)
Рис. 1. Популяция дрожжей: (а) - первичное изображение; (Ь) - изображение, полученное в результате обработки. Увеличение х!50
а) Ь) с)
Рис. 2. Внешний вид популяции дрожжей Зассіїаготусея сегеуіяіае, выращиваемых на мелассном сусле: (а) - через час после начала эксперимента, (Ь) - через 8 часов, (с) - через 12 часов после начала эксперимента. Увеличение х 140
Осо .млДмлмин)
Рис. 3. Изменение основных показателей в процессе роста дрожжей: а) - концентрация клеток; Ь) - морфологический параметр изображений; с) - концентрация сухих веществ питательной среды; с!) - концентрация этанола; е) - скорость продукции СОг
В ходе эксперимента наблюдение вели за следующими показателями процесса: концентрациями клеток микроорганизмов, сухих веществ питательной среды, этанола и объемом двуокиси углерода, выделяемой культурой за промежуток времени (10 мин) (рис. 3). По ним рассчитывали физиологические показатели: скорости роста биомассы, потребления субстратов, выделения метаболитов и дыхания. На рис. 4 приведены графики временных зависимостей этих показателей, которые хорошо коррелируют с морфологическим параметром изображений (рис. ЗЬ, 4). Коэффициент корреляции во всех случаях более 0,8.
На рис. 5 показаны графики зависимостей морфологического параметра от основных показателей, наблюдаемых в процессе культивирования дрожжей. Значения коэффициентов корреляции с морфологическим параметром изображений в этом случае составляет от 0,7 до 0,9.
Существование подобных зависимостей имеет значительную практическую ценность, так как некоторые физиологические характеристики популяции (концентрация клеток и связанные с ней скорости роста биомассы, потребления субстратов, дыхания и выделения метаболитов; распределение клеток по размерам и связанный с этим показателем средний биологический возраст культуры) можно вычислить из анализа их изображений.
Рис. 4. Изменение основных физиологических показателей в процессе роста дрожжей, а) - скорость роста биомассы (кл/(млч)); Ь) - скорость потребления субстратов (%мас/ч); с) - скорость выделения метаболитов (%мас/ч)
•2 и
а) Ь)
уО о
1 / 1 • ^
о/ О
Уи
0.8 □ □ 0.8
. □ о
— -СГ □
□
0.6 0.6 *
2 4 6 N х 107, кл/мл 0.4 0.6 асс,2,мл/(мл мин)
1.2
С) d)
А о
А 1 9
А'ч. /о
о/
\ А * А 0.8 О о, у
А о
А 0.6
0.5 0.75 1 S. %масс. 2 3 4 Р, %масс
Рис. 5. Графики зависимостей морфологического параметра изображений от основных показателей, наблюдаемых в процессе роста дрожжей: а) - от концентрации клеток; Ь) - скорости выделения СОг; с) - концентрации сухих веществ питательной среды; d) - концентрации этанола
Поскольку подсчет морфологического параметра размерности цифрового изображения занимает незначительное время, такой способ позволяет использовать указанную технологию для управления процессами культивирования микроорганизмов.
ЛИТЕРАТУРА
1. Mandelbrot В.В. The Fractal Geometry of Nature. San Francisco: Freeman, 1982.
2. Пайтген X.-O., Рихтер П.Х. Красота фракталов. Образы комплексных динамических систем: Пер. с англ. М.: Мир, 1993.
3. Фракталы и прикладная синергетика / Тез. докл. Первого междис- 8.
циплинарного семинара. М., 1999.
4. Федср Е. Фракталы. М.: Мир. 1991. 9.
5. Рудакова Г.З., Божокин С.В. Диагностика состояния дрожжей по
измерению фрактальной симметрии структур газоразрядного све- 10.
чения воздуха, образующихся при воздействии высокого напряжения на жидкофазный биологический // Биофизика. 1991.
Вып. 2. С. 298-303.
6. Божокин С.В. Математическая модель морфологического строения фибов // Биофизика. 1996. Вып. 6. С. 1298-1300.
7. Инструкция по технологическому контролю спиртового производства. М.: Пищевая пром-сть. 1967.
Voss R.F. Random Fractals: Characterization and Measurement, Scaling Phenomena in Disordered Systems. N. Y.: Plenum Press, 1985. Кроновер P.M. Фракталы и хаос в динамических системах. Основы теории. М.: Постмаркет, 2000.
Слетков Д.В.. Арзамасцев А.А. Сравнение различных алгоритмов вычисления фрактальной размерности // Вестн. Тамбов, ун-та. Сер. Естеств. и техн. науки. Тамбов, 2003. Т. 8. Вып. 2. С. 282-285.
Поступила в редакцию 24 сентября 2003 г.
