1
|
Zhang Z, Yang W, Zhu T, Wang L, Zhao X, Zhao G, Qu L, Jia Y. Genetic Parameter Estimation and Whole Sequencing Analysis of the Genetic Architecture of Chicken Keel Bending. Front Genet 2022; 13:833132. [PMID: 35401685 PMCID: PMC8984200 DOI: 10.3389/fgene.2022.833132] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2021] [Accepted: 02/24/2022] [Indexed: 11/18/2022] Open
Abstract
Bone health is particularly important for high-yielding commercial layer chickens. The keel of poultry is an extension of the abdomen side of the sternum along the sagittal plane and is one of the most important bones. In this study, the keel phenotype of White Leghorns laying hen flocks showed significant individual differences. To clarify its genetic mechanism, we first estimated the heritability of keel bend (KB) in White Leghorn, recorded the production performance of the chicken flock, examined the blood biochemical indexes and bone quality in KB and keel normal (KN) chickens, and performed whole-genome pooled sequencing in KB and KN chickens. We then performed selection elimination analysis to determine the genomic regions that may affect the keel phenotypes. The results show that KB is a medium heritability trait. We found that cage height had a significant effect on the KB (p < 0.01). At 48 weeks, there were significant differences in the number of eggs, the number of normal eggs, and eggshell strength (p < 0.05). The content of parathyroid hormone was lower (p < 0.01) and that of calcitonin was higher (p < 0.01) in KB chickens than in KN chickens. The differences in bone mineral density, bone strength, and bone cortical thickness of the humerus and femur were extremely significant (p < 0.01), with all being lower in KB chickens than in KN chickens. In addition, the bones of KB chickens contained more fat organization. A total of 128 genes were identified in selective sweep regions. We identified 10 important candidate genes: ACP5, WNT1, NFIX, CNN1, CALR, FKBP11, TRAPPC5, MAP2K7, RELA, and ENSGALG00000047166. Among the significantly enriched Kyoto Encyclopedia of Genes and Genomes pathways found, we identifed two bone-related pathways, one involving “osteoclast differentiation” and the other the “MAPK signaling pathway.” These results may help us better understand the molecular mechanism of bone traits in chickens and other birds and provide new insights for the genetic breeding of chickens.
Collapse
Affiliation(s)
- Zhihao Zhang
- Institute of Animal Sciences, Chinese Academy of Agricultural Science, Beijing, China
| | - Weifang Yang
- Beijing General Station of Animal Husbandry, Beijing, China
| | - Tao Zhu
- State Key Laboratory of Animal Nutrition, Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Liang Wang
- Beijing General Station of Animal Husbandry, Beijing, China
| | - Xiaoyu Zhao
- Hebei Dawu Poultry Breeding Co., Ltd., Hebei, China
| | | | - Lujiang Qu
- State Key Laboratory of Animal Nutrition, Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
- *Correspondence: Lujiang Qu, ; Yaxiong Jia,
| | - Yaxiong Jia
- Institute of Animal Sciences, Chinese Academy of Agricultural Science, Beijing, China
- *Correspondence: Lujiang Qu, ; Yaxiong Jia,
| |
Collapse
|
2
|
Magris G, Marroni F, D’Agaro E, Vischi M, Chiabà C, Scaglione D, Kijas J, Messina M, Tibaldi E, Morgante M. ddRAD-seq reveals the genetic structure and detects signals of selection in Italian brown trout. Genet Sel Evol 2022; 54:8. [PMID: 35100964 PMCID: PMC8805291 DOI: 10.1186/s12711-022-00698-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2021] [Accepted: 01/14/2022] [Indexed: 01/29/2023] Open
Abstract
Background Brown trout is one of the most widespread fresh-water fish species in Europe. The evolutionary history of and phylogenetic relationships between brown trout populations are complex, and this is especially true for Italian populations, which are heavily influenced in different ways by stocking practices. The characterization of the genetic structure of Italian brown trout populations may give information on the risk of losing endemic Italian populations due to lack of genetic diversity or to admixture with stocking populations. The identification of signatures of selection, and the information deriving from dense genotyping data will help genotype-informed breeding programs. We used a ddRAD-seq approach to obtain more than 100,000 single nucleotide polymorphisms (SNPs), and to characterize the population structure and signatures of selection in 90 brown trout samples. Results Italian brown trout populations are genetically differentiated, although the stocking practices have introduced strong admixture in endemic Italian trout, especially with the Atlantic lineage. Most of the analysed populations showed high levels of kinship and inbreeding. We detected putative signatures of selection using different approaches, and investigated if the regions were enriched for functional categories. Several regions putatively under selection and characterized by a reduction in heterozygosity across all the studied populations are enriched for genes involved in the response to viral infections. Conclusions Our results, which show evidence of admixture with the Atlantic lineage (commonly used for stocking), confirm the need for controlling stocking practices, in order to avoid the erosion of the endemic gene pool; given the apparently high levels of kinship and inbreeding in local populations, our results also show the need to take action for increasing gene diversity. In addition, we used the genetically-distinct lineages to detect signatures of selection and we identified putative signatures of selection in several regions associated with resistance to infectious diseases. These constitute candidate regions for the study of resistance to infections in wild and farmed trout. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-022-00698-7.
Collapse
|
3
|
Runs of homozygosity analysis reveals consensus homozygous regions affecting production traits in Chinese Simmental beef cattle. BMC Genomics 2021; 22:678. [PMID: 34548021 PMCID: PMC8454143 DOI: 10.1186/s12864-021-07992-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 09/07/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Genomic regions with a high frequency of runs of homozygosity (ROH) are related to important traits in farm animals. We carried out a comprehensive analysis of ROH and evaluated their association with production traits using the BovineHD (770 K) SNP array in Chinese Simmental beef cattle. RESULTS We detected a total of 116,953 homozygous segments with 2.47Gb across the genome in the studied population. The average number of ROH per individual was 99.03 and the average length was 117.29 Mb. Notably, we detected 42 regions with a frequency of more than 0.2. We obtained 17 candidate genes related to body size, meat quality, and reproductive traits. Furthermore, using Fisher's exact test, we found 101 regions were associated with production traits by comparing high groups with low groups in terms of production traits. Of those, we identified several significant regions for production traits (P < 0.05) by association analysis, within which candidate genes including ECT2, GABRA4, and GABRB1 have been previously reported for those traits in beef cattle. CONCLUSIONS Our study explored ROH patterns and their potential associations with production traits in beef cattle. These results may help to better understand the association between production traits and genome homozygosity and offer valuable insights into managing inbreeding by designing reasonable breeding programs in farm animals.
Collapse
|
4
|
de Sousa MAP, de Athayde FRF, Maldonado MBC, de Lima AO, Fortes MRS, Lopes FL. Single nucleotide polymorphisms affect miRNA target prediction in bovine. PLoS One 2021; 16:e0249406. [PMID: 33882076 PMCID: PMC8059806 DOI: 10.1371/journal.pone.0249406] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Accepted: 03/17/2021] [Indexed: 02/06/2023] Open
Abstract
Single nucleotide polymorphisms (SNPs) can have significant effects on phenotypic characteristics in cattle. MicroRNAs (miRNAs) are small, non-coding RNAs that act as post-transcriptional regulators by binding them to target mRNAs. In the present study, we scanned ~56 million SNPs against 1,064 bovine miRNA sequences and analyzed, in silico, their possible effects on target binding prediction, primary miRNA formation, association with QTL regions and the evolutionary conservation for each SNP locus. Following target prediction, we show that 71.6% of miRNA predicted targets were altered as a consequence of SNPs located within the seed region of the mature miRNAs. Next, we identified variations in the Minimum Free Energy (MFE), which represents the capacity to alter molecule stability and, consequently, miRNA maturation. A total of 48.6% of the sequences analyzed showed values within those previously reported as sufficient to alter miRNA maturation. We have also found 131 SNPs in 46 miRNAs, with altered target prediction, occurring in QTL regions. Lastly, analysis of evolutionary conservation scores for each SNP locus suggested that they have a conserved biological function through the evolutionary process. Our results suggest that SNPs in microRNAs have the potential to affect bovine phenotypes and could be of great value for genetic improvement studies, as well as production.
Collapse
Affiliation(s)
- Marco Antônio Perpétuo de Sousa
- Department of Production and Animal Health, São Paulo State University (Unesp), School of Veterinary Medicine, Araçatuba, São Paulo, Brazil
| | - Flavia Regina Florêncio de Athayde
- Department of Production and Animal Health, São Paulo State University (Unesp), School of Veterinary Medicine, Araçatuba, São Paulo, Brazil
| | | | - Andressa Oliveira de Lima
- Department of Production and Animal Health, São Paulo State University (Unesp), School of Veterinary Medicine, Araçatuba, São Paulo, Brazil
| | - Marina Rufino S. Fortes
- School of Chemistry and Molecular Bioscience, The University of Queensland, Brisbane, Queensland, Australia
| | - Flavia Lombardi Lopes
- Department of Production and Animal Health, São Paulo State University (Unesp), School of Veterinary Medicine, Araçatuba, São Paulo, Brazil
| |
Collapse
|
5
|
Ghoreishifar SM, Rochus CM, Moghaddaszadeh-Ahrabi S, Davoudi P, Salek Ardestani S, Zinovieva NA, Deniskova TE, Johansson AM. Shared Ancestry and Signatures of Recent Selection in Gotland Sheep. Genes (Basel) 2021; 12:genes12030433. [PMID: 33802939 PMCID: PMC8002741 DOI: 10.3390/genes12030433] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 03/01/2021] [Accepted: 03/10/2021] [Indexed: 12/13/2022] Open
Abstract
Gotland sheep, a breed native to Gotland, Sweden (an island in the Baltic Sea), split from the Gute sheep breed approximately 100 years ago, and since, has probably been crossed with other breeds. This breed has recently gained popularity, due to its pelt quality. This study estimates the shared ancestors and identifies recent selection signatures in Gotland sheep using 600 K single nucleotide polymorphism (SNP) genotype data. Admixture analysis shows that the Gotland sheep is a distinct breed, but also has shared ancestral genomic components with Gute (~50%), Karakul (~30%), Romanov (~20%), and Fjällnäs (~10%) sheep breeds. Two complementary methods were applied to detect selection signatures: A Bayesian population differentiation FST and an integrated haplotype homozygosity score (iHS). Our results find that seven significant SNPs (q-value < 0.05) using the FST analysis and 55 significant SNPs (p-value < 0.0001) using the iHS analysis. Of the candidate genes that contain significant markers, or are in proximity to them, we identify several belongings to the keratin genes, RXFP2, ADCY1, ENOX1, USF2, COX7A1, ARHGAP28, CRYBB2, CAPNS1, FMO3, and GREB1. These genes are involved in wool quality, polled and horned phenotypes, fertility, twining rate, meat quality, and growth traits. In summary, our results provide shared founders of Gotland sheep and insight into genomic regions maintained under selection after the breed was formed. These results contribute to the detection of candidate genes and QTLs underlying economic traits in sheep.
Collapse
Affiliation(s)
- Seyed Mohammad Ghoreishifar
- Department of Animal Science, University College of Agriculture and Natural Resources, University of Tehran, Karaj 31587-11167, Iran;
| | - Christina Marie Rochus
- Animal Breeding and Genomics, Wageningen University and Research, P.O. Box 338, 6700 AH Wageningen, The Netherlands;
| | - Sima Moghaddaszadeh-Ahrabi
- Department of Animal Science, Faculty of Agriculture and Natural Resources, Islamic Azad University, Tabriz Branch, Tabriz 5157944533, Iran;
| | - Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS B2N5E3, Canada; (P.D.); (S.S.A.)
| | - Siavash Salek Ardestani
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS B2N5E3, Canada; (P.D.); (S.S.A.)
| | - Natalia A. Zinovieva
- L.K. Ernst Federal Research Center for Animal Husbandry, 142132 Podolsk, Russia; (N.A.Z.); (T.E.D.)
| | - Tatiana E. Deniskova
- L.K. Ernst Federal Research Center for Animal Husbandry, 142132 Podolsk, Russia; (N.A.Z.); (T.E.D.)
| | - Anna M. Johansson
- Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, SE-75007 Uppsala, Sweden
- Correspondence:
| |
Collapse
|
6
|
Genomic Scan for Selection Signature Reveals Fat Deposition in Chinese Indigenous Sheep with Extreme Tail Types. Animals (Basel) 2020; 10:ani10050773. [PMID: 32365604 PMCID: PMC7278473 DOI: 10.3390/ani10050773] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2020] [Revised: 04/21/2020] [Accepted: 04/22/2020] [Indexed: 02/03/2023] Open
Abstract
Simple Summary According to the tail types, sheep can be briefly classified into three groups (fat-tailed, fat-rumped, and thin-tailed sheep). In this study, we used these three typical breeds from Chinese indigenous sheep breeds to perform a genome scan for selective sweeps using Ovine Infinium HD SNP BeadChip genotype data. Results showed that 25 genomic regions exhibited selection signals and harbored 73 positional candidate genes. These genes were documented not only to be associated with tail fat formation, but also be related to reproduction, body conformation, and appearance. Our findings contributed to understanding genetic basis of fat deposition in sheep tail and provide a reference for developing new sheep breeds with an ideal tail type. Abstract It is a unique feature that fat can be deposited in sheep tails and rumps. To elucidate the genetic mechanism underlying this trait, we collected 120 individuals from three Chinese indigenous sheep breeds with extreme tail types, namely large fat-tailed sheep (n = 40), Altay sheep (n = 40), and Tibetan sheep (n = 40), and genotyped them using the Ovine Infinium HD SNP BeadChip. Then genomic scan for selection signatures was performed using the hapFLK. In total, we identified 25 genomic regions exhibiting evidence of having been under selection. Bioinformatic analysis of the genomic regions showed that selection signatures related to multiple candidate genes had a demonstrated role in phenotypic variation. Nine genes have documented association with sheep tail types, including WDR92, TBX12, WARS2, BMP2, VEGFA, PDGFD, HOXA10, ALX4, and ETAA1. Moreover, a number of genes were of particular interest, including RXFP2 associated with the presence/absence and morphology of horns; MITF involved in coat color; LIN52 and SYNDIG1L related to the number of teats; MSRB3 gene associated with ear sizes; LTBP2 considered as a positional candidate genes for number of ribs; JAZF1 regulating lipid metabolism; PGRMC2, SPAG17, TSHR, GTF2A1, and LARP1B implicated with reproductive traits. Our findings provide insights into fat tail formation and a reference for carrying out molecular breeding and conservation in sheep.
Collapse
|
7
|
Dolebo AT, Khayatzadeh N, Melesse A, Wragg D, Rekik M, Haile A, Rischkowsky B, Rothschild MF, Mwacharo JM. Genome-wide scans identify known and novel regions associated with prolificacy and reproduction traits in a sub-Saharan African indigenous sheep (Ovis aries). Mamm Genome 2019; 30:339-352. [PMID: 31758253 PMCID: PMC6884434 DOI: 10.1007/s00335-019-09820-5] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2019] [Accepted: 11/12/2019] [Indexed: 02/05/2023]
Abstract
Maximizing the number of offspring born per female is a key functionality trait in commercial- and/or subsistence-oriented livestock enterprises. Although the number of offspring born is closely associated with female fertility and reproductive success, the genetic control of these traits remains poorly understood in sub-Saharan Africa livestock. Using selection signature analysis performed on Ovine HD BeadChip data from the prolific Bonga sheep in Ethiopia, 41 candidate regions under selection were identified. The analysis revealed one strong selection signature on a candidate region on chromosome X spanning BMP15, suggesting this to be the primary candidate prolificacy gene in the breed. The analysis also identified several candidate regions spanning genes not reported before in prolific sheep but underlying fertility and reproduction in other species. The genes associated with female reproduction traits included SPOCK1 (age at first oestrus), GPR173 (mediator of ovarian cyclicity), HB-EGF (signalling early pregnancy success) and SMARCAL1 and HMGN3a (regulate gene expression during embryogenesis). The genes involved in male reproduction were FOXJ1 (sperm function and successful fertilization) and NME5 (spermatogenesis). We also observed genes such as PKD2L2, MAGED1 and KDM3B, which have been associated with diverse fertility traits in both sexes of other species. The results confirm the complexity of the genetic mechanisms underlying reproduction while suggesting that prolificacy in the Bonga sheep, and possibly African indigenous sheep is partly under the control of BMP15 while other genes that enhance male and female fertility are essential for reproductive fitness.
Collapse
Affiliation(s)
- Asrat Tera Dolebo
- Southern Agricultural Research Institute (SARI), P.O. Box 06, Hawassa, Ethiopia
- Department of Animal and Range Sciences, Hawassa University, P.O Box 5, Hawassa, Ethiopia
| | - Negar Khayatzadeh
- Department of Sustainable Agricultural Systems, Division of Livestock Sciences, University of Natural Resources and Life Sciences (BOKU), Gregor-Mendel-Strasse, 1180, Vienna, Austria
| | - Aberra Melesse
- Department of Animal and Range Sciences, Hawassa University, P.O Box 5, Hawassa, Ethiopia
| | - David Wragg
- Centre for Tropical Livestock Genetics and Health, The Roslin Institute, Edinburgh, UK
| | - Mourad Rekik
- Small Ruminant Genomics, International Centre for Agricultural Research in the Dry Areas (ICARDA), P.O. Box 5689, Addis Ababa, Ethiopia
| | - Aynalem Haile
- Small Ruminant Genomics, International Centre for Agricultural Research in the Dry Areas (ICARDA), P.O. Box 5689, Addis Ababa, Ethiopia
| | - Barbara Rischkowsky
- Small Ruminant Genomics, International Centre for Agricultural Research in the Dry Areas (ICARDA), P.O. Box 5689, Addis Ababa, Ethiopia
| | - Max F Rothschild
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011-3150, USA
| | - Joram M Mwacharo
- Small Ruminant Genomics, International Centre for Agricultural Research in the Dry Areas (ICARDA), P.O. Box 5689, Addis Ababa, Ethiopia.
| |
Collapse
|
8
|
Yuan Z, Li W, Li F, Yue X. Selection signature analysis reveals genes underlying sheep milking performance. Arch Anim Breed 2019; 62:501-508. [PMID: 31807661 PMCID: PMC6859915 DOI: 10.5194/aab-62-501-2019] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2019] [Accepted: 07/15/2019] [Indexed: 11/30/2022] Open
Abstract
Sheep milk is the most important feed resource for newborn lambs and an important food resource for humans. Sheep milk production and ingredients are influenced by genetic and environmental factors. In this study, we implemented selection signature analysis using Illumina Ovine SNP50 BeadChip data of 78 meat Lacaune and 103 milk Lacaune sheep, which have similar genetic backgrounds, from the Sheep HapMap project to identify candidate genes related to ovine milk traits. Since different methods can detect different variation types and complement each other, we used a haplotype-based method (hapFLK) to implement selection signature analysis. The results revealed six selection signature regions showing signs of being selected ( P < 0.001 ): chromosomes 1, 2, 3, 6, 13 and 18. In addition, 38 quantitative trait loci (QTLs) related to sheep milk performance were identified in selection signature regions, which contain 334 candidate genes. Of those, SUCNR1 (succinate receptor 1) and PPARGC1A (PPARG coactivator 1 alpha) may be the most significant genes that affect sheep milking performance, which supply a significant indication for future studies to investigate candidate genes that play an important role in milk production and quality.
Collapse
Affiliation(s)
- Zehu Yuan
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
| | - Wanhong Li
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
| | - Fadi Li
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
- Engineering Laboratory of Sheep Breeding and Reproduction
Biotechnology in Gansu Province, Minqin, 733300, P. R. China
| | - Xiangpeng Yue
- State Key Laboratory of Grassland Agro-Ecosystems; Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural
Affairs; College of Pastoral Agriculture Science and Technology, Lanzhou
University, Lanzhou, 730020, P. R. China
| |
Collapse
|
9
|
Ahbara A, Bahbahani H, Almathen F, Al Abri M, Agoub MO, Abeba A, Kebede A, Musa HH, Mastrangelo S, Pilla F, Ciani E, Hanotte O, Mwacharo JM. Genome-Wide Variation, Candidate Regions and Genes Associated With Fat Deposition and Tail Morphology in Ethiopian Indigenous Sheep. Front Genet 2019; 9:699. [PMID: 30687385 PMCID: PMC6334744 DOI: 10.3389/fgene.2018.00699] [Citation(s) in RCA: 46] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2018] [Accepted: 12/13/2018] [Indexed: 12/03/2022] Open
Abstract
Variations in body weight and in the distribution of body fat are associated with feed availability, thermoregulation, and energy reserve. Ethiopia is characterized by distinct agro-ecological and human ethnic farmer diversity of ancient origin, which have impacted on the variation of its indigenous livestock. Here, we investigate autosomal genome-wide profiles of 11 Ethiopian indigenous sheep populations using the Illumina Ovine 50 K SNP BeadChip assay. Sheep from the Caribbean, Europe, Middle East, China, and western, northern and southern Africa were included to address globally, the genetic variation and history of Ethiopian populations. Population relationship and structure analysis separated Ethiopian indigenous fat-tail sheep from their North African and Middle Eastern counterparts. It indicates two main genetic backgrounds and supports two distinct genetic histories for African fat-tail sheep. Within Ethiopian sheep, our results show that the short fat-tail sheep do not represent a monophyletic group. Four genetic backgrounds are present in Ethiopian indigenous sheep but at different proportions among the fat-rump and the long fat-tail sheep from western and southern Ethiopia. The Ethiopian fat-rump sheep share a genetic background with Sudanese thin-tail sheep. Genome-wide selection signature analysis identified eight putative candidate regions spanning genes influencing growth traits and fat deposition (NPR2, HINT2, SPAG8, INSR), development of limbs and skeleton, and tail formation (ALX4, HOXB13, BMP4), embryonic development of tendons, bones and cartilages (EYA2, SULF2), regulation of body temperature (TRPM8), body weight and height variation (DIS3L2), control of lipogenesis and intracellular transport of long-chain fatty acids (FABP3), the occurrence and morphology of horns (RXFP2), and response to heat stress (DNAJC18). Our findings suggest that Ethiopian fat-tail sheep represent a uniquely admixed but distinct genepool that presents an important resource for understanding the genetic control of skeletal growth, fat metabolism and associated physiological processes.
Collapse
Affiliation(s)
- Abulgasim Ahbara
- School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
- Department of Zoology, Faculty of Sciences, Misurata University, Misurata, Libya
| | - Hussain Bahbahani
- Department of Biological Sciences, Faculty of Science, Kuwait University, Safat, Kuwait
| | - Faisal Almathen
- Department of Veterinary Public Health and Animal Husbandry, College of Veterinary Medicine, King Faisal University, Al-Ahsa, Saudi Arabia
| | - Mohammed Al Abri
- Department of Animal and Veterinary Sciences, College of Agriculture and Marine Sciences, Sultan Qaboos University, Muscat, Oman
| | | | - Ayelle Abeba
- Debre Berhan Research Centre, Debre Berhan, Ethiopia
| | - Adebabay Kebede
- Amhara Regional Agricultural Research Institute, Bahir Dar, Ethiopia
- LiveGene, International Livestock Research Institute, Addis Ababa, Ethiopia
| | - Hassan Hussein Musa
- Faculty of Medical Laboratory Sciences, University of Khartoum, Khartoum, Sudan
| | - Salvatore Mastrangelo
- Dipartimento di Scienze Agrarie e Forestali, Viale delle Scienze, Università Palermo, Palermo, Italy
| | - Fabio Pilla
- Dipartimento Agricoltura, Ambiente e Alimenti, Università degli Studi del Molise, Campobasso, Italy
| | - Elena Ciani
- Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, Università degli Studi di Bari “Aldo Moro ”, Bari, Italy
| | - Olivier Hanotte
- School of Life Sciences, University of Nottingham, Nottingham, United Kingdom
- LiveGene, International Livestock Research Institute, Addis Ababa, Ethiopia
| | - Joram M. Mwacharo
- Small Ruminant Genomics, International Center for Agricultural Research in the Dry Areas (ICARDA), Addis Ababa, Ethiopia
| |
Collapse
|
10
|
Avila F, Mickelson JR, Schaefer RJ, McCue ME. Genome-Wide Signatures of Selection Reveal Genes Associated With Performance in American Quarter Horse Subpopulations. Front Genet 2018; 9:249. [PMID: 30105047 PMCID: PMC6060370 DOI: 10.3389/fgene.2018.00249] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2018] [Accepted: 06/22/2018] [Indexed: 11/13/2022] Open
Abstract
Selective breeding for athletic performance in various disciplines has resulted in population stratification within the American Quarter Horse (QH) breed. The goals of this study were to utilize high density genotype data to: (1) identify genomic regions undergoing positive selection within and among QH subpopulations; (2) investigate haplotype structure within each QH subpopulation; and (3) identify candidate genes within genomic regions of interest (ROI), as well as biological pathways, predicted to play a role in elite performance in each group. For that, 65K SNP genotyping data on 143 elite individuals from 6 QH subpopulations (cutting, halter, racing, reining, western pleasure, and working cow) were imputed to 2M SNPs. Signatures of selection were identified using FST-based (di ) and haplotype-based (hapFLK) analyses, accompanied by identification of local haplotype structure and sharing within subpopulations (hapQTL). Regions undergoing positive selection were identified on all 31 autosomes, and ROI on 2 chromosomes were identified by all 3 methods combined. Genes within each ROI were retrieved and used to identify pathways and genes that might contribute to performance in each subpopulation. These included, among others, candidate genes associated with skeletal muscle development, metabolism, and central nervous system development. This work improves our understanding of equine breed development, and provides breeders with a better understanding of how selective breeding impacts the performance of QH populations.
Collapse
Affiliation(s)
- Felipe Avila
- Department of Veterinary Population Medicine, College of Veterinary Medicine, University of Minnesota, St. Paul, MN, United States
| | - James R Mickelson
- Department of Veterinary and Biomedical Sciences, College of Veterinary Medicine, University of Minnesota, St. Paul, MN, United States
| | - Robert J Schaefer
- Department of Veterinary Population Medicine, College of Veterinary Medicine, University of Minnesota, St. Paul, MN, United States
| | - Molly E McCue
- Department of Veterinary Population Medicine, College of Veterinary Medicine, University of Minnesota, St. Paul, MN, United States
| |
Collapse
|
11
|
Pendleton AL, Shen F, Taravella AM, Emery S, Veeramah KR, Boyko AR, Kidd JM. Comparison of village dog and wolf genomes highlights the role of the neural crest in dog domestication. BMC Biol 2018; 16:64. [PMID: 29950181 PMCID: PMC6022502 DOI: 10.1186/s12915-018-0535-2] [Citation(s) in RCA: 89] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2018] [Accepted: 05/23/2018] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND Domesticated from gray wolves between 10 and 40 kya in Eurasia, dogs display a vast array of phenotypes that differ from their ancestors, yet mirror other domesticated animal species, a phenomenon known as the domestication syndrome. Here, we use signatures persisting in dog genomes to identify genes and pathways possibly altered by the selective pressures of domestication. RESULTS Whole-genome SNP analyses of 43 globally distributed village dogs and 10 wolves differentiated signatures resulting from domestication rather than breed formation. We identified 246 candidate domestication regions containing 10.8 Mb of genome sequence and 429 genes. The regions share haplotypes with ancient dogs, suggesting that the detected signals are not the result of recent selection. Gene enrichments highlight numerous genes linked to neural crest and central nervous system development as well as neurological function. Read depth analysis suggests that copy number variation played a minor role in dog domestication. CONCLUSIONS Our results identify genes that act early in embryogenesis and can confer phenotypes distinguishing domesticated dogs from wolves, such as tameness, smaller jaws, floppy ears, and diminished craniofacial development as the targets of selection during domestication. These differences reflect the phenotypes of the domestication syndrome, which can be explained by alterations in the migration or activity of neural crest cells during development. We propose that initial selection during early dog domestication was for behavior, a trait influenced by genes which act in the neural crest, which secondarily gave rise to the phenotypes of modern dogs.
Collapse
Affiliation(s)
- Amanda L Pendleton
- Department of Human Genetics, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Feichen Shen
- Department of Human Genetics, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Angela M Taravella
- Department of Human Genetics, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Sarah Emery
- Department of Human Genetics, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Krishna R Veeramah
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY, 11794, USA
| | - Adam R Boyko
- Department of Biomedical Sciences, Cornell University, Ithaca, New York, 14853, USA
| | - Jeffrey M Kidd
- Department of Human Genetics, University of Michigan, Ann Arbor, MI, 48109, USA.
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, MI, 48109, USA.
| |
Collapse
|
12
|
Gutiérrez-Gil B, Esteban-Blanco C, Wiener P, Chitneedi PK, Suarez-Vega A, Arranz JJ. High-resolution analysis of selection sweeps identified between fine-wool Merino and coarse-wool Churra sheep breeds. Genet Sel Evol 2017; 49:81. [PMID: 29115919 PMCID: PMC5674817 DOI: 10.1186/s12711-017-0354-x] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Accepted: 10/19/2017] [Indexed: 12/20/2022] Open
Abstract
BACKGROUND With the aim of identifying selection signals in three Merino sheep lines that are highly specialized for fine wool production (Australian Industry Merino, Australian Merino and Australian Poll Merino) and considering that these lines have been subjected to selection not only for wool traits but also for growth and carcass traits and parasite resistance, we contrasted the OvineSNP50 BeadChip (50 K-chip) pooled genotypes of these Merino lines with the genotypes of a coarse-wool breed, phylogenetically related breed, Spanish Churra dairy sheep. Genome re-sequencing datasets of the two breeds were analyzed to further explore the genetic variation of the regions initially identified as putative selection signals. RESULTS Based on the 50 K-chip genotypes, we used the overlapping selection signals (SS) identified by four selection sweep mapping analyses (that detect genetic differentiation, reduced heterozygosity and patterns of haplotype diversity) to define 18 convergence candidate regions (CCR), five associated with positive selection in Australian Merino and the remainder indicating positive selection in Churra. Subsequent analysis of whole-genome sequences from 15 Churra and 13 Merino samples identified 142,400 genetic variants (139,745 bi-allelic SNPs and 2655 indels) within the 18 defined CCR. Annotation of 1291 variants that were significantly associated with breed identity between Churra and Merino samples identified 257 intragenic variants that caused 296 functional annotation variants, 275 of which were located across 31 coding genes. Among these, four synonymous and four missense variants (NPR2_His847Arg, NCAPG_Ser585Phe, LCORL_Asp1214Glu and LCORL_Ile1441Leu) were included. CONCLUSIONS Here, we report the mapping and genetic variation of 18 selection signatures that were identified between Australian Merino and Spanish Churra sheep breeds, which were validated by an additional contrast between Spanish Merino and Churra genotypes. Analysis of whole-genome sequencing datasets allowed us to identify divergent variants that may be viewed as candidates involved in the phenotypic differences for wool, growth and meat production/quality traits between the breeds analyzed. The four missense variants located in the NPR2, NCAPG and LCORL genes may be related to selection sweep regions previously identified and various QTL reported in sheep in relation to growth traits and carcass composition.
Collapse
Affiliation(s)
- Beatriz Gutiérrez-Gil
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| | - Cristina Esteban-Blanco
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
- Fundación Centro Supercomputación de Castilla y León, Campus de Vegazana, León, 24071 Spain
| | - Pamela Wiener
- Roslin Institute and R(D)SVS, University of Edinburgh, Easter Bush Campus, Midlothian, EH25 9RG UK
| | - Praveen Krishna Chitneedi
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| | - Aroa Suarez-Vega
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| | - Juan-Jose Arranz
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| |
Collapse
|
13
|
Brito LF, Kijas JW, Ventura RV, Sargolzaei M, Porto-Neto LR, Cánovas A, Feng Z, Jafarikia M, Schenkel FS. Genetic diversity and signatures of selection in various goat breeds revealed by genome-wide SNP markers. BMC Genomics 2017; 18:229. [PMID: 28288562 PMCID: PMC5348779 DOI: 10.1186/s12864-017-3610-0] [Citation(s) in RCA: 89] [Impact Index Per Article: 12.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2016] [Accepted: 03/07/2017] [Indexed: 01/08/2023] Open
Abstract
Background The detection of signatures of selection has the potential to elucidate the identities of genes and mutations associated with phenotypic traits important for livestock species. It is also very relevant to investigate the levels of genetic diversity of a population, as genetic diversity represents the raw material essential for breeding and has practical implications for implementation of genomic selection. A total of 1151 animals from nine goat populations selected for different breeding goals and genotyped with the Illumina Goat 50K single nucleotide polymorphisms (SNP) Beadchip were included in this investigation. Results The proportion of polymorphic SNPs ranged from 0.902 (Nubian) to 0.995 (Rangeland). The overall mean HO and HE was 0.374 ± 0.021 and 0.369 ± 0.023, respectively. The average pairwise genetic distance (D) ranged from 0.263 (Toggenburg) to 0.323 (Rangeland). The overall average for the inbreeding measures FEH, FVR, FLEUT, FROH and FPED was 0.129, −0.012, −0.010, 0.038 and 0.030, respectively. Several regions located on 19 chromosomes were potentially under selection in at least one of the goat breeds. The genomic population tree constructed using all SNPs differentiated breeds based on selection purpose, while genomic population tree built using only SNPs in the most significant region showed a great differentiation between LaMancha and the other breeds. We hypothesized that this region is related to ear morphogenesis. Furthermore, we identified genes potentially related to reproduction traits, adult body mass, efficiency of food conversion, abdominal fat deposition, conformation traits, liver fat metabolism, milk fatty acids, somatic cells score, milk protein, thermo-tolerance and ear morphogenesis. Conclusions In general, moderate to high levels of genetic variability were observed for all the breeds and a characterization of runs of homozygosity gave insights into the breeds’ development history. The information reported here will be useful for the implementation of genomic selection and other genomic studies in goats. We also identified various genome regions under positive selection using smoothed FST and hapFLK statistics and suggested genes, which are potentially under selection. These results can now provide a foundation to formulate biological hypotheses related to selection processes in goats. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3610-0) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Luiz F Brito
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.
| | - James W Kijas
- CSIRO Agriculture & Food, Brisbane, Queensland, Australia
| | - Ricardo V Ventura
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.,Beef Improvement Opportunities, Guelph, Ontario, Canada
| | - Mehdi Sargolzaei
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.,The Semex Alliance, Guelph, Ontario, Canada
| | | | - Angela Cánovas
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada
| | - Zeny Feng
- Department of Mathematics and Statistics, University of Guelph, Guelph, Ontario, Canada
| | - Mohsen Jafarikia
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada.,Canadian Centre for Swine Improvement Inc., Ottawa, Ontario, Canada
| | - Flávio S Schenkel
- Centre for Genetic Improvement of Livestock, University of Guelph, Guelph, Ontario, Canada
| |
Collapse
|
14
|
HAN YJ, CHEN Y, LIU Y, LIU XL. Sequence variants of the LCORL gene and its association with growth and carcass traits in Qinchuan cattle in China. J Genet 2017; 96:9-17. [DOI: 10.1007/s12041-016-0732-0] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
15
|
Population structure of eleven Spanish ovine breeds and detection of selective sweeps with BayeScan and hapFLK. Sci Rep 2016; 6:27296. [PMID: 27272025 PMCID: PMC4895181 DOI: 10.1038/srep27296] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2015] [Accepted: 04/26/2016] [Indexed: 11/09/2022] Open
Abstract
The goals of the current work were to analyse the population structure of 11 Spanish ovine breeds and to detect genomic regions that may have been targeted by selection. A total of 141 individuals were genotyped with the Infinium 50 K Ovine SNP BeadChip (Illumina). We combined this dataset with Spanish ovine data previously reported by the International Sheep Genomics Consortium (N = 229). Multidimensional scaling and Admixture analyses revealed that Canaria de Pelo and, to a lesser extent, Roja Mallorquina, Latxa and Churra are clearly differentiated populations, while the remaining seven breeds (Ojalada, Castellana, Gallega, Xisqueta, Ripollesa, Rasa Aragonesa and Segureña) share a similar genetic background. Performance of a genome scan with BayeScan and hapFLK allowed us identifying three genomic regions that are consistently detected with both methods i.e. Oar3 (150–154 Mb), Oar6 (4–49 Mb) and Oar13 (68–74 Mb). Neighbor-joining trees based on polymorphisms mapping to these three selective sweeps did not show a clustering of breeds according to their predominant productive specialization (except the local tree based on Oar13 SNPs). Such cryptic signatures of selection have been also found in the bovine genome, posing a considerable challenge to understand the biological consequences of artificial selection.
Collapse
|
16
|
Uncovering Adaptation from Sequence Data: Lessons from Genome Resequencing of Four Cattle Breeds. Genetics 2016; 203:433-50. [PMID: 27017625 DOI: 10.1534/genetics.115.181594] [Citation(s) in RCA: 60] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2015] [Accepted: 03/03/2016] [Indexed: 01/25/2023] Open
Abstract
Detecting the molecular basis of adaptation is one of the major questions in population genetics. With the advance in sequencing technologies, nearly complete interrogation of genome-wide polymorphisms in multiple populations is becoming feasible in some species, with the expectation that it will extend quickly to new ones. Here, we investigate the advantages of sequencing for the detection of adaptive loci in multiple populations, exploiting a recently published data set in cattle (Bos taurus). We used two different approaches to detect statistically significant signals of positive selection: a within-population approach aimed at identifying hard selective sweeps and a population-differentiation approach that can capture other selection events such as soft or incomplete sweeps. We show that the two methods are complementary in that they indeed capture different kinds of selection signatures. Our study confirmed some of the well-known adaptive loci in cattle (e.g., MC1R, KIT, GHR, PLAG1, NCAPG/LCORL) and detected some new ones (e.g., ARL15, PRLR, CYP19A1, PPM1L). Compared to genome scans based on medium- or high-density SNP data, we found that sequencing offered an increased detection power and a higher resolution in the localization of selection signatures. In several cases, we could even pinpoint the underlying causal adaptive mutation or at least a very small number of possible candidates (e.g., MC1R, PLAG1). Our results on these candidates suggest that a vast majority of adaptive mutations are likely to be regulatory rather than protein-coding variants.
Collapse
|
17
|
Detection of Selection Signatures on the X Chromosome in Three Sheep Breeds. Int J Mol Sci 2015; 16:20360-74. [PMID: 26343642 PMCID: PMC4613208 DOI: 10.3390/ijms160920360] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2015] [Revised: 08/18/2015] [Accepted: 08/19/2015] [Indexed: 12/24/2022] Open
Abstract
Artificial selection has played a critical role in animal breeding. Detection of artificial selection footprints in genomic regions can provide insights for understanding the function of specific phenotypic traits and better guide animal breeding. To more fully understand the relationship between genomic composition and phenotypic diversity arising from breed development, a genome-wide scan was conducted using an OvineSNP50 BeadChip and integrated haplotype score and fixation index analyses to detect selection signatures on the X chromosome in three sheep breeds. We identified 49, 34, and 55 candidate selection regions with lengths of 27.49, 16.47, and 25.42 Mb in German Mutton, Dorper, and Sunit sheep, respectively. Bioinformatics analysis showed that some of the genes in these regions with selection signatures, such as BMP15, were relevant to reproduction. We also identified some selection regions harboring genes that had human orthologs, including BKT, CENPI, GUCY2F, MSN, PCDH11X, PLP1, VSIG4, PAK3, WAS, PCDH19, PDHA1, and SRPX2. The VSIG4 and PCDH11X genes are associated with the immune system and disease, PDHA1 is associated with biosynthetic related pathways, and PCDH19 is expressed in the nervous system and skin. These genes may be useful as candidate genes for molecular breeding.
Collapse
|
18
|
Wang H, Zhang L, Cao J, Wu M, Ma X, Liu Z, Liu R, Zhao F, Wei C, Du L. Genome-Wide Specific Selection in Three Domestic Sheep Breeds. PLoS One 2015; 10:e0128688. [PMID: 26083354 PMCID: PMC4471085 DOI: 10.1371/journal.pone.0128688] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2015] [Accepted: 04/29/2015] [Indexed: 01/01/2023] Open
Abstract
BACKGROUND Commercial sheep raised for mutton grow faster than traditional Chinese sheep breeds. Here, we aimed to evaluate genetic selection among three different types of sheep breed: two well-known commercial mutton breeds and one indigenous Chinese breed. RESULTS We first combined locus-specific branch lengths and di statistical methods to detect candidate regions targeted by selection in the three different populations. The results showed that the genetic distances reached at least medium divergence for each pairwise combination. We found these two methods were highly correlated, and identified many growth-related candidate genes undergoing artificial selection. For production traits, APOBR and FTO are associated with body mass index. For meat traits, ALDOA, STK32B and FAM190A are related to marbling. For reproduction traits, CCNB2 and SLC8A3 affect oocyte development. We also found two well-known genes, GHR (which affects meat production and quality) and EDAR (associated with hair thickness) were associated with German mutton merino sheep. Furthermore, four genes (POL, RPL7, MSL1 and SHISA9) were associated with pre-weaning gain in our previous genome-wide association study. CONCLUSIONS Our results indicated that combine locus-specific branch lengths and di statistical approaches can reduce the searching ranges for specific selection. And we got many credible candidate genes which not only confirm the results of previous reports, but also provide a suite of novel candidate genes in defined breeds to guide hybridization breeding.
Collapse
Affiliation(s)
- Huihua Wang
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Li Zhang
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Jiaxve Cao
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Mingming Wu
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Xiaomeng Ma
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Zhen Liu
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Ruizao Liu
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Fuping Zhao
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Caihong Wei
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| | - Lixin Du
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, People’s Republic of China, Beijing, China
| |
Collapse
|