Комплексная оценка генотипа романовской породы овец в управлении селекционным процессом племенных стад
Автор: Костылев Михаил Николаевич, Косяченко Николай Михайлович, Абрамова Марина Владимировна, Барышева Мария Сергеевна
Журнал: АгроЗооТехника @azt-journal
Рубрика: Разведение, селекция и генетика сельскохозяйственных животных
Статья в выпуске: 4 т.2, 2019 года.
Бесплатный доступ
Цель настоящей работы - оценка генотипов овец романовской породы по комплексу признаков для прогнозирования племенных качеств животных, ожидаемого генетического улучшения по экономически значимым признакам в микропопуляциях, определения прогнозируемого экономического эффекта от использования селекционного индекса. Впервые предложен усовершенствованный метод для оценки генотипа овец романовской породы, учитывающий показатели многоплодия и мясной продуктивности. При выполнении работы были использованы вариационно-статистические, коррелятивно-регрессионные и дисперсионно-ковариансные методы обработки данных. Установлены существенные достоверные различия между стадами по живой массе маток по первому окоту и плодовитости, что обуславливает необходимость индивидуального подхода к конструированию индекса в каждом конкретном стаде. Уровень наследуемости по изучаемым показателям входит в низкие и средние границы, по хозяйствам варьирует в пределах 0,293-0,335 по живой массе и 0,121-0,184 по плодовитости. По изучаемым признакам в подконтрольных стадах овец романовской породы выявлена довольно высокая положительная фенотипическая и генотипическая зависимость между показателями живой массы и плодовитости. Показатели находились в пределах 0,487-0,549 для фенотипического и 0,491-0,581 для генотипического коэффициента корреляции. Это говорит о том, что использование данных признаков в одной модели селекционного индекса позволит вести достоверную оценку и эффективный отбор одновременно по обоим признакам. Точность оценки генотипа с применением предложенной модели по всем стадам превышала 0,400. Следовательно, индекс может быть использован для оценки генотипа овец и племенного отбора. Установлено, что при использовании предложенного селекционного индекса ожидаемый генетический эффект по живой массе при первом окоте составит от 0,591 до 1,356 кг на голову в год, ожидаемый выход ягнят увеличится на 0,036-0,045 ягненка на овцематку в год. Экономический эффект от использования модели составит до+51,08 тыс. руб. от реализации мяса баранины и +18,60 тыс. руб. от племенной продажи молодняка.
Овцы, романовская порода, плодовитость, живая масса, селекционно-генетические параметры, оценка генотипа, селекционный индекс
Короткий адрес: https://sciup.org/147225596
IDR: 147225596 | DOI: 10.15838/alt.2019.2.4.2
Текст научной статьи Комплексная оценка генотипа романовской породы овец в управлении селекционным процессом племенных стад
Цель настоящей работы – оценка генотипов овец романовской породы по комплексу признаков для прогнозирования племенных качеств животных, ожидаемого генетического улучшения по экономически значимым признакам в микропопуляциях, определения прогнозируемого экономического эффекта от использования селекционного индекса. Впервые предложен усовершенствованный метод для оценки генотипа овец романовской породы, учитывающий показатели многоплодия и мясной продуктивности. При выполнении работы были использованы вариационно-статистические, коррелятивно-регрессионные и дисперсионно-кова-риансные методы обработки данных. Установлены существенные достоверные различия между стадами по живой массе маток по первому окоту и плодовитости, что обуславливает необходимость индивидуального подхода к конструированию индекса в каждом конкретном стаде. Уровень наследуемости по изучаемым показателям входит в низкие и средние границы, по хозяйствам варьирует в пределах 0,293–0,335 по живой массе и 0,121–0,184 по плодовитости. По изучаемым признакам в подконтрольных стадах овец романовской породы выявлена довольно высокая положительная фенотипическая и генотипиче ская зависимость между показателями живой массы и плодовитости. Показатели находились в пределах 0,487–0,549 для фенотипического и 0,491–0,581 для генотипи ческого коэффициента корреляции. Это говорит о том, что использование данных признаков в одной модели селекционного индекса позволит вести достоверную оценку и эффективный отбор одновременно по обоим признакам. Точность оценки генотипа с применением предложенной модели по всем стадам превышала 0,400. Следователь но, индекс может быть использован для оценки генотипа овец и племенного отбора Установлено, что при использовании предложенного селекционного индекса ожидае мый генетический эффект по живой массе при первом окоте составит от 0,591 до 1,356 кг на голову в год, ожидаемый выход ягнят увеличится на 0,036–0,045 ягненка на овцематку в год. Экономический эффект от использования модели составит до +51,08 тыс. руб. от реализации мяса баранины и +18,60 тыс. руб. от племенной про дажи молодняка.
Овцы, романовская порода, плодовитость, живая масса, селекционно-генетические параметры, оценка генотипа, селекционный индекс.
На современном этапе развития отрасли овцеводства изменяются направления экономических и селекционных предпочтений, происходит смена ориентиров и продуктивной направленности. Все чаще для улучшения желательных в селекционном плане признаков используется межпородное скрещивание со специализированными породами. В связи с этим проблема сохранения генофонда отечественных пород сельскохозяйственных животных в Российской Федерации стоит как никогда остро. Романовская порода является старейшей отечественной породой и занимает особое положение среди многочисленных грубошерстных северных овец. Романовские овцы обладают комплексом ценных свойств, определяющих их важное народнохозяйственное значение. Совершенствование породы осуществляется при чистопородном разведении с применением современных научно обоснованных методов [1, с. 3–4; 2, с. 6–7; 3, с. 37].
Современный селекционный процесс представляет собой сложную систему, включающую множество факторов как генетической, так и паратипической природы, влияющих на окончательный результат. В связи с этим для анализа практической деятельности и прогноза важнейших характеристик селекционного процесса необходимо использовать современный математический аппарат, базирующийся на последних достижениях популяционной генетики и использующий возможности быстродействующих вычислительных средств. Основная задача селекции, в т. ч. и в романовском овцеводстве, заключается в том, чтобы из большого количества животных популяции отобрать для репродукции только тех, кто гарантированно даст потомство лучшего качества желательного направления продуктивности.
Рост продуктивных показателей может быть достигнут, в первую очередь, путем применения передовых методов селекционно-племенной работы. В связи с этим возникает необходимость создания селекционных программ, основным элементом которых является разработка и внедрение эффективных методов оценки племенных качеств животных для их дальнейшего практического применения в отдельных стадах и популяции романовских овец в целом [4, с. 274–275].
Оценка животных, в том числе и в овцеводстве, на практике проводится на различных этапах селекционно-племенной работы: предварительная – по происхождению, по собственному фенотипу и окончательная – по качеству потомства. При этом большое внимание уделяется экономической значимости признаков, по которым происходят оценка и отбор желательных особей. Классические методы селекции в племенном животноводстве по нескольким признакам имеют различные подходы и варианты отбора.
Тандемная селекция представляет собой метод, предусматривающий отбор сначала по одному признаку, далее включают второй, третий и так далее, последовательно доводя значения включенных в программу селекции признаков до желаемого уровня.
Отбор по независимым границам предполагает выставление так называемых «границ отбора» для выделенных признаков, представляющих собой селекционный оптимум фенотипического проявления. Из последующего разведения исключаются особи, не достигшие минимального порога по одному признаку, при прочих показателях, удовлетворяющих границам отбора. Недостаток метода в том, что выбраковываются животные, хотя бы по одному признаку не входящие в выделенные границы, по отдельным признакам превышающие средний уровень популяции. Существует угроза превращения стада в безлидерную группу с ограниченными селекционными возможностями.
Симультативная селекция (селекция по независимым уровням отбора) основывается на единовременном улучшении комплекса признаков, входящих в программу селекции с применением селекционных индексов.
В практике племенного дела индексная селекция – один из наиболее точных инструментов оценки и отбора. Метод селекционных индексов позволяет ускорить оценку племенных животных и повысить ее достоверность, однако он наиболее тру- доемкий и требует создания информационных баз больших объемов для оценок популяционно-генетических характеристик [5, с. 59–62].
Селекционные индексы, рассчитанные на основе комплекса генетических и фенотипических корреляционных связей, а также наследуемости и экономических значений признаков, дают наиболее полную оценку суммарного генотипа особей. Их использование позволяет более полно рассмотреть племенную ценность пробанда по совокупности селекционируемых признаков как с учетом их корреляционных связей, так и экономической значимости [6, с. 23–24; 7, с. 9; 8, с. 72–79].
В современных условиях ведения романовского овцеводства именно комплексное изучение признаков живой массы (в аспекте производства мяса баранины) и плодовитости (для получения большего числа племенного и товарного молодняка) является управляющим фактором в экономическом благополучии сельхозтоваропроизводителей [9, c. 32].
Научно обоснованный подход к оценке и отбору животных, включающий характеристику фенотипических и генетических связей, в практической селекционной работе дает возможность более точного прогнозирования племенных качеств животных следующих генераций.
Целью настоящей работы стала оценка генотипов овец романовской породы по комплексу признаков для прогнозирования племенных качеств животных, оценки ожидаемого генетического улучшения по экономически значимым признакам в микропопуляциях (стадах), определения прогнозируемого экономического эффекта от использования селекционных индексов.
В задачи исследований включены следующие направления:
-
- проведение ковариансного анализа для определения селекционно-генетических параметров, необходимых для вы-
- числения генетико-экономических весов в конструировании индексов;
-
- оценка взаимосвязей фенотипических показателей продуктивности у исследуемых животных с их генетическими оценками;
-
- проведение комплексной оценки генотипа племенных овец романовской породы с применением селекционного индекса;
-
- определение точности предложенного метода и ожидаемого генетического эффекта от его применения.
Новизна проводимых исследований заключается в отсутствии комплексных методов для оценки генотипа овец романовской породы, учитывающих показатели многоплодия и мясной продуктивности.
Материал и методы
Материалом для исследований послужили сведения по продуктивности, генеалогической структуре стад, сохранности молодняка и бонитировке стад овец романовской породы в племенных хозяйствах Ярославской области КХ Абдулатипова С.М. (Гаврилов-Ямс-кий район), ООО «Атис СХ» (Борисоглебский район), ООО «Агрофирма «Авангард» (Угличский район). В обработку вошла информация по овцематкам и баранам-производителям общей численностью 1176 голов.
Были рассчитаны основные селекционно-генетические параметры продуктивных признаков: средняя арифметическая (М), стандартная ошибка (m), стандартное отклонение признака ( σ ), коэффициент вариации (Cv), наследуемость (h2), повторяемость (r S ), фенотипическая (rP) и генетическая (r g ) корреляция [10, с. 161–175; 11, с. 267–294].
Оценка генотипа овец проводилась в соответствии с «Методическими рекомендациями по использованию селекционных индексов в племенной работе и анализу селекционно-генетических параметров признаков с альтернативной изменчивостью» и авторскими методиками [12, с. 25–26; 13, с. 66–91].
Отбор по селекционному индексу наиболее рационален, так как в нем учитывается комплекс генотипических и фенотипических связей, наследуемость и изменчивость признаков, а также их экономическая значимость для селекции. Рабочая модель имеет вид:
^^^^^^^^^^^^^^^ ^^^^^^^^^^^^^^^— —
-
/ = в 1 • (Х1 - Х1) + в 2 • (Х 2 - Х2) , (1)
где: _
-
Х1 – Х1 – отклонение по первому признаку;
Х2 – Х2 – отклонение по второму признаку;
в1 и в2 – коэффициенты линейной регрессии индекса.
Наиболее трудоемкий этап в построении индекса – расчет коэффициента линейной регрессии. Его вычисляют однократно (для конкретного стада) по уравнению:
P · в = G · a, (2)
где:
Р – матрица фенотипических варианс и кова-рианс;
G – матрица генетических варианс и коварианс; а – вектор-столбец относительных экономических весовых коэффициентов.
Это выражение представлено в виде матрицы:
[ вв 1 ] I" " ^"[с11 с12]х[“ 1 ] , (3)
L в 1-1 ["21 "22.1 1с21 с22-1 1-1111
где:
Р11 = σ 1 ; Р22 = σ 22; Р21 = Р12 = σ 2 σ 1rp1/2; G11 = h12 σ 12;
G22 = h22 σ 22; G12 = G21 = rg1/2 h1 σ 1h2 σ 2;
σ – стандартное отклонение признака;
h2 – коэффициент наследуемости признака.
Для признаков, по которым невозможно определить фактический экономический вес, определяется относительный экономический вес по формуле:
ai = a0^ Гр1 р0^ — , (4)
О1
где:
-
а0 – экономическое значение для основного признака (для которого есть стоимостное выражение);
-
σ 0 – стандартное фенотипическое отклонение основного признака;
-
σ 1 – стандартное отклонение для признака с неизвестным экономическим значением;
-
rp1p0 – коэффициент фенотипической корреляции.
Точность оценки генотипа (RJH) рассчитывали по формуле:
R jh = — , (5)
Он где:
-
О] =Jв 2 1 • Р1 1 + в 22 • ^2 2 + 2 • в 1 1 • в 12 • Р12 , (6)
Он = ДО2 • Си + 02 • G22 + 2 • Oi • «2 • ^12 (7)
Ожидаемый эффект селекции (Δg) за счет проводимых мероприятий рассчитывается по формуле:
в 11 • G11 + в 12 • 1^12 ./О\
^01 =--i,
°)
в 11 • 612 + в 12 • G22 ,
^92 --I ,
°)
где:
С ∙ П
Э = Ц ∙ ∙ Л ∙ К , (10)
где:
Э – стоимость дополнительной основной продукции, руб.;
Ц – закупочная цена единицы продукции в масштабах цен, действовавших в области, руб.;
С – средняя продуктивность животных данного вида или породы;
П – средняя прибавка основной продукции, выраженная в процентах на 1 голову животного;
Л – постоянный коэффициент уменьшения результата, связанного с дополнительными затратами на прибавочную продукцию, равный 0,75;
К – численность поголовья сельскохозяйственных животных улучшенного селекционного достижения.
Биометрическая обработка данных проводилась с использованием встроенных функций Microsoft Office Excel 2016.
Результаты исследований
При построении селекционных индексов необходимо оценить степень разнородности популяции и целесообразность проведения селекции по избранному направлению. Предложенная селекционная модель должна удовлетворять следующим требованиям:
-
1) оптимальный селекционный индекс максимально гарантирует отбор особей с лучшим генотипом, а значит и более высокой племенной ценностью вне зависимости от фенотипического проявления признака;
-
2) с экономической точки зрения взаимозависимость между селекционным индексом и суммарным генотипом и ожидаемое генетическое улучшение имеют максимальное значение;
-
3) ожидаемое генетическое улучшение должно составлять не менее 5 процентов от фенотипической характеристики признака.
. Для решения этих задач в подконтрольных стадах была проведена оценка селекционногенетических параметров продуктивных
Таблица 1. Популяционно-генетические параметры стад овец романовской породы
Параметр |
Показатель |
|||||
живая масса по 1 окоту |
плодовитость по 1 окоту |
|||||
КХ Абдулатипова С.М. |
ООО «Атис СХ» |
ООО «Агрофирма «Авангард» |
КХ Абдулатипова С.М. |
ООО «Атис СХ» |
ООО «Агрофирма «Авангард» |
|
Подконтрольное поголовье (n) |
268 |
43 |
259 |
268 |
43 |
259 |
Среднее значение признаков (М ± m) |
44,3±0,54 |
47,4±0,70*** |
53,7±0,31*** |
1,96±0,035 |
1,84±0,093*** |
2,17±0,046*** |
Стандартное отклонение (σ) |
8,86 |
4,58 |
4,93 |
0,57 |
0,61 |
0,74 |
Фенотипическая изменчивость (Сv) |
20,02 |
9,67 |
9,66 |
29,06 |
33,22 |
10,08 |
Генетическая изменчивость (h2) |
0,318 |
0,293 |
0,335 |
0,184 |
0,121 |
0,142 |
Повторяемость (rs) |
0,292 |
0,311 |
0,311 |
0,198 |
0,155 |
0,153 |
Источник: результаты собственных исследований.
Достоверность разности: * P ≥ 0,95; ** Р ≥ 0,99; *** Р ≥ 0,999.
признаков, являющихся приоритетными на современном этапе разведения романовской породы (табл. 1) .
Полученные результаты выявили существенные различия между стадами по живой массе маток по первому окоту и плодовитости, достоверность разницы между крайними значениями относится к третьему порогу критериев, что, в свою очередь обуславливает необходимость индивидуального подхода к конструированию индекса в каждом конкретном стаде.
По мнению ученых [15, с. 16–19; 16, с. 226], эффективность селекции зависит от величины генетической изменчивости признака. Уровень наследуемости по изучаемым показателям входит в низкие и средние границы, по хозяйствам варьирует в пределах 0,293–0,335 по живой массе и 0,121– 0,184 по плодовитости (рис. 1) .
Полученные результаты указывают на то, что по живой массе овец в КХ Абдулатипова С.М. фенотипическая и генетическая изменчивость входит в средние границы. Следовательно, популяция характеризуется как высокопродуктивная, стабильная по уровню продуктивности. Существует

■ живая масса ■ плодовитость
Рис. 1. Коэффициенты генетической изменчивости по экономически значимым признакам возможность повышения фенотипических величин признака при углубленной селекции и отборе лучших генотипов, оцененных по селекционному индексу. Повышение продуктивности может происходить сравнительно быстро при хороших условиях кормления и содержания.
При среднем значении наследуемости и низкой фенотипической изменчивости улучшение показателей живой массы в ООО «Атис СХ» и ООО «Агрофирма «Авангард» возможно с отбором высокопродук- тивных генотипов путем комплексной индексной селекции (семейная по материнской основе и индивидуальная по рангу индекса). Как обязательный элемент, должна учитываться оценка наследственных качеств особей.
Оценка данных по плодовитости, приведенных в таблице 1, позволяет сделать вывод о том, что в КХ Абдулатипова С.М. и ООО «Атис СХ» при низкой наследуемости и высокой степени изменчивости для улучшения признака целесообразно закрепление внутригрупповых наследственных особенностей [17, с. 63–64]. Существует возможность регулирования, в частности уменьшения, фенотипической изменчивости за счет изменения условий среды. В ООО «Агрофирма «Авангард» при низкой наследуемости отмечен и низкий уровень изменчивости, поэтому в данном стаде наряду с индексной селекцией рекомендуется проводить кратковременный отбор по независимым границам с целью ускоренного повышения продуктивности племенного ядра.
По изучаемым признакам в подконтрольных стадах овец романовской породы выявлена довольно высокая положительная фенотипическая и генотипическая зависимость между показателями живой массы и плодовитости. Показатели находились в пределах 0,487–0,549 для фенотипиче- ского и 0,491–0,581 для генотипического коэффициента корреляции. Это говорит о том, что использование указанных признаков в одной модели селекционного индекса позволит вести достоверную оценку и эффективный отбор одновременно по обоим признакам.
На рис. 2 приведена характеристика распределения оцененных животных в разрезе стад.
Из приведенного графика видно, что в пределы М+1 σ по стаду КХ Абдулатипова С.М. входит 77% оцененных животных, по стадам ООО «Атис СХ» и ООО «Агрофирма «Авангард» – 76 и 78% соответственно, что удовлетворяет требованиям отбора в селекционную часть (уровень выбраковки 30%).
В современных условиях программы и системы селекции разрабатываются с использованием передовых методов оценки племенной ценности и отбора животных для дальнейшего разведения. В связи с этим повышение точности оценки генотипа животных становится одной из важнейших задач. Об эффективности того или иного метода позволяет судить показатель точности оценки ( R JH ), который представляет собой коэффициент регрессии индекса на генотип (рис. 3) .
Выявлена высокая точность оценки генотипа с применением предложенной модели, показатель превысил 0,400. Следо-

КХ Абдулатипова С.М.
ООО «Атис СХ» ООО «Агрофирма «Авангард»
Рис. 2. Распределение оцененных по индексу овец в разрезе подконтрольных стад

Рис. 4. Степень ожидаемого генетического улучшения по живой массе
Рис. 3. Точность комплексной оценки генотипа пробанда вательно, индекс может быть использован для оценки генотипа овец и племенного отбора.
Основные цели селекции – улучшение хозяйственно-ценных признаков и увеличение продуктивности животных, соответствующих новым экономическим условиям. Следующим этапом работы стало выявление ожидаемого эффекта селекции за счет проводимых мероприятий.
Установлено, что при использовании предложенного селекционного индекса ожидаемый генетический эффект по живой массе при первом окоте составит от 0,591 до 1,356 кг на голову в год (рис. 4) .
При использовании усовершенствованного метода ожидаемый выход ягнят увеличится на 0,036–0,045 ягненка на овцематку в год (рис. 5) .
Полученные показатели были использованы при расчете экономической эффективности предлагаемой модели оценки генотипа (табл. 2) .
Из полученных результатов экономических оценок видно, что эффективность
0.039


0.045

КХ Абдулатипова С.М.
ООО «Атис СХ»
ООО «Агрофирма «Авангард»
Рис. 5. Степень ожидаемого генетического улучшения по плодовитости отбора по индексу в существенной степени зависит от уровня селекционной работы, что подтверждается ранее сделанными оценками селекционно-генетических ситуаций в стадах.
Заключение
Оценка селекционно-генетических параметров в племенных стадах романовских овец Ярославской области выявила неравномерную среднюю степень фенотипической изменчивости показателя живой мас-
Таблица 2. Экономическая эффективность индексной оценки генотипа, руб.
Хозяйство |
Получено дополнительной продукции |
|||
от увеличения живой массы |
от племпродажи молодняка |
|||
всего по стаду |
на 1 голову |
всего по стаду |
на 1 голову |
|
КХ Абдулатипова С.М. |
78328,48 |
292,27 |
13362,28 |
49,86 |
ООО «Атис СХ» |
6957,74 |
161,81 |
2158,81 |
50,20 |
ООО «Агрофирма «Авангард» |
51079,31 |
197,22 |
18603,31 |
71,83 |
Источник: результаты собственных исследований.
сы и плодовитости овцематок при первом окоте. Фенотипическая и генотипическая связь между признаками прямая умеренная, что в совокупности со средними показателями наследуемости позволяет использовать эти признаки в комплексной оценке генотипа овец.
Предлагаемый усовершенствованный метод обладает высокой степенью точности. Ожидаемый генетический эффект от его внедрения 5,67–7,83%, что в экономическом выражении составляет 161,81– 292,27 руб. на голову по живой массе и 49,86–71,83 руб. по племпродаже на одну голову. Полученные в ходе исследований результаты позволяют реализовать его в наставлении по оценке генотипа с последующим использованием при разработке селекционных программ и планов племенной работы.
Список литературы Комплексная оценка генотипа романовской породы овец в управлении селекционным процессом племенных стад
- Ерохин А.И., Карасев Е.А., Ерохин С.А. Состояние, динамика и тенденции в развитии овцеводства в мире и в России // Овцы, козы, шерстяное дело. 2019. № 3. С. 3–7.
- Абонеев В.В., Марченко В.В., Абонеева Е.В. Повышение эффективности научного обеспечения современного состояния овцеводства России // Овцы, козы, шерстяное дело. 2019. № 2. С. 5–9.
- Костылев М.Н., Барышева М.С. Продуктивность овец романовской породы в племенных хозяйствах Ярославской области // Овцы, козы, шерстяное дело. 2019. № 2. С. 37–39.
- Филинская О.В. Анализ влияния различных факторов на продуктивные признаки овцематок // Инновационный путь развития предприятий АПК: сб. науч. тр. по мат-лам XL Междунар. науч.-практ. конф. проф.-препод. состава. Ярославль, 2017. С. 274–278.
- Furlotte N.A., Eskin E. Efficient Multiple-Trait Association and Estimation of Genetic Correlation Using the Matrix-Variate Linear Mixed Model. Genetic, 2015, vol. 200, pp. 59–68.
- Комплексная оценка овец методом селекционных индексов / А.И. Горлов [и др.] // Сб. науч. тр. Ставроп. науч.-исслед. ин-та животноводства и кормопроизводства. 2009. Т. 2. № 2-2. С. 23–25.
- Герман Ю.И., Садыков Е.В., Селионова М.И. Индексная оценка племенной ценности овец мясо-шерстного направления продуктивности // Овцы, козы, шерстяное дело. 2019. № 3. С. 18–21.
- GebreMichael S.G. Sheep Resources of Ethiopia. Genetic diversity and breeding strategy. PhD thesis. Wageningen University. The Netherlands, 2008. 149 р. Available at: https://edepot.wur.nl/122006 (accessed 25.03.2019).
- Костылев М.Н., Абрамова М.В., Барышева М.С. Селекционная оценка продуктивных качеств племенных овец романовской породы // Вестн. АПК Ставрополья. 2018. № 3 (31). С. 30–33. DOI: 10.31279/2222-9345-2018-7-31-30-33
- Теоретические основы селекции животных / З.С. Никоро [и др.]. М.: Колос, 1968. 439 с.
- Урбах В.Ю. Биометрические методы. М.: Наука, 1964. 415 с.
- Косяченко Н.М. Анализ и оценка генетического потенциала ярославской породы крупного рогатого скота с разработкой методов по его контролю и управлению: автореф. дис. … д-ра биол. наук: 06.02.01. СПб., 1998. 35 с.
- Басовский Н.З., Попов В.П., Рудаков А.Н. Методические рекомендации по использованию селекционных индексов в племенной работе и анализу селекционно-генетических параметров признаков с альтернативной изменчивостью. Л.: ВНИИРГЖ, 1976. 121 с.
- Методика определения экономической эффективности использования в сельском хозяйстве результатов научно-исследовательских и опытно-конструкторских работ, новой техники, изобретений и рационализаторских предложений / Г.М. Лоза [и др.]. М.: Колос, 1980. 111 с.
- Зубенко Э.В. Оценка быков-производителей молочных пород по пожизненной продуктивности потомства: автореф. дис. … д-ра с.-х. наук: 06.02.07. Лесные Поляны Московской обл., 2015. 43 с.
- Популяционная генетика для животноводов-селекционеров / В. Шталь [и др.]. М.: Колос, 1973. 439 с.
- Коваленко В.Н. Математическое моделирование в селекции животных. Киев: Урожай, 1980. 162 с.