1
|
Silva Neto JB, Mota LFM, Londoño-Gil M, Schmidt PI, Rodrigues GRD, Ligori VA, Arikawa LM, Magnabosco CU, Brito LF, Baldi F. Genotype-by-environment interactions in beef and dairy cattle populations: A review of methodologies and perspectives on research and applications. Anim Genet 2024. [PMID: 39377556 DOI: 10.1111/age.13483] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2024] [Revised: 09/18/2024] [Accepted: 09/23/2024] [Indexed: 10/09/2024]
Abstract
Modern livestock production systems are characterized by a greater focus on intensification, involving managing larger numbers of animals to achieve higher productive efficiency and animal health and welfare within herds. Therefore, animal breeding programs need to be strategically designed to select animals that can effectively enhance production performance and animal welfare across a range of environmental conditions. Thus, this review summarizes the main methodologies used for assessing the levels of genotype-by-environment interaction (G × E) in cattle populations. In addition, we explored the importance of integrating genomic and phenotypic information to quantify and account for G × E in breeding programs. An overview of the structure of cattle breeding programs is provided to give insights into the potential outcomes and challenges faced when considering G × E to optimize genetic gains in breeding programs. The role of nutrigenomics and its impact on gene expression related to metabolism in cattle are also discussed, along with an examination of current research findings and their potential implications for future research and practical applications. Out of the 116 studies examined, 60 and 56 focused on beef and dairy cattle, respectively. A total of 83.62% of these studies reported genetic correlations across environmental gradients below 0.80, indicating the presence of G × E. For beef cattle, 69.33%, 24%, 2.67%, 2.67%, and 1.33% of the studies evaluated growth, reproduction, carcass and meat quality, survival, and feed efficiency traits, respectively. By contrast, G × E research in dairy cattle populations predominantly focused on milk yield and milk composition (79.36% of the studies), followed by reproduction and fertility (19.05%), and survival (1.59%) traits. The importance of G × E becomes particularly evident when considering complex traits such as heat tolerance, disease resistance, reproductive performance, and feed efficiency, as highlighted in this review. Genomic models provide a valuable avenue for studying these traits in greater depth, allowing for the identification of candidate genes and metabolic pathways associated with animal fitness, adaptation, and environmental efficiency. Nutrigenetics and nutrigenomics are emerging fields that require extensive investigation to maximize our understanding of gene-nutrient interactions. By studying various transcription factors, we can potentially improve animal metabolism, improving performance, health, and quality of products such as meat and milk.
Collapse
Affiliation(s)
- João B Silva Neto
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
- Department of Animal Sciences, Purdue University, West Lafayette, Indiana, USA
| | - Lucio F M Mota
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
| | - Marisol Londoño-Gil
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
| | - Patrícia I Schmidt
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
| | - Gustavo R D Rodrigues
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
- Beef Cattle Research Center, Institute of Animal Science, Sertãozinho, Brazil
| | - Viviane A Ligori
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
- Beef Cattle Research Center, Institute of Animal Science, Sertãozinho, Brazil
| | - Leonardo M Arikawa
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
| | | | - Luiz F Brito
- Department of Animal Sciences, Purdue University, West Lafayette, Indiana, USA
| | - Fernando Baldi
- Department of Animal Science, School of Agricultural and Veterinarian Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Brazil
| |
Collapse
|
2
|
Teodoro M, Maiorano AM, Campos GS, de Albuquerque LG, de Oliveira HN. Genetic parameters, genomic prediction, and identification of regulatory regions located on chromosome 14 for weight traits in Nellore cattle. J Anim Breed Genet 2024. [PMID: 39189106 DOI: 10.1111/jbg.12895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Revised: 08/06/2024] [Accepted: 08/11/2024] [Indexed: 08/28/2024]
Abstract
This study aimed to investigate functional variants in chromosome 14 (BTA14) and its impact in genomic selection for birth weight (BW), weaning weight (WW), and yearling weight (YW) in Nellore cattle. Genetic parameter estimation and the weighted single-step genomic best linear unbiased prediction (WssGBLUP) analyses were performed. Direct additive heritability estimates were high for WW and YW, and moderate for BW. Trait-associated variants distributed across multiple regions on BTA14 were observed in the weighted single-step genome-wide association studies (WssGWAS) results, implying a polygenic genetic architecture for weight in different ages. Several genes have been found in association with the weight traits, including the CUB And Sushi multiple domains 3 (CSMD3), thyroglobulin (TG), and diacylglycerol O-acyltransferase 1 (DGAT1) genes. The variance explained per SNP was higher in six functional classes of gene regulatory regions (5UTR, CpG islands, downstream, upstream, long non-coding RNA, and transcription factor binding sites (TFBS)), highlighting their importance for weight traits in Nellore cattle. A marginal increase in accuracy was observed when the selected functional variants (SV) information was considered in the WssGBLUP method, probably because of the small number of SV available on BTA14. The identified genes, pathways, and functions contribute to a better understanding of the genetic and physiological mechanisms regulating weight traits in the Nellore breed.
Collapse
Affiliation(s)
- Miller Teodoro
- Department of Animal Science, São Paulo State University, Jaboticabal, Brazil
| | | | | | | | | |
Collapse
|
3
|
Mota LFM, Giannuzzi D, Pegolo S, Toledo-Alvarado H, Schiavon S, Gallo L, Trevisi E, Arazi A, Katz G, Rosa GJM, Cecchinato A. Combining genetic markers, on-farm information and infrared data for the in-line prediction of blood biomarkers of metabolic disorders in Holstein cattle. J Anim Sci Biotechnol 2024; 15:83. [PMID: 38851729 PMCID: PMC11162571 DOI: 10.1186/s40104-024-01042-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2024] [Accepted: 04/28/2024] [Indexed: 06/10/2024] Open
Abstract
BACKGROUND Various blood metabolites are known to be useful indicators of health status in dairy cattle, but their routine assessment is time-consuming, expensive, and stressful for the cows at the herd level. Thus, we evaluated the effectiveness of combining in-line near infrared (NIR) milk spectra with on-farm (days in milk [DIM] and parity) and genetic markers for predicting blood metabolites in Holstein cattle. Data were obtained from 388 Holstein cows from a farm with an AfiLab system. NIR spectra, on-farm information, and single nucleotide polymorphisms (SNP) markers were blended to develop calibration equations for blood metabolites using the elastic net (ENet) approach, considering 3 models: (1) Model 1 (M1) including only NIR information, (2) Model 2 (M2) with both NIR and on-farm information, and (3) Model 3 (M3) combining NIR, on-farm and genomic information. Dimension reduction was considered for M3 by preselecting SNP markers from genome-wide association study (GWAS) results. RESULTS Results indicate that M2 improved the predictive ability by an average of 19% for energy-related metabolites (glucose, cholesterol, NEFA, BHB, urea, and creatinine), 20% for liver function/hepatic damage, 7% for inflammation/innate immunity, 24% for oxidative stress metabolites, and 23% for minerals compared to M1. Meanwhile, M3 further enhanced the predictive ability by 34% for energy-related metabolites, 32% for liver function/hepatic damage, 22% for inflammation/innate immunity, 42.1% for oxidative stress metabolites, and 41% for minerals, compared to M1. We found improved predictive ability of M3 using selected SNP markers from GWAS results using a threshold of > 2.0 by 5% for energy-related metabolites, 9% for liver function/hepatic damage, 8% for inflammation/innate immunity, 22% for oxidative stress metabolites, and 9% for minerals. Slight reductions were observed for phosphorus (2%), ferric-reducing antioxidant power (1%), and glucose (3%). Furthermore, it was found that prediction accuracies are influenced by using more restrictive thresholds (-log10(P-value) > 2.5 and 3.0), with a lower increase in the predictive ability. CONCLUSION Our results highlighted the potential of combining several sources of information, such as genetic markers, on-farm information, and in-line NIR infrared data improves the predictive ability of blood metabolites in dairy cattle, representing an effective strategy for large-scale in-line health monitoring in commercial herds.
Collapse
Affiliation(s)
- Lucio F M Mota
- Department of Agronomy, Food, Natural resources, Animals and Environment (DAFNAE), University of Padova, Legnaro, Padova, 35020, Italy
| | - Diana Giannuzzi
- Department of Agronomy, Food, Natural resources, Animals and Environment (DAFNAE), University of Padova, Legnaro, Padova, 35020, Italy.
| | - Sara Pegolo
- Department of Agronomy, Food, Natural resources, Animals and Environment (DAFNAE), University of Padova, Legnaro, Padova, 35020, Italy
| | - Hugo Toledo-Alvarado
- Department of Genetics and Biostatistics, School of Veterinary Medicine and Zootechnics, National Autonomous University of Mexico, Ciudad Universitaria, Mexico City, 04510, Mexico
| | - Stefano Schiavon
- Department of Agronomy, Food, Natural resources, Animals and Environment (DAFNAE), University of Padova, Legnaro, Padova, 35020, Italy
| | - Luigi Gallo
- Department of Agronomy, Food, Natural resources, Animals and Environment (DAFNAE), University of Padova, Legnaro, Padova, 35020, Italy
| | - Erminio Trevisi
- Department of Animal Science, Food and Nutrition (DIANA) and the Romeo and Enrica Invernizzi Research Center for Sustainable Dairy Production (CREI), Faculty of Agricultural, Food, and Environmental Sciences, Università Cattolica del Sacro Cuore, Piacenza, 29122, Italy
| | | | - Gil Katz
- Afimilk LTD, Afikim, 15148, Israel
| | - Guilherme J M Rosa
- Department of Animal and Dairy Sciences, University of Wisconsin, Madison, WI, 53706, USA
| | - Alessio Cecchinato
- Department of Agronomy, Food, Natural resources, Animals and Environment (DAFNAE), University of Padova, Legnaro, Padova, 35020, Italy
| |
Collapse
|
4
|
Khazaei-Koohpar H, Gholizadeh M, Hafezian SH, Esmaeili-Fard SM. Weighted single-step genome-wide association study for direct and maternal genetic effects associated with birth and weaning weights in sheep. Sci Rep 2024; 14:13120. [PMID: 38849438 PMCID: PMC11161479 DOI: 10.1038/s41598-024-63974-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Accepted: 06/04/2024] [Indexed: 06/09/2024] Open
Abstract
Body weight is an important economic trait for sheep meat production, and its genetic improvement is considered one of the main goals in the sheep breeding program. Identifying genomic regions that are associated with growth-related traits accelerates the process of animal breeding through marker-assisted selection, which leads to increased response to selection. In this study, we conducted a weighted single-step genome-wide association study (WssGWAS) to identify potential candidate genes for direct and maternal genetic effects associated with birth weight (BW) and weaning weight (WW) in Baluchi sheep. The data used in this research included 13,408 birth and 13,170 weaning records collected at Abbas-Abad Baluchi Sheep Breeding Station, Mashhad-Iran. Genotypic data of 94 lambs genotyped by Illumina 50K SNP BeadChip for 54,241 markers were used. The proportion of variance explained by genomic windows was calculated by summing the variance of SNPs within 1 megabase (Mb). The top 10 window genomic regions explaining the highest percentages of additive and maternal genetic variances were selected as candidate window genomic regions associated with body weights. Our findings showed that for BW, the top-ranked genomic regions (1 Mb windows) explained 4.30 and 4.92% of the direct additive and maternal genetic variances, respectively. The direct additive genetic variance explained by the genomic window regions varied from 0.31 on chromosome 1 to 0.59 on chromosome 8. The highest (0.84%) and lowest (0.32%) maternal genetic variances were explained by genomic windows on chromosome 10 and 17, respectively. For WW, the top 10 genomic regions explained 6.38 and 5.76% of the direct additive and maternal genetic variances, respectively. The highest and lowest contribution of direct additive genetic variances were 1.37% and 0.42%, respectively, both explained by genomic regions on chromosome 2. For maternal effects on WW, the highest (1.38%) and lowest (0.41%) genetic variances were explained by genomic windows on chromosome 2. Further investigation of these regions identified several possible candidate genes associated with body weight. Gene ontology analysis using the DAVID database identified several functional terms, such as translation repressor activity, nucleic acid binding, dehydroascorbic acid transporter activity, growth factor activity and SH2 domain binding.
Collapse
Affiliation(s)
- Hava Khazaei-Koohpar
- Department of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University (SANRU), Sari, Iran
| | - Mohsen Gholizadeh
- Department of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University (SANRU), Sari, Iran.
| | - Seyed Hasan Hafezian
- Department of Animal Science and Fisheries, Sari Agricultural Sciences and Natural Resources University (SANRU), Sari, Iran
| | | |
Collapse
|
5
|
Londoño-Gil M, López-Correa R, Aguilar I, Magnabosco CU, Hidalgo J, Bussiman F, Baldi F, Lourenco D. Strategies for genomic predictions of an indicine multi-breed population using single-step GBLUP. J Anim Breed Genet 2024. [PMID: 38812461 DOI: 10.1111/jbg.12882] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Revised: 05/10/2024] [Accepted: 05/15/2024] [Indexed: 05/31/2024]
Abstract
Brazilian livestock breeding programmes strive to enhance the genetics of beef cattle, with a strong emphasis on the Nellore breed, which has an extensive database and has achieved significant genetic progress in the last years. There are other indicine breeds that are economically important in Brazil; however, these breeds have more modest sets of phenotypes, pedigree and genotypes, slowing down their genetic progress as their predictions are less accurate. Combining several breeds in a multi-breed evaluation could help enhance predictions for those breeds with less information available. This study aimed to evaluate the feasibility of multi-breed, single-step genomic best linear unbiased predictor genomic evaluations for Nellore, Brahman, Guzerat and Tabapua. Multi-breed evaluations were contrasted to the single-breed ones. Data were sourced from the National Association of Breeders and Researchers of Brazil and included pedigree (4,207,516), phenotypic (328,748), and genomic (63,492) information across all breeds. Phenotypes were available for adjusted weight at 210 and 450 days of age, and scrotal circumference at 365 days of age. Various scenarios were evaluated to ensure pedigree and genomic information compatibility when combining different breeds, including metafounders (MF) or building the genomic relationship matrix with breed-specific allele frequencies. Scenarios were compared using the linear regression method for bias, dispersion and accuracy. The results showed that using multi-breed evaluations significantly improved accuracy, especially for smaller breeds like Guzerat and Tabapua. The validation statistics indicated that the MF approach provided accurate predictions, albeit with some bias. While single-breed evaluations tended to have lower accuracy, merging all breeds in multi-breed evaluations increased accuracy and reduced dispersion. This study demonstrates that multi-breed genomic evaluations are proper for indicine beef cattle breeds. The MF approach may be particularly beneficial for less-represented breeds, addressing limitations related to small reference populations and incompatibilities between G and A22. By leveraging genomic information across breeds, breeders and producers can make more informed selection decisions, ultimately improving genetic gain in these cattle populations.
Collapse
Affiliation(s)
- Marisol Londoño-Gil
- Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, Brazil
- Department of Animal and Dairy Science, University of Georgia, Athens, Georgia, USA
| | - Rodrigo López-Correa
- Departamento de Genética y Mejoramiento Animal, Facultad de Veterinaria, Universidad de la República, Montevideo, Uruguay
| | - Ignacio Aguilar
- Instituto Nacional de Investigación Agropecuaria (INIA), Montevideo, Uruguay
| | | | - Jorge Hidalgo
- Department of Animal and Dairy Science, University of Georgia, Athens, Georgia, USA
| | - Fernando Bussiman
- Department of Animal and Dairy Science, University of Georgia, Athens, Georgia, USA
| | - Fernando Baldi
- Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista Júlio de Mesquita Filho, Jaboticabal, SP, Brazil
| | - Daniela Lourenco
- Department of Animal and Dairy Science, University of Georgia, Athens, Georgia, USA
| |
Collapse
|
6
|
Pan R, Qi L, Xu Z, Zhang D, Nie Q, Zhang X, Luo W. Weighted single-step GWAS identified candidate genes associated with carcass traits in a Chinese yellow-feathered chicken population. Poult Sci 2024; 103:103341. [PMID: 38134459 PMCID: PMC10776626 DOI: 10.1016/j.psj.2023.103341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2023] [Revised: 11/26/2023] [Accepted: 11/28/2023] [Indexed: 12/24/2023] Open
Abstract
Carcass traits in broiler chickens are complex traits that are influenced by multiple genes. To gain deeper insights into the genetic mechanisms underlying carcass traits, here we conducted a weighted single-step genome-wide association study (wssGWAS) in a population of Chinese yellow-feathered chicken. The objective was to identify genomic regions and candidate genes associated with carcass weight (CW), eviscerated weight with giblets (EWG), eviscerated weight (EW), breast muscle weight (BMW), drumstick weight (DW), abdominal fat weight (AFW), abdominal fat percentage (AFP), gizzard weight (GW), and intestine length (IL). A total of 1,338 broiler chickens with phenotypic and pedigree information were included in this study. Of these, 435 chickens were genotyped using a 600K single nucleotide polymorphism chip for association analysis. The results indicate that the most significant regions for 9 traits explained 2.38% to 5.09% of the phenotypic variation, from which the region of 194.53 to 194.63Mb on chromosome 1 with the gene RELT and FAM168A identified on it was significantly associated with CW, EWG, EW, BMW, and DW. Meanwhile, the 5 traits have a strong genetic correlation, indicating that the region and the genes can be used for further research. In addition, some candidate genes associated with skeletal muscle development, fat deposition regulation, intestinal repair, and protection were identified. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses suggested that the genes are involved in processes such as vascular development (CD34, FGF7, FGFR3, ITGB1BP1, SEMA5A, LOXL2), bone formation (FGFR3, MATN1, MEF2D, DHRS3, SKI, STC1, HOXB1, HOXB3, TIPARP), and anatomical size regulation (ADD2, AKT1, CFTR, EDN3, FLII, HCLS1, ITGB1BP1, SEMA5A, SHC1, ULK1, DSTN, GSK3B, BORCS8, GRIP2). In conclusion, the integration of phenotype, genotype, and pedigree information without creating pseudo-phenotype will facilitate the genetic improvement of carcass traits in chickens, providing valuable insights into the genetic architecture and potential candidate genes underlying carcass traits, enriching our understanding and contributing to the breeding of high-quality broiler chickens.
Collapse
Affiliation(s)
- Rongyang Pan
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Xugang Yellow Poultry Seed Industry Group Co., Ltd, Jiangmen City, Guangdong Province, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, Guangzhou 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Lin Qi
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, Guangzhou 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Zhenqiang Xu
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, Guangzhou 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Dexiang Zhang
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, Guangzhou 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Qinghua Nie
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, Guangzhou 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Xiquan Zhang
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, Guangzhou 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Wen Luo
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, Guangzhou 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou 510642, China.
| |
Collapse
|