Введение. Глобальная тенденция к индустриализации сельского хозяйства сопряжена с рядом рисков. Одним из таких рисков является уменьшение национальных генетических ресурсов, или генофондов животных [1].Анализ полиморфных вариантов генов, отвечающих за молочную продуктивность, предоставляет возможность оценить уровень вариации аллельных форм в стадах крупного рогатого скота и определить степень гетерозиготности. Сохранение этих генов и их низкая степень эволюции у крупного рогатого скота по сравнению с другими плацентарными видами предоставляет возможность использовать их в качестве ДНК-маркеров в селекционном процессе [2–4]. Это особенно актуально при отборе быков-производителей и маточного поголовья по показателям продуктивности. В дальнейшем это позволит сформировать микропопуляции, способные производить продукцию с необходимыми технологическими свойствами [5, 6]. Изучение генетического полиморфизма открывает перспективы для вычисления генетических расстояний и составления генетических карт. Анализ аллельных вариантов генов, ответственных за белковую молочность, также способствует решению задачи сохранения чистоты отечественных молочных пород, особенно на фоне активного процесса голштинизации скота, наблюдаемого в последние десятилетия.Улучшение отечественных пород сельскохозяйственных животных с целью увеличения объемов производства и удовлетворения потребностей населения зависит от наличия генетического разнообразия популяций, а также от применения современных методов селекции. Использование метода ДНК-маркирования способствует выявлению генетических характеристик и ускорению селекционного процесса за счет отбора животных с предпочтительными генотипами [7–9].Рост объемов валового производства молока достигается благодаря увеличению молочной продуктивности коров, что представляет собой главную тенденцию развития отрасли молочного скотоводства в Российской Федерации [10]. Сокращение поголовья пород отечественной селекции вызывает опасения, так как ценные генетические ресурсы этих животных заключают в себе высокую адаптивность к разным климатическим условиям, равно как и связь с историческими и культурными традициями российских регионов [11–13].Цель исследования – изучение влияния полиморфизма генов на молочную продуктивность и производственную типичность крупного рогатого скота ярославской породы.Задачи: изучить распределение генотипов и генетическое разнообразие по генам молочной продуктивности в популяции крупного рогатого скота ярославской породы; оценить молочную продуктивность коров в зависимости от генотипа по генам: синтаза жирных кислот (FASN), тиреоглобулин 5 (TG5), гипофизарный фактор транскрипции 1 (Pit 1), диацилглицерол О-ацилтрансфераза 1 (DGAT 1).Объекты и методы. Научно-производственный опыт проводился в племенных стадах Ярославской области в период 2024–2025 гг. Объект исследования – чистопородные коровы-первотелки ярославской породы в количестве 229 голов. Материалом для исследования послужили показатели молочной продуктивности (удой за 305 дней, содержание жира и белка, количество молочного жира и белка), промеры экстерьера (2–4-й месяц первой лактации), биологический материал (цельная кровь).Исследования по изучению полиморфизма генов проводились в лаборатории генетики и биотехнологии Ярославского НИИЖК – филиала ФНЦ «ВИК им. В.Р. Вильямса».Для анализа генетического разнообразия крупного рогатого скота изучен полиморфизм генов молочной продуктивности: синтаза жирных кислот (FASN), тиреоглобулин 5 (TG5), гипофизарный фактор транскрипции 1 (Pit 1), диацилглицерол О-ацилтрансфераза 1 (DGAT 1).Геномную ДНК выделяли из цельной крови, используя наборы «ДНК-Экстран-1» (ООО «Синтол», Россия) согласно инструкциям производителя. Полимеразно-цепную реакцию (ПЦР) проводили на термоциклере C1000 Touch (ООО «Био-Рад Лаборатории», Россия). При проведении ПЦР использовали готовую инкубационную смесь 5X ScreenMix-HS (ЗАО «Евроген») с добавлением 0,2 мкМ каждого праймера (ООО «Синтол», Россия).Праймерные пары для проведения амплификации указаны в таблице 1. Таблица 1Праймеры для амплификации участков геновPrimers for amplification of gene regions Наименование генаПоследовательность праймераСинтаза жирных кислот (FASN)F: 5′– AGAGCTGACGGACTCCACAC – 3′R: 5′– GCCGATGCACTCGATGTAG – 3′Тиреоглобулин (TG5)F: 5′ – GGGGATGACTACGAGTATGACTG – 3′R: 5′ – GTGAAAATCTTGTGGAGGCTGTA – 3′.Гипофизарный фактор транскрипции, Pit1 / HinfIF: 5′ – CAATGAGAAAGTTGGTGC – 3′R: 5′ – TCTGCATTCGAGATGCTC – 3′АцилСоАдиацилглицеролацилтрансфераза 1 (DGAT1)F: 5' – TGCTGGCCCTGATGGTCTACAC – 3'R: 5' – GAAGGAAGCAAGCGGACAGT – 3'  Температура отжига праймеров, концентрация ДНК подбирались эмпирически на основе теоретических расчетов. Для проверки и поиска оптимальных температур и времени этапов амплификации проводилась множественная градиентная ПЦР. Диапазон температур для подбора температуры на этапе отжига праймеров составил от 58 до 68 °C, разведения ДНК были кратны 0,5; 1; 2.Полученные продукты амплификации обрабатывали при необходимости эндонуклеазами рестрикции (ООО «СибЭнзайм», Россия). Срок инкубации составлял не менее 16 ч при температуре 37 ºС в сушильном шкафу Binder ED 115 (BINDER GmbH, Германия). Разделение рестрикционных фрагментов проводили с помощью электрофореза в 2 %-м агарозном геле (ООО «Компания Хеликон», Россия) окрашенным бромистым этидием (ООО «Компания Хеликон», Россия) в электрофорезной горизонтальной камере SubCell GT (ООО «Био-Рад Лаборатории», Россия) при напряжении 120 В с помощью источника питания «Эльф-4» («ООО ДНК-Технология», Россия) в течение 80 мин.Результаты электрофореза анализировали с помощью гель-документирующей системы MolecularImagerGelDocXR (ООО «Био-Рад Лаборатории», Россия).Для оценки генетического разнообразия рассчитаны показатели на основе формул, представленных в работах Л.А. Животовского (1991) [14] и В.М. Кузнецова (2014) [15].Количество аллелей – Na.Количество эффективных аллелей:                              Ne=1Σpi2  ,                         (1) где pi – частота i-й аллели.  Информационный индекс Шеннона                    I= -1∙ pi∙ lnpi  ,                (2) где pi – частота i-й аллели.Наблюдаемая гетерозиготность                              Ho=NгетN  ,                          (3) где N – общее количество особей; Nгет – количество особей-гетерозигот.Ожидаемая гетерозиготность                    He=1- pi2  ,                    (4) где pi– частота i-й аллели.Индекс фиксации                            Fis=(He-Ho)He  ,                       (5) где Ho – наблюдаемая гетерозиготность; He – ожидаемая гетерозиготность.Число степеней свободы                  DF=r-1∙c-1  ,              (6) Критерий Хи-квадрат                  X2= i=1rj=1c(Oij- Eij)2Eij  ,             (7) где i – номер строки; j – номер столбца; r – общее количество строк; c – общее количество столбцов; Oij – фактическое количество наблюдений в ячейке ij; Eji – ожидаемое число наблюдений в ячейке ij.Для оценки статистической значимости различий между эмпирическими и теоретическими частотами встречаемости генотипов применяли критерий Хи-квадрат (χ2) Пирсона [16]. Теоретические частоты встречаемости рассчитывали на основе закона Харди – Вайнберга [14]. Для выявления частных различий эмпирического и теоретического распределений применялся z-тест. Критический уровень значимости принимали равным 0,05 (α = 0,05). Для обработки и статистического анализа данных использовали программное обеспечение: Python 3.9.7, MS Excel, GenAlEx 6.51 [16].Количество молока, полученного на 100 кг живой массы подопытных коров, – коэффициент молочности (КМ) – определяли по формуле                         КМ = УЖМ∙ 100  ,                     (8) где У – удой коровы за 305 дней первой лактации, кг; ЖМ – живая масса по первой лактации, кг.Коэффициент производственной типичности (КПТ) определяли по формуле, предложенной Ничиком Б. (1987) [17]:                              КПТ = У∙ИДВ∙ИС  ,                         (9) где У – удой за 305 дней первой лактации, кг; ИД – индекс длинноногости, %; В – живая масса по первой лактации, кг; ИС – индекс сбитости, %.Производственную типичность определяли по значению коэффициента следующим образом: КПТ больше 3 – молочное направление; КПТ от 2 до 3 – молочно-мясное направление; КПТ от 1 до 2 – мясо-молочное направление.Промеры экстерьера определялись по стандартным методикам, биометрическая обработка данных проводилась с использованием «пакета анализа», встроенного в MS Excel, и программы Statistica 10 по методикам Н.А. Плохинского [18].Результаты и их обсуждение. В результате исследований установлено, что во всех генах наблюдается аллельный полиморфизм, кроме диацилглицерол О-ацилтрансферазы 1, для которого выявлен только один аллель и гомозиготный генотип DGATKK. Аналогичные результаты по отсутствию полиморфизма в гене (DGAT1) установили на коровах белорусской черно-пестрой породы А.Н. Михалюк и Л.А. Танана (2023) [19].Для гена синтазы жирных кислот наблюдается значительное отклонение от ожидаемых значений: гомозиготный генотип FASNAA встречается реже (4,21 %), чем ожидаемо в равновесном распределении (7,93 %), в то время как гетерозиготный генотип FASNAG и гомозиготный генотип FASNGG встречаются чаще ожидаемого (47,89 % к 40,46 % и 47,89 % к 51,61 % соответственно).Для гена гипофизарного фактора транскрипции отмечается преобладание гомозиготного генотипа PIT 1/Hinf 1AA (60,56 % к 55,72 %), при этом гетерозиготный генотип PIT 1/Hinf 1AB встречается реже ожидаемого (28,17 % к 37,85 %). Гомозиготный генотип PIT 1/Hinf 1BB также демонстрирует значительное отклонение, встречаясь чаще ожидаемого (11,27 % к 6,43 %).В случае гена тиреоглобулина наблюдается преобладание гетерозиготного генотипа TG5CT относительно ожидаемого распределения (61,57 % против 47,62 %). Гомозиготный генотип TG5CC встречается реже ожидаемого (30,13 % против 37,11 %), а гомозиготный генотип TG5TT также демонстрирует значительное отклонение, встречаясь реже (8,30 % против 15,27 %).По-видимому, такое распределение частот генотипов характерно для ярославской породы, поскольку, изучая распределение полиморфных генотипов в популяции татарстанского типа холмогорской породы, Ю.Р. Юльметьевой, Ш.К. Шакировым (2020) было установлено преобладание гомозиготного генотипа TG5CС, частота которого составила 63,5 %. Количество гетерозигот составило 32,1 %, а гомозиготный генотип TG5ТT имел наименьшую частоту встречаемости – 4,4 % [20].Анализ данных показывает, что наиболее выраженные отклонения от ожидаемых частот наблюдаются для генотипов TG5, PIT 1/Hinf 1 и FASN, это может свидетельствовать о наличии селективного давления на данные аллели в исследуемой популяции. Полученные данные подтверждаются наличием статистически значимых различий между теоретически ожидаемым по закону Харди – Вайнберга и наблюдаемым распределением частот генотипов (табл. 2).Ген диацилглицерол О-ацилтрансферазы 1 (DAGAT1) характеризуется полным отсутствием генетического разнообразия, что подтверждается значениями числа аллелей (Na = 1,00), эффективного числа аллелей (Ne = 1,00), информационным индексом Шеннона (I = 0,00) и наблюдаемой гетерозиготностью (Ho = 0,00), указывающими на наличие только одной аллели и мономорфность гена в исследуемой популяции (табл. 3).  Таблица 2Анализ соответствия распределения частот генотипов теоретически ожидаемомуравновесному распределению Харди – ВайнбергаAnalysis of the correspondence of the genotype frequency distributionto the theoretically expected Hardy-Weinberg equilibrium distribution ГенЧисло степеней свободы DFЗначение критерия Хи-квадрат χ2Достигнутый уровень значимости p-valueГен диацилглицеролО-ацилтрансферазы 1–––Синтаза жирных кислот16,4190,011Гипофизарный фактор транскрипции19,2890,002Тиреоглобулин119,6710,000  Ген синтазы жирных кислот (FASN) имеет умеренное генетическое разнообразие с двумя аллелями, эффективным числом аллелей 1,68 и информационным индексом Шеннона, равным 0,59. Наблюдаемая гетерозиготность (Ho = 0,48) незначительно превышает ожидаемую (He = 0,40), что подтверждается отрицательным значением индекса фиксации (Fis = –0,18) (табл. 3).Гипофизарный фактор транскрипции (PIT1) также имеет два аллеля, но более низкие эффективное число аллелей (Ne = 1,61) и информационный индекс Шеннона (I = 0,57). При этом наблюдается существенное отклонение наблюдаемой гетерозиготности (Ho = 0,28) от ожидаемой (He = 0,38), что подтверждается положительным значением индекса фиксации (Fis = 0,26) и говорит о дефиците гетерозигот по этому гену (см. табл. 3).Наиболее высоким генетическим разнообразием характеризуется ген тиреоглобулина (TG5), имеющий два аллеля, эффективное число аллелей 1,91 и информационный индекс Шеннона I = 0,67. Наблюдаемая гетерозиготность (Ho = 0,62) значительно превышает ожидаемую (He = 0,48), что подтверждается отрицательным значением индекса фиксации (Fis = –0,29) (см. табл. 3).  Таблица 3Показатели генетического разнообразия популяции крупного рогатого скотаярославской породыIndicators of the genetic diversity of the Yaroslavl cattle population ПоказательГен диацилглицеролО-ацилтрансферазы 1Синтазажирных кислотГипофизарныйфактор транскрипцииТиреоглобулинОбъем выборки N70190142229Число аллелей Na1,002,002,002,00Эффективное числоаллелей Ne1,001,681,611,91Информационный индекс Шеннона I0,000,590,570,67Наблюдаемаягетерозиготность Ho0,000,480,280,62Ожидаемаягетерозиготность He0,000,400,380,48Индекс фиксации Fis–-0,180,26-0,29  Оценка взаимосвязи генотипов по изученным генам с молочной продуктивностью коров показала: животные с генотипом FASNAG достоверно превосходили гомозиготных животных с генотипом FASNGG по надою за 305 дней первой лактации на 639 кг (p ≤ 0,001) молока (табл. 4). По содержанию молочного жира достоверное превосходство также выявлено у первотелок с генотипом FASNAG над животными с генотипом FASNAA 0,29 % (p ≤ 0,05): оно составило 4,54 %.  Таблица 4Молочная продуктивность коров по первой лактации в зависимости от генотипа (M±m)Milk productivity of cows after the first lactation, depending on the genotype (M±m) ПризнакFASNTG5Pit 1/HinfAAAGGGTTCTCCAAABBBУдой за 305 дней, кг5700±517,35699±144,9***5060±115,05758±332,75549±114,7*5172±103,55941±157,95816±269,75437±343,9МДЖ, %4,25±0,144,54±0,06*4,48±0,064,36±0,104,46±0,044,37±0,074,41±0,044,34±0,064,48±0,08Количествомолочного жира, кг243,4±24,7257,7±6,7***226,2±24,7250,5±14,6247,1±5,4*226,3±6,0262,0±7,4255,5±12,7243,2±15,5МДБ, %3,40±0,06*3,27±0,023,27±0,023,31±0,063,32±0,02***3,20±0,023,26±0,023,33±0,043,32±0,05Количествомолочного белка, кг194,3±18,2187,4±5,4**166,4±4,3191,1±11,7*184,5±4,2***165,7±3,6194,3±5,7197,9±10,9181,6±12,7Живая масса, кг519±21,5499±5,4***476±3,4477±8,7494±4,1***474±4,3501±5,3509±10,7487±15,7Коэффициентмолочности, кг1094,0±92,11219±32,7**1104,0±28,21376,0±133,71169,0±25,01175,0±32,61268±39,2*1170,0±54,71123,0±58,4  По комплексному выходу молочного жира и белка на 100 кг живой массы также особи с генотипом FASNAG превосходили показатели гомозигот с генотипом FASNGG на 7,5 кг (p ≤ 0,01) (см. табл. 4).Особи с генотипом FASNAG произвели больше молока на 100 кг живой массы по сравнению с животными с генотипами FASNАА и FASNGG на 125 кг и 115 кг соответственно, и КМ составил 1219,0 кг молока.По производственной типичности коровы всех генотипов FASN относились к молочному типу (КПТ > 3,0), однако особи с генотипом FASNAG отличались более молочным типом (p ≤ 0,01). Показатель КПТ составил у FASNAG 4,3, FASNAА – 4,3; FASNGG – 3,9 (рис. 1).    Рис. 1. Расчет различий средних методом LSMeans по гену FASNCalculation of differences in averages using the LS Means method for the FASN gene  Установлена статистически значимая разница между показателями коэффициента производственной типичности у особей с генотипами FASNAG и FASNGG (p = 0,0018) (рис. 1).При оценке продуктивных показателей коров, генотипированных по гену TG5, была выявлена тенденция превосходства особей с гетерозиготным генотипом TG5СТ по качественным показателям молока (см. табл. 4). Так, по содержанию жира показатель гетерозигот равнялся 4,46 %, превосходство над гомозиготами TG5ТТ и TG5СС составило 0,10 и 0,09 абс. % соответственно. По содержанию массовой доли белка животные с генотипом TG5СТ превосходили показатели гомозигот TG5СС на 0,12 абс. % (p ≤ 0,01).Изучая взаимосвязь продуктивных признаков с полиморфизмом гена TG5, Е.О. Крупин и Ш.К. Шакиров (2019) выявили существенные различия по содержанию массовой доли жира в молоке у коров с генотипами CC и TC (на 0,29 %) [21].По удою за 305 дней первой лактации, количеству молока на 100 кг живой массы гомозиготы TG5ТТ превосходили особей с генотипами TG5СТ и TG5СС. Это подтверждается и значением коэффициента производственной типичности, первотелки с генотипом TG5ТТ имели более молочный тип (КПТ = 4,6), однако разница не была статистически значима (р = 0,3472) (рис. 2).У животных, генотипированных по гену Pit 1/Hinf, достоверных различий по показателям молочной продуктивности установить не удалось. Однако проведенный анализ показал, что особи с генотипом Pit 1АА превосходили показатели животных с генотипами Pit 1АВ и Pit 1ВВ по удою за первую лактацию на 125 и 504 кг молока соответственно, по общему выходу молочного жира и белка – на 2,9 и 31,6 кг соответственно, по общему выходу молочного жира и белка на 100 кг живой массы – на 3,4 и 8,5 кг соответственно. Также от коров с гомозиготным генотипом Pit 1АА получено больше молока на 100 кг живой массы (КМ = 1268,0 кг) по сравнению с особями с генотипом Pit 1АВ (КМ = 1170,0 кг) и Pit 1ВВ (КМ = 1123,0 кг), т. е. они имели лучшую конверсию корма.   Рис. 2. Расчет различий средних методом LSMeans по гену TG5Calculation of differences in averages by the LS Means method for the TG5 gene  Аналогичную зависимость молочной продуктивности от генотипа Pit 1АА установили М.В. Позовникова и Г.Н. Сердюк (2017) у помесных коров черно-пестрой и голштинской пород Ленинградской популяции [22].Также гомозиготы Pit 1АА имели более молочный тип (КПТ = 4,5) по сравнению со сверстницами с генотипами Pit 1АВ (КПТ = 4,3) и Pit 1ВВ (КПТ = 4,1) (рАА–АВ = 0,1534; рАА–ВВ = 0,1231) (рис. 3).Т.Ф. Лефлер, И.В. Крашенинникова (2022) сообщают о влиянии генотипа коров на молочную продуктивность и коэффициент молочности, который находился в пределах 1064,6–1261,7, что выше общепринятого стандарта (800) на 133,1–157,7 %, это свидетельствует о соответствии коров молочному типу [23].    Рис. 3. Расчет различий средних методом LSMeans по гену Pit 1Calculation of differences in averages by the LS Means method for the Pit 1 gene Заключение. Исследование генетического разнообразия крупного рогатого скота ярославской породы выявило существенные различия между анализируемыми генами. Наиболее важным результатом является полное отсутствие полиморфизма в гене DAGAT1, что может свидетельствовать о его высокой консервативности или действии направленного стабилизирующего селекционного отбора. Ген тиреоглобулина проявляет максимальное генетическое разнообразие (Fis = –0,29), что указывает на преобладание гетерозигот и аутбредный способ подбора родительских пар. Умеренное разнообразие FASN (Fis = –0,18) указывает на его значимость в метаболических процессах, а особенности распределения PIT1 (Fis = 0,26) могут быть связаны с его функцией в регуляции транскрипции. Полученные данные имеют важное значение для понимания генетической структуры породы и могут быть использованы при планировании селекционной работы.На основании оценки молочной продуктивности коров по первой лактации установлено, что наибольшее статистически значимое влияние на удой молока оказывает аллельный полиморфизм гена FASNAG. От них получено на 639 кг молока больше (p ≤ 0,001) относительно гомозигот FASNGG. Коровы с генотипом TG5СТ имели больший удой по сравнению со сверстницами с генотипом TG5СС на 377 кг молока (p ≤ 0,05).Качественные показатели молока в исследуемых стадах были в наибольшей степени обусловлены полиморфизмом генов FASN и TG5: первотелки с генотипом FASNAG превосходили по массовой доле жира животных с генотипом FASNAA на 0,29 % (p ≤ 0,05); животные с генотипом TG5СТ превосходили показатели гомозигот TG5СС по массовой доле белка на 0,12 % (p ≤ 0,001).Наибольший коэффициент молочности установлен у особей с генотипом FASNAG и Pit 1АА – 1219 и 1268 кг молока на 100 кг живой массыЪ соответственно. Согласно значению коэффициента производственной типичности коров с генотипами FASNAGTG5ТТ и Pit 1АА можно отнести к высокомолочному типу.