АППРОКСИМАЦИЯ ГЕНЕТИЧЕСКИХ КОРРЕЛЯЦИЙ ПРИ СЕЛЕКЦИИ СЕВЕРНЫХ ОЛЕНЕЙ (RANGIFER TARANDUS) ПО ЖИВОЙ МАССЕ И ПРОМЕРАМ С ПОМОЩЬЮ ФЕНОТИПИЧЕСКИХ КОРРЕЛЯЦИЙ
АППРОКСИМАЦИЯ ГЕНЕТИЧЕСКИХ КОРРЕЛЯЦИЙ ПРИ СЕЛЕКЦИИ СЕВЕРНЫХ ОЛЕНЕЙ (RANGIFER TARANDUS) ПО ЖИВОЙ МАССЕ И ПРОМЕРАМ С ПОМОЩЬЮ ФЕНОТИПИЧЕСКИХ КОРРЕЛЯЦИЙ
Аннотация
Цель исследований — оценка возможности использования фенотипических корреляций (r) между живой массой, промерами и индексами телосложения в качестве аппроксимации генетических корреляций у домашних северных оленей. Исследована выборка из 98 важенок (самок) уральского экотипа ненецкой породы, отобранных случайным образом, без пороков телосложения и имевших теленка. В результате проведенных исследований в выборке наблюдались положительные фенотипические корреляции между живой массой и «обхватом груди» (+0,63), между живой массой и «глубиной груди» (+0,60). Средние значимые величины фенотипических корреляций с живой массой показали промеры: «высота в холке» (+0,59) «косая длина туловища» (+0,58), «длина головы» (+0,56), «ширина в маклоках» (+0,55), «ширина груди» (+0,42), «высота в локте» (+0,37). Сопоставление отклонений генетических корреляционных коэффициентов от фенотипических показало, что значения статистически совпадают в 26 случаях из 36 при разности, не превышающей 0,05, максимальная разница составила 0,18 между коэффициентами промеров «глубина груди» и «длина туловища». Такое высокое совпадение величин генетических и фенотипических корреляций позволяет делать вывод о высокой степени аппроксимации показателей, что важно для практической селекции по живой массе и линейным размерам животных. Статистически значимые величины обнаружены у 30 фенотипических корреляций, у генетических — 26 из 36 исследованных. Сравнительный анализ между фенотипическими и генетическими корреляционными матрицами показал совпадение значимых величин 87%, что дает основание считать аппроксимацию убедительной. Фенотипические корреляции между индексами телосложения и живой массой показали положительную, но низкую степень корреляций +0,1–0,2. Предлагается использовать фенотипические корреляции в качестве приближённой оценки генетических корреляций при селекции северных оленей.
1. Введение
Северные олени — животные универсальной продуктивности: от них получают мясо, субпродукты, кожевенно-меховую продукцию, панты и рога, эндокринно-ферментное сырье. Олени используются для передвижения коренного населения в тундре и тайге, для перевозок грузов. В ряде районов практикуется доение важенок (отелившихся самок). Но главной продукцией северного оленеводства по-прежнему остается мясная, она дает основной доход хозяйствам всех форм собственности.
Наибольшая по численности ненецкая порода северных оленей выпасается в трех федеральных округах: Северо-Западном, Уральском и Сибирском и насчитывает 1 млн. 222 тыс. голов, что составляет 78% от всего российского поголовья домашних северных оленей . По данным отчета «Анализ рынка оленины в России», подготовленного BusinesStat в 2022 г, в 2021 г., производство мяса оленей в РФ выросло на 6,7%: с 9 до 9,6 тыс. т (в убойном весе). Хозяйства, выпасающие ненецкую породу оленей, произвели 77,1% товарной оленины .
Чтобы поддерживать и совершенствовать показатели мясной продуктивности северных оленей, в отрасли разрабатываются и внедряются инновационные приемы селекции, включая молекулярно-генетический анализ .
Микроэволюционные изменения в популяциях определяются взаимодействием наследственности (генов) и условий окружающей среды. Сочетание этих факторов, а также взаимное влияние формирует внешний вид и характеристики организма (фенотип). В животноводстве для понимания этих процессов и повышения эффективности селекции используют корреляционный анализ. Этот статистический метод позволяет оценить, как связаны между собой различные признаки. Фенотипическая корреляция, то есть наблюдаемая связь между признаками, складывается из генетической (наследственной) и негенетической (связанной с условиями среды) составляющих. Анализ генетической корреляции позволяет понять, как селекция влияет на разные признаки одновременно. Если признаки связаны положительно, то отбор по одному из них улучшает и другой. Однако отрицательная корреляция затрудняет одновременное улучшение обоих признаков.
В силу уникальности северного оленеводства и суровых условий среды, точно определить генетическую взаимосвязь между признаками у оленей, крайне затруднительно. Тем не менее, в научной литературе распространено мнение, что при соблюдении условий, наблюдаемые фенотипические корреляции могут служить индикатором генетических корреляций . Эта концепция, известная как гипотеза Шеверуда, предполагает тесную связь между генетическими и фенотипическими корреляциями . Эмпирические исследования на растениях, насекомых и животных подтверждают эту гипотезу , , , что дает возможность использовать фенотипические корреляции для оценки генетических.
Живая масса — показатель, определяющий мясную продуктивность сельскохозяйственных животных. Ранее, при изучении полезных хозяйственных характеристик пород и популяций северных оленей, уже применялся корреляционный анализ , .
В предыдущих исследованиях нами было установлено, что у северных оленей ненецкой породы наблюдается статистически значимая изменчивость, связанная с географическим положением. В Ямало-Ненецком автономном округе было проведено исследование, направленное на оценку возможности улучшения продуктивности домашних северных оленей путем селекции по живой массе. Анализ показал, что вес матери оказывает статистически значимое влияние на вес потомства в первые полгода жизни. Кроме того, вес оленей в 6 месяцев (r ≤ 0,6) коррелирует с весом в 1,5 года, что может быть использовано для раннего отбора перспективных животных. Однако эта связь ослабевает к возрасту 2,5 лет (r ≤ 0,5) , .
Внедрение мультиплексных панелей микросателлитов , , стало прогрессом в анализе генетической структуры популяций северного оленя, обеспечив точную оценку происхождения и генетического разнообразия. Фенотипическая оценка, несмотря на значимость в селекции, требует валидации генетическими методами. Подтверждение соответствия фенотипических корреляций генетическим позволит оптимизировать стратегию селекции, сократив число анализируемых признаков или скорректировав их выбор, без ущерба для селекционной эффективности.
Цель исследования — оценка возможности использования фенотипических корреляций (r) между живой массой, промерами и индексами телосложения в качестве аппроксимации генетических корреляций у домашних северных оленей.
2. Методы и принципы исследования
В 2020 году проведено зоотехническое обследование оленей, принадлежащих Ямальской опытной станции ТюмНЦ СО РАН. В рамках этого обследования для дальнейшего изучения случайным образом выбрали 98 здоровых важенок от 3 до 9 лет, без видимых отклонений в телосложении и имевших потомство. Оценка живой массы и измерения тела оленей проводились в соответствии с методическими рекомендациями, разработанными отделом оленеводства НИИСХ Крайнего Севера .
С точностью до 1 см измерялись следующие параметры:
– Мерной палкой: высота в холке и в локте, глубина и ширина груди (за лопатками).
– Мерной лентой: косая длина туловища и обхват груди.
– Мерным циркулем: ширина в маклоках и длина головы.
На основе полученных данных рассчитаны индексы телосложения: массивности, растянутости, сбитости и эйрисомии. Эти индексы позволяют оценить пропорции тела животных, мясную продуктивность.
Дополнительно у всех обмеренных животных был отобран биологический материал (пробы хрящевой ткани ушной раковины) для дальнейших генетических исследований.
Генетические исследования и математическая обработка проведена в лаборатории генетики животных Института общей генетики им. Н.И. Вавилова по 16 микросателлитным маркерам, в которых выявлено 119 аллелей. Число аллелей на локус варьировалось от 2 до 13 . Были вычислены фенотипические и генотипические линейные корреляции промеров и индексов по Пирсону.
Полученные данные обработаны статистически с использованием программного обеспечения MS Statistica 27.
3. Основные результаты
Фенотипические корреляции. Полученные значения корреляций были положительными, но различались по величине (силе). Положительные показатели обнаружены между «живой массой» и «обхватом груди» (+0,63), между «живой массой» и «глубиной груди» (+0,60), а также между «высотой в холке» (+0,60) и «высотой в локте» (+0,84), между «высотой в холке» и «глубиной груди» (+0,71), (табл. 1).
Таблица 1 - Фенотипические корреляции между живой массой и размерами тела, а также статистическая значимость (p-значения), были определены у домашних северных оленей
Фенотипические корреляции | Живая масса | Высота в холке | Высота в локте | Глубина груди | Обхват груди | Ширина груди | Ширина в маклоках | Косая дл. туловища | Длина головы |
Живая масса | 1 |
| |||||||
Высота в холке | 0,59* | 1 |
| ||||||
Высота в локте | 0,37* | 0,84** | 1 |
| |||||
Глубина груди | 0,6** | 0,71** | 0,27* | 1 |
| ||||
Обхват груди | 0,63** | 0,46* | 0,24 | 0,55* | 1 |
| |||
Ширина груди | 0,42* | 0,27* | 0,13 | 0,31* | 0,55* | 1 |
| ||
Ширина в маклоках | 0,55* | 0,36* | 0,24 | 0,31* | 0,37* | 0,46* | 1 |
| |
Косая дл. туловища | 0,58* | 0,45* | 0,31* | 0,39* | 0,23 | 0,1 | 0,31* | 1 |
|
Длина головы | 0,56* | 0,37* | 0,26 | 0,35* | 0,37* | 0,29* | 0,37* | 0,39* | 1 |
Примечание: ** p<0,01; *p <0,05
Статистически значимыми показатели были у 30 фенотипических корреляций. Максимально значимые величины фенотипических корреляций с живой массой показали промеры: «высота в холке» (+0,59), «косая длина туловища» (+0,58), «длина головы» (+0,56), «ширина в маклоках» (+0,55), «ширина груди» (+0,42), «высота в локте» (+0,37).
Генетические корреляции. Для определения генетических корреляций между признаками использовалась линейная регрессионная модель, в которой каждый анализируемый признак разбивается на две компоненты: генетическую (Gβ) и средовую (e). В модели Y представляет собой вектор значений признака, X — матрица ковариат, состоящая из двух столбцов, где первый столбец содержит вектор единиц, а второй — вектор возраста. G — это матрица генотипов с элементами 0, 1 и 2, отражающими количество определённого аллеля в конкретном локусе у данного животного. Вектор фиксированных эффектов ковариат (интерсепт и возраст) обозначен как β — вектор фиксированных генетических эффектов на признак, а e — вектор случайных остатков модели. Предполагается, что регрессионные остатки имеют нормальное распределение с параметрами, где 2 — это дисперсия признака.
Генетические корреляции (Gβ) между генетическими компонентами каждой пары изучаемых признаков рассчитаны и представлены в таблице 2. Согласованность между фенотипическими и генетическими корреляционными матрицами оценивалась посредством ранговой корреляции Спирмена. Полученные результаты свидетельствуют о соответствии между фенотипическими и генотипическими корреляциями (ρ = 0,91, p <2,2 × 10^-16).
Таблица 2 - Оценка генетических корреляций между живой массой и морфометрическими признаками у домашних северных оленей
Генетические корреляции | Живая масса | Высота в холке | Высота в локте | Глубина груди | Обхват груди | Ширина груди | Ширина в маклоках | Косая дл. туловища | Длина головы |
Живая масса | 1 |
| |||||||
Высота в холке | 0,54* | 1 |
| ||||||
Высота в локте | 0,31* | 0,83** | 1 |
| |||||
Глубина груди | 0,52* | 0,71** | 0,26 | 1 |
| ||||
Обхват груди | 0,57* | 0,47* | 0,23 | 0,57* | 1 |
| |||
Ширина груди | 0,45* | 0,31* | 0,13 | 0,38* | 0,61* | 1 |
| ||
Ширина в маклоках | 0,46* | 0,36* | 0,23 | 0,25 | 0,27* | 0,52* | 1 |
| |
Косая дл. туловища | 0,54* | 0,41* | 0,32* | 0,23 | 0,23 | 0,11 | 0,21 | 1 |
|
Длина головы | 0,58* | 0,32* | 0,23 | 0,28* | 0,43* | 0,36* | 0,38* | 0,38* | 1 |
Примечание: **p <0,01; *p <0,05
Статистически значимыми были 26 генетических корреляций.
Сопоставление отклонений генетических корреляционных коэффициентов от фенотипических (Таблица 3) показало, что значения статистически совпадают в 26 случаях из 36 при разности, не превышающей 0,05, разница составила 0,18 между коэффициентами промеров «глубина груди» и «длина туловища». Такое совпадение величин генетических и фенотипических корреляций позволяет делать вывод о значительной величине аппроксимации показателей, что важно для практической селекции по живой массе и линейным размерам животных.
Таблица 3 - Значения отклонений генетических корреляционных коэффициентов от фенотипических
Признак | Живая масса | Высота в холке | Высота в локте | Глубина груди | Обхват груди | Ширина груди | Ширина в маклоках | Косая длина туловища | Длина головы |
Живая масса | 1 | 0,06 | –0,06 | –0,08 | –0,06 | 0,04 | –0,09 | –0,04 | 0,02 |
Высота в холке |
| 1 | –0,01 | 0,00 | 0,02 | 0,04 | 0,00 | –0,04 | –0,06 |
Высота в локте |
|
| 1 | –0,01 | –0,02 | 0,00 | –0,01 | 0,01 | –0,02 |
Глубина груди |
|
|
| 1 | 0,02 | 0,06 | –0,06 | –0,16 | –0,07 |
Обхват груди |
|
|
|
| 1 | 0,06 | –0,09 | 0,00 | 0,06 |
Ширина груди |
|
|
|
|
| 1 | 0,06 | 0,02 | 0,07 |
Ширина в маклоках |
|
|
|
|
|
| 1 | –0,10 | 0,01 |
Косая длина туловища |
|
|
|
|
|
|
| 1 | –0,02 |
Длина головы |
|
|
|
|
|
|
|
| 1 |
Фенотипические корреляции индексов, как и промеров, дают многомерную цифровую модель домашнего животного, которая может быть использована в селекции и планировании хозяйственного использования. В таблице 4 представлены средние показатели индексов важенок ненецкой породы, а также фенотипические корреляции индексов с живой массой и между собой. Были получены положительные корреляции исследованных индексов с живой массой, однако величина их была на низком уровне: от +0,10 до +0,20.
Таблица 4 - Живая масса, индексы телосложения и их фенотипические корреляции у важенок ненецкой породы
Показатель | Живая масса, кг | Индексы телосложения, % | |||
Массивности | Растянутости | Сбитости | Эйрисомии | ||
82,7±0,78 | 121,7±0,49 | 107,7±0,42 | 113,0±0,57 | 24,6±0,16 | |
+0,12 | +0,12 | +0,1 | +0,2 | ||
И. массивности | +0,12 | 1 | +0,21* | +0,64** | +0,43** |
И. растянутости | +0,12 | +0,21* | 1 | -0,62** | -0,17 |
И. сбитости | +0,10 | +0,64** | -0,62** | 1 | +0,48** |
И. эйрисомии | +0,20 | +0,43** | -0,17 | +0,48** | 1 |
Примечание: **p <0,01; *p <0,05
Индекс массивности связан с другими индексами более тесно: связь с и. растянутости составила +0,21 при первом пороге достоверности, с И. сбитости +0,64 (р <0,01), с И. эйрисомии +0,43 ((р <0,01). И. сбитости показал достоверную фенотипическую связь с И. эйрисомии +0,48 (р <0,01). И. растянутости и И. сбитости показали достоверно отрицательную связь – -0,62 ((р <0,01), что закономерно, поскольку индексы отражают противоположные характеристики телосложения. Отрицательной в невысокой степени -0,17 оказалась связь И. эйрисомии и И. растянутости.
4. Заключение
Несмотря на развитие генетических методов, оценка животных по внешним признакам (фенотипу) остается частью селекции. Однако для повышения точности, результаты фенотипической оценки необходимо подтверждать генетическим анализом.
У домашних северных оленей обнаружена сильная взаимосвязь между живой массой и обхватом груди (r = +0,63), а также глубиной груди (r = +0,60). Также выявлена умеренная, но значимая связь живой массы с высотой в холке (+0,59), косой длиной туловища (+0,58), длиной головы (+0,56), шириной в маклоках (+0,55) и шириной груди (+0,42).
Генетический анализ с использованием микросателлитных маркеров подтвердил эти взаимосвязи, показав сходные результаты с фенотипическими корреляциями в 87% случаев. Сравнение фенотипических и генетических корреляционных матриц выявило высокую степень соответствия между ними. Взаимосвязь между индексами телосложения и живой массой оказалась положительной, но не выраженной.
В селекции северных оленей, направленной на увеличение живой массы и размеров тела, рекомендуется использовать фенотипические корреляции в качестве приближенной оценки генетических взаимосвязей.
