1
|
Kafi Z, Masoudi AA, Torshizi RV, Ehsani A. Copy number variations affecting growth curve parameters in a crossbred chicken population. Gene 2024; 927:148710. [PMID: 38901536 DOI: 10.1016/j.gene.2024.148710] [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: 12/31/2023] [Revised: 06/01/2024] [Accepted: 06/17/2024] [Indexed: 06/22/2024]
Abstract
Copy number variations (CNVs) are key structural variations in the genome and may contribute to phenotypic differences. In this study, we used a F2 chicken population created from reciprocal crossing between fast-growing Arian broiler line and Urmia native chickens. The chickens were genotyped by 60 K SNP BeadChip, and PennCNV algorithm was used to detect genome-wide CNVs. The growth curve parameters of W0, k, L, Wf, Wi, ti and average GR were used as phenotypic data. The association between CNV and growth curve parameters was carried out using the CNVRanger R/Bioconductor package. Five CNV regions (CNVRs) were chosen for the validation experiment using qPCR. Gene enrichment analysis was done using WebGestalt. The STRING database was used to search for significant pathways. The results identified 966 CNVs and 600 CNVRs including 468 gains, 67 losses, and 65 both events on autosomal chromosomes. Validation of the CNVRs obtained from the qPCR assay were 79 % consistent with the prediction by PennCNV. A total of 43 significant CNVs were obtained for the seven growth curve parameters. The 416 genes annotated for significant CNVs. Six genes out of 416 genes were most related to growth curve parameters. These genes were LCP2, Dock2, CD80, CYFIP1, NIPA1 and NIPA2. Some of these genes in their biological process were associated with the growth, reproduction and development of cells or organs that ultimately lead to the growth of the body. The results of the study could pave the way for better understanding the molecular process of CNVs and growth curve parameters in birds.
Collapse
Affiliation(s)
- Zeinab Kafi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| | - Ali Akbar Masoudi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| | - Rasoul Vaez Torshizi
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| | - Alireza Ehsani
- Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran.
| |
Collapse
|
2
|
Tu TC, Lin CJ, Liu MC, Hsu ZT, Chen CF. Comparison of genomic prediction accuracy using different models for egg production traits in Taiwan country chicken. Poult Sci 2024; 103:104063. [PMID: 39098301 DOI: 10.1016/j.psj.2024.104063] [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: 12/15/2023] [Revised: 06/20/2024] [Accepted: 07/01/2024] [Indexed: 08/06/2024] Open
Abstract
In local chickens targeted for niche markets, genotyping costs are relatively high due to the small population size and diverse breeding goals. The single-step genomic best linear unbiased prediction (ssGBLUP) model, which combines pedigree and genomic information, has been introduced to increase the accuracy of genomic estimated breeding value (GEBV). Therefore, this model may be more beneficial than the genomic BLUP (GBLUP) model for genomic selection in local chickens. Additionally, the single-step genome-wide association study (ssGWAS) can be used to extend the ssGBLUP model results to animals with available phenotypic information but without genotypic data. In this study, we compared the accuracy of (G)EBVs using the pedigree-based BLUP (PBLUP), GBLUP, and ssGBLUP models. Moreover, we conducted single-SNP GWAS (SNP-GWAS), GBLUP-GWAS, and ssGWAS methods to identify genes associated with egg production traits in the NCHU-G101 chicken to understand the feasibility of using genomic selection in a small population. The average prediction accuracy of (G)EBV for egg production traits using the PBLUP, GBLUP, and ssGBLUP models is 0.536, 0.531, and 0.555, respectively. In total, 22 suggestive- and 5% Bonferroni genome-wide significant-level SNPs for total egg number (EN), average laying rate (LR), average clutch length, and total clutch number are detected using 3 GWAS methods. These SNPs are mapped onto Gallus gallus chromosomes (GGA) 4, 6, 10, 18, and 25 in NCHU-G101 chicken. Furthermore, through SNP-GWAS and ssGWAS methods, we identify 2 genes on GGA4 associated with EN and LR: ENSGALG00000023172 and PPARGC1A. In conclusion, the ssGBLUP model demonstrates superior prediction accuracy, performing on average 3.41% than the PBLUP model. The implications of our gene results may guide future selection strategies for Taiwan Country chickens. Our results highlight the applicability of the ssGBLUP model for egg production traits selection in a small population, specifically NCHU-G101 chicken in Taiwan.
Collapse
Affiliation(s)
- Tsung-Che Tu
- Department of Animal Science, National Chung Hsing University, Taichung 402, Taiwan; Ray Hsing Agricultural Biotechnology Co. Ltd., Yunlin 633, Taiwan
| | - Chen-Jyuan Lin
- Department of Animal Science, National Chung Hsing University, Taichung 402, Taiwan
| | - Ming-Che Liu
- Ray Hsing Agricultural Biotechnology Co. Ltd., Yunlin 633, Taiwan
| | - Zhi-Ting Hsu
- Ray Hsing Agricultural Biotechnology Co. Ltd., Yunlin 633, Taiwan
| | - Chih-Feng Chen
- Department of Animal Science, National Chung Hsing University, Taichung 402, Taiwan; The iEGG and Animal Biotechnology Center, National Chung Hsing University, Taichung 402, Taiwan.
| |
Collapse
|
3
|
Yang J, Wang DF, Huang JH, Zhu QH, Luo LY, Lu R, Xie XL, Salehian-Dehkordi H, Esmailizadeh A, Liu GE, Li MH. Structural variant landscapes reveal convergent signatures of evolution in sheep and goats. Genome Biol 2024; 25:148. [PMID: 38845023 PMCID: PMC11155191 DOI: 10.1186/s13059-024-03288-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 05/21/2024] [Indexed: 06/10/2024] Open
Abstract
BACKGROUND Sheep and goats have undergone domestication and improvement to produce similar phenotypes, which have been greatly impacted by structural variants (SVs). Here, we report a high-quality chromosome-level reference genome of Asiatic mouflon, and implement a comprehensive analysis of SVs in 897 genomes of worldwide wild and domestic populations of sheep and goats to reveal genetic signatures underlying convergent evolution. RESULTS We characterize the SV landscapes in terms of genetic diversity, chromosomal distribution and their links with genes, QTLs and transposable elements, and examine their impacts on regulatory elements. We identify several novel SVs and annotate corresponding genes (e.g., BMPR1B, BMPR2, RALYL, COL21A1, and LRP1B) associated with important production traits such as fertility, meat and milk production, and wool/hair fineness. We detect signatures of selection involving the parallel evolution of orthologous SV-associated genes during domestication, local environmental adaptation, and improvement. In particular, we find that fecundity traits experienced convergent selection targeting the gene BMPR1B, with the DEL00067921 deletion explaining ~10.4% of the phenotypic variation observed in goats. CONCLUSIONS Our results provide new insights into the convergent evolution of SVs and serve as a rich resource for the future improvement of sheep, goats, and related livestock.
Collapse
Affiliation(s)
- Ji Yang
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Dong-Feng Wang
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, 100101, China
- College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, 100049, China
| | - Jia-Hui Huang
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Qiang-Hui Zhu
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, 100101, China
- College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, 100049, China
| | - Ling-Yun Luo
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Ran Lu
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China
| | - Xing-Long Xie
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, 100101, China
- College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, 100049, China
| | - Hosein Salehian-Dehkordi
- CAS Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences (CAS), Beijing, 100101, China
- College of Life Sciences, University of Chinese Academy of Sciences (UCAS), Beijing, 100049, China
| | - Ali Esmailizadeh
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, Kerman, 76169-133, Iran
| | - George E Liu
- Animal Genomics and Improvement Laboratory, BARC, USDA-ARS, Beltsville, MD, 20705, USA
| | - Meng-Hua Li
- State Key Laboratory of Animal Biotech Breeding, China Agricultural University, Beijing, 100193, China.
- College of Animal Science and Technology, China Agricultural University, Beijing, 100193, China.
| |
Collapse
|
4
|
Xiao L, Qi L, Fu R, Nie Q, Zhang X, Luo W. A large-scale comparison of the meat quality characteristics of different chicken breeds in South China. Poult Sci 2024; 103:103740. [PMID: 38701629 PMCID: PMC11087722 DOI: 10.1016/j.psj.2024.103740] [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: 01/19/2024] [Revised: 03/26/2024] [Accepted: 04/03/2024] [Indexed: 05/05/2024] Open
Abstract
Meat quality traits are essential for producing high-quality broilers, but the genetic improvement has been limited by the complexity of measurement methods and the numerous traits involved. To systematically understand the meat quality characteristics of different broiler breeds, this study collected data on slaughter performance, skin color, fat deposition, and meat quality traits of 434 broilers from 12 different breeds in South China. The results showed that there was no significant difference in the live weight and slaughter weight of various broiler breeds at their respective market ages. Commercial broiler breeds such as Xiaobai and Huangma chickens had higher breast muscle and leg muscle rates. The skin and abdominal fat of Huangma chickens cultivated in the consumer market in South China exhibited significantly higher levels of yellowness compared to other varieties. Concerning fat traits, we observed that Wenchang chickens exhibited a strong ability to fat deposition, while the younger breeds showed lower fat deposition. Additionally, there were significant positive correlations found among different traits, including traits related to weight, traits related to fat, and skin color of different parts. Hierarchical clustering analysis revealed that fast-growing and large broiler Xiaobai chickens formed a distinct cluster based on carcass characteristics, skin color, and meat quality traits. Principal component analysis (PCA) was used to extract multiple principal components as substitutes for complex meat quality indicators, establishing a chicken meat quality evaluation model to differentiate between different breeds of chickens. At the same time, we identified 46, 22, and 20 SNP loci and their adjacent genes that were significantly associated with muscle mass traits, fat deposition, and skin color through genome-wide association studies (GWAS). The above results are helpful for systematically understanding the differences and characteristics of meat quality traits among different breeds.
Collapse
Affiliation(s)
- Liangchao Xiao
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, 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
| | - Lin Qi
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, 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
| | - Rong Fu
- Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, 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
| | - Qinghua Nie
- State Key Laboratory of Livestock and Poultry Breeding, and Lingnan Guangdong Laboratory of Agriculture, 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; 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
| | - Xiquan Zhang
- State Key Laboratory of Livestock and Poultry Breeding, and Lingnan Guangdong Laboratory of Agriculture, 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; 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
| | - Wen Luo
- State Key Laboratory of Livestock and Poultry Breeding, and Lingnan Guangdong Laboratory of Agriculture, 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; 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.
| |
Collapse
|
5
|
Wang SZ, Wang MD, Wang JY, Yuan M, Li YD, Luo PT, Xiao F, Li H. Genome-wide association study of growth curve parameters reveals novel genomic regions and candidate genes associated with metatarsal bone traits in chickens. Animal 2024; 18:101129. [PMID: 38574453 DOI: 10.1016/j.animal.2024.101129] [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: 10/26/2023] [Revised: 03/02/2024] [Accepted: 03/05/2024] [Indexed: 04/06/2024] Open
Abstract
The growth and development of chicken bones have an enormous impact on the health and production performance of chickens. However, the development pattern and genetic regulation of the chicken skeleton are poorly understood. This study aimed to evaluate metatarsal bone growth and development patterns in chickens via non-linear models, and to identify the genetic determinants of metatarsal bone traits using a genome-wide association study (GWAS) based on growth curve parameters. Data on metatarsal length (MeL) and metatarsal circumference (MeC) were obtained from 471 F2 chickens (generated by crossing broiler sires, derived from a line selected for high abdominal fat, with Baier layer dams) at 4, 6, 8, 10, and 12 weeks of age. Four non-linear models (Gompertz, Logistic, von Bertalanffy, and Brody) were used to fit the MeL and MeC growth curves. Subsequently, the estimated growth curve parameters of the mature MeL or MeC (A), time-scale parameter (b), and maturity rate (K) from the non-linear models were utilized as substitutes for the original bone data in GWAS. The Logistic and Brody models displayed the best goodness-of-fit for MeL and MeC, respectively. Single-trait and multi-trait GWASs based on the growth curve parameters of the Logistic and Brody models revealed 4 618 significant single nucleotide polymorphisms (SNPs), annotated to 332 genes, associated with metatarsal bone traits. The majority of these significant SNPs were located on Gallus gallus chromosome (GGA) 1 (167.433-176.318 Mb), GGA2 (96.791-103.543 Mb), GGA4 (65.003-83.104 Mb) and GGA6 (64.685-95.285 Mb). Notably, we identified 12 novel GWAS loci associated with chicken metatarsal bone traits, encompassing 35 candidate genes. In summary, the combination of single-trait and multi-trait GWASs based on growth curve parameters uncovered numerous genomic regions and candidate genes associated with chicken bone traits. The findings benefit an in-depth understanding of the genetic architecture underlying metatarsal growth and development in chickens.
Collapse
Affiliation(s)
- S Z Wang
- Key Laboratory of Chicken Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Harbin 150030, PR China; Key Laboratory of Animal Genetics, Breeding and Reproduction, Education Department of Heilongjiang Province, Harbin 150030, PR China; College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, PR China
| | - M D Wang
- Key Laboratory of Chicken Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Harbin 150030, PR China; Key Laboratory of Animal Genetics, Breeding and Reproduction, Education Department of Heilongjiang Province, Harbin 150030, PR China; College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, PR China
| | - J Y Wang
- Key Laboratory of Chicken Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Harbin 150030, PR China; Key Laboratory of Animal Genetics, Breeding and Reproduction, Education Department of Heilongjiang Province, Harbin 150030, PR China; College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, PR China
| | - M Yuan
- Key Laboratory of Chicken Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Harbin 150030, PR China; Key Laboratory of Animal Genetics, Breeding and Reproduction, Education Department of Heilongjiang Province, Harbin 150030, PR China; College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, PR China
| | - Y D Li
- Key Laboratory of Chicken Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Harbin 150030, PR China; Key Laboratory of Animal Genetics, Breeding and Reproduction, Education Department of Heilongjiang Province, Harbin 150030, PR China; College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, PR China
| | - P T Luo
- Fujian Sunnzer Biotechnology Development Co. Ltd, Guangze, Fujian Province 354100, PR China
| | - F Xiao
- Fujian Sunnzer Biotechnology Development Co. Ltd, Guangze, Fujian Province 354100, PR China
| | - H Li
- Key Laboratory of Chicken Genetics and Breeding, Ministry of Agriculture and Rural Affairs, Harbin 150030, PR China; Key Laboratory of Animal Genetics, Breeding and Reproduction, Education Department of Heilongjiang Province, Harbin 150030, PR China; College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, PR China.
| |
Collapse
|
6
|
Qiu Z, Cai W, Liu Q, Liu K, Liu C, Yang H, Huang R, Li P, Zhao Q. Unravelling novel and pleiotropic genes for cannon bone circumference and bone mineral density in Yorkshire pigs. J Anim Sci 2024; 102:skae036. [PMID: 38330300 PMCID: PMC10914368 DOI: 10.1093/jas/skae036] [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: 10/07/2023] [Accepted: 02/03/2024] [Indexed: 02/10/2024] Open
Abstract
Leg weakness is a prevalent health condition in pig farms. The augmentation of cannon bone circumference and bone mineral density can effectively improve limb strength in pigs and alleviate leg weakness. This study measured forelimb cannon bone circumference (fCBC) and rear limb cannon bone circumference (rCBC) using an inelastic tapeline and rear limb metatarsal area bone mineral density (raBMD) using a dual-energy X-ray absorptiometry bone density scanner. The samples of Yorkshire castrated boars were genotyped using a 50K single-nucleotide polymorphism (SNP) array. The SNP-chip data were imputed to the level of whole-genome sequencing data (iWGS). This study used iWGS data to perform genome-wide association studies and identified novel significant SNPs associated with fCBC on SSC6, SSC12, and SSC13, rCBC on SSC12 and SSC14, and raBMD on SSC7. Based on the high phenotypic and genetic correlations between CBC and raBMD, multi-trait meta-analysis was performed to identify pleiotropic SNPs. A significant potential pleiotropic quantitative trait locus (QTL) regulating both CBC and raBMD was identified on SSC15. Bayes fine mapping was used to establish the confidence intervals for these novel QTLs with the most refined confidence interval narrowed down to 56 kb (15.11 to 15.17 Mb on SSC12 for fCBC). Furthermore, the confidence interval for the potential pleiotropic QTL on SSC15 in the meta-analysis was narrowed down to 7.45 kb (137.55 to137.56 Mb on SSC15). Based on the biological functions of genes, the following genes were identified as novel regulatory candidates for different phenotypes: DDX42, MYSM1, FTSJ3, and MECOM for fCBC; SMURF2, and STC1 for rCBC; RGMA for raBMD. Additionally, RAMP1, which was determined to be located 23.68 kb upstream of the confidence interval of the QTL on SSC15 in the meta-analysis, was identified as a potential pleiotropic candidate gene regulating both CBC and raBMD. These findings offered valuable insights for identifying pathogenic genes and elucidating the genetic mechanisms underlying CBC and BMD.
Collapse
Affiliation(s)
- Zijian Qiu
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
| | - Wenwu Cai
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
| | - Qian Liu
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
| | - Kaiyue Liu
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
| | - Chenxi Liu
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
| | - Huilong Yang
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
| | - Ruihua Huang
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
- Huaian Academy, Nanjing Agricultural University, Huaian 223005, China
| | - Pinghua Li
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
- Huaian Academy, Nanjing Agricultural University, Huaian 223005, China
| | - Qingbo Zhao
- Key Laboratory in Nanjing for Evaluation and Utilization of Pigs Resources, Ministry of Agriculture and Rural Areas of China, Institute of Swine Science, Nanjing Agricultural University, Nanjing 210095, China
| |
Collapse
|
7
|
Cendron F, Cassandro M, Penasa M. Genome-wide investigation to assess copy number variants in the Italian local chicken population. J Anim Sci Biotechnol 2024; 15:2. [PMID: 38167097 PMCID: PMC10763469 DOI: 10.1186/s40104-023-00965-7] [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: 07/17/2023] [Accepted: 12/01/2023] [Indexed: 01/05/2024] Open
Abstract
BACKGROUND Copy number variants (CNV) hold significant functional and evolutionary importance. Numerous ongoing CNV studies aim to elucidate the etiology of human diseases and gain insights into the population structure of livestock. High-density chips have enabled the detection of CNV with increased resolution, leading to the identification of even small CNV. This study aimed to identify CNV in local Italian chicken breeds and investigate their distribution across the genome. RESULTS Copy number variants were mainly distributed across the first six chromosomes and primarily associated with loss type CNV. The majority of CNV in the investigated breeds were of types 0 and 1, and the minimum length of CNV was significantly larger than that reported in previous studies. Interestingly, a high proportion of the length of chromosome 16 was covered by copy number variation regions (CNVR), with the major histocompatibility complex being the likely cause. Among the genes identified within CNVR, only those present in at least five animals across breeds (n = 95) were discussed to reduce the focus on redundant CNV. Some of these genes have been associated to functional traits in chickens. Notably, several CNVR on different chromosomes harbor genes related to muscle development, tissue-specific biological processes, heat stress resistance, and immune response. Quantitative trait loci (QTL) were also analyzed to investigate potential overlapping with the identified CNVR: 54 out of the 95 gene-containing regions overlapped with 428 QTL associated to body weight and size, carcass characteristics, egg production, egg components, fat deposition, and feed intake. CONCLUSIONS The genomic phenomena reported in this study that can cause changes in the distribution of CNV within the genome over time and the comparison of these differences in CNVR of the local chicken breeds could help in preserving these genetic resources.
Collapse
Affiliation(s)
- Filippo Cendron
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale Dell'Università 16, 35020, Legnaro, PD, Italy.
| | - Martino Cassandro
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale Dell'Università 16, 35020, Legnaro, PD, Italy
- Federazione Delle Associazioni Nazionali Di Razza E Specie, Via XXIV Maggio 43, 00187, Rome, Italy
| | - Mauro Penasa
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale Dell'Università 16, 35020, Legnaro, PD, Italy
| |
Collapse
|
8
|
Rice ES, Alberdi A, Alfieri J, Athrey G, Balacco JR, Bardou P, Blackmon H, Charles M, Cheng HH, Fedrigo O, Fiddaman SR, Formenti G, Frantz LAF, Gilbert MTP, Hearn CJ, Jarvis ED, Klopp C, Marcos S, Mason AS, Velez-Irizarry D, Xu L, Warren WC. A pangenome graph reference of 30 chicken genomes allows genotyping of large and complex structural variants. BMC Biol 2023; 21:267. [PMID: 37993882 PMCID: PMC10664547 DOI: 10.1186/s12915-023-01758-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2023] [Accepted: 11/02/2023] [Indexed: 11/24/2023] Open
Abstract
BACKGROUND The red junglefowl, the wild outgroup of domestic chickens, has historically served as a reference for genomic studies of domestic chickens. These studies have provided insight into the etiology of traits of commercial importance. However, the use of a single reference genome does not capture diversity present among modern breeds, many of which have accumulated molecular changes due to drift and selection. While reference-based resequencing is well-suited to cataloging simple variants such as single-nucleotide changes and short insertions and deletions, it is mostly inadequate to discover more complex structural variation in the genome. METHODS We present a pangenome for the domestic chicken consisting of thirty assemblies of chickens from different breeds and research lines. RESULTS We demonstrate how this pangenome can be used to catalog structural variants present in modern breeds and untangle complex nested variation. We show that alignment of short reads from 100 diverse wild and domestic chickens to this pangenome reduces reference bias by 38%, which affects downstream genotyping results. This approach also allows for the accurate genotyping of a large and complex pair of structural variants at the K feathering locus using short reads, which would not be possible using a linear reference. CONCLUSIONS We expect that this new paradigm of genomic reference will allow better pinpointing of exact mutations responsible for specific phenotypes, which will in turn be necessary for breeding chickens that meet new sustainability criteria and are resilient to quickly evolving pathogen threats.
Collapse
Affiliation(s)
- Edward S Rice
- Bond Life Sciences Center, University of Missouri, Columbia, MO, USA
- Faculty of Veterinary Medicine, Ludwig-Maximilians-Universität, Munich, Germany
| | - Antton Alberdi
- Center for Evolutionary Hologenomics, Globe Institute, University of Copenhagen (UCPH), Copenhagen, Denmark
| | - James Alfieri
- Department of Ecology & Evolutionary Biology, Texas A&M University, College Station, TX, USA
| | - Giridhar Athrey
- Department of Poultry Science, Texas A&M University, College Station, TX, USA
| | - Jennifer R Balacco
- Vertebrate Genome Laboratory, The Rockefeller University, New York, NY, USA
| | - Philippe Bardou
- Sigenae, GenPhySE, Université de Toulouse, INRAE, ENVT, Castanet Tolosan, 31326, France
| | - Heath Blackmon
- Department of Biology, Texas A&M University, College Station, TX, USA
| | - Mathieu Charles
- University Paris-Saclay, INRAE, AgroParisTech, GABI, Sigenae, Jouy-en-Josas, France
| | - Hans H Cheng
- Avian Disease and Oncology Laboratory, USDA, ARS, USNPRC, East Lansing, MI, USA
| | - Olivier Fedrigo
- Vertebrate Genome Laboratory, The Rockefeller University, New York, NY, USA
| | | | - Giulio Formenti
- Vertebrate Genome Laboratory, The Rockefeller University, New York, NY, USA
| | - Laurent A F Frantz
- Faculty of Veterinary Medicine, Ludwig-Maximilians-Universität, Munich, Germany
- School of Biological and Behavioural Sciences, Queen Mary University of London, London, E1 4DQ, UK
| | - M Thomas P Gilbert
- Center for Evolutionary Hologenomics, Globe Institute, University of Copenhagen (UCPH), Copenhagen, Denmark
| | - Cari J Hearn
- Avian Disease and Oncology Laboratory, USDA, ARS, USNPRC, East Lansing, MI, USA
| | - Erich D Jarvis
- Vertebrate Genome Laboratory, The Rockefeller University, New York, NY, USA
- The Howard Hughes Medical Institute, Chevy Chase, MD, USA
| | - Christophe Klopp
- Sigenae, Genotoul Bioinfo, MIAT UR875, INRAE, Castanet Tolosan, France
| | - Sofia Marcos
- Center for Evolutionary Hologenomics, Globe Institute, University of Copenhagen (UCPH), Copenhagen, Denmark
- Applied Genomics and Bioinformatics, University of the Basque Country (UPV/EHU), Leioa, Bilbao, Spain
| | | | | | - Luohao Xu
- Key Laboratory of Freshwater Fish Reproduction and Development (Ministry of Education), Key Laboratory of Aquatic Science of Chongqing, School of Life Sciences, Southwest University, Chongqing, 400715, China
| | - Wesley C Warren
- Department of Animal Sciences, University of Missouri, Columbia, MO, USA.
| |
Collapse
|
9
|
Shi H, Li T, Su M, Wang H, Li Q, Lang X, Ma Y. Identification of copy number variation in Tibetan sheep using whole genome resequencing reveals evidence of genomic selection. BMC Genomics 2023; 24:555. [PMID: 37726692 PMCID: PMC10510117 DOI: 10.1186/s12864-023-09672-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2022] [Accepted: 09/12/2023] [Indexed: 09/21/2023] Open
Abstract
BACKGROUND Copy number variation (CNV) is an important source of structural variation in the mammalian genome. CNV assays present a new method to explore the genomic diversity of environmental adaptations in animals and plants and genes associated with complex traits. In this study, the genome-wide CNV distribution characteristics of 20 Tibetan sheep from two breeds (10 Oula sheep and 10 Panou sheep) were analysed using whole-genome resequencing to investigate the variation in the genomic structure of Tibetan sheep during breeding. RESULTS CNVs were detected using CNVnator, and the overlapping regions of CNVs between individual sheep were combined. Among them, a total of 60,429 CNV events were detected between the indigenous sheep breed (Oula) and the synthetic sheep breed (Panou). After merging the overlapping CNVs, 4927 CNV regions (CNVRs) were finally obtained. Of these, 4559 CNVRs were shared by two breeds, and there were 368 differential CNVRs. Deletion events have a higher percentage of occurrences than duplication events. Functional enrichment analysis showed that the shared CNVRs were significantly enriched in 163 GO terms and 62 KEGG pathways, which were mainly associated with organ development, neural regulation, immune regulation, digestion and metabolism. In addition, 140 QTLs overlapped with some of the CNVRs at more than 1 kb, such as average daily gain QTL, body weight QTL, and total lambs born QTL. Many of the CNV-overlapping genes such as PPP3CA, SSTR1 and FASN, overlap with the average daily weight gain and carcass weight QTL regions. Moreover, VST analysis showed that XIRP2, ABCB1, CA1, ASPA and EEF2 differed significantly between the synthetic breed and local sheep breed. The duplication of the ABCB1 gene may be closely related to adaptation to the plateau environment in Panou sheep, which deserves further study. Additionally, cluster analysis, based on all individuals, showed that the CNV clustering could be divided into two origins, indicating that some Tibetan sheep CNVs are likely to arise independently in different populations and contribute to population differences. CONCLUSIONS Collectively, we demonstrated the genome-wide distribution characteristics of CNVs in Panou sheep by whole genome resequencing. The results provides a valuable genetic variation resource and help to understand the genetic characteristics of Tibetan sheep. This study also provides useful information for the improvement and breeding of Tibetan sheep in the future.
Collapse
Affiliation(s)
- Huibin Shi
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China
- Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070, China
- College of Animal Science & Technology, Henan University of Animal Husbandry and Economy, Zhengzhou, 450046, China
| | - Taotao Li
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China
- Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070, China
| | - Manchun Su
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China
- Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070, China
| | - Huihui Wang
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China
- Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070, China
| | - Qiao Li
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China
- Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070, China
| | - Xia Lang
- Institute of Animal & Pasture Science and Green Agriculture, Gansu Academy of Agricultural Science, Lanzhou, 730070, China
| | - Youji Ma
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070, China.
- Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070, China.
| |
Collapse
|
10
|
Meng G, La Y, Bao Q, Wu X, Ma X, Huang C, Chu M, Liang C, Yan P. Early Growth and Development and Nonlinear Model Fitting Analysis of Ashidan Yak. Animals (Basel) 2023; 13:ani13091545. [PMID: 37174583 PMCID: PMC10177478 DOI: 10.3390/ani13091545] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Revised: 04/27/2023] [Accepted: 05/03/2023] [Indexed: 05/15/2023] Open
Abstract
Understanding animal growth plays an important role in improving animal genetics and breeding. In order to explore the early growth and development law of Ashidan yak, the body weight (BW), wither height (WH), body oblique length (BL) and chest girth (CG) of 260 female Ashidan yaks were measured. These individuals grew under grazing conditions, and growth traits were measured at 6, 12, 18 and 30 months of age. Then the absolute growth and relative growth of Ashidan yak were calculated, and five nonlinear models (Logistic model, Gompertz model, Brody model, von Bertalanffy model and Richards model) were used to fit the growth curve of Ashidan yak. The fitting effect of the model was evaluated according to MSE, AIC and BIC. The results showed that the growth rate of Ashidan yak was the fastest from 12 to 18 months old, and the growth was slow or even stagnant from 6 to 12 months old. The AIC and BIC values of the Richards model were the lowest among the five models, with an AIC value of 4543.98 and a BIC value of 4563.19. The Richards model estimated body weight at 155.642 kg. In summary, the growth rate of female Ashidan yak changes with the seasons, growing faster in warm seasons and slower in cold seasons. Richards model is the best model to describe the growth curve of female Ashidan yak in five nonlinear models.
Collapse
Affiliation(s)
- Guangyao Meng
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Yongfu La
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Qi Bao
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Xiaoyun Wu
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Xiaoming Ma
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Chun Huang
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Min Chu
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Chunnian Liang
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| | - Ping Yan
- Lanzhou Institute of Husbandry and Pharmaceutical Sciences, Chinese Academy of Agricultural Sciences, Lanzhou 730050, China
- Key Laboratory of Animal Genetics and Breeding on Tibetan Plateau, Ministry of Agriculture and Rural Affairs, Lanzhou 730050, China
| |
Collapse
|
11
|
Shi H, Li T, Su M, Wang H, Li Q, Lang X, Ma Y. Whole genome sequencing revealed genetic diversity, population structure, and selective signature of Panou Tibetan sheep. BMC Genomics 2023; 24:50. [PMID: 36707771 PMCID: PMC9883975 DOI: 10.1186/s12864-023-09146-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Accepted: 01/20/2023] [Indexed: 01/29/2023] Open
Abstract
BACKGROUND The detection of selective traits in different populations can not only reveal current mechanisms of artificial selection for breeding, but also provide new insights into phenotypic variation in new varieties and the search for genes associated with important traits. Panou sheep is a cultivated breed of Tibetan sheep in China with stable genetic performance, consistent appearance and fast growth and development after decades of artificial selection and cultivation. Due to long-term adaptation to the high altitude, cold and hypoxic environment in the plateau area, they may have formed a unique gene pool that is different from other Tibetan sheep breeds. To explore the genetic resources of Panou sheep, we used next-generation sequencing technology for the first time to investigate the genome-wide population structure, genetic diversity, and candidate signatures of positive selection in Panou sheep. RESULTS Comparative genomic analysis with the closely related species Oula sheep (a native breed of Tibetan sheep in China) was used to screen the population selection signal of Panou sheep. Principal component analysis and neighbor joining tree showed that Panou sheep and Oula sheep had differences in population differentiation. Furthermore, analyses of population structure, they came from the same ancestor, and when K = 2, the two populations could be distinguished. Panou sheep exhibit genetic diversity comparable to Oula sheep, as shown by observed heterozygosity, expected heterozygosity and runs of homozygosity. Genome-wide scanning using the Fst and π ratio methods revealed a list of potentially selected related genes in Panou sheep compared to Oula sheep, including histone deacetylase 9 (HDAC9), protein tyrosine kinase 2 (PTK2), microphthalmia-related transcription factor (MITF), vesicular amine transporter 1 (VAT1), trichohyalin-like 1 (TCHHL1), amine oxidase, copper containing 3 (AOC3), interferon-inducible protein 35 (IFI35). CONCLUSIONS The results suggest that traits related to growth and development and plateau adaptation may be selection targets for the domestication and breeding improvement of Tibetan sheep. This study provides the fundamental footprints for Panou sheep breeding and management.
Collapse
Affiliation(s)
- Huibin Shi
- grid.411734.40000 0004 1798 5176College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070 China ,Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070 China
| | - Taotao Li
- grid.411734.40000 0004 1798 5176College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070 China ,Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070 China
| | - Manchun Su
- grid.411734.40000 0004 1798 5176College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070 China ,Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070 China
| | - Huihui Wang
- grid.411734.40000 0004 1798 5176College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070 China ,Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070 China
| | - Qiao Li
- grid.411734.40000 0004 1798 5176College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070 China ,Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070 China
| | - Xia Lang
- grid.464277.40000 0004 0646 9133Institute of Animal & Pasture Science and Green Agriculture, Gansu Academy of Agricultural Science, Lanzhou, 730070 China
| | - Youji Ma
- grid.411734.40000 0004 1798 5176College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, 730070 China ,Gansu Key Laboratory of Animal Generational Physiology and Reproductive Regulation, Lanzhou, 730070 China
| |
Collapse
|
12
|
Shibak A, Maghsoudi A, Rokouei M, Farhangfar H, Faraji-Arough H. Investigation of egg production curve in ostrich using nonlinear functions. Poult Sci 2022; 102:102333. [PMID: 36463766 PMCID: PMC9719868 DOI: 10.1016/j.psj.2022.102333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Revised: 10/16/2022] [Accepted: 10/18/2022] [Indexed: 11/21/2022] Open
Abstract
In most countries, ostrich farming is considered a developing branch of the efficient poultry industry. The profitability of ostrich farm requires specific consideration of productions features such as the female fertility, egg production, hatchability, and growth performance. Hence, this study aimed to fit nonlinear functions to describe the ostrich egg production pattern to achieve the most appropriate and recommendable mathematical function for future studies. For this purpose, 14,507 daily records of 184 female ostriches in 5 production seasons (periods) during 2016 to 2021 were used. Five nonlinear functions including Incomplete gamma (Wood function), Corrected gamma (McNally), nonlinear Logistic (Yang), Logistic (Nelder), and Lokhorst were fitted for modeling the egg production curve in ostrich. The goodness of fit criteria's including Mean Square Error (MSE), Likelihood Ratio Test (LRT), Akaike Information Criterion (AIC), and Bayesian Information Criterion (BIC) were used to evaluate and selection of the best function. The results indicated that the Wood and the McNally functions with a slight difference in all fitting criteria were the best-fitted functions and the Yang function with the highest values of MSE, LRT, AIC, BIC, were the most inappropriate function to describe the ostrich egg production curve. The McNally and the Wood can be recommended as appropriate functions to describe egg production during 5 production seasons in the studied ostrich flock.
Collapse
Affiliation(s)
- Abbas Shibak
- Department of Animal Science, Faculty of Agriculture, University of Zabol, Zabol, Iran
| | - Ali Maghsoudi
- Department of Animal Science, Faculty of Agriculture, University of Zabol, Zabol, Iran,Department of Animal Science, Faculty of Agriculture, Tarbiat Modares University, Tehran, Iran,Corresponding author:
| | - Mohammad Rokouei
- Department of Animal Science, Faculty of Agriculture, University of Zabol, Zabol, Iran
| | - Homayoun Farhangfar
- Department of Animal Science, Faculty of Agriculture, University of Birjand, Birjand, Iran
| | - Hadi Faraji-Arough
- Department of Ostrich, Special Domestic Animals Institute, Research Institute of Zabol, Zabol, Iran
| |
Collapse
|