doi: 10.15389/agrobiology.2022.6.1083rus
УДК 636.2:637.043:577.21
Исследование выполнено при финансовой поддержке Российского научного фонда (проект № 21-76-20046).
ИСПОЛЬЗОВАНИЕ ВЫСОКОПРОИЗВОДИТЕЛЬНОЙ ИНФРАКРАСНОЙ СПЕКТРОМЕТРИИ ПРИ ИЗУЧЕНИИ ПОЛНОГЕНОМНЫХ АССОЦИАЦИЙ С ЖИРНОКИСЛОТНЫМ И КОМПОНЕНТНЫМ СОСТАВОМ МОЛОКА У КОРОВ (Bos taurus)
А.А. СЕРМЯГИН✉, Л.П. ИГНАТЬЕВА, И.А. ЛАШНЕВА, А.А. КОСИЦИН,
О.В. КОСИЦИНА, А.С. АБДЕЛЬМАНОВА, Н.А. ЗИНОВЬЕВА
Количество жира в молоке коров относится к признакам, наиболее подверженным высокой изменчивости, и зависит как от условий среды (кормление, технология содержания), так и от генетических факторов (порода, генотип). Особый интерес представляет содержание жирных кислот (ЖК) как биомаркера контроля физиологического состояния животных и критерия оценки показателей качества сырого молока, в части его пригодности к переработке (выход сыра, масла и сливок). Соотношение ЖК в молоке по числу атомов углерода, а также длине цепи, степени ее насыщения различается между особями и на популяционном уровне. Поэтому изучение генетической и геномной изменчивости признаков молочной продуктивности для повышения эффективности управления отбором животных остается актуальной задачей. Цель наших исследований заключалась в поиске полногеномных ассоциаций и полиморфизмов в генах, детерминирующих жирнокислотный состав молока, на основе инфракрасной спектрометрии как одного из наиболее быстрых и точных экспресс-методов физико-химического анализа состава молока. Популяционно-генетические параметры и изменчивость содержания жирных кислот в молоке изучали на популяции голштинизированного черно-пестрого и голштинского скота 14 племенных стад из Московской области (2017-2018 годы). Всего для оценки суточных показателей молочной продуктивности использовали 36982 образца. Расчет коэффициентов наследуемости (h2) и корреляции (rg) показателей состава молока коров проводили на основе метода REML (residual maximumlikelihood) с использованием семейства программ BLUPF90. Поиск SNP проводили в выборке коров из экспериментального стада голштинизированного черно-пестрого скота (ПЗ «Ладожский» — филиал ФГБНУ ФИЦ ВИЖ им. Л.К. Эрнста, 2020-2021 годы). Фракционный состав молока определяли с помощью автоматического анализатора MilkoScan 7 DC («FOSS», Дания), принцип действия которого основан на экспресс-оценке методом инфракрасной спектроскопии. При индивидуальной оценке была сформирована группа из 144 генотипированных с помощью биочипа Bovine GGP 150K («Neogen», США) коров с полным фенотипическим описанием спектра жирных кислот и компонентов молока. Контроль качества генотипирования (110884 SNP), анализ полногеномных ассоциаций (GWAS, genome-wide association study) и многомерное шкалирование (MDS, multidimensional scaling) выполняли с помощью программы Plink 1.9. Поиск генов по выявленным значимым полиморфизмам проводили в браузере Ensembl по сборке генома крупного рогатого скота Bos taurus UMD 3.1.1 (https://www.ncbi.nlm.nih.gov/assembly/). Аннотацию генов для определения локусов количественных признаков на хромосомах животных осуществляли по международной базе данных Animal QTLdb. Наследуемость показателей жирных кислот молока варьировала от низкой для полиненасыщенных ЖК (h2 = 0,018) до умеренных для средне- (h2 = 0,125), длинноцепочечных (h2 = 0,155)и миристиновой кислоты (h2 = 0,155), мононенасыщенных ЖК (h2 = 0,176) и олеиновой кислоты (h2 = 0,196). Методом многомерного шкалирования оценили генетическую структуру экспериментальной группы животных, где мы наблюдали умеренный размах вариабельности по первой (PC1 = 7,82 %) и по второй (PC2 = 4,65 %) компонентам изменчивости. Для миристиновой и пальмитиновой ЖК выявлены общие кластеры (QTL) на хромосомах BTA5, BTA10, BTA14, BTA18 и BTA27; для стеариновой и олеиновой ЖК (как входящих в группу длинноцепочечные ЖК) показана схожая локализация QTL на хромосомах BTA9, BTA10, BTA11, BTA14, BTA17, BTA18, BTA19, BTA20 и BTA29. Для коротко- и среднецепочечных ЖК обнаружены ассоциации на хромосомах BTA1, BTA5, BTA10, BTA11, BTA14, BTA18, BTA19 и BTA24, для длинноцепочечных ЖК детектированы QTL на BTA6, BTA7, BTA9, BTA10, BTA11, BTA17, BTA18 и BTA29. Для коротко- и среднецепочечных ЖК, насыщенных ЖК, C14:0, C16:0, C18:0 и C18:1 установлены гены, образующие QTL на хромосомах BTA10, BTA11 и BTA14, — CACNA1C, GCH1,ATG14, KCNH5, PRKCE, CTNNA2, CYHR1, VPS28, DGAT1, ZC3H3, RHPN1, TSNARE1. Коротко- и среднецепочечные ЖК, миристиновая и пальмитиновая ЖК, насыщенные ЖК показали связь с полиморфизмами в генах MED12L, EPHB1, GRIN2B, PRMT8, ERC1, PELI2, ARHGAP39, MROH1, MAF1, GSDMD, LY6D. Для длинноцепочечных и мононенасыщенных ЖК, стеариновой и олеиновой ЖК в результате аннотации получены селекционно значимые гены RPS6KA2, CPQ, CPE, FTO, FAT3, LUZP2. Продолжение изучения генетических механизмов наследования содержания жирных кислот и компонентов молока необходимо для формирования стратегии селекции молочного скота с лучшим жирнокислотным профилем и составом компонентов.
Ключевые слова: корова, жирные кислоты, компоненты молока, коэффициент наследуемости, GWAS, SNP, QTL, гены.
A.A. Sermyagin✉, L.P. Ignatieva, I.A. Lashneva, A.A. Kositsin,
O.V. Kositsina, A.S. Abdelmanova, N.A. Zinovieva
Milk fat percentage is highly variabile and depends on environmental conditions which include feeding and farm technology, and on genetic factors such as breed and genotype features. The content of fatty acids (FA) is a biomarker for the physiological state of animals and a parameter of raw milk suitability for processing (yield of cheese, butter and cream). FA profile of mild in terms of C number, the chain length and saturation degree differs between individuals and at the population level. Therefore, the study of genetic and genomic variability of milk production traits to improve the efficiency of animal selection remains relevant. This study aimed at searching for genome-wide associations and polymorphisms in genes involved in milk fatty acid production. In the study, infrared spectrometry was used as an accurate and rapid method to analyze milk composition. Population variability of milk fatty acid profiles was studied using 36982 milk samples from Holsteinized Black-and-White and Holstein cows of 14 breeding herds from the Moscow region in 2017-2018. The heritability (h2) and correlation (rg) coefficients for cows’ milk components were calculated using REML (residual maximum likelihood) method with BLUPF90 family software. SNPs were detected for a dataset of Holsteinized Black-and-White cows from an experimental herd (the breeding farm Ladozhsky, branch of Ernst Federal Research Center for Animal Husbandry, Krasnodar Territory, 2020-2021). Milk composition was determined using an infrared spectroscopy-based automatic MilkoScan 7 DC analyzer (FOSS, Denmark). A group of 144 cows subjected to phenotyping for fatty acids and milk components were individually genotyped (Bovine GGP 150K biochip, Neogen, USA). Plink 1.9 software was applied to control genotyping quality (110884 SNPs) and to perform GWAS (genome-wide association study) analysis and multidimensional scaling (MDS). Searching genes by identified significant polymorphisms was performed using the bovine genome assembly Bos taurus UMD 3.1.1 (https://www.ncbi.nlm.nih.gov/assembly/) and the Ensembl genome browser. QTL annotation was carried out using the Animal QTLdb database. In general, milk fatty acids showed a heritability level that ranged from low to moderate, varying from h2 = 0.018 for polyunsaturated fatty acids to h2 = 0.125 for medium-chain FAs, h2 = 0.155 for long-chain FAs, h2 = 0.155 for myristic acid, h2 = 0.176 for monounsaturated FAs, and h2 = 0.196 for oleic acid. Visualizing experimental cows’ population structure by multidimensional scaling showed a moderate range of variability (PC1 = 7.82 %, PC2 = 4.65 %). For myristic and palmitic acids, common QTL clusters are identified on BTA5, BTA10, BTA14, BTA18, and BTA27. For stearic and oleic acids (as members of the long-chain FA family), similar location of QTLs is found on BTA9, BTA10, BTA11, BTA14, BTA17, BTA18, BTA19, BTA20, and BTA29. For short- and medium-chain FAs, there are associations revealed on BTA1, BTA5, BTA10, BTA11, BTA14, BTA18, BTA19, and BTA24. For long-chain FAs, QTLs are detected on BTA6, BTA7, BTA9, BTA10, BTA11, BTA17, BTA18, and BTA29. For short- and medium-chain FAs, saturated FAs, C14:0, C16:0, C18:0 and C18:1, the genes CACNA1C, GCH1, ATG14, KCNH5, PRKCE, CTNNA2, CYHR1, VPS28, DGAT1, ZC3H3, RHPN1, TSNARE1 are identified which form QTLs on BTA10, BTA11 and BTA14. Short- and medium-chain FAs, myristic and palmitic acids and saturated FAs show associations with polymorphisms in the MED12L, EPHB1, GRIN2B, PRMT8, ERC1, PELI2,ARHGAP39, MROH1, MAF1, GSDMD, and LY6D genes. For long-chain, monounsaturated fatty acids, stearic and oleic acids, there are significant associations with genesRPS6KA2, CPQ, CPE, FTO, FAT3, and LUZP2 which may be valuable for genetic improvement of dairy cattle. Continued study of the inheritance of cows’ milk fatty acids and other components is necessary to develop a strategy for breeding dairy cattle with a better fatty acid profile and milk composition.
Keywords: cow, fatty acids, milk components, heritability, GWAS, SNP, QTL, genes.
ФГБНУ Федеральный исследовательский центр |
Поступила в редакцию |
