Введение. При выборе сырья в системе переработки сои на пищевые цели, как правило, возникает вопрос о качестве зерна, включающем различные его показатели, начиная от физических параметров, определяющих пригодность зерна сои к переработке, и заканчивая химическим составом, обуславливающим пищевую ценность того или иного продукта. Одними из важнейших критериев оценки пищевой ценности соевого зерна являются: общее содержание белка (36–48 %) и аминокислот, особенно незаменимых, таких как лизин, треонин, триптофан; масла (от 16 до 27 %) и жироподобных веществ, обеспечивающих выполнение в организме человека различных физиологических и биохимических функций [1–3]. Ценность соевого масла определяют входящие в его состав жирные кислоты, среди которых на линолевую кислоту приходится более 50 %; олеиновую – 20–36; стеариновую – 4,5–7,0; линоленовую – 3–5; пальмитиновую – 2,5–6,0 %. Семена сои также содержат значительные количества фосфолипидов, витаминов (А, D, C, E, группы В), широкий спектр макро- и микроэлементов [4, 5].Ввиду большого сортового разнообразия и широкого спектра изучаемых признаков для их анализа и систематизации актуальное значение имеет использование многомерных статистических методов, в частности кластерного анализа. Методы кластеризации применяются в различных научно-технических областях с целью решения задач сегментации для объединения объектов в группы со схожими характеристиками, а также сжатия большого объема информации и построения научно обоснованных классификаций [6–10].Цель исследования – систематизация сортов сои амурской селекции по биохимическим показателям качества на основе кластерного анализа для рационального использования в производстве пищевых продуктов и добавок.Задачи: провести кластерный анализ сортов сои по критерию наибольшего сходства биохимических признаков; дать сравнительную оценку выделенных кластеров по степени выраженности признаков; определить группы сортов с комплексом качественных показателей.Объекты и методы. Исследование проводилось в 2018–2020 гг. в лаборатории переработки сельскохозяйственной продукции ФНЦ ВНИИ сои. В качестве объектов исследования были взяты 18 сортов сои селекции Всероссийского НИИ сои, выращенных в условиях Амурской области и включенных в Государственный реестр селекционных достижений РФ: Персона, Бонус, Сентябринка, Интрига, Чародейка, Топаз, Нега 1, Куханна, Невеста, Кружевница, Китросса, Статная, Умка, Пепелина, Золушка, Алена, Лебедушка, Журавушка [11]. Образцы для исследования получены из коллекционных питомников опытного поля лаборатории селекции и генетики сои ФНЦ ВНИИ сои. Отбор проб исследуемых образцов проводили в фазу полной спелости зерна в соответствии с ГОСТ 13586.3-2015.В ходе биохимических исследований зерна сои была сформирована база данных и определена совокупность основных признаков, характеризующих качество зерна, по которым осуществлялась группировка сортов. В качестве признаков для кластерного анализа было выбрано 13 показателей: белок, незаменимые аминокислоты (лизин, лейцин, изолейцин, валин, треонин, фенилаланин, метионин + цистеин), масло, жирные кислоты (олеиновая, линолевая, линоленовая), активность ингибиторов трипсина.Содержание белка, масла, аминокислотный состав и жирные кислоты определяли методом диффузного отражения света в инфракрасной области спектра на ИК-сканере модели FOSS NIR Systems 5000 (Дания). Для определения белка и масла в зерне сои использовали стандартные калибровочные уравнения фирмы FOSS Analytical A/S в соответствии с ГОСТ 32749-2014. Калибровочные уравнения для определения состава аминокислот и жирных кислот разрабатывались в 2009–2010 гг. в аналитической группе ВНИИ сои [5, 12]. Трипсинингибирующую активность (ТИА) измеряли с помощью спектрофотометра Cary 50 (VARIAN) в соответствии с ГОСТ 33427-2015.Статистическую обработку данных проводили с использованием пакета программ MS Excel и STATISTICA 10, используя метод кластерного анализа [13]. В качестве алгоритма кластеризации использовали метод Варда (мера близости объектов – Манхэттенское расстояние) и метод k-средних (5 кластеров, 10 итераций). Показатели с различиями в размерных единицах были приведены в единую систему методом стандартизации средних значений каждого признака в диапазоне от –2,7 до +2,7 ед. При выборе числа кластеров уровень сходства 15 % считали оптимальным. Проверку различий между кластерами по каждому признаку рассчитывали методом параметрического дисперсионного анализа.Результаты и их обсуждение. На основании средних стандартизированных значений биохимических признаков был проведен кластерный анализ методом Варда, результатом которого являлось построение дендрограммы в виде графической последовательности объединения и разделения сортов сои на кластеры (рис. 1).    Рис. 1. Дендрограмма сходства-различия сортов сои по комплексу биохимических признаков Анализ дендрограммы позволил сгруппировать сорта сои по принципу минимального взаимодействия внутри группы и максимального – между кластерами и сформировать пять основных кластеров, наиболее схожих по комплексу изучаемых признаков:1-й кластер – Лебедушка, Журавушка, Пепелина, Золушка, Алена;2-й кластер – Сентябринка, Интрига;3-й кластер – Куханна, Невеста, Кружевница, Китросса, Статная, Умка;4-й кластер – Персона, Бонус;5-й кластер – Чародейка, Топаз, Нега 1.Для получения более детальной информации о выделенных кластерах и оценки результатов биохимических параметров исследуемых сортов сои на следующем этапе исследований применяли метод k-средних. Оценку качества кластеризации наблюдений и проверку значимости различий кластеров по исследуемым характеристикам проводили посредством параметрического дисперсионного анализа, используя процедуру, сравнивающую внутрикластерные и межкластерные дисперсии по каждому признаку (табл. 1). Таблица 1Результаты дисперсионного анализа кластеризации наблюдений по биохимическим признакам (при р = 0,05) Признак Межгрупповая дисперсияЧисло степеней свободыВнутригрупповая дисперсияЧисло степеней свободыFУровень значимости pБелок13,7097643,290241313,542090,000144Лизин14,3859642,614041317,885830,000034Лейцин12,1247244,87528138,082690,001679Изолейцин13,5646543,435351312,832790,000189Валин13,8480343,151971314,278720,000110Фенилаланин14,4061142,593891318,050090,000032Треонин14,7862642,213741321,707750,000012Метионин + цистеин14,9457242,054281323,645040,000007Масло13,7574243,242581313,788890,000132Олеиновая 13,4896243,510381312,489040,000217Линолевая15,4461641,553841332,307100,000001Линоленовая14,6584642,341541320,345620,000017ТИА13,6534241,246581320,491100,000152  Результаты дисперсионного анализа указывают на высокое качество разбиения 18 наблюдений на 5 кластеров по каждому показателю, поскольку признаком качественной классификации является минимум изменчивости внутри кластера и максимум между кластерами.Для визуализации полученных результатов была реализована процедура построения усредненных графиков, демонстрирующих значимые отличия кластеров и степень выраженности биохимических признаков (рис. 2). Анализ графиков позволил провести градацию наблюдений в порядковой шкале стандартизированных значений, где в качестве уровня «средний +» – стандартизированный интервал (0,0; +0,7), «выше среднего» – (> +0,7), «средний –» – стандартизированный интервал (0,0; –0,7), «ниже среднего» –(< –0,7). Результаты, полученные методом k-средних, в полной мере соответствуют дендрограмме, построенной методом иерархической кластеризации.Основываясь на полученных результатах кластеризации, проведен сравнительный анализ выделенных кластеров по средним значениям биохимических признаков (табл. 2).    Рис. 2. Графики средних стандартизированных значений для каждого кластера Таблица 2Сравнительный анализ биохимических признаков сортов сои амурской селекции по кластерам (M±∆ при р = 0,05) ПризнакСреднее значение за три года, %Кластер 1Кластер 2Кластер 3Кластер 4Кластер 5Белок37,21±0,6638,10±0,4839,11±0,3841,10±0,1936,83±0,37Лизин6,07±0,296,04±0,225,92±0,605,73±0,826,27±0,11Лейцин8,28±0,657,91±0,078,14±0,787,83±0,218,20±0,36Изолейцин5,69±0,515,67±0,315,81±0,556,01±0,845,34±0,21Валин7,14±0,666,78±0,256,72±0,386,16±0,367,48±0,34Фенилаланин4,27±0,564,27±0,664,33±0,174,39±0,664,20±0,09Треонин3,44±0,413,32±0,343,35±0,443,20±0,683,40±0,07Метионин + цистеин1,73±0,221,75±0,021,61±0,521,46±0,481,94±0,50Масло18,89±0,6119,12±0,3918,06±0,3017,46±0,8620,23±0,22Олеиновая (ω9)18,74±0,3920,61±0,1121,65±0,3123,79±1,4518,5±0,36Линолевая (ω6)50,83±0,3650,91±0,0150,36±0,3350,05±0,4951,46±0,34Линоленовая (ω3)6,58±0,497,23±0,328,01±0,339,85±1,715,41±0,12ТИА, мг/г22,84±0,6223,57±1,1123,17±1,1623,14±0,9324,29±1,07  Первый кластер представлен группой сортов сои, отличающихся средними показателями белка (37,2 %) и масла (18,9 %), высокой концентрацией незаменимых аминокислот (лейцина, треонина) и минимальным уровнем трипсинигибирующей активности, являющимся одним из важных показателей в зерне сои пищевого направления.Сорта второго и третьего кластеров по изучаемым признакам значимо не отличались друг от друга и характеризовались средними значениями по всем показателям. При этом третий кластер образуют сорта со значениями признаков выше среднего по содержанию фенилаланина (4,3 %) в белке и олеиновой кислоты (21,7 %) в масле, а также значениями ниже среднего по признакам содержания метионина+цистеина (1,6 %) в белке, масла (18,1 %) и линолевой кислоты (50,4 %).Наибольшие различия по большинству признаков выявлены в четвертом и пятом кластерах. Сорта, сформировавшие четвертый кластер, отличались максимальными значениями признаков белка (41,1 %), изолейцина (6,0 %), фенилаланина (4,4 %), а также олеиновой и линоленовой кислот (соответственно 23,8 и 9,9 %). Помимо этого в сортах четвертого кластера отмечено наибольшее значение суммы полиненасыщенных жирных кислот в масле – 83,8 %, доля данных кислот в сортах сои других кластеров составила 75,4–80,1 %.Пятый кластер представлен высокомасличными сортами с повышенным содержанием линолевой (ω6) и минимальным – линоленовой (ω3) кислот, что делает эти сорта перспективными для производства масла на пищевые цели. Кроме того, представители этого кластера отличались наибольшим содержанием лизина, валина и метионина + цистеина. Поскольку сорта в этой группе являются высоколизиновыми, то их можно рекомендовать для использования в произ­водстве лизина в качестве пищевой добавки.Следует отметить, что, несмотря на межсортовые различия, представители всех кластеров имеют сбалансированный аминокислотный состав и оптимальное соотношение ω6 и ω3 жирных кислот (5,1–9,5:1), что соответствует нормам физиологической потребности взрослого человека [14].По белково-масличному составу зерна выделенные кластеры сортов сои условно классифицированы на три группы: универсальные (с содержанием белка – не менее 38 %, масла – не менее 18 %), масличные (с содержанием белка – менее 38 %, масла – более 18 %), высокобелковые (с содержанием белка – более 40 %, масла – менее 18 %) [15].Универсальную группу сформировали второй и третий кластеры, включающие сорта Сентябринка, Интрига, Куханна, Невеста, Кружевница, Китросса, Статная, Умка, которые могут быть использованы как для производства масла, так и белковых продуктов. Масличная группа объединила первый и пятый кластеры, включающие сорта: Лебедушка, Журавушка, Пепелина, Золушка, Алена, Чародейка, Топаз, Нега 1, – соответственно их лучше использовать в маслоэкстракционном производстве. В высокобелковую группу вошли сорта четвертого кластера (Чародейка, Топаз, Нега 1), они наиболее пригодны для производства соевой муки, белковых концентратов, изолятов, текстуратов, соевого молока, пищевых добавок и т. д.Заключение. Применение методов кластерного анализа позволило систематизировать сорта сои по степени выраженности биохимических параметров зерна, сформировать пять основных кластеров, наиболее схожих по комплексу изучаемых признаков, и провести их сравнительную оценку. Определены три группы сортов по критериям белково-масличного состава (универсальные, масличные, высокобелковые), что позволит обеспечить целенаправленный выбор сырья в производстве пищевых продуктов (добавок) и прогнозирование качества готовой продукции.