1
|
Zhao L, Yuan L, Li F, Zhang X, Tian H, Ma Z, Zhang D, Zhang Y, Zhao Y, Huang K, Li X, Cheng J, Xu D, Yang X, Han K, Weng X, Wang W. Whole-genome resequencing of Hu sheep identifies candidate genes associated with agronomic traits. J Genet Genomics 2024:S1673-8527(24)00068-7. [PMID: 38582298 DOI: 10.1016/j.jgg.2024.03.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2023] [Revised: 03/30/2024] [Accepted: 03/30/2024] [Indexed: 04/08/2024]
Abstract
The phenotypic diversity resulting from artificial or natural selection of sheep has made a significant contribution to human civilization. Hu sheep are a local sheep breed unique to China with high reproductive rates and rapid growth. Genomic selection signatures have been widely used to investigate the genetic mechanisms underlying phenotypic variation in livestock. Here, we conduct whole-genome sequencing of 207 Hu sheep and compare them with the wild ancestors of domestic sheep (Asiatic mouflon) to investigate the genetic characteristics and selection signatures of Hu sheep. Based on six signatures of selection approaches, we detect genomic regions containing genes related to reproduction (BMPR1B, BMP2, PGFS, CYP19, CAMK4, GGT5, and GNAQ), vision (ALDH1A2, SAG, and PDE6B), nervous system (NAV1), and immune response (GPR35, SH2B2, PIK3R3, and HRAS). Association analysis with a population of 1299 Hu sheep reveal those missense mutations in the GPR35 (GPR35 g.952651 A>G; GPR35 g.952496 C>T) and NAV1 (NAV1 g.84216190 C>T; NAV1 g.84227412 G>A) genes are significantly associated (P < 0.05) with immune and growth traits in Hu sheep, respectively. This research offers unique insights into the selection characteristics of Hu sheep and facilitates further genetic improvement and molecular investigations.
Collapse
Affiliation(s)
- Liming Zhao
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Lvfeng Yuan
- Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences, Lanzhou, Gansu 730046, China
| | - Fadi Li
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Xiaoxue Zhang
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, Gansu 730070, China
| | - Huibin Tian
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Zongwu Ma
- College of Animal Science and Technology, Gansu Agricultural University, Lanzhou, Gansu 730070, China
| | - Deyin Zhang
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Yukun Zhang
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Yuan Zhao
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Kai Huang
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Xiaolong Li
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Jiangbo Cheng
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Dan Xu
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Xiaobin Yang
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Kunchao Han
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Xiuxiu Weng
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China
| | - Weimin Wang
- State Key Laboratory of Herbage Improvement and Grassland Agro-Ecosystems, Key Laboratory of Grassland Livestock Industry Innovation, Ministry of Agriculture and Rural Affairs, Engineering Research Center of Grassland Industry, Ministry of Education, College of Pastoral Agriculture Science and Technology, Lanzhou University, Lanzhou, Gansu 730020, China.
| |
Collapse
|
2
|
Nayak SS, Panigrahi M, Rajawat D, Ghildiyal K, Sharma A, Jain K, Bhushan B, Dutt T. Deciphering climate resilience in Indian cattle breeds by selection signature analyses. Trop Anim Health Prod 2024; 56:46. [PMID: 38233536 DOI: 10.1007/s11250-023-03879-8] [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: 08/11/2023] [Accepted: 12/07/2023] [Indexed: 01/19/2024]
Abstract
The signature of selection is a crucial concept in evolutionary biology that refers to the pattern of genetic variation which arises in a population due to natural selection. In the context of climate adaptation, the signature of selection can reveal the genetic basis of adaptive traits that enable organisms to survive and thrive in changing environmental conditions. Breeds living in diverse agroecological zones exhibit genetic "footprints" within their genomes that mirror the influence of climate-induced selective pressures, subsequently impacting phenotypic variance. It is assumed that the genomes of animals residing in these regions have been altered through selection for various climatic adaptations. These regions are known as signatures of selection and can be identified using various summary statistics. We examined genotypic data from eight different cattle breeds (Gir, Hariana, Kankrej, Nelore, Ongole, Red Sindhi, Sahiwal, and Tharparkar) that are adapted to diverse regional climates. To identify selection signature regions in this investigation, we used four intra-population statistics: Tajima's D, CLR, iHS, and ROH. In this study, we utilized Bovine 50 K chip data and four genome scan techniques to assess the genetic regions of positive selection for high-temperature adaptation. We have also performed a genome-wide investigation of genetic diversity, inbreeding, and effective population size in our target dataset. We identified potential regions for selection that are likely to be caused by adverse climatic conditions. We observed many adaptation genes in several potential selection signature areas. These include genes like HSPB2, HSPB3, HSP20, HSP90AB1, HSF4, HSPA1B, CLPB, GAP43, MITF, and MCHR1 which have been reported in the cattle populations that live in varied climatic regions. The findings demonstrated that genes involved in disease resistance and thermotolerance were subjected to intense selection. The findings have implications for marker-assisted breeding, understanding the genetic landscape of climate-induced adaptation, putting breeding and conservation programs into action.
Collapse
Affiliation(s)
- Sonali Sonejita Nayak
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India.
| | - Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India
| | - Kanika Ghildiyal
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India
| | - Anurodh Sharma
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India
| | - Karan Jain
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India
| | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India
| | - Triveni Dutt
- Livestock Production and Management Section, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, 243122, UP, India
| |
Collapse
|
3
|
Persichilli C, Senczuk G, Mastrangelo S, Marusi M, van Kaam JT, Finocchiaro R, Di Civita M, Cassandro M, Pilla F. Exploring genome-wide differentiation and signatures of selection in Italian and North American Holstein populations. J Dairy Sci 2023; 106:5537-5553. [PMID: 37291034 DOI: 10.3168/jds.2022-22159] [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: 04/05/2022] [Accepted: 02/07/2023] [Indexed: 06/10/2023]
Abstract
Among Italian dairy cattle, the Holstein is the most reared breed for the production of Parmigiano Reggiano protected designation of origin cheese, which represents one of the most renowned products in the entire Italian dairy industry. In this work, we used a medium-density genome-wide data set consisting of 79,464 imputed SNPs to study the genetic structure of Italian Holstein breed, including the population reared in the area of Parmigiano Reggiano cheese production, and assessing its distinctiveness from the North American population. Multidimensional scaling and ADMIXTURE approaches were used to explore the genetic structure among populations. We also investigated putative genomic regions under selection among these 3 populations by combining 4 different statistical methods based either on allele frequencies (single marker and window-based) or extended haplotype homozygosity (EHH; standardized log-ratio of integrated EHH and cross-population EHH). The genetic structure results allowed us to clearly distinguish the 3 Holstein populations; however, the most remarkable difference was observed between Italian and North American stock. Selection signature analyses identified several significant SNPs falling within or closer to genes with known roles in several traits such as milk quality, resistance to disease, and fertility. In particular, a total of 22 genes related to milk production have been identified using the 2 allele frequency approaches. Among these, a convergent signal has been found in the VPS8 gene which resulted to be involved in milk traits, whereas other genes (CYP7B1, KSR2, C4A, LIPE, DCDC1, GPR20, and ST3GAL1) resulted to be associated with quantitative trait loci related to milk yield and composition in terms of fat and protein percentage. In contrast, a total of 7 genomic regions were identified combining the results of standardized log-ratio of integrated EHH and cross-population EHH. In these regions candidate genes for milk traits were also identified. Moreover, this was also confirmed by the enrichment analyses in which we found that the majority of the significantly enriched quantitative trait loci were linked to milk traits, whereas the gene ontology and pathway enrichment analysis pointed to molecular functions and biological processes involved in AA transmembrane transport and methane metabolism pathway. This study provides information on the genetic structure of the examined populations, showing that they are distinguishable from each other. Furthermore, the selection signature analyses can be considered as a starting point for future studies in the identification of causal mutations and consequent implementation of more practical application.
Collapse
Affiliation(s)
- Christian Persichilli
- Department of Agricultural, Environmental and Food Sciences, University of Molise, Via A. De sanctis, 86100 Campobasso (CB), Italy
| | - Gabriele Senczuk
- Department of Agricultural, Environmental and Food Sciences, University of Molise, Via A. De sanctis, 86100 Campobasso (CB), Italy.
| | - Salvatore Mastrangelo
- Department of Agricultural, Food and Forest Sciences, University of Palermo, Viale delle Scienze, 90128 Palermo (PA), Italy
| | - Maurizio Marusi
- National Association of Italian Holstein, Brown and Jersey Breeders, Via Bergamo, 292, 26100 Cremona (CR), Italy
| | - Jan-Thijs van Kaam
- National Association of Italian Holstein, Brown and Jersey Breeders, Via Bergamo, 292, 26100 Cremona (CR), Italy
| | - Raffaella Finocchiaro
- National Association of Italian Holstein, Brown and Jersey Breeders, Via Bergamo, 292, 26100 Cremona (CR), Italy
| | - Marika Di Civita
- Department of Agricultural, Environmental and Food Sciences, University of Molise, Via A. De sanctis, 86100 Campobasso (CB), Italy
| | - Martino Cassandro
- National Association of Italian Holstein, Brown and Jersey Breeders, Via Bergamo, 292, 26100 Cremona (CR), Italy; Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale dell'Università 16, 35020 Legnaro (PD), Italy
| | - Fabio Pilla
- Department of Agricultural, Environmental and Food Sciences, University of Molise, Via A. De sanctis, 86100 Campobasso (CB), Italy
| |
Collapse
|
4
|
Revealing Genetic Diversity and Population Structure of Endangered Altay White-Headed Cattle Population Using 100 k SNP Markers. Animals (Basel) 2022; 12:ani12223214. [PMID: 36428441 PMCID: PMC9686749 DOI: 10.3390/ani12223214] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 10/17/2022] [Accepted: 11/18/2022] [Indexed: 11/22/2022] Open
Abstract
Understanding the genetic basis of native cattle populations that have adapted to the local environment is of great significance for formulating appropriate strategies and programs for genetic improvement and protection. Therefore, it is necessary to understand the genetic diversity and population structure of Altay white-headed cattle so as to meet the current production needs under various environments, carry out continuous genetic improvement, and promote rapid adaptation to changing environments and breeding objectives. A total of 46 individual samples of endangered Xinjiang Altay white-headed cattle were collected in this study, including nine bulls and 37 cows. To collect genotype data, 100 k SNP markers were used, and then studies of genetic diversity, genetic structure, inbreeding degree, and family analysis were carried out. A total of 101,220 SNP loci were detected, and the genotype detection rate for individuals was ≥90%. There were 85,993 SNP loci that passed quality control, of which 93.5% were polymorphic. The average effective allele number was 0.036, the Polymorphism Information Content was 0.304 and the minimum allele frequency was 0.309, the average observed heterozygosity was 0.413, and the average expected heterozygosity was 0.403. The average genetic distance of Idengtical By State (IBS) was 0.3090, there were 461 ROH (genome-length homozygous fragments), 76.1% of which were between 1 and 5 MB in length, and the average inbreeding coefficient was 0.016. The 46 Altay white-headed cattle were divided into their families, and the individual numbers of each family were obviously different. To sum up, the Altay white-headed cattle conservation population had low heterozygosity, a high inbreeding degree, few families, and large differences in the number of individuals in each family, which can easily cause a loss of genetic diversity. In the follow-up seed conservation process, seed selection and matching should be carried out according to the divided families to ensure the long-term protection of Altay white-headed cattle genetic resources.
Collapse
|
5
|
Duarte INH, Bessa AFDO, Rola LD, Genuíno MVH, Rocha IM, Marcondes CR, Regitano LCDA, Munari DP, Berry DP, Buzanskas ME. Cross-population selection signatures in Canchim composite beef cattle. PLoS One 2022; 17:e0264279. [PMID: 35363779 PMCID: PMC8975110 DOI: 10.1371/journal.pone.0264279] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2021] [Accepted: 02/07/2022] [Indexed: 12/15/2022] Open
Abstract
Analyses of livestock genomes have been used to detect selection signatures, which are genomic regions associated with traits under selection leading to a change in allele frequency. The objective of the present study was to characterize selection signatures in Canchim composite beef cattle using cross-population analyses with the founder Nelore and Charolais breeds. High-density single nucleotide polymorphism genotypes were available on 395 Canchim representing the target population, along with genotypes from 809 Nelore and 897 Charolais animals representing the reference populations. Most of the selection signatures were co-located with genes whose functions agree with the expectations of the breeding programs; these genes have previously been reported to associate with meat quality, as well as reproductive traits. Identified genes were related to immunity, adaptation, morphology, as well as behavior, could give new perspectives for understanding the genetic architecture of Canchim. Some selection signatures identified genes that were recently introduced in Canchim, such as the loci related to the polled trait.
Collapse
Affiliation(s)
| | | | - Luciana Diniz Rola
- Departamento de Zootecnia, Universidade Federal da Paraíba, Areia, Paraíba, Brazil
| | | | - Iasmin Marques Rocha
- Departamento de Zootecnia, Universidade Federal da Paraíba, Areia, Paraíba, Brazil
| | | | | | - Danísio Prado Munari
- Departamento de Engenharia e Ciências Exatas, Universidade Estadual Paulista, Jaboticabal, São Paulo, Brazil
| | - Donagh Pearse Berry
- Teagasc, Animal & Grassland Research and Innovation Centre, Moorepark, Fermoy Co. Cork., Ireland
| | - Marcos Eli Buzanskas
- Departamento de Zootecnia, Universidade Federal da Paraíba, Areia, Paraíba, Brazil
- * E-mail:
| |
Collapse
|
6
|
Rekik E, Ahbara AM, Abate Z, Goshme S, Getachew T, Haile A, Rischkowsky B, Mwacharo JM. Genomic analysis of 10 years of artificial selection in community‐based breeding programs in two Ethiopian indigenous sheep breeds. Anim Genet 2022; 53:447-451. [PMID: 35428998 PMCID: PMC10138745 DOI: 10.1111/age.13190] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Revised: 02/25/2022] [Accepted: 03/01/2022] [Indexed: 11/30/2022]
Abstract
In recent times, community-based breeding programs (CBBPs) have been advocated as the best strategy for genetic improvement of local breeds in smallholder farms in developing countries. Since 2009, CBBPs have been implemented for Ethiopian Bonga and Menz sheep to improve growth rates resulting in significant genetic gains in 6-month weights. With the hypothesis that selection could be impacting their genomes, we systematically screened for possible genome changes in the two breeds by analyzing 600K BeadChip genotype data of 151 individuals (with the highest breeding values for 6-month weights) from CBBP flocks against 98 individuals from non-CBBP flocks. We observed no differences in genetic diversity and demographic dynamics between CBBP and non-CBBP flocks. Selection signature analysis employing ROH, logistic regression genome-wide association study , FST , XP-EHH and iHS revealed 5 (Bonga) and 11 (Menz) overlapping regions under selection, that co-localized with QTLs for production (body size/weight, growth, milk yield), meat/milk quality, and health/parasite resistance, suggesting that the decade-long selection has likely started to impact their genomes. However, genome-wide genetic differentiation between the CBBP and non-CBBP flocks is not yet clearly evident.
Collapse
Affiliation(s)
- Emna Rekik
- Small Ruminant Genomics International Centre for Agricultural Research in the Dry areas (ICARDA) Addis Ababa Ethiopia
| | - Abulgasim M. Ahbara
- Small Ruminant Genomics International Centre for Agricultural Research in the Dry areas (ICARDA) Addis Ababa Ethiopia
- Department of Zoology Faculty of Sciences Misurata University Misurata Libya
| | - Zelalem Abate
- Animal Sciences Case Team Bonga Agricultural Research Center Bonga Ethiopia
| | - Shenkute Goshme
- Debre‐Birhan Agricultural Research Center Debre‐Birhan Ethiopia
| | - Tesfaye Getachew
- Small Ruminant Genomics International Centre for Agricultural Research in the Dry areas (ICARDA) Addis Ababa Ethiopia
| | - Aynalem Haile
- Small Ruminant Genomics International Centre for Agricultural Research in the Dry areas (ICARDA) Addis Ababa Ethiopia
| | - Barbara Rischkowsky
- Small Ruminant Genomics International Centre for Agricultural Research in the Dry areas (ICARDA) Addis Ababa Ethiopia
| | - Joram M. Mwacharo
- Small Ruminant Genomics International Centre for Agricultural Research in the Dry areas (ICARDA) Addis Ababa Ethiopia
- Animal and Veterinary Sciences Scotland Rural College and Centre for Tropical Livestock Genetics and Health (CTLGH) The Roslin Institute Building Easter Bush Midlothian UK
| |
Collapse
|
7
|
Hartfield M, Poulsen NA, Guldbrandtsen B, Bataillon T. Using singleton densities to detect recent selection in Bos taurus. Evol Lett 2021; 5:595-606. [PMID: 34917399 PMCID: PMC8645200 DOI: 10.1002/evl3.263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2020] [Revised: 10/07/2021] [Accepted: 10/08/2021] [Indexed: 11/05/2022] Open
Abstract
Many quantitative traits are subject to polygenic selection, where several genomic regions undergo small, simultaneous changes in allele frequency that collectively alter a phenotype. The widespread availability of genome data, along with novel statistical techniques, has made it easier to detect these changes. We apply one such method, the "Singleton Density Score" (SDS), to the Holstein breed of Bos taurus to detect recent selection (arising up to around 740 years ago). We identify several genes as candidates for targets of recent selection, including some relating to cell regulation, catabolic processes, neural-cell adhesion and immunity. We do not find strong evidence that three traits that are important to humans-milk protein content, milk fat content, and stature-have been subject to directional selection. Simulations demonstrate that because B. taurus recently experienced a population bottleneck, singletons are depleted so the power of SDS methods is reduced. These results inform on which genes underlie recent genetic change in B. taurus, while providing information on how polygenic selection can be best investigated in future studies.
Collapse
Affiliation(s)
- Matthew Hartfield
- Bioinformatics Research CentreAarhus UniversityAarhusDK‐8000Denmark
- Institute of Evolutionary BiologyUniversity of EdinburghEdinburghEH9 3FLUnited Kingdom
| | | | - Bernt Guldbrandtsen
- Center for Quantitative Genetics and Genomics, Department of Molecular Biology and GeneticsAarhus UniversityTjeleDK‐8830Denmark
- Rheinische Friedrich‐Wilhelms‐Universität BonnInstitut für TierwissenschaftenBonnDE‐53115Germany
- Department of Veterinary SciencesCopenhagen UniversityFrederiksberg CDK‐1870Denmark
| | - Thomas Bataillon
- Bioinformatics Research CentreAarhus UniversityAarhusDK‐8000Denmark
| |
Collapse
|
8
|
Makanjuola BO, Maltecca C, Miglior F, Marras G, Abdalla EA, Schenkel FS, Baes CF. Identification of unique ROH regions with unfavorable effects on production and fertility traits in Canadian Holsteins. Genet Sel Evol 2021; 53:68. [PMID: 34461820 PMCID: PMC8406729 DOI: 10.1186/s12711-021-00660-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Accepted: 08/19/2021] [Indexed: 12/02/2022] Open
Abstract
Background The advent of genomic information and the reduction in the cost of genotyping have led to the use of genomic information to estimate genomic inbreeding as an alternative to pedigree inbreeding. Using genomic measures, effects of genomic inbreeding on production and fertility traits have been observed. However, there have been limited studies on the specific genomic regions causing the observed negative association with the trait of interest. Our aim was to identify unique run of homozygosity (ROH) genotypes present within a given genomic window that display negative associations with production and fertility traits and to quantify the effects of these identified ROH genotypes. Methods In total, 50,575 genotypes based on a 50K single nucleotide polymorphism (SNP) array and 259,871 pedigree records were available. Of these 50,575 genotypes, 46,430 cows with phenotypic records for production and fertility traits and having a first calving date between 2008 and 2018 were available. Unique ROH genotypes identified using a sliding-window approach were fitted into an animal mixed model as fixed effects to determine their effect on production and fertility traits. Results In total, 133 and 34 unique ROH genotypes with unfavorable effects were identified for production and fertility traits, respectively, at a 1% genome-wise false discovery rate. Most of these ROH regions were located on bovine chromosomes 8, 13, 14 and 19 for both production and fertility traits. For production traits, the average of all the unfavorably identified unique ROH genotypes effects were estimated to decrease milk yield by 247.30 kg, fat yield by 11.46 kg and protein yield by 8.11 kg. Similarly, for fertility traits, an average 4.81-day extension in first service to conception, a 0.16 increase in number of services, and a − 0.07 incidence in 56-day non-return rate were observed. Furthermore, a ROH region located on bovine chromosome 19 was identified that, when homozygous, had a negative effect on production traits. Signatures of selection proximate to this region have implicated GH1 as a potential candidate gene, which encodes the growth hormone that binds the growth hormone receptor. This observed negative effect could be a consequence of unfavorable alleles in linkage disequilibrium with favorable alleles. Conclusions ROH genotypes with unfavorable effects on production and fertility traits were identified within and across multiple traits on most chromosomes. These identified ROH genotypes could be included in mate selection programs to minimize their frequency in future generations. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00660-z.
Collapse
Affiliation(s)
- Bayode O Makanjuola
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| | - Christian Maltecca
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.,Department of Animal Science and Genetics Program, North Carolina State University, Raleigh, NC, 27607, USA
| | - Filippo Miglior
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | | | - Emhimad A Abdalla
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Flavio S Schenkel
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Christine F Baes
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.,Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| |
Collapse
|
9
|
Peng Q, Ehlers CL. Long tracks of homozygosity predict the severity of alcohol use disorders in an American Indian population. Mol Psychiatry 2021; 26:2200-2211. [PMID: 33398086 PMCID: PMC8254832 DOI: 10.1038/s41380-020-00989-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/13/2020] [Revised: 11/30/2020] [Accepted: 12/07/2020] [Indexed: 11/20/2022]
Abstract
Runs of homozygosity (ROH) arise when an individual inherits two copies of the same haplotype segment. While ROH are ubiquitous across human populations, Native populations-with shared parental ancestry arising from isolation and endogamy-can carry a substantial enrichment for ROH. We have been investigating genetic and environmental risk factors for alcohol use disorders (AUD) in a group of American Indians (AI) who have higher rates of AUD than the general U. S. population. Here we explore whether ROH might be associated with incidence and severity of AUD in this admixed AI population (n = 742) that live on geographically contiguous reservations, using low-coverage whole genome sequences. We have found that the genomic regions in the ROH that were identified in this population had significantly elevated American Indian heritage compared with the rest of the genome. Increased ROH abundance and ROH burden are likely risk factors for AUD severity in this AI population, especially in those diagnosed with severe and moderate AUD. The association between ROH and AUD was mostly driven by ROH of moderate lengths between 1 and 2 Mb. An ROH island on chromosome 1p32.3 and a rare ROH pool on chromosome 3p12.3 were found to be significantly associated with AUD severity. They contain genes involved in lipid metabolism, oxidative stress and inflammatory responses; and OSBPL9 was found to reside on the consensus part of the ROH island. These data demonstrate that ROH are associated with risk for AUD severity in this AI population.
Collapse
Affiliation(s)
- Qian Peng
- Department of Neuroscience, The Scripps Research Institute, La Jolla, CA, 92037, USA.
| | - Cindy L Ehlers
- Department of Neuroscience, The Scripps Research Institute, La Jolla, CA, 92037, USA.
| |
Collapse
|
10
|
Geibel J, Reimer C, Weigend S, Weigend A, Pook T, Simianer H. How array design creates SNP ascertainment bias. PLoS One 2021; 16:e0245178. [PMID: 33784304 PMCID: PMC8009414 DOI: 10.1371/journal.pone.0245178] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2020] [Accepted: 12/22/2020] [Indexed: 12/30/2022] Open
Abstract
Single nucleotide polymorphisms (SNPs), genotyped with arrays, have become a widely used marker type in population genetic analyses over the last 10 years. However, compared to whole genome re-sequencing data, arrays are known to lack a substantial proportion of globally rare variants and tend to be biased towards variants present in populations involved in the development process of the respective array. This affects population genetic estimators and is known as SNP ascertainment bias. We investigated factors contributing to ascertainment bias in array development by redesigning the Axiom™ Genome-Wide Chicken Array in silico and evaluating changes in allele frequency spectra and heterozygosity estimates in a stepwise manner. A sequential reduction of rare alleles during the development process was shown. This was mainly caused by the identification of SNPs in a limited set of populations and a within-population selection of common SNPs when aiming for equidistant spacing. These effects were shown to be less severe with a larger discovery panel. Additionally, a generally massive overestimation of expected heterozygosity for the ascertained SNP sets was shown. This overestimation was 24% higher for populations involved in the discovery process than not involved populations in case of the original array. The same was observed after the SNP discovery step in the redesign. However, an unequal contribution of populations during the SNP selection can mask this effect but also adds uncertainty. Finally, we make suggestions for the design of specialized arrays for large scale projects where whole genome re-sequencing techniques are still too expensive.
Collapse
Affiliation(s)
- Johannes Geibel
- Department of Animal Sciences, Animal Breeding and Genetics Group, University of Goettingen, Göttingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Göttingen, Germany
- * E-mail:
| | - Christian Reimer
- Department of Animal Sciences, Animal Breeding and Genetics Group, University of Goettingen, Göttingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Göttingen, Germany
| | - Steffen Weigend
- Center for Integrated Breeding Research, University of Goettingen, Göttingen, Germany
- Institute of Farm Animal Genetics, Friedrich-Loeffler-Institut, Neustadt-Mariensee, Germany
| | - Annett Weigend
- Institute of Farm Animal Genetics, Friedrich-Loeffler-Institut, Neustadt-Mariensee, Germany
| | - Torsten Pook
- Department of Animal Sciences, Animal Breeding and Genetics Group, University of Goettingen, Göttingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Göttingen, Germany
| | - Henner Simianer
- Department of Animal Sciences, Animal Breeding and Genetics Group, University of Goettingen, Göttingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Göttingen, Germany
| |
Collapse
|
11
|
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
|
12
|
Guo L, Sun H, Zhao Q, Xu Z, Zhang Z, Liu D, Qadri QR, Ma P, Wang Q, Pan Y. Positive selection signatures in Anqing six-end-white pig population based on reduced-representation genome sequencing data. Anim Genet 2021; 52:143-154. [PMID: 33458851 DOI: 10.1111/age.13034] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/10/2020] [Indexed: 12/26/2022]
Abstract
Anqing six-end-white (AQ) pig performs well on resistance to coarse fodder and disease, reproduction and meat quality, offering high potential for exploitation. Environmental conditions and strict selections from local farmers have cultivated the AQ pig to be an outstanding and unique local pig breed. Thus we aim to detect genetic positive selection signatures within the AQ pig population to explore underlying genetic mechanisms. A relative extended haplotype homozygosity (REHH) test was performed in the population of 79 AQ pigs to seek evidence demonstrating that selective actions have left an imprint on the whole genome. In total, 430 500 REHH tests were performed on 53 067 core regions with average REHH tests of 8.11, average lengths of 11.50 kb and an overall length of 610.38 Mb which accounted for 26.94% of the whole genome. Finally, a total of 1819 core haplotypes (P < 0.01) and 586 candidate genes were obtained. These genes were mainly related to meat quality (MYOG, SNX19), resistance to disease (CRISPLD2, CD14) and reproduction traits (ERBB2, NRP2). A panel of genes within the 30 top significant REHH tests was mainly categorized to traits of meat quality and disease resistance. Among 13 KEGG pathways, MAPK, GnRH and Oxytocin signaling pathways, associated with the biological processes of crucial economic traits, were noteworthy. The excellent characteristics of the AQ pig benefited from the combination of natural and human factors. We provide a sketch map that shows the distribution of selection footprints on the whole genome of AQ pig and found potential genes for future studies.
Collapse
Affiliation(s)
- L Guo
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - H Sun
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - Q Zhao
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - Z Xu
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - Z Zhang
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - D Liu
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - Q R Qadri
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - P Ma
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Dongchuan Road, Shanghai, East, 200240, China
| | - Q Wang
- Department of Animal Science, College of Animal Sciences, Zhejiang University, Yuhangtang Road, Hangzhou, East, 310058, China
| | - Y Pan
- Department of Animal Science, College of Animal Sciences, Zhejiang University, Yuhangtang Road, Hangzhou, East, 310058, China
| |
Collapse
|
13
|
Kumar S, Deng CH, Hunt M, Kirk C, Wiedow C, Rowan D, Wu J, Brewer L. Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear ( Pyrus spp.). FRONTIERS IN PLANT SCIENCE 2021; 11:590846. [PMID: 33469460 PMCID: PMC7813798 DOI: 10.3389/fpls.2020.590846] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/03/2020] [Accepted: 12/04/2020] [Indexed: 06/12/2023]
Abstract
Runs of homozygosity (ROH) have been widely used to study population history and trait architecture in humans and livestock species, but their application in self-incompatible plants has not been reported. The distributions of ROH in 199 accessions representing Asian pears (45), European pears (109), and interspecific hybrids (45) were investigated using genotyping-by-sequencing in this study. Fruit phenotypes including fruit weight, firmness, Brix, titratable acidity, and flavor volatiles were measured for genotype-phenotype analyses. The average number of ROH and the average total genomic length of ROH were 6 and 11 Mb, respectively, in Asian accessions, and 13 and 30 Mb, respectively, in European accessions. Significant associations between genomic inbreeding coefficients (FROH) and phenotypes were observed for 23 out of 32 traits analyzed. An overlap between ROH islands and significant markers from genome-wide association analyses was observed. Previously published quantitative trait loci for fruit traits and disease resistances also overlapped with some of the ROH islands. A prominent ROH island at the bottom of linkage group 17 overlapped with a recombination-supressed genomic region harboring the self-incompatibility locus. The observed ROH patterns suggested that systematic breeding of European pears would have started earlier than of Asian pears. Our research suggest that FROH would serve as a novel tool for managing inbreeding in gene-banks of self-incompatible plant species. ROH mapping provides a complementary strategy to unravel the genetic architecture of complex traits, and to evaluate differential selection in outbred plants. This seminal work would provide foundation for the ROH research in self-incompatible plants.
Collapse
Affiliation(s)
- Satish Kumar
- Hawke’s Bay Research Centre, The New Zealand Institute for Plant and Food Research Limited, Havelock North, New Zealand
| | - Cecilia Hong Deng
- Mount Albert Research Centre, The New Zealand Institute for Plant and Food Research Limited, Auckland, New Zealand
| | - Martin Hunt
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Chris Kirk
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Claudia Wiedow
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Daryl Rowan
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Jun Wu
- Centre of Pear Engineering Technology Research, Nanjing Agricultural University, Nanjing, China
| | - Lester Brewer
- Motueka Research Centre, The New Zealand Institute for Plant and Food Research Limited, Motueka, New Zealand
| |
Collapse
|
14
|
Liu Y, Xu L, Yang L, Zhao G, Li J, Liu D, Li Y. Discovery of Genomic Characteristics and Selection Signatures in Southern Chinese Local Cattle. Front Genet 2020; 11:533052. [PMID: 33391332 PMCID: PMC7775540 DOI: 10.3389/fgene.2020.533052] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2020] [Accepted: 11/27/2020] [Indexed: 01/20/2023] Open
Abstract
Chinese local cattle with a high level of genetic diversity mainly originate from two subspecies; the cattle in northern China are primarily Bos Taurus, and the cattle in southern China are primarily Bos indicus. Cattle from southern China are characterized by a specific phenotype and adapted to the local environment. This study explored the genetic diversity, degree of admixture, and selection signature in eight local cattle breeds in southern China. The lowest level of heterozygosity was found in Hainan and Nandan cattle from Hainan and Guangxi province, respectively, whereas the highest level of heterozygosity was detected in Zhaotong cattle from Yunnan province. A neighbor-joining phylogenetic tree analysis clearly separated Lufeng cattle from other breeds, whereas Leiqiong and Hainan cattle have some crossover. Based on linkage disequilibrium-filtered single nucleotide polymorphisms (SNPs), the admixture analysis revealed two clusters corresponding to the taurine and indicine cattle lineages, and the local cattle breeds from southern China showed a certain degree of admixture. When K = 4 and 9, we found a slight separation among Leiqiong, Lufeng, and Hainan cattle. Meanwhile, we performed a selection signature analysis in Hainan, Leiqiong, and Lufeng cattle distributed in the extreme south of China, using the integrated haplotype score (iHS), Rsb statistic, and BayeScan software. Using the iHS approach, we identified 251, 270, and 256 candidate regions in Lufeng, Leiqiong, and Hainan cattle, respectively. Moreover, we identified 184, 174, and 146 candidate regions in pairwise comparisons of Leiqiong vs. Lufeng, Leiqiong vs. Hainan, and Hainan vs. Lufeng cattle using the Rsb approach. In addition, we identified 76 loci with a total of 48 genes under selection, based on the FST approach. Several candidate genes under selection were found to be related to meat quality, immunity, and adaptation to the local environment in southern China. Our results provide significant information about the genetic differences among the cattle breeds from southern China and the possible cause of difference in breed-specific characteristics. Selection signature analysis identified a few candidate SNPs and genes related to certain important traits of these cattle. In general, our results provide valuable insights into the genetic basis of specific traits under selection in certain local cattle breeds.
Collapse
Affiliation(s)
- Yuqiang Liu
- College of Animal Science, South China Agricultural University, Guangzhou, China.,Innovation Team of Cattle Genetic Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China.,Lingnan Guangdong Laboratory of Modern Agriculture, Guangzhou, China
| | - Lingyang Xu
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Liu Yang
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, China
| | - Guoyao Zhao
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Junya Li
- Innovation Team of Cattle Genetic Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Dewu Liu
- College of Animal Science, South China Agricultural University, Guangzhou, China.,Lingnan Guangdong Laboratory of Modern Agriculture, Guangzhou, China
| | - Yaokun Li
- College of Animal Science, South China Agricultural University, Guangzhou, China.,Lingnan Guangdong Laboratory of Modern Agriculture, Guangzhou, China
| |
Collapse
|
15
|
Lee YL, Bosse M, Mullaart E, Groenen MAM, Veerkamp RF, Bouwman AC. Functional and population genetic features of copy number variations in two dairy cattle populations. BMC Genomics 2020; 21:89. [PMID: 31992181 PMCID: PMC6988284 DOI: 10.1186/s12864-020-6496-1] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2019] [Accepted: 01/14/2020] [Indexed: 12/27/2022] Open
Abstract
BACKGROUND Copy Number Variations (CNVs) are gain or loss of DNA segments that are known to play a role in shaping a wide range of phenotypes. In this study, we used two dairy cattle populations, Holstein Friesian and Jersey, to discover CNVs using the Illumina BovineHD Genotyping BeadChip aligned to the ARS-UCD1.2 assembly. The discovered CNVs were investigated for their functional impact and their population genetics features. RESULTS We discovered 14,272 autosomal CNVs, which were aggregated into 1755 CNV regions (CNVR) from 451 animals. These CNVRs together cover 2.8% of the bovine autosomes. The assessment of the functional impact of CNVRs showed that rare CNVRs (MAF < 0.01) are more likely to overlap with genes, than common CNVRs (MAF ≥ 0.05). The Population differentiation index (Fst) based on CNVRs revealed multiple highly diverged CNVRs between the two breeds. Some of these CNVRs overlapped with candidate genes such as MGAM and ADAMTS17 genes, which are related to starch digestion and body size, respectively. Lastly, linkage disequilibrium (LD) between CNVRs and BovineHD BeadChip SNPs was generally low, close to 0, although common deletions (MAF ≥ 0.05) showed slightly higher LD (r2 = ~ 0.1 at 10 kb distance) than the rest. Nevertheless, this LD is still lower than SNP-SNP LD (r2 = ~ 0.5 at 10 kb distance). CONCLUSIONS Our analyses showed that CNVRs detected using BovineHD BeadChip arrays are likely to be functional. This finding indicates that CNVs can potentially disrupt the function of genes and thus might alter phenotypes. Also, the population differentiation index revealed two candidate genes, MGAM and ADAMTS17, which hint at adaptive evolution between the two populations. Lastly, low CNVR-SNP LD implies that genetic variation from CNVs might not be fully captured in routine animal genetic evaluation, which relies solely on SNP markers.
Collapse
Affiliation(s)
- Young-Lim Lee
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands.
| | - Mirte Bosse
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| | | | - Martien A M Groenen
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| | - Roel F Veerkamp
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| | - Aniek C Bouwman
- Wageningen University & Research, Animal Breeding and Genomics, P.O. Box 338, Wageningen, AH, 6700, the Netherlands
| |
Collapse
|
16
|
Qanbari S. On the Extent of Linkage Disequilibrium in the Genome of Farm Animals. Front Genet 2020; 10:1304. [PMID: 32010183 PMCID: PMC6978288 DOI: 10.3389/fgene.2019.01304] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2019] [Accepted: 11/26/2019] [Indexed: 11/13/2022] Open
Abstract
Given the importance of linkage disequilibrium (LD) in gene mapping and evolutionary inferences, I characterize in this review the pattern of LD and discuss the influence of human intervention during domestication, breed establishment, and subsequent genetic improvement on shaping the genome of livestock species. To this end, I summarize data on the profile of LD based on array genotypes vs. sequencing data in cattle and chicken, two major livestock species, and compare to the human case. This comparison provides insights into the real dimension of the pairwise allelic correlation and haplo-block structuring. The dependency of LD on allelic frequency is pictured and a recently introduced metric for moderating it is outlined. In the context of the contact farm animals had with human, the impact of genetic forces including admixture, mutation, recombination rate, selection, and effective population size on LD is discussed. The review further highlights the interplay of LD with runs of homozygosity and concludes with the operational implications of the widely used association and selection mapping studies in relation to LD.
Collapse
Affiliation(s)
- Saber Qanbari
- Leibniz Institute for Farm Animal Biology (FBN), Institute of Genetics and Biometry, Dummerstorf, Germany.,Animal Breeding and Genetics Group, Department of Animal Sciences, Center for Integrated Breeding Research, University of Göttingen, Göttingen, Germany
| |
Collapse
|
17
|
Whole genome detection of recent selection signatures in Sarabi cattle: a unique Iranian taurine breed. Genes Genomics 2019; 42:203-215. [DOI: 10.1007/s13258-019-00888-6] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2019] [Accepted: 11/12/2019] [Indexed: 12/13/2022]
|
18
|
Kim JY, Jeong S, Kim KH, Lim WJ, Lee HY, Jeong N, Moon JK, Kim N. Dissection of soybean populations according to selection signatures based on whole-genome sequences. Gigascience 2019; 8:giz151. [PMID: 31869408 PMCID: PMC6927394 DOI: 10.1093/gigascience/giz151] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2019] [Revised: 08/21/2019] [Accepted: 12/05/2019] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND Domestication and improvement processes, accompanied by selections and adaptations, have generated genome-wide divergence and stratification in soybean populations. Simultaneously, soybean populations, which comprise diverse subpopulations, have developed their own adaptive characteristics enhancing fitness, resistance, agronomic traits, and morphological features. The genetic traits underlying these characteristics play a fundamental role in improving other soybean populations. RESULTS This study focused on identifying the selection signatures and adaptive characteristics in soybean populations. A core set of 245 accessions (112 wild-type, 79 landrace, and 54 improvement soybeans) selected from 4,234 soybean accessions was re-sequenced. Their genomic architectures were examined according to the domestication and improvement, and accessions were then classified into 3 wild-type, 2 landrace, and 2 improvement subgroups based on various population analyses. Selection and gene set enrichment analyses revealed that the landrace subgroups have selection signals for soybean-cyst nematode HG type 0 and seed development with germination, and that the improvement subgroups have selection signals for plant development with viability and seed development with embryo development, respectively. The adaptive characteristic for soybean-cyst nematode was partially underpinned by multiple resistance accessions, and the characteristics related to seed development were supported by our phenotypic findings for seed weights. Furthermore, their adaptive characteristics were also confirmed as genome-based evidence, and unique genomic regions that exhibit distinct selection and selective sweep patterns were revealed for 13 candidate genes. CONCLUSIONS Although our findings require further biological validation, they provide valuable information about soybean breeding strategies and present new options for breeders seeking donor lines to improve soybean populations.
Collapse
Affiliation(s)
- Jae-Yoon Kim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Gwahak-ro 125, Yuseong-gu, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Gajeong-ro 217, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Seongmun Jeong
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Gwahak-ro 125, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Kyoung Hyoun Kim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Gwahak-ro 125, Yuseong-gu, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Gajeong-ro 217, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Won-Jun Lim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Gwahak-ro 125, Yuseong-gu, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Gajeong-ro 217, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Ho-Yeon Lee
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Gwahak-ro 125, Yuseong-gu, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Gajeong-ro 217, Yuseong-gu, Daejeon 34141, Republic of Korea
| | - Namhee Jeong
- National Institute of Crop Science, Rural Development Administration, Nongsaengmyeong-ro 370, Deokjin-gu, Jeon-Ju 54874, Republic of Korea
| | - Jung-Kyung Moon
- National Institute of Crop Science, Rural Development Administration, Nongsaengmyeong-ro 370, Deokjin-gu, Jeon-Ju 54874, Republic of Korea
| | - Namshin Kim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Gwahak-ro 125, Yuseong-gu, Daejeon 34141, Republic of Korea
- Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Gajeong-ro 217, Yuseong-gu, Daejeon 34141, Republic of Korea
| |
Collapse
|
19
|
Yin T, Jaeger M, Scheper C, Grodkowski G, Sakowski T, Klopčič M, Bapst B, König S. Multi-breed genome-wide association studies across countries for electronically recorded behavior traits in local dual-purpose cows. PLoS One 2019; 14:e0221973. [PMID: 31665138 PMCID: PMC6821105 DOI: 10.1371/journal.pone.0221973] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2019] [Accepted: 08/16/2019] [Indexed: 12/20/2022] Open
Abstract
Basic bovine behavior is a crucial parameter influencing cattle domestication. In addition, behavior has an impact on cattle productivity, welfare and adaptation. The aim of the present study was to infer quantitative genetic and genomic mechanisms contributing to natural dual-purpose cow behavior in grazing systems. In this regard, we genotyped five dual-purpose breeds for a dense SNP marker panel from four different European countries. All cows from the across-country study were equipped with the same electronic recording devices. In this regard, we analyzed 97,049 longitudinal sensor behavior observations from 319 local dual-purpose cows for rumination, feeding, basic activity, high active, not active and ear temperature. According to the specific sensor behaviors and following a welfare protocol, we computed two different welfare indices. For genomic breed characterizations and multi-breed genome-wide association studies, sensor traits and test-day production records were merged with 35,826 SNP markers per cow. For the estimation of variance components, we used the pedigree relationship matrix and a combined similarity matrix that simultaneously included both pedigree and genotypes. Heritabilities for feeding, high active and not active were in a moderate range from 0.16 to 0.20. Estimates were very similar from both relationship matrix-modeling approaches and had quite small standard errors. Heritabilities for the remaining sensor traits (feeding, basic activity, ear temperature) and welfare indices were lower than 0.09. Five significant SNPs on chromosomes 11, 17, 27 and 29 were associated with rumination, and two different SNPs significantly influenced the sensor traits “not active” (chromosome 13) and “feeding” (chromosome 23). Gene annotation analyses inferred 22 potential candidate genes with a false discovery rate lower than 20%, mostly associated with rumination (13 genes) and feeding (8 genes). Mendelian randomization based on genomic variants (i.e., the instrumental variables) was used to infer causal inference between an exposure and an outcome. Significant regression coefficients among behavior traits indicate that all specific behavioral mechanisms contribute to similar physiological processes. The regression coefficients of rumination and feeding on milk yield were 0.10 kg/% and 0.12 kg/%, respectively, indicating their positive influence on dual-purpose cow productivity. Genomically, an improved welfare behavior of grazing cattle, i.e., a higher score for welfare indices, was significantly associated with increased fat and protein percentages.
Collapse
Affiliation(s)
- Tong Yin
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, Gießen, Germany
| | - Maria Jaeger
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, Gießen, Germany
| | - Carsten Scheper
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, Gießen, Germany
| | - Gregorz Grodkowski
- Institute of Genetics and Animal Breeding of the Polish Academy of Sciences, Jastrzębiec, Poland
| | - Tomasz Sakowski
- Institute of Genetics and Animal Breeding of the Polish Academy of Sciences, Jastrzębiec, Poland
| | - Marija Klopčič
- University of Ljubljana, Biotechnical Faculty, Department of Animal Science, Domzale, Slovenia
| | - Beat Bapst
- Genetic evaluation center, Qualitas AG, Switzerland
| | - Sven König
- Institute of Animal Breeding and Genetics, Justus-Liebig-University of Gießen, Gießen, Germany
- * E-mail:
| |
Collapse
|
20
|
Ablondi M, Viklund Å, Lindgren G, Eriksson S, Mikko S. Signatures of selection in the genome of Swedish warmblood horses selected for sport performance. BMC Genomics 2019; 20:717. [PMID: 31533613 PMCID: PMC6751828 DOI: 10.1186/s12864-019-6079-1] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2019] [Accepted: 09/04/2019] [Indexed: 01/09/2023] Open
Abstract
Background A growing demand for improved physical skills and mental attitude in modern sport horses has led to strong selection for performance in many warmblood studbooks. The aim of this study was to detect genomic regions with low diversity, and therefore potentially under selection, in Swedish Warmblood horses (SWB) by analysing high-density SNP data. To investigate if such signatures could be the result of selection for equestrian sport performance, we compared our SWB SNP data with those from Exmoor ponies, a horse breed not selected for sport performance traits. Results The genomic scan for homozygous regions identified long runs of homozygosity (ROH) shared by more than 85% of the genotyped SWB individuals. Such ROH were located on ECA4, ECA6, ECA7, ECA10 and ECA17. Long ROH were instead distributed evenly across the genome of Exmoor ponies in 77% of the chromosomes. Two population differentiation tests (FST and XP-EHH) revealed signatures of selection on ECA1, ECA4, and ECA6 in SWB horses. Conclusions Genes related to behaviour, physical abilities and fertility, appear to be targets of selection in the SWB breed. This study provides a genome-wide map of selection signatures in SWB horses, and ground for further functional studies to unravel the biological mechanisms behind complex traits in horses.
Collapse
Affiliation(s)
- Michela Ablondi
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden.,Department of Veterinary Science, Università degli Studi di Parma, 43126, Parma, Italy
| | - Åsa Viklund
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden
| | - Gabriella Lindgren
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden.,Livestock Genetics, Department of Biosystems, Leuven, KU, Belgium
| | - Susanne Eriksson
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden
| | - Sofia Mikko
- Dept. of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, PO Box 7023, S-750 07, Uppsala, Sweden.
| |
Collapse
|
21
|
Kim JY, Jeong S, Kim KH, Lim WJ, Lee HY, Kim N. Discovery of Genomic Characteristics and Selection Signatures in Korean Indigenous Goats Through Comparison of 10 Goat Breeds. Front Genet 2019; 10:699. [PMID: 31440273 PMCID: PMC6694180 DOI: 10.3389/fgene.2019.00699] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2018] [Accepted: 07/03/2019] [Indexed: 12/16/2022] Open
Abstract
Indigenous breeds develop their own genomic characteristics by adapting to local environments or cultures over long periods of time. Most of them are not particularly productive in commercial terms, but they have abilities to survive in harsh environments or tolerate to specific diseases. Their adaptive characteristics play an important role as genetic materials for improving commercial breeds. As a step toward this goal, we analyzed the genome of Korean indigenous goats within 10 goat breeds. We collected 136 goat individuals by sequencing 46 new goats and employing 90 publicly available goats. Our whole-genome data was comprised of three indigenous breeds (Korean indigenous goat, Iranian indigenous goat, and Moroccan indigenous goat; n = 29, 18, 20), six commercial breeds (Saanen, Boer, Anglo-Nubian, British Alpine, Alpine, and Korean crossbred; n = 16, 11, 5, 5, 2, 13), and their ancestral species (Capra aegagrus; n = 17). We identified that the Iranian indigenous goat and the Moroccan indigenous goat have relatively similar genomic characteristics within a large category of genomic diversity but found that the Korean indigenous goat has unique genomic characteristics distinguished from the other nine breeds. Through population analysis, we confirmed that these characteristics have resulted from a near-isolated environment with strong genetic drift. The Korean indigenous goat experienced a severe genetic bottleneck upon entering the Korean Peninsula about 2,000 years ago, and has subsequently rarely experienced genetic interactions with other goat breeds. From selection analysis and gene-set enrichment analysis, we revealed selection signals for Salmonella infection and cardiomyopathy in the genome of the Korean indigenous goat. These adaptive characteristics were further identified with genomic-based evidence. We uncovered genomic regions of selective sweeps in the LBP and BPI genes (Salmonella infection) and the TTN and ITGB6 genes (cardiomyopathy), among several candidate genes. Our research presents unique genomic characteristics and distinctive selection signals of the Korean indigenous goat based on the extensive comparison. Although the adaptive traits require further validation through biological experiments, our findings are expected to provide a direction for future biodiversity conservation strategies and to contribute another option to genomic-based breeding programmes for improving the viability of Capra hircus.
Collapse
Affiliation(s)
- Jae-Yoon Kim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, South Korea.,Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Daejeon, South Korea
| | - Seongmun Jeong
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, South Korea
| | - Kyoung Hyoun Kim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, South Korea.,Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Daejeon, South Korea
| | - Won-Jun Lim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, South Korea.,Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Daejeon, South Korea
| | - Ho-Yeon Lee
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, South Korea.,Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Daejeon, South Korea
| | - Namshin Kim
- Genome Editing Research Center, Korea Research Institute of Bioscience and Biotechnology (KRIBB), Daejeon, South Korea.,Department of Bioinformatics, KRIBB School of Bioscience, University of Science and Technology (UST), Daejeon, South Korea
| |
Collapse
|
22
|
Genome-wide scan reveals genetic divergence and diverse adaptive selection in Chinese local cattle. BMC Genomics 2019; 20:494. [PMID: 31200634 PMCID: PMC6570941 DOI: 10.1186/s12864-019-5822-y] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 05/21/2019] [Indexed: 01/18/2023] Open
Abstract
Background Understanding the population structure and genetic bases of well-adapted cattle breeds to local environments is one of the most essential tasks to develop appropriate genetic improvement programs. Results We performed a comprehensive study to investigate the population structure, divergence and selection signatures at genome-wide level in diverse Chinese local cattle using Bovine HD SNPs array, including two breeds from North China, one breed from Northwest China, three breeds from Southwest China and two breeds from South China. Population genetic analyses revealed the genetic structures of these populations were mostly related to the geographic locations. Notably, we detected 294 and 1263 candidate regions under selection using the di and iHS approaches, respectively. A series of group-specific and breed-specific candidate genes were identified, which are involved in immune response, sexual maturation, stature related, birth and bone weight, embryonic development, coat colors and adaptation. Furthermore, haplotype diversity and network pattern for candidate genes, including LPGAT1, LCORL, PPP1R8, RXFP2 and FANCA, suggest that these genes have been under differential selection pressures in various environmental conditions. Conclusions Our results shed insights into diverse selection during breed formation in Chinese local cattle. These findings may promote the application of genome-assisted breeding for well-adapted local breeds with economic and ecological importance. Electronic supplementary material The online version of this article (10.1186/s12864-019-5822-y) contains supplementary material, which is available to authorized users.
Collapse
|
23
|
Wang Z, Ma H, Xu L, Zhu B, Liu Y, Bordbar F, Chen Y, Zhang L, Gao X, Gao H, Zhang S, Xu L, Li J. Genome-Wide Scan Identifies Selection Signatures in Chinese Wagyu Cattle Using a High-Density SNP Array. Animals (Basel) 2019; 9:ani9060296. [PMID: 31151238 PMCID: PMC6617538 DOI: 10.3390/ani9060296] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2019] [Revised: 05/23/2019] [Accepted: 05/24/2019] [Indexed: 12/31/2022] Open
Abstract
Selective breeding can lead to genetic diversity and diverse phenotypes in farm animals. Analysis of the genomic regions under selection can provide important insights into the genetic basis of complex traits. In this study, a high-density SNP array was used for analysis of genome selection signatures in Chinese Wagyu cattle. In total, we obtained 478,903 SNPs and 24,820 no-overlap regions for |iHS| (integrated haplotype score) estimations. Under the threshold of the top 1%, 239 regions were finally identified as candidate selected regions and 162 candidate genes were found based on the UMD3.1 genome assembly. These genes were reported to be associated with fatty acids, such as Bos taurus nitric oxide synthase 1 adaptor protein (NOS1AP), Bos taurus hydroxysteroid 17-beta dehydrogenase 7 (HSD17B7), Bos taurus WD repeat domain 7 (WDR7), Bos taurus ELOVL fatty acid elongase 2 (ELOVL2), Bos taurus calpain 1 (CAPN1), Bos taurus parkin RBR E3 ubiquitin protein ligase (PRKN, also known as PARK2), Bos taurus mitogen-activated protein kinase kinase 6 (MAP2K6), meat quality, including Bos taurus ADAM metallopeptidase domain 12 (ADAM12), Bos taurus 5'-aminolevulinate synthase 1 (ALAS1), Bos taurus small integral membrane protein 13 (SMIM13) and Bos taurus potassium two pore domain channel subfamily K member 2 (KCNK2), growth, and developmental traits, such as Bos taurus insulin like growth factor 2 receptor (IGF2R), Bos taurus RAR related orphan receptor A (RORA), Bos taurus fibroblast growth factor 14 (FGF14), Bos taurus paired box 6 (PAX6) and Bos taurus LIM homeobox 6 (LHX6). In addition, we identified several genes that are associated with body size and weight, including Bos taurus sorting nexin 29 (SNX29), Bos taurus zinc finger imprinted 2 (ZIM2), Bos taurus family with sequence similarity 110 member A (FAM110A), immune system, including Bos taurus toll like receptor 9 (TLR9), Bos taurus TAFA chemokine like family member 1 (TAFA1), Bos taurus glutathione peroxidase 8 (putative) (GPX8), Bos taurus interleukin 5 (IL5), Bos taurus PR domain containing 9 (PRDM9), Bos taurus glutamate ionotropic receptor kainate type subunit 2 (GRIK2) and feed intake efficiency, Bos taurus sodium voltage-gated channel alpha subunit 9 (SCN9A), Bos taurus relaxin family peptide/INSL5 receptor 4 (RXFP4), Bos taurus RNA polymerase II associated protein 3 (RPAP3). Moreover, four GO terms of biological regulation (GO:0009987, GO:0008152) and metabolic process (GO:0003824, GO:0005488) were found based on these genes. In addition, we found that 232 candidate regions (~18 Mb) overlapped with the Quantitative trait loci (QTL)regions extracted from cattle QTLdb. Our findings imply that many genes were selected for important traits in Chinese Wagyu cattle. Moreover, these results can contribute to the understanding of the genetic basis of the studied traits during the formation of this population.
Collapse
Affiliation(s)
- Zezhao Wang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China.
| | - Haoran Ma
- College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, China.
| | - Lei Xu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China.
- Institute of Animal Husbandry and Veterinary Research, Anhui Academy of Agricultural Sciences, Hefei 230031, China.
| | - Bo Zhu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Ying Liu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Farhad Bordbar
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Yan Chen
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Lupei Zhang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Xue Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Huijiang Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Shengli Zhang
- National Engineering Laboratory for Animal Breeding, Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China.
| | - Lingyang Xu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| | - Junya Li
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing 100193, China.
| |
Collapse
|
24
|
Ma L, Sonstegard TS, Cole JB, VanTassell CP, Wiggans GR, Crooker BA, Tan C, Prakapenka D, Liu GE, Da Y. Genome changes due to artificial selection in U.S. Holstein cattle. BMC Genomics 2019; 20:128. [PMID: 30744549 PMCID: PMC6371544 DOI: 10.1186/s12864-019-5459-x] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2018] [Accepted: 01/17/2019] [Indexed: 01/13/2023] Open
Abstract
Background The availability of a unique unselected Holstein line since 1964 provided a direct comparison between selected and unselected Holstein genomes whereas large Holstein samples provided unprecedented statistical power for identifying high-confidence SNP effects. Utilizing these unique resources, we aimed to identify genome changes affected by selection since 1964. Results Direct comparison of genome-wide SNP markers between a Holstein line unselected since 1964 and contemporary Holsteins showed that the 40 years of artificial selection since 1964 resulted in genome landscape changes. Among the regions affected by selection, the regions containing 198 genes with fertility functions had a larger negative correlation than that of all SNPs between the SNP effects on milk yield and daughter pregnancy rate. These results supported the hypothesis that hitchhiking of genetic selection for milk production by negative effects of fertility genes contributed to the unintended declines in fertility since 1964. The genome regions subjected to selection also contained 67 immunity genes, the bovine MHC region of Chr23 with significantly decreased heterozygosity in contemporary Holsteins, and large gene clusters including T-cell receptor and immunoglobulin genes. Conclusions This study for the first time provided direct evidence that genetic selection for milk production affected fertility and immunity genes and that the hitchhiking of genetic selection for milk production by negative fertility effects contributed to the fertility declines since 1964, and identified a large number of candidate fertility and immunity genes affected by selection. The results provided novel understanding about genome changes due to artificial selection and their impact on fertility and immunity genes and could facilitate developing genetic methods to reverse the declines in fertility and immunity in Holstein cattle. Electronic supplementary material The online version of this article (10.1186/s12864-019-5459-x) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Li Ma
- Department of Animal and Avian Sciences, University of Maryland, College Park, MD, USA
| | | | - John B Cole
- Animal Genomics and Improvement Laboratory, USDA-ARS, Beltsville, MD, USA
| | | | - George R Wiggans
- Animal Genomics and Improvement Laboratory, USDA-ARS, Beltsville, MD, USA
| | - Brian A Crooker
- Department of Animal Science, University of Minnesota, Saint Paul, MN, USA
| | - Cheng Tan
- College of Animal Science and National Engineering Research Center for Breeding Swine Industry, South China Agricultural University, Guangzhou, People's Republic of China
| | - Dzianis Prakapenka
- Department of Animal Science, University of Minnesota, Saint Paul, MN, USA
| | - George E Liu
- Animal Genomics and Improvement Laboratory, USDA-ARS, Beltsville, MD, USA
| | - Yang Da
- Department of Animal Science, University of Minnesota, Saint Paul, MN, USA.
| |
Collapse
|
25
|
Johnsson M. Integrating Selection Mapping With Genetic Mapping and Functional Genomics. Front Genet 2018; 9:603. [PMID: 30619447 PMCID: PMC6295561 DOI: 10.3389/fgene.2018.00603] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2018] [Accepted: 11/19/2018] [Indexed: 01/23/2023] Open
Abstract
Genomic scans for signatures of selection allow us to, in principle, detect variants and genes that underlie recent adaptations. By combining selection mapping with genetic mapping of traits known to be relevant to adaptation, we can simultaneously investigate whether genes and variants show signals of recent selection and whether they impact traits that have likely been selected. There are three ways to integrate selection mapping with genetic mapping or functional genomics: (1) To use genetic mapping data from other populations as a form of genome annotation. (2) To perform experimental evolution or artificial selection to be able to study selected variants when they segregate, either by performing genetic mapping before selection or by crossing the selected individuals to some reference population. (3) To perform a comparative study of related populations facing different selection regimes. This short review discusses these different ways of integrating selection mapping with genetic mapping and functional genomics, with examples of how each has been done.
Collapse
Affiliation(s)
- Martin Johnsson
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, The University of Edinburgh, Edinburgh, United Kingdom.,Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
26
|
A genome-wide detection of selection signatures in conserved and commercial pig breeds maintained in Poland. BMC Genet 2018; 19:95. [PMID: 30348079 PMCID: PMC6198424 DOI: 10.1186/s12863-018-0681-0] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2018] [Accepted: 10/03/2018] [Indexed: 01/19/2023] Open
Abstract
BACKGROUND Identification of selection signatures can provide a direct insight into the mechanism of artificial selection and allow further disclosure of the candidate genes related to the animals' phenotypic variation. Domestication and subsequent long-time selection have resulted in extensive phenotypic changes in domestic pigs, involving a number of traits, like behavior, body composition, disease resistance, reproduction and coat color. In this study, based on genotypes obtained from PorcineSNP60 Illumina assay we attempt to detect both diversifying and within-breed selection signatures in 530 pigs belonging to four breeds: Polish Landrace, Puławska, Złotnicka White and Złotnicka Spotted, of which the last three are a subject of conservative breeding and substantially represent the native populations. RESULTS A two largely complementary statistical methods were used for signatures detection, including: pairwise FST and relative extended haplotype homozygosity (REHH) test. Breed-specific diversifying selection signals included several genes involved in processes connected with fertility, growth and metabolism which are potentially responsible for different phenotypes of the studied breeds. The diversifying selection signals also comprised PPARD gene that was previously found to have a large effect on the shape of the external ear in pigs or two genes encoding neuropeptide Y receptors (Y2 and Y5) involved in fat deposition and stress response which are important features differentiating the studied breeds. REHH statistics allowed detecting several within-breed selection signatures overlapping with genes connected with a range of functions including, among others: metabolic pathways, immune system response or implantation and development of the embryo. CONCLUSIONS The study provides many potential candidate genes with implication for traits selected in the individual breeds and gives strong basis for further studies aiming at identification of sources of variation among the studied pig breeds.
Collapse
|
27
|
Hanuš O, Samková E, Křížová L, Hasoňová L, Kala R. Role of Fatty Acids in Milk Fat and the Influence of Selected Factors on Their Variability-A Review. Molecules 2018; 23:E1636. [PMID: 29973572 PMCID: PMC6100482 DOI: 10.3390/molecules23071636] [Citation(s) in RCA: 95] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2018] [Revised: 06/29/2018] [Accepted: 07/02/2018] [Indexed: 11/16/2022] Open
Abstract
Fatty acids (FAs) of milk fat are considered to be important nutritional components of the diets of a significant portion of the human population and substantially affect human health. With regard to dairy farming, the FA profile is also seen as an important factor in the technological quality of raw milk. In this sense, making targeted modifications to the FA profile has the potential to significantly contribute to the production of dairy products with higher added value. Thus, FAs also have economic importance. Current developments in analytical methods and their increasing efficiency enable the study of FA profiles not only for scientific purposes but also in terms of practical technological applications. It is important to study the sources of variability of FAs in milk, which include population genetics, type of farming, and targeted animal nutrition. It is equally important to study the health and technological impacts of FAs. This review summarizes current knowledge in the field regarding sources of FA variability, including the impact of factors such as: animal nutrition, seasonal feed changes, type of animal farming (conventional and organic), genetic parameters (influence of breed), animal individuality, lactation, and milk yield. Potential practical applications (to improve food technology and consumer health) of FA profile information are also reviewed.
Collapse
Affiliation(s)
- Oto Hanuš
- Dairy Research Institute Ltd., 16000 Prague, Czech Republic.
| | - Eva Samková
- Department of Food Biotechnologies and Agricultural Products´ Quality, Faculty of Agriculture, University of South Bohemia, 37005 České Budějovice, Czech Republic.
| | - Ludmila Křížová
- Department of Animal Nutrition, Faculty of Veterinary Hygiene and Ecology, University of Veterinary and Pharmaceutical Sciences Brno, 61242 Brno, Czech Republic.
| | - Lucie Hasoňová
- Department of Food Biotechnologies and Agricultural Products´ Quality, Faculty of Agriculture, University of South Bohemia, 37005 České Budějovice, Czech Republic.
| | - Robert Kala
- Department of Food Biotechnologies and Agricultural Products´ Quality, Faculty of Agriculture, University of South Bohemia, 37005 České Budějovice, Czech Republic.
| |
Collapse
|
28
|
Weaver SR, Hernandez LL. Could use of Selective Serotonin Reuptake Inhibitors During Lactation Cause Persistent Effects on Maternal Bone? J Mammary Gland Biol Neoplasia 2018; 23:5-25. [PMID: 29603039 DOI: 10.1007/s10911-018-9390-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/14/2017] [Accepted: 03/13/2018] [Indexed: 02/07/2023] Open
Abstract
The lactating mammary gland elegantly coordinates maternal homeostasis to provide calcium for milk. During lactation, the monoamine serotonin regulates the synthesis and release of various mammary gland-derived factors, such as parathyroid hormone-related protein (PTHrP), to stimulate bone resorption. Recent evidence suggests that bone mineral lost during prolonged lactation is not fully recovered following weaning, possibly putting women at increased risk of fracture or osteoporosis. Selective Serotonin Reuptake Inhibitor (SSRI) antidepressants have also been associated with reduced bone mineral density and increased fracture risk. Therefore, SSRI exposure while breastfeeding may exacerbate lactational bone loss, compromising long-term bone health. Through an examination of serotonin and calcium homeostasis during lactation, lactational bone turnover and post-weaning recovery of bone mineral, and the effect of peripartum depression and SSRI on the mammary gland and bone, this review will discuss the hypothesis that peripartum SSRI exposure causes persistent reductions in bone mineral density through mammary-derived PTHrP signaling with bone.
Collapse
Affiliation(s)
- Samantha R Weaver
- Endocrine and Reproductive Physiology Program, University of Wisconsin-Madison, Madison, WI, USA
| | - Laura L Hernandez
- Department of Dairy Science, University of Wisconsin-Madison, Madison, WI, USA.
| |
Collapse
|
29
|
Environmental selection during the last ice age on the mother-to-infant transmission of vitamin D and fatty acids through breast milk. Proc Natl Acad Sci U S A 2018; 115:E4426-E4432. [PMID: 29686092 PMCID: PMC5948952 DOI: 10.1073/pnas.1711788115] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
The frequency of the human-specific EDAR V370A isoform is highly elevated in North and East Asian populations. The gene is known to have several pleiotropic effects, among which are sweat gland density and ductal branching in the mammary gland. The former has led some geneticists to argue that the near-fixation of this allele was caused by selection for modulation of thermoregulatory sweating. We provide an alternative hypothesis, that selection instead acted on the allele’s effect of increasing ductal branching in the mammary gland, thereby amplifying the transfer of critical nutrients to infants via mother’s milk. This is likely to have occurred during the Last Glacial Maximum when a human population was genetically isolated in the high-latitude environment of the Beringia. Because of the ubiquitous adaptability of our material culture, some human populations have occupied extreme environments that intensified selection on existing genomic variation. By 32,000 years ago, people were living in Arctic Beringia, and during the Last Glacial Maximum (LGM; 28,000–18,000 y ago), they likely persisted in the Beringian refugium. Such high latitudes provide only very low levels of UV radiation, and can thereby lead to dangerously low levels of biosynthesized vitamin D. The physiological effects of vitamin D deficiency range from reduced dietary absorption of calcium to a compromised immune system and modified adipose tissue function. The ectodysplasin A receptor (EDAR) gene has a range of pleiotropic effects, including sweat gland density, incisor shoveling, and mammary gland ductal branching. The frequency of the human-specific EDAR V370A allele appears to be uniquely elevated in North and East Asian and New World populations due to a bout of positive selection likely to have occurred circa 20,000 y ago. The dental pleiotropic effects of this allele suggest an even higher occurrence among indigenous people in the Western Hemisphere before European colonization. We hypothesize that selection on EDAR V370A occurred in the Beringian refugium because it increases mammary ductal branching, and thereby may amplify the transfer of critical nutrients in vitamin D-deficient conditions to infants via mothers’ milk. This hypothesized selective context for EDAR V370A was likely intertwined with selection on the fatty acid desaturase (FADS) gene cluster because it is known to modulate lipid profiles transmitted to milk from a vitamin D-rich diet high in omega-3 fatty acids.
Collapse
|
30
|
Doekes HP, Veerkamp RF, Bijma P, Hiemstra SJ, Windig JJ. Trends in genome-wide and region-specific genetic diversity in the Dutch-Flemish Holstein-Friesian breeding program from 1986 to 2015. Genet Sel Evol 2018; 50:15. [PMID: 29642838 PMCID: PMC5896142 DOI: 10.1186/s12711-018-0385-y] [Citation(s) in RCA: 42] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2017] [Accepted: 03/27/2018] [Indexed: 12/19/2022] Open
Abstract
Background In recent decades, Holstein–Friesian (HF) selection schemes have undergone profound changes, including the introduction of optimal contribution selection (OCS; around 2000), a major shift in breeding goal composition (around 2000) and the implementation of genomic selection (GS; around 2010). These changes are expected to have influenced genetic diversity trends. Our aim was to evaluate genome-wide and region-specific diversity in HF artificial insemination (AI) bulls in the Dutch-Flemish breeding program from 1986 to 2015. Methods Pedigree and genotype data (~ 75.5 k) of 6280 AI-bulls were used to estimate rates of genome-wide inbreeding and kinship and corresponding effective population sizes. Region-specific inbreeding trends were evaluated using regions of homozygosity (ROH). Changes in observed allele frequencies were compared to those expected under pure drift to identify putative regions under selection. We also investigated the direction of changes in allele frequency over time. Results Effective population size estimates for the 1986–2015 period ranged from 69 to 102. Two major breakpoints were observed in genome-wide inbreeding and kinship trends. Around 2000, inbreeding and kinship levels temporarily dropped. From 2010 onwards, they steeply increased, with pedigree-based, ROH-based and marker-based inbreeding rates as high as 1.8, 2.1 and 2.8% per generation, respectively. Accumulation of inbreeding varied substantially across the genome. A considerable fraction of markers showed changes in allele frequency that were greater than expected under pure drift. Putative selected regions harboured many quantitative trait loci (QTL) associated to a wide range of traits. In consecutive 5-year periods, allele frequencies changed more often in the same direction than in opposite directions, except when comparing the 1996–2000 and 2001–2005 periods. Conclusions Genome-wide and region-specific diversity trends reflect major changes in the Dutch-Flemish HF breeding program. Introduction of OCS and the shift in breeding goal were followed by a drop in inbreeding and kinship and a shift in the direction of changes in allele frequency. After introduction of GS, rates of inbreeding and kinship increased substantially while allele frequencies continued to change in the same direction as before GS. These results provide insight in the effect of breeding practices on genomic diversity and emphasize the need for efficient management of genetic diversity in GS schemes. Electronic supplementary material The online version of this article (10.1186/s12711-018-0385-y) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Harmen P Doekes
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands. .,Centre for Genetic Resources the Netherlands, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands.
| | - Roel F Veerkamp
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Piter Bijma
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Sipke J Hiemstra
- Centre for Genetic Resources the Netherlands, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| | - Jack J Windig
- Animal Breeding and Genomics, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands.,Centre for Genetic Resources the Netherlands, Wageningen University & Research, P.O. Box 338, 6700 AH, Wageningen, The Netherlands
| |
Collapse
|
31
|
Signer-Hasler H, Burren A, Neuditschko M, Frischknecht M, Garrick D, Stricker C, Gredler B, Bapst B, Flury C. Population structure and genomic inbreeding in nine Swiss dairy cattle populations. Genet Sel Evol 2017; 49:83. [PMID: 29115934 PMCID: PMC5674839 DOI: 10.1186/s12711-017-0358-6] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2016] [Accepted: 10/26/2017] [Indexed: 12/25/2022] Open
Abstract
BACKGROUND Domestication, breed formation and intensive selection have resulted in divergent cattle breeds that likely exhibit their own genomic signatures. In this study, we used genotypes from 27,612 autosomal single nucleotide polymorphisms to characterize population structure based on 9214 sires representing nine Swiss dairy cattle populations: Brown Swiss (BS), Braunvieh (BV), Original Braunvieh (OB), Holstein (HO), Red Holstein (RH), Swiss Fleckvieh (SF), Simmental (SI), Eringer (ER) and Evolèner (EV). Genomic inbreeding (F ROH) and signatures of selection were determined by calculating runs of homozygosity (ROH). The results build the basis for a better understanding of the genetic development of Swiss dairy cattle populations and highlight differences between the original populations (i.e. OB, SI, ER and EV) and those that have become more popular in Switzerland as currently reflected by their larger populations (i.e. BS, BV, HO, RH and SF). RESULTS The levels of genetic diversity were highest and lowest in the SF and BS breeds, respectively. Based on F ST values, we conclude that, among all pairwise comparisons, BS and HO (0.156) differ more than the other pairs of populations. The original Swiss cattle populations OB, SI, ER, and EV are clearly genetically separated from the Swiss cattle populations that are now more common and represented by larger numbers of cows. Mean levels of F ROH ranged from 0.027 (ER) to 0.091 (BS). Three of the original Swiss cattle populations, ER (F ROH: 0.027), OB (F ROH: 0.029), and SI (F ROH: 0.039), showed low levels of genomic inbreeding, whereas it was much higher in EV (F ROH: 0.074). Private signatures of selection for the original Swiss cattle populations are reported for BTA4, 5, 11 and 26. CONCLUSIONS The low levels of genomic inbreeding observed in the original Swiss cattle populations ER, OB and SI compared to the other breeds are explained by a lesser use of artificial insemination and greater use of natural service. Natural service results in more sires having progeny at each generation and thus this breeding practice is likely the major reason for the remarkable levels of genetic diversity retained within these populations. The fact that the EV population is regionally restricted and its small census size of herd-book cows explain its high level of genomic inbreeding.
Collapse
Affiliation(s)
- Heidi Signer-Hasler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | - Alexander Burren
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| | | | - Mirjam Frischknecht
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
- Qualitas AG, Zug, Switzerland
| | | | | | | | | | - Christine Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, Zollikofen, Switzerland
| |
Collapse
|
32
|
Kasarapu P, Porto-Neto LR, Fortes MRS, Lehnert SA, Mudadu MA, Coutinho L, Regitano L, George A, Reverter A. The Bos taurus-Bos indicus balance in fertility and milk related genes. PLoS One 2017; 12:e0181930. [PMID: 28763475 PMCID: PMC5538644 DOI: 10.1371/journal.pone.0181930] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2017] [Accepted: 07/10/2017] [Indexed: 12/16/2022] Open
Abstract
Numerical approaches to high-density single nucleotide polymorphism (SNP) data are often employed independently to address individual questions. We linked independent approaches in a bioinformatics pipeline for further insight. The pipeline driven by heterozygosity and Hardy-Weinberg equilibrium (HWE) analyses was applied to characterize Bos taurus and Bos indicus ancestry. We infer a gene co-heterozygosity network that regulates bovine fertility, from data on 18,363 cattle with genotypes for 729,068 SNP. Hierarchical clustering separated populations according to Bos taurus and Bos indicus ancestry. The weights of the first principal component were subjected to Normal mixture modelling allowing the estimation of a gene’s contribution to the Bos taurus-Bos indicus axis. We used deviation from HWE, contribution to Bos indicus content and association to fertility traits to select 1,284 genes. With this set, we developed a co-heterozygosity network where the group of genes annotated as fertility-related had significantly higher Bos indicus content compared to other functional classes of genes, while the group of genes associated with milk production had significantly higher Bos taurus content. The network analysis resulted in capturing novel gene associations of relevance to bovine domestication events. We report transcription factors that are likely to regulate genes associated with cattle domestication and tropical adaptation. Our pipeline can be generalized to any scenarios where population structure requires scrutiny at the molecular level, particularly in the presence of a priori set of genes known to impact a phenotype of evolutionary interest such as fertility.
Collapse
Affiliation(s)
- Parthan Kasarapu
- CSIRO Agriculture and Food, Queensland Bioscience Precinct, St. Lucia, Brisbane, Queensland, Australia
| | - Laercio R. Porto-Neto
- CSIRO Agriculture and Food, Queensland Bioscience Precinct, St. Lucia, Brisbane, Queensland, Australia
| | - Marina R. S. Fortes
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
| | - Sigrid A. Lehnert
- CSIRO Agriculture and Food, Queensland Bioscience Precinct, St. Lucia, Brisbane, Queensland, Australia
| | | | - Luiz Coutinho
- Centro de Genomica Funcional ESALQ, University of São Paulo, Piracicaba, Sao Paulo, Brazil
| | - Luciana Regitano
- Embrapa Southeast Livestock, Rodovia Washington Luiz, São Carlos, Sao Paulo, Brazil
| | - Andrew George
- CSIRO, DATA61, Ecosciences Precinct Brisbane, Brisbane, Queensland, Australia
| | - Antonio Reverter
- CSIRO Agriculture and Food, Queensland Bioscience Precinct, St. Lucia, Brisbane, Queensland, Australia
- * E-mail:
| |
Collapse
|
33
|
Genome-Wide Detection of Selective Signatures in Chicken through High Density SNPs. PLoS One 2016; 11:e0166146. [PMID: 27820849 PMCID: PMC5098818 DOI: 10.1371/journal.pone.0166146] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2016] [Accepted: 10/24/2016] [Indexed: 11/25/2022] Open
Abstract
Chicken is recognized as an excellent model for studies of genetic mechanism of phenotypic and genomic evolution, with large effective population size and strong human-driven selection. In the present study, we performed Extended Haplotype Homozygosity (EHH) tests to identify significant core regions employing 600K SNP Chicken chip in an F2 population of 1,534 hens, which was derived from reciprocal crosses between White Leghorn and Dongxiang chicken. Results indicated that a total of 49,151 core regions with an average length of 9.79 Kb were identified, which occupied approximately 52.15% of genome across all autosomes, and 806 significant core regions attracted us mostly. Genes in candidate regions may experience positive selection and were considered to have possible influence on beneficial economic traits. A panel of genes including AASDHPPT, GDPD5, PAR3, SOX6, GPC1 and a signal pathway of AKT1 were detected with the most extreme P-values. Further enrichment analyses indicated that these genes were associated with immune function, sensory organ development and neurogenesis, and may have experienced positive selection in chicken. Moreover, some of core regions exactly overlapped with genes excavated in our previous GWAS, suggesting that these genes have undergone positive selection may affect egg production. Findings in our study could draw a comparatively integrate genome-wide map of selection signature in the chicken genome, and would be worthy for explicating the genetic mechanisms of phenotypic diversity in poultry breeding.
Collapse
|
34
|
Identification of selective sweeps reveals divergent selection between Chinese Holstein and Simmental cattle populations. Genet Sel Evol 2016; 48:76. [PMID: 27716022 PMCID: PMC5054554 DOI: 10.1186/s12711-016-0254-5] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2015] [Accepted: 09/26/2016] [Indexed: 12/23/2022] Open
Abstract
Background The identification of signals left by recent positive selection provides a feasible approach for targeting genomic variants that underlie complex traits and fitness. A better understanding of the selection mechanisms that occurred during the evolution of species can also be gained. In this study, we simultaneously detected the genome-wide footprints of recent positive selection that occurred within and between Chinese Holstein and Simmental populations, which have been subjected to artificial selection for distinct purposes. We conducted analyses using various complementary approaches, including LRH, XP-EHH and FST, based on the Illumina 770K high-density single nucleotide polymorphism (SNP) array, to enable more comprehensive detection. Results We successfully constructed profiles of selective signals in both cattle populations. To further annotate these regions, we identified a set of novel functional genes related to growth, reproduction, immune response and milk production. There were no overlapping candidate windows between the two breeds. Finally, we investigated the distribution of SNPs that had low FST values across five distinct functional regions in the genome. In the low-minor allele frequency bin, we found a higher proportion of low-FST SNPs in the exons of the bovine genome, which indicates strong purifying selection of the exons. Conclusions The selection signatures identified in these two populations demonstrated positive selection pressure on a set of important genes with potential functions that are involved in many biological processes. We also demonstrated that in the bovine genome, exons were under strong purifying selection. Our findings provide insight into the mechanisms of artificial selection and will facilitate follow-up functional studies of potential candidate genes that are related to various economically important traits in cattle. Electronic supplementary material The online version of this article (doi:10.1186/s12711-016-0254-5) contains supplementary material, which is available to authorized users.
Collapse
|
35
|
Fu W, Lee WR, Abasht B. Detection of genomic signatures of recent selection in commercial broiler chickens. BMC Genet 2016; 17:122. [PMID: 27565946 PMCID: PMC5002100 DOI: 10.1186/s12863-016-0430-1] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2016] [Accepted: 08/22/2016] [Indexed: 02/06/2023] Open
Abstract
Background Identification of the genomic signatures of recent selection may help uncover causal polymorphisms controlling traits relevant to recent decades of selective breeding in livestock. In this study, we aimed at detecting signatures of recent selection in commercial broiler chickens using genotype information from single nucleotide polymorphisms (SNPs). A total of 565 chickens from five commercial purebred lines, including three broiler sire (male) lines and two broiler dam (female) lines, were genotyped using the 60K SNP Illumina iSelect chicken array. To detect genomic signatures of recent selection, we applied two methods based on population comparison, cross-population extended haplotype homozygosity (XP-EHH) and cross-population composite likelihood ratio (XP-CLR), and further analyzed the results to find genomic regions under recent selection in multiple purebred lines. Results A total of 321 candidate selection regions spanning approximately 1.45 % of the chicken genome in each line were detected by consensus of results of both XP-EHH and XP-CLR methods. To minimize false discovery due to genetic drift, only 42 of the candidate selection regions that were shared by 2 or more purebred lines were considered as high-confidence selection regions in the study. Of these 42 regions, 20 were 50 kb or less while 4 regions were larger than 0.5 Mb. In total, 91 genes could be found in the 42 regions, among which 19 regions contained only 1 or 2 genes, and 9 regions were located at gene deserts. Conclusions Our results provide a genome-wide scan of recent selection signatures in five purebred lines of commercial broiler chickens. We found several candidate genes for recent selection in multiple lines, such as SOX6 (Sex Determining Region Y-Box 6) and cTR (Thyroid hormone receptor beta). These genes may have been under recent selection due to their essential roles in growth, development and reproduction in chickens. Furthermore, our results suggest that in some candidate regions, the same or opposite alleles have been under recent selection in multiple lines. Most of the candidate genes in the selection regions are novel, and as such they should be of great interest for future research into the genetic architecture of traits relevant to modern broiler breeding. Electronic supplementary material The online version of this article (doi:10.1186/s12863-016-0430-1) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Weixuan Fu
- Department of Animal and Food Sciences, University of Delaware, Newark, DE, 19716, USA
| | | | - Behnam Abasht
- Department of Animal and Food Sciences, University of Delaware, Newark, DE, 19716, USA.
| |
Collapse
|
36
|
Sorbolini S, Gaspa G, Steri R, Dimauro C, Cellesi M, Stella A, Marras G, Marsan PA, Valentini A, Macciotta NPP. Use of canonical discriminant analysis to study signatures of selection in cattle. Genet Sel Evol 2016; 48:58. [PMID: 27521154 PMCID: PMC4983034 DOI: 10.1186/s12711-016-0236-7] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2016] [Accepted: 08/01/2016] [Indexed: 01/08/2023] Open
Abstract
BACKGROUND Cattle include a large number of breeds that are characterized by marked phenotypic differences and thus constitute a valuable model to study genome evolution in response to processes such as selection and domestication. Detection of "signatures of selection" is a useful approach to study the evolutionary pressures experienced throughout history. In the present study, signatures of selection were investigated in five cattle breeds farmed in Italy using a multivariate approach. METHODS A total of 4094 bulls from five breeds with different production aptitudes (two dairy breeds: Italian Holstein and Italian Brown Swiss; two beef breeds: Piemontese and Marchigiana; and one dual purpose breed: Italian Simmental) were genotyped using the Illumina BovineSNP50 v.1 beadchip. Canonical discriminant analysis was carried out on the matrix of single nucleotide polymorphisms (SNP) genotyping data, separately for each chromosome. Scores for each canonical variable were calculated and then plotted in the canonical space to quantify the distance between breeds. SNPs for which the correlation with the canonical variable was in the 99th percentile for a specific chromosome were considered to be significantly associated with that variable. Results were compared with those obtained using an FST-based approach. RESULTS Based on the results of the canonical discriminant analysis, a large number of signatures of selection were detected, among which several had strong signals in genomic regions that harbour genes known to have an impact on production and morphological bovine traits, including MSTN, LCT, GHR, SCD, NCAPG, KIT, and ASIP. Moreover, new putative candidate genes were identified, such as GCK, B3GALNT1, MGAT1, GALNTL1, PRNP, and PRND. Similar results were obtained with the FST-based approach. CONCLUSIONS The use of canonical discriminant analysis on 50 K SNP genotypes allowed the extraction of new variables that maximize the separation between breeds. This approach is quite straightforward, it can compare more than two groups simultaneously, and relative distances between breeds can be visualized. The genes that were highlighted in the canonical discriminant analysis were in concordance with those obtained using the FST index.
Collapse
Affiliation(s)
- Silvia Sorbolini
- Dipartimento di Agraria, Sezione di Scienze Zootecniche, Università degli Studi di Sassari, V. le Italia, 9, 07100, Sassari, Italy
| | - Giustino Gaspa
- Dipartimento di Agraria, Sezione di Scienze Zootecniche, Università degli Studi di Sassari, V. le Italia, 9, 07100, Sassari, Italy
| | - Roberto Steri
- Consiglio per la Ricerca e la Sperimentazione in Agricoltura, via Salaria 31, 00015, Monterotondo, Italy
| | - Corrado Dimauro
- Dipartimento di Agraria, Sezione di Scienze Zootecniche, Università degli Studi di Sassari, V. le Italia, 9, 07100, Sassari, Italy
| | - Massimo Cellesi
- Dipartimento di Agraria, Sezione di Scienze Zootecniche, Università degli Studi di Sassari, V. le Italia, 9, 07100, Sassari, Italy
| | | | | | - Paolo Ajmone Marsan
- Istituto di Zootecnica, Università Cattolica del Sacro Cuore, Piacenza, Italy
| | - Alessio Valentini
- Dipartimento per l'Innovazione dei Sistemi Biologici Agroalimentari e Forestali DIBAF, Università della Tuscia, Viterbo, Italy
| | - Nicolò Pietro Paolo Macciotta
- Dipartimento di Agraria, Sezione di Scienze Zootecniche, Università degli Studi di Sassari, V. le Italia, 9, 07100, Sassari, Italy.
| |
Collapse
|
37
|
Sánchez-Molano E, Tsiokos D, Chatziplis D, Jorjani H, Degano L, Diaz C, Rossoni A, Schwarzenbacher H, Seefried F, Varona L, Vicario D, Nicolazzi EL, Banos G. A practical approach to detect ancestral haplotypes in livestock populations. BMC Genet 2016; 17:91. [PMID: 27342071 PMCID: PMC4921009 DOI: 10.1186/s12863-016-0405-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2016] [Accepted: 06/21/2016] [Indexed: 12/04/2022] Open
Abstract
Background The effects of different evolutionary forces are expected to lead to the conservation, over many generations, of particular genomic regions (haplotypes) due to the development of linkage disequilibrium (LD). The detection and identification of early (ancestral) haplotypes can be used to clarify the evolutionary dynamics of different populations as well as identify selection signatures and genomic regions of interest to be used both in conservation and breeding programs. The aims of this study were to develop a simple procedure to identify ancestral haplotypes segregating across several generations both within and between populations with genetic links based on whole-genome scanning. This procedure was tested with simulated and then applied to real data from different genotyped populations of Spanish, Fleckvieh, Simmental and Brown-Swiss cattle. Results The identification of ancestral haplotypes has shown coincident patterns of selection across different breeds, allowing the detection of common regions of interest on different bovine chromosomes and mirroring the evolutionary dynamics of the studied populations. These regions, mainly located on chromosomes BTA5, BTA6, BTA7 and BTA21 are related with certain animal traits such as coat colour and milk protein and fat content. Conclusion In agreement with previous studies, the detection of ancestral haplotypes provides useful information for the development and comparison of breeding and conservation programs both through the identification of selection signatures and other regions of interest, and as indicator of the general genetic status of the populations. Electronic supplementary material The online version of this article (doi:10.1186/s12863-016-0405-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Enrique Sánchez-Molano
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, Scotland, UK.
| | - Dimitrios Tsiokos
- Laboratory of Agrobiotechnology and Inspection of Agricultural Products, Department of Agricultural Technology, School of Agricultural Technology, Food Technology and Nutrition, Alexander Technological Educational Institute of Thessaloniki, Thessaloniki, Greece
| | - Dimitrios Chatziplis
- Laboratory of Agrobiotechnology and Inspection of Agricultural Products, Department of Agricultural Technology, School of Agricultural Technology, Food Technology and Nutrition, Alexander Technological Educational Institute of Thessaloniki, Thessaloniki, Greece
| | | | - Lorenzo Degano
- Associazione Nazionale Allevatori Bovini di razza Pezzata Rossa Italiana, Udine, Italy
| | - Clara Diaz
- Departamento de Mejora Genética Animal, INIA, Madrid, 28040, Spain
| | - Attilio Rossoni
- Associazione Nazionale Allevatori Bovini della Razza Bruna, Verona, Italy
| | | | | | - Luis Varona
- Departamento de Anatomía, Embriología y Genética, Universidad de Zaragoza, Zaragoza, 50013, Spain.,Instituto Agroalimentario de Aragón (IA2), Zaragoza, 50013, Spain
| | - Daniele Vicario
- Associazione Nazionale Allevatori Bovini di razza Pezzata Rossa Italiana, Udine, Italy
| | - Ezequiel L Nicolazzi
- Bioinformatics core facility, Fondazione Parco Tecnologico Padano, Via Einstein, Loc. CascinaCodazza, Lodi, 26900, Italy
| | - Georgios Banos
- The Roslin Institute and Royal (Dick) School of Veterinary Studies, University of Edinburgh, Easter Bush, Midlothian, EH25 9RG, Scotland, UK.,SRUC,The Roslin Institute Building, Easter Bush, Midlothian, EH25 9RG, Edinburgh, UK.,School of Veterinary Medicine, Aristotle University of Thessaloniki, Thessaloniki, Greece
| |
Collapse
|
38
|
Asif AR, Qadri S, Ijaz N, Javed R, Ansari AR, Awais M, Younus M, Riaz H, Du X. Genetic signature of strong recent positive selection at interleukin-32 gene in goat. ASIAN-AUSTRALASIAN JOURNAL OF ANIMAL SCIENCES 2016; 30:912-919. [PMID: 27165029 PMCID: PMC5495668 DOI: 10.5713/ajas.15.0941] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/19/2015] [Revised: 01/06/2016] [Accepted: 03/25/2016] [Indexed: 02/08/2023]
Abstract
OBJECTIVE Identification of the candidate genes that play key roles in phenotypic variations can provide new information about evolution and positive selection. Interleukin (IL)-32 is involved in many biological processes, however, its role for the immune response against various diseases in mammals is poorly understood. Therefore, the current investigation was performed for the better understanding of the molecular evolution and the positive selection of single nucleotide polymorphisms in IL-32 gene. METHODS By using fixation index (FST ) based method, IL-32 (9375) gene was found to be outlier and under significant positive selection with the provisional combined allocation of mean heterozygosity and FST . Using nucleotide sequences of 11 mammalian species from National Center for Biotechnology Information database, the evolutionary selection of IL-32 gene was determined using Maximum likelihood model method, through four models (M1a, M2a, M7, and M8) in Codeml program of phylogenetic analysis by maximum liklihood. RESULTS IL-32 is detected under positive selection using the FST simulations method. The phylogenetic tree revealed that goat IL-32 was in close resemblance with sheep IL-32. The coding nucleotide sequences were compared among 11 species and it was found that the goat IL-32 gene shared identity with sheep (96.54%), bison (91.97%), camel (58.39%), cat (56.59%), buffalo (56.50%), human (56.13%), dog (50.97%), horse (54.04%), and rabbit (53.41%) respectively. CONCLUSION This study provides evidence for IL-32 gene as under significant positive selection in goat.
Collapse
Affiliation(s)
- Akhtar Rasool Asif
- Key Lab of Animal Genetics, Breeding and Reproduction of Ministry Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, China.,The Cooperative Innovation Center for Sustainable Pig Production, Wuhan 430070, China.,Theriogenology Department, College of Veterinary and Animal Science, Jhang, Sub campus of University of Veterinary and Animal Sciences, Lahore 54000, Pakistan
| | - Sumayyah Qadri
- Theriogenology Department, College of Veterinary and Animal Science, Jhang, Sub campus of University of Veterinary and Animal Sciences, Lahore 54000, Pakistan
| | - Nabeel Ijaz
- Key Lab of Animal Genetics, Breeding and Reproduction of Ministry Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
| | - Ruheena Javed
- Key Lab of Animal Genetics, Breeding and Reproduction of Ministry Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, China.,The Cooperative Innovation Center for Sustainable Pig Production, Wuhan 430070, China
| | - Abdur Rahman Ansari
- Theriogenology Department, College of Veterinary and Animal Science, Jhang, Sub campus of University of Veterinary and Animal Sciences, Lahore 54000, Pakistan.,Department of Anatomy, Histology and Embryology, College of Animal and Veterinary Medicine, Huazhong Agricultural University, Wuhan 430070, China
| | - Muhammd Awais
- Key Lab of Animal Genetics, Breeding and Reproduction of Ministry Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, China
| | - Muhammad Younus
- Theriogenology Department, College of Veterinary and Animal Science, Jhang, Sub campus of University of Veterinary and Animal Sciences, Lahore 54000, Pakistan
| | - Hasan Riaz
- Department of Biosciences, COMSATS Institute of Information Technology, Sahiwal 57000, Pakistan
| | - Xiaoyong Du
- Key Lab of Animal Genetics, Breeding and Reproduction of Ministry Education, College of Animal Science and Technology, Huazhong Agricultural University, Wuhan 430070, China.,The Cooperative Innovation Center for Sustainable Pig Production, Wuhan 430070, China.,Hubei Key Laboratory of Agricultural Bioinformatics, College of Informatics, Huazhong Agricultural University, Wuhan China
| |
Collapse
|
39
|
Randhawa IAS, Khatkar MS, Thomson PC, Raadsma HW. A Meta-Assembly of Selection Signatures in Cattle. PLoS One 2016; 11:e0153013. [PMID: 27045296 PMCID: PMC4821596 DOI: 10.1371/journal.pone.0153013] [Citation(s) in RCA: 49] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2015] [Accepted: 03/22/2016] [Indexed: 12/31/2022] Open
Abstract
Since domestication, significant genetic improvement has been achieved for many traits of commercial importance in cattle, including adaptation, appearance and production. In response to such intense selection pressures, the bovine genome has undergone changes at the underlying regions of functional genetic variants, which are termed “selection signatures”. This article reviews 64 recent (2009–2015) investigations testing genomic diversity for departure from neutrality in worldwide cattle populations. In particular, we constructed a meta-assembly of 16,158 selection signatures for individual breeds and their archetype groups (European, African, Zebu and composite) from 56 genome-wide scans representing 70,743 animals of 90 pure and crossbred cattle breeds. Meta-selection-scores (MSS) were computed by combining published results at every given locus, within a sliding window span. MSS were adjusted for common samples across studies and were weighted for significance thresholds across and within studies. Published selection signatures show extensive coverage across the bovine genome, however, the meta-assembly provides a consensus profile of 263 genomic regions of which 141 were unique (113 were breed-specific) and 122 were shared across cattle archetypes. The most prominent peaks of MSS represent regions under selection across multiple populations and harboured genes of known major effects (coat color, polledness and muscle hypertrophy) and genes known to influence polygenic traits (stature, adaptation, feed efficiency, immunity, behaviour, reproduction, beef and dairy production). As the first meta-assembly of selection signatures, it offers novel insights about the hotspots of selective sweeps in the bovine genome, and this method could equally be applied to other species.
Collapse
Affiliation(s)
- Imtiaz A. S. Randhawa
- Reprogen - Animal Bioscience Group, Faculty of Veterinary Science, The University of Sydney, 425 Werombi Road, Camden, 2570, NSW, Australia
- * E-mail:
| | - Mehar S. Khatkar
- Reprogen - Animal Bioscience Group, Faculty of Veterinary Science, The University of Sydney, 425 Werombi Road, Camden, 2570, NSW, Australia
| | - Peter C. Thomson
- Reprogen - Animal Bioscience Group, Faculty of Veterinary Science, The University of Sydney, 425 Werombi Road, Camden, 2570, NSW, Australia
| | - Herman W. Raadsma
- Reprogen - Animal Bioscience Group, Faculty of Veterinary Science, The University of Sydney, 425 Werombi Road, Camden, 2570, NSW, Australia
| |
Collapse
|
40
|
Metzger J, Karwath M, Tonda R, Beltran S, Águeda L, Gut M, Gut IG, Distl O. Runs of homozygosity reveal signatures of positive selection for reproduction traits in breed and non-breed horses. BMC Genomics 2015; 16:764. [PMID: 26452642 PMCID: PMC4600213 DOI: 10.1186/s12864-015-1977-3] [Citation(s) in RCA: 83] [Impact Index Per Article: 9.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2015] [Accepted: 10/03/2015] [Indexed: 11/24/2022] Open
Abstract
Background Modern horses represent heterogeneous populations specifically selected for appearance and performance. Genomic regions under high selective pressure show characteristic runs of homozygosity (ROH) which represent a low genetic diversity. This study aims at detecting the number and functional distribution of ROHs in different horse populations using next generation sequencing data. Methods Next generation sequencing was performed for two Sorraia, one Dülmen Horse, one Arabian, one Saxon-Thuringian Heavy Warmblood, one Thoroughbred and four Hanoverian. After quality control reads were mapped to the reference genome EquCab2.70. ROH detection was performed using PLINK, version 1.07 for a trimmed dataset with 11,325,777 SNPs and a mean read depth of 12. Stretches with homozygous genotypes of >40 kb as well as >400 kb were defined as ROHs. SNPs within consensus ROHs were tested for neutrality. Functional classification was done for genes annotated within ROHs using PANTHER gene list analysis and functional variants were tested for their distribution among breed or non-breed groups. Results ROH detection was performed using whole genome sequences of ten horses of six populations representing various breed types and non-breed horses. In total, an average number of 3492 ROHs were detected in windows of a minimum of 50 consecutive homozygous SNPs and an average number of 292 ROHs in windows of 500 consecutive homozygous SNPs. Functional analyses of private ROHs in each horse revealed a high frequency of genes affecting cellular, metabolic, developmental, immune system and reproduction processes. In non-breed horses, 198 ROHs in 50-SNP windows and seven ROHs in 500-SNP windows showed an enrichment of genes involved in reproduction, embryonic development, energy metabolism, muscle and cardiac development whereas all seven breed horses revealed only three common ROHs in 50-SNP windows harboring the fertility-related gene YES1. In the Hanoverian, a total of 18 private ROHs could be shown to be located in the region of genes potentially involved in neurologic control, signaling, glycogen balance and reproduction. Comparative analysis of homozygous stretches common in all ten horses displayed three ROHs which were all located in the region of KITLG, the ligand of KIT known to be involved in melanogenesis, haematopoiesis and gametogenesis. Conclusions The results of this study give a comprehensive insight into the frequency and number of ROHs in various horses and their potential influence on population diversity and selection pressures. Comparisons of breed and non-breed horses suggest a significant artificial as well as natural selection pressure on reproduction performance in all types of horse populations. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1977-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Julia Metzger
- Institute for Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Bünteweg 17p, 30559, Hannover, Germany.
| | - Matthias Karwath
- Lower Saxony State Office for the Environment, Agriculture and Geology, Unit 74, Animal Breeding and Hygiene, Schlossallee 1, 01468, Moritzburg, Germany.
| | - Raul Tonda
- Centro Nacional de Análisis Genómico, Parc Científic de Barcelona, Torre I Baldiri Reixac, 4, 08028, Barcelona, Spain.
| | - Sergi Beltran
- Centro Nacional de Análisis Genómico, Parc Científic de Barcelona, Torre I Baldiri Reixac, 4, 08028, Barcelona, Spain.
| | - Lídia Águeda
- Centro Nacional de Análisis Genómico, Parc Científic de Barcelona, Torre I Baldiri Reixac, 4, 08028, Barcelona, Spain.
| | - Marta Gut
- Centro Nacional de Análisis Genómico, Parc Científic de Barcelona, Torre I Baldiri Reixac, 4, 08028, Barcelona, Spain.
| | - Ivo Glynne Gut
- Centro Nacional de Análisis Genómico, Parc Científic de Barcelona, Torre I Baldiri Reixac, 4, 08028, Barcelona, Spain.
| | - Ottmar Distl
- Institute for Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Bünteweg 17p, 30559, Hannover, Germany.
| |
Collapse
|
41
|
Jeong H, Song KD, Seo M, Caetano-Anollés K, Kim J, Kwak W, Oh JD, Kim E, Jeong DK, Cho S, Kim H, Lee HK. Exploring evidence of positive selection reveals genetic basis of meat quality traits in Berkshire pigs through whole genome sequencing. BMC Genet 2015; 16:104. [PMID: 26289667 PMCID: PMC4545873 DOI: 10.1186/s12863-015-0265-1] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2015] [Accepted: 08/13/2015] [Indexed: 01/04/2023] Open
Abstract
BACKGROUND Natural and artificial selection following domestication has led to the existence of more than a hundred pig breeds, as well as incredible variation in phenotypic traits. Berkshire pigs are regarded as having superior meat quality compared to other breeds. As the meat production industry seeks selective breeding approaches to improve profitable traits such as meat quality, information about genetic determinants of these traits is in high demand. However, most of the studies have been performed using trained sensory panel analysis without investigating the underlying genetic factors. Here we investigate the relationship between genomic composition and this phenotypic trait by scanning for signatures of positive selection in whole-genome sequencing data. RESULTS We generated genomes of 10 Berkshire pigs at a total of 100.6 coverage depth, using the Illumina Hiseq2000 platform. Along with the genomes of 11 Landrace and 13 Yorkshire pigs, we identified genomic variants of 18.9 million SNVs and 3.4 million Indels in the mapped regions. We identified several associated genes related to lipid metabolism, intramuscular fatty acid deposition, and muscle fiber type which attribute to pork quality (TG, FABP1, AKIRIN2, GLP2R, TGFBR3, JPH3, ICAM2, and ERN1) by applying between population statistical tests (XP-EHH and XP-CLR). A statistical enrichment test was also conducted to detect breed specific genetic variation. In addition, de novo short sequence read assembly strategy identified several candidate genes (SLC25A14, IGF1, PI4KA, CACNA1A) as also contributing to lipid metabolism. CONCLUSIONS Results revealed several candidate genes involved in Berkshire meat quality; most of these genes are involved in lipid metabolism and intramuscular fat deposition. These results can provide a basis for future research on the genomic characteristics of Berkshire pigs.
Collapse
Affiliation(s)
- Hyeonsoo Jeong
- Interdisciplinary Program in Bioinformatics, Seoul National University, Kwan-ak St. 599, Seoul, Kwan-ak Gu, 151-741, Republic of Korea.
| | - Ki-Duk Song
- Department of Animal Biotechnology, Chonbuk National University, Jeonju, 561-756, Republic of Korea.
| | - Minseok Seo
- Interdisciplinary Program in Bioinformatics, Seoul National University, Kwan-ak St. 599, Seoul, Kwan-ak Gu, 151-741, Republic of Korea.
| | | | - Jaemin Kim
- Interdisciplinary Program in Bioinformatics, Seoul National University, Kwan-ak St. 599, Seoul, Kwan-ak Gu, 151-741, Republic of Korea.
| | - Woori Kwak
- Interdisciplinary Program in Bioinformatics, Seoul National University, Kwan-ak St. 599, Seoul, Kwan-ak Gu, 151-741, Republic of Korea.
- C&K genomics, Main Bldg. #514, SNU Research Park, Seoul, 151-919, Republic of Korea.
| | - Jae-Don Oh
- Department of Animal Biotechnology, Chonbuk National University, Jeonju, 561-756, Republic of Korea.
| | - EuiSoo Kim
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA.
| | - Dong Kee Jeong
- Department of Animal Biotechnology, Faculty of Biotechnology, Jeju National University, Ara-1 Dong, Jeju-Do, Jeju, 690-756, Republic of Korea.
| | - Seoae Cho
- C&K genomics, Main Bldg. #514, SNU Research Park, Seoul, 151-919, Republic of Korea.
| | - Heebal Kim
- Interdisciplinary Program in Bioinformatics, Seoul National University, Kwan-ak St. 599, Seoul, Kwan-ak Gu, 151-741, Republic of Korea.
- C&K genomics, Main Bldg. #514, SNU Research Park, Seoul, 151-919, Republic of Korea.
- Department of Agricultural Biotechnology, Seoul National University, Seoul, 151-742, South Korea.
| | - Hak-Kyo Lee
- Department of Animal Biotechnology, Chonbuk National University, Jeonju, 561-756, Republic of Korea.
| |
Collapse
|
42
|
Kasarda R, Moravčíková N, Trakovická A, Mészáros G, Kadlečík O. Genome-wide selection signatures in Pinzgau cattle. POTRAVINARSTVO 2015. [DOI: 10.5219/478] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The aim of this study was to identify the evidence of recent selection based on estimation of the integrated Haplotype Score (iHS), population differentiation index (FST) and characterize affected regions near QTL associated with traits under strong selection in Pinzgau cattle. In total 21 Austrian and 19 Slovak purebreed bulls genotyped with Illumina bovineHD and bovineSNP50 BeadChip were used to identify genomic regions under selection. Only autosomal loci with call rate higher than 90%, minor allele frequency higher than 0.01 and Hardy-Weinberg equlibrium limit of 0.001 were included in the subsequent analyses of selection sweeps presence. The final dataset was consisted from 30538 SNPs with 81.86 kb average adjacent SNPs spacing. The iHS score were averaged into non-overlapping 500 kb segments across the genome. The FST values were also plotted against genome position based on sliding windows approach and averaged over 8 consecutive SNPs. Based on integrated Haplotype Score evaluation only 7 regions with iHS score higher than 1.7 was found. The average iHS score observed for each adjacent syntenic regions indicated slight effect of recent selection in analysed group of Pinzgau bulls. The level of genetic differentiation between Austrian and Slovak bulls estimated based on FST index was low. Only 24% of FST values calculated for each SNP was greather than 0.01. By using sliding windows approach was found that 5% of analysed windows had higher value than 0.01. Our results indicated use of similar selection scheme in breeding programs of Slovak and Austrian Pinzgau bulls. The evidence for genome-wide association between signatures of selection and regions affecting complex traits such as milk production was insignificant, because the loci in segments identified as affected by selection were very distant from each other. Identification of genomic regions that may be under pressure of selection for phenotypic traits to better understanding of the relationship between genotype and phenotype is one of the challenges for livestock genetics.
Collapse
|
43
|
Sorbolini S, Marras G, Gaspa G, Dimauro C, Cellesi M, Valentini A, Macciotta NP. Detection of selection signatures in Piemontese and Marchigiana cattle, two breeds with similar production aptitudes but different selection histories. Genet Sel Evol 2015; 47:52. [PMID: 26100250 PMCID: PMC4476081 DOI: 10.1186/s12711-015-0128-2] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2014] [Accepted: 05/20/2015] [Indexed: 12/21/2022] Open
Abstract
BACKGROUND Domestication and selection are processes that alter the pattern of within- and between-population genetic variability. They can be investigated at the genomic level by tracing the so-called selection signatures. Recently, sequence polymorphisms at the genome-wide level have been investigated in a wide range of animals. A common approach to detect selection signatures is to compare breeds that have been selected for different breeding goals (i.e. dairy and beef cattle). However, genetic variations in different breeds with similar production aptitudes and similar phenotypes can be related to differences in their selection history. METHODS In this study, we investigated selection signatures between two Italian beef cattle breeds, Piemontese and Marchigiana, using genotyping data that was obtained with the Illumina BovineSNP50 BeadChip. The comparison was based on the fixation index (Fst), combined with a locally weighted scatterplot smoothing (LOWESS) regression and a control chart approach. In addition, analyses of Fst were carried out to confirm candidate genes. In particular, data were processed using the varLD method, which compares the regional variation of linkage disequilibrium between populations. RESULTS Genome scans confirmed the presence of selective sweeps in the genomic regions that harbour candidate genes that are known to affect productive traits in cattle such as DGAT1, ABCG2, CAPN3, MSTN and FTO. In addition, several new putative candidate genes (for example ALAS1, ABCB8, ACADS and SOD1) were detected. CONCLUSIONS This study provided evidence on the different selection histories of two cattle breeds and the usefulness of genomic scans to detect selective sweeps even in cattle breeds that are bred for similar production aptitudes.
Collapse
Affiliation(s)
- Silvia Sorbolini
- Dipartimento di Agraria, Sezione di Scienze Zootecniche Università degli Studi di Sassari, 07100, Sassari, Italy.
| | - Gabriele Marras
- Dipartimento di Agraria, Sezione di Scienze Zootecniche Università degli Studi di Sassari, 07100, Sassari, Italy.
| | - Giustino Gaspa
- Dipartimento di Agraria, Sezione di Scienze Zootecniche Università degli Studi di Sassari, 07100, Sassari, Italy.
| | - Corrado Dimauro
- Dipartimento di Agraria, Sezione di Scienze Zootecniche Università degli Studi di Sassari, 07100, Sassari, Italy.
| | - Massimo Cellesi
- Dipartimento di Agraria, Sezione di Scienze Zootecniche Università degli Studi di Sassari, 07100, Sassari, Italy.
| | - Alessio Valentini
- Dipartimento per l'Innovazione dei Sistemi Biologici Agroalimentari e Forestali DIBAF, Università della Tuscia, Viterbo, Italy.
| | - Nicolò Pp Macciotta
- Dipartimento di Agraria, Sezione di Scienze Zootecniche Università degli Studi di Sassari, 07100, Sassari, Italy.
| |
Collapse
|
44
|
Zhao F, McParland S, Kearney F, Du L, Berry DP. Detection of selection signatures in dairy and beef cattle using high-density genomic information. Genet Sel Evol 2015; 47:49. [PMID: 26089079 PMCID: PMC4472243 DOI: 10.1186/s12711-015-0127-3] [Citation(s) in RCA: 120] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2014] [Accepted: 05/19/2015] [Indexed: 01/05/2023] Open
Abstract
Background Artificial selection for economically important traits in cattle is expected to have left distinctive selection signatures on the genome. Access to high-density genotypes facilitates the accurate identification of genomic regions that have undergone positive selection. These findings help to better elucidate the mechanisms of selection and to identify candidate genes of interest to breeding programs. Results Information on 705 243 autosomal single nucleotide polymorphisms (SNPs) in 3122 dairy and beef male animals from seven cattle breeds (Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental) were used to detect selection signatures by applying two complementary methods, integrated haplotype score (iHS) and global fixation index (FST). To control for false positive results, we used false discovery rate (FDR) adjustment to calculate adjusted iHS within each breed and the genome-wide significance level was about 0.003. Using the iHS method, 83, 92, 91, 101, 85, 101 and 86 significant genomic regions were detected for Angus, Belgian Blue, Charolais, Hereford, Holstein-Friesian, Limousin and Simmental cattle, respectively. None of these regions was common to all seven breeds. Using the FST approach, 704 individual SNPs were detected across breeds. Annotation of the regions of the genome that showed selection signatures revealed several interesting candidate genes i.e. DGAT1, ABCG2, MSTN, CAPN3, FABP3, CHCHD7, PLAG1, JAZF1, PRKG2, ACTC1, TBC1D1, GHR, BMP2, TSG1, LYN, KIT and MC1R that play a role in milk production, reproduction, body size, muscle formation or coat color. Fifty-seven common candidate genes were found by both the iHS and global FST methods across the seven breeds. Moreover, many novel genomic regions and genes were detected within the regions that showed selection signatures; for some candidate genes, signatures of positive selection exist in the human genome. Multilevel bioinformatic analyses of the detected candidate genes suggested that the PPAR pathway may have been subjected to positive selection. Conclusions This study provides a high-resolution bovine genomic map of positive selection signatures that are either specific to one breed or common to a subset of the seven breeds analyzed. Our results will contribute to the detection of functional candidate genes that have undergone positive selection in future studies. Electronic supplementary material The online version of this article (doi:10.1186/s12711-015-0127-3) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Fuping Zhao
- National Center for Molecular Genetics and Breeding of Animal, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| | - Sinead McParland
- Animal and Grassland Research and Innovation Centre, Teagasc, Moorpark, Co., Cork, Ireland.
| | - Francis Kearney
- Irish Cattle Breeding Federation, Highfield House, Bandon, Co., Cork, Ireland.
| | - Lixin Du
- National Center for Molecular Genetics and Breeding of Animal, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China.
| | - Donagh P Berry
- Animal and Grassland Research and Innovation Centre, Teagasc, Moorpark, Co., Cork, Ireland.
| |
Collapse
|
45
|
Gutiérrez-Gil B, Arranz JJ, Wiener P. An interpretive review of selective sweep studies in Bos taurus cattle populations: identification of unique and shared selection signals across breeds. Front Genet 2015; 6:167. [PMID: 26029239 PMCID: PMC4429627 DOI: 10.3389/fgene.2015.00167] [Citation(s) in RCA: 57] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/31/2014] [Accepted: 04/13/2015] [Indexed: 12/11/2022] Open
Abstract
This review compiles the results of 21 genomic studies of European Bos taurus breeds and thus provides a general picture of the selection signatures in taurine cattle identified by genome-wide selection-mapping scans. By performing a comprehensive summary of the results reported in the literature, we compiled a list of 1049 selection sweeps described across 37 cattle breeds (17 beef breeds, 14 dairy breeds, and 6 dual-purpose breeds), and four different beef-vs.-dairy comparisons, which we subsequently grouped into core selective sweep (CSS) regions, defined as consecutive signals within 1 Mb of each other. We defined a total of 409 CSSs across the 29 bovine autosomes, 232 (57%) of which were associated with a single-breed (Single-breed CSSs), 134 CSSs (33%) were associated with a limited number of breeds (Two-to-Four-breed CSSs) and 39 CSSs (9%) were associated with five or more breeds (Multi-breed CSSs). For each CSS, we performed a candidate gene survey that identified 291 genes within the CSS intervals (from the total list of 5183 BioMart-extracted genes) linked to dairy and meat production, stature, and coat color traits. A complementary functional enrichment analysis of the CSS positional candidates highlighted other genes related to pathways underlying behavior, immune response, and reproductive traits. The Single-breed CSSs revealed an over-representation of genes related to dairy and beef production, this was further supported by over-representation of production-related pathway terms in these regions based on a functional enrichment analysis. Overall, this review provides a comparative map of the selection sweeps reported in European cattle breeds and presents for the first time a characterization of the selection sweeps that are found in individual breeds. Based on their uniqueness, these breed-specific signals could be considered as “divergence signals,” which may be useful in characterizing and protecting livestock genetic diversity.
Collapse
Affiliation(s)
| | - Juan J Arranz
- Departamento de Producción Animal, Universidad de León León, Spain
| | - Pamela Wiener
- Division of Genetics and Genomics, Roslin Institute and R(D)SVS, University of Edinburgh Midlothian, UK
| |
Collapse
|
46
|
Kim ES, Sonstegard TS, Rothschild MF. Recent artificial selection in U.S. Jersey cattle impacts autozygosity levels of specific genomic regions. BMC Genomics 2015; 16:302. [PMID: 25887761 PMCID: PMC4409734 DOI: 10.1186/s12864-015-1500-x] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2014] [Accepted: 03/30/2015] [Indexed: 01/01/2023] Open
Abstract
BACKGROUND Genome signatures of artificial selection in U.S. Jersey cattle were identified by examining changes in haplotype homozygosity for a resource population of animals born between 1953 and 2007. Genetic merit of this population changed dramatically during this period for a number of traits, especially milk yield. The intense selection underlying these changes was achieved through extensive use of artificial insemination (AI), which also increased consanguinity of the population to a few superior Jersey bulls. As a result, allele frequencies are shifted for many contemporary animals, and in numerous cases to a homozygous state for specific genomic regions. The goal of this study was to identify those selection signatures that occurred after extensive use of AI since the 1960, using analyses of shared haplotype segments or Runs of Homozygosity. When combined with animal birth year information, signatures of selection associated with economically important traits were identified and compared to results from an extended haplotype homozygosity analysis. RESULTS Overall, our results reveal that more recent selection increased autozygosity across the entire genome, but some specific regions increased more than others. A genome-wide scan identified more than 15 regions with a substantial change in autozygosity. Haplotypes found to be associated with increased milk, fat and protein yield in U.S. Jersey cattle also consistently increased in frequency. CONCLUSIONS The analyses used in this study was able to detect directional selection over the last few decades when individual production records for Jersey animals were available.
Collapse
Affiliation(s)
- Eui-Soo Kim
- United States Department of Agriculture, Animal Genomics & Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, Beltsville, MD, 20705, USA.
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA.
| | - Tad S Sonstegard
- United States Department of Agriculture, Animal Genomics & Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, Beltsville, MD, 20705, USA.
| | - Max F Rothschild
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA.
| |
Collapse
|
47
|
Bomba L, Nicolazzi EL, Milanesi M, Negrini R, Mancini G, Biscarini F, Stella A, Valentini A, Ajmone-Marsan P. Relative extended haplotype homozygosity signals across breeds reveal dairy and beef specific signatures of selection. Genet Sel Evol 2015; 47:25. [PMID: 25888030 PMCID: PMC4383072 DOI: 10.1186/s12711-015-0113-9] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2014] [Accepted: 03/19/2015] [Indexed: 02/02/2023] Open
Abstract
BACKGROUND A number of methods are available to scan a genome for selection signatures by evaluating patterns of diversity within and between breeds. Among these, "extended haplotype homozygosity" (EHH) is a reliable approach to detect genome regions under recent selective pressure. The objective of this study was to use this approach to identify regions that are under recent positive selection and shared by the most representative Italian dairy and beef cattle breeds. RESULTS A total of 3220 animals from Italian Holstein (2179), Italian Brown (775), Simmental (493), Marchigiana (485) and Piedmontese (379) breeds were genotyped with the Illumina BovineSNP50 BeadChip v.1. After standard quality control procedures, genotypes were phased and core haplotypes were identified. The decay of linkage disequilibrium (LD) for each core haplotype was assessed by measuring the EHH. Since accurate estimates of local recombination rates were not available, relative EHH (rEHH) was calculated for each core haplotype. Genomic regions that carry frequent core haplotypes and with significant rEHH values were considered as candidates for recent positive selection. Candidate regions were aligned across to identify signals shared by dairy or beef cattle breeds. Overall, 82 and 87 common regions were detected among dairy and beef cattle breeds, respectively. Bioinformatic analysis identified 244 and 232 genes in these common genomic regions. Gene annotation and pathway analysis showed that these genes are involved in molecular functions that are biologically related to milk or meat production. CONCLUSIONS Our results suggest that a multi-breed approach can lead to the identification of genomic signatures in breeds of cattle that are selected for the same production goal and thus to the localisation of genomic regions of interest in dairy and beef production.
Collapse
Affiliation(s)
- Lorenzo Bomba
- Istituto di Zootecnica, UCSC, via Emilia Parmense 84, Piacenza, 29122, Italy.
| | - Ezequiel L Nicolazzi
- Fondazione Parco Tecnologico Padano, Via Einstein, Loc. Cascina Codazza, Lodi, 26900, Italy.
| | - Marco Milanesi
- Istituto di Zootecnica, UCSC, via Emilia Parmense 84, Piacenza, 29122, Italy.
| | - Riccardo Negrini
- Associazione Italiana Allevatori (AIA), Via Tomassetti 9, Rome, 00161, Italy.
| | - Giordano Mancini
- Center for Computational Chemistry and Cosmology, Scuola Normale Superiore, Via Consoli del Mare 2, Pisa, 56126, Italy.
| | - Filippo Biscarini
- Fondazione Parco Tecnologico Padano, Via Einstein, Loc. Cascina Codazza, Lodi, 26900, Italy.
| | - Alessandra Stella
- Fondazione Parco Tecnologico Padano, Via Einstein, Loc. Cascina Codazza, Lodi, 26900, Italy. .,Istituto di biologia e biotecnologia Agraria (IBBA-CNR), Consiglio Nazionale delle Ricerche, Via Einstein, Cascina Codazza, Lodi, 26900, Italy.
| | - Alessio Valentini
- Dipartimento per l'Innovazione nei Sistemi Biologici, Agroalimentari e Forestali (DIBAF), via de Lellis, Viterbo, 01100, Italy.
| | - Paolo Ajmone-Marsan
- Istituto di Zootecnica, UCSC, via Emilia Parmense 84, Piacenza, 29122, Italy.
| |
Collapse
|
48
|
Howard JT, Maltecca C, Haile-Mariam M, Hayes BJ, Pryce JE. Characterizing homozygosity across United States, New Zealand and Australian Jersey cow and bull populations. BMC Genomics 2015; 16:187. [PMID: 25879195 PMCID: PMC4460752 DOI: 10.1186/s12864-015-1352-4] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2014] [Accepted: 02/18/2015] [Indexed: 11/19/2022] Open
Abstract
Background Dairy cattle breeding objectives are in general similar across countries, but environment and management conditions may vary, giving rise to slightly different selection pressures applied to a given trait. This potentially leads to different selection pressures to loci across the genome that, if large enough, may give rise to differential regions with high levels of homozygosity. The objective of this study was to characterize differences and similarities in the location and frequency of homozygosity related measures of Jersey dairy cows and bulls from the United States (US), Australia (AU) and New Zealand (NZ). Results The populations consisted of a subset of genotyped Jersey cows born in US (n = 1047) and AU (n = 886) and Jersey bulls progeny tested from the US (n = 736), AU (n = 306) and NZ (n = 768). Differences and similarities across populations were characterized using a principal component analysis (PCA) and a run of homozygosity (ROH) statistic (ROH45), which counts the frequency of a single nucleotide polymorphism (SNP) being in a ROH of at least 45 SNP. Regions that exhibited high frequencies of ROH45 and those that had significantly different ROH45 frequencies between populations were investigated for their association with milk yield traits. Within sex, the PCA revealed slight differentiation between the populations, with the greatest occurring between the US and NZ bulls. Regions with high levels of ROH45 for all populations were detected on BTA3 and BTA7 while several other regions differed in ROH45 frequency across populations, the largest number occurring for the US and NZ bull contrast. In addition, multiple regions with different ROH45 frequencies across populations were found to be associated with milk yield traits. Conclusion Multiple regions exhibited differential ROH45 across AU, NZ and US cow and bull populations, an interpretation is that locations of the genome are undergoing differential directional selection. Two regions on BTA3 and BTA7 had high ROH45 frequencies across all populations and will be investigated further to determine the gene(s) undergoing directional selection. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1352-4) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Jeremy T Howard
- Department of Animal Science, North Carolina State University, Raleigh, NC, 27695-7627, USA.
| | - Christian Maltecca
- Department of Animal Science, North Carolina State University, Raleigh, NC, 27695-7627, USA.
| | - Mekonnen Haile-Mariam
- Dairy Futures Cooperative Research Centre, 5 Ring Road, Bundoora, Victoria, 3083, Australia. .,Biosciences Research Division, Department of Environment and Primary Industries Victoria, 5 Ring Road, Bundoora, 3083, Australia.
| | - Ben J Hayes
- Dairy Futures Cooperative Research Centre, 5 Ring Road, Bundoora, Victoria, 3083, Australia. .,La Trobe University, Bundoora, Victoria, 3086, Australia. .,Biosciences Research Division, Department of Environment and Primary Industries Victoria, 5 Ring Road, Bundoora, 3083, Australia.
| | - Jennie E Pryce
- Dairy Futures Cooperative Research Centre, 5 Ring Road, Bundoora, Victoria, 3083, Australia. .,La Trobe University, Bundoora, Victoria, 3086, Australia. .,Biosciences Research Division, Department of Environment and Primary Industries Victoria, 5 Ring Road, Bundoora, 3083, Australia.
| |
Collapse
|
49
|
Ma Y, Wei J, Zhang Q, Chen L, Wang J, Liu J, Ding X. A genome scan for selection signatures in pigs. PLoS One 2015; 10:e0116850. [PMID: 25756180 PMCID: PMC4355907 DOI: 10.1371/journal.pone.0116850] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2014] [Accepted: 12/15/2014] [Indexed: 11/24/2022] Open
Abstract
Identifying signatures of selection can provide a straightforward insight into the mechanism of artificial selection and further uncover the causal genes related to the phenotypic variation. Based on Illumina Porcine60KSNP chip data, four complementary methods, Long-Range Haplotype (LRH), Tajima’s D, Cross Population Extend Haplotype Homozygosity Test (XPEHH) and FST, were implemented in this study to detect the selection signatures in the whole genome of one typical Chinese indigenous breed, Rongchang, one Chinese cultivated breed, Songliao, and two western breeds, Landrace and Yorkshire. False Discovery Rate (FDR) was implemented to control the false positive rates. In our study, a total of 159, 127, 179 and 159 candidate selection regions with average length of 0.80 Mb, 0.73 Mb, 0.78 Mb and 0.73 Mb were identified in Landrace, Rongchang, Songliao and Yorkshire, respectively, that span approximately 128.00 Mb, 92.38 Mb, 130.30 Mb and 115.40 Mb and account for approximately 3.74–5.33% of genome across all autosomes. The selection regions of 11.52 Mb shared by Landrace and Yorkshire were the longest when chosen pairs from the pool of the four breeds were examined. The overlaps between Yorkshire and Songliao, approximately 9.20 Mb, were greater than those of Yorkshire and Rongchang. Meanwhile, the overlaps between Landrace and Songliao were greater than those of Landrace and Rongchang but less than those of Songliao and Ronchang. Bioinformatics analysis showed that the genes/QTLs relevant to fertility, coat color, and ear morphology were found in candidate selection regions. Some genes, such as LEMD3, MC1R, KIT, TRHR etc. that were reported under selection, were confirmed in our study, and this analysis also demonstrated the diversity of breeds.
Collapse
Affiliation(s)
- Yunlong Ma
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, National Engineering Laboratory for Animal Breeding, College Animal Science and Technology, China Agricultural University, Beijing, P.R. China
| | - Julong Wei
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, National Engineering Laboratory for Animal Breeding, College Animal Science and Technology, China Agricultural University, Beijing, P.R. China
| | - Qin Zhang
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, National Engineering Laboratory for Animal Breeding, College Animal Science and Technology, China Agricultural University, Beijing, P.R. China
| | - Lei Chen
- Chongqing Academy of Animal Science, Chongqing, P.R. China
| | - Jinyong Wang
- Chongqing Academy of Animal Science, Chongqing, P.R. China
| | - Jianfeng Liu
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, National Engineering Laboratory for Animal Breeding, College Animal Science and Technology, China Agricultural University, Beijing, P.R. China
- * E-mail: (JFL); (XDD)
| | - Xiangdong Ding
- Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture, National Engineering Laboratory for Animal Breeding, College Animal Science and Technology, China Agricultural University, Beijing, P.R. China
- * E-mail: (JFL); (XDD)
| |
Collapse
|
50
|
García-Ruiz A, Ruiz-López FDJ, Van Tassell CP, Montaldo HH, Huson HJ. Genetic differentiation of Mexican Holstein cattle and its relationship with Canadian and U.S. Holsteins. Front Genet 2015; 6:7. [PMID: 25709615 PMCID: PMC4321612 DOI: 10.3389/fgene.2015.00007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2014] [Accepted: 01/09/2015] [Indexed: 11/22/2022] Open
Abstract
The Mexican Holstein (HO) industry has imported Canadian and US (CAN + USA) HO germplasm for use in two different production systems, the conventional (Conv) and the low income (Lowi) system. The objective of this work was to study the genetic composition and differentiation of the Mexican HO cattle, considering the production system in which they perform and their relationship with the Canadian and US HO populations. The analysis included information from 149, 303, and 173 unrelated or with unknown pedigree HO animals from the Conv, Lowi, and CAN + USA populations, respectively. Canadian and US Jersey (JE) and Brown Swiss (BS) genotypes (162 and 86, respectively) were used to determine if Mexican HOs were hybridized with either of these breeds. After quality control filtering, a total of 6,617 out of 6,836 single nucleotide polymorphism markers were used. To describe the genetic diversity across the populations, principal component (PC), admixture composition, and linkage disequilibrium (LD; r2) analyses were performed. Through the PC analysis, HO × JE and HO × BS crossbreeding was detected in the Lowi system. The Conv system appeared to be in between Lowi and CAN + USA populations. Admixture analysis differentiated between the genetic composition of the Conv and Lowi systems, and five ancestry groups associated to sire’s country of origin were identified. The minimum distance between markers to estimate a useful LD was found to be 54.5 kb for the Mexican HO populations. At this average distance, the persistence of phase across autosomes of Conv and Lowi systems was 0.94, for Conv and CAN + USA was 0.92 and for the Lowi and CAN + USA was 0.91. Results supported the flow of germplasm among populations being Conv a source for Lowi, and dependent on migration from CAN + USA. Mexican HO cattle in Conv and Lowi populations share common ancestry with CAN + USA but have different genetic signatures.
Collapse
Affiliation(s)
- Adriana García-Ruiz
- Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México Ajuchitlán, Mexico
| | - Felipe de J Ruiz-López
- Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México Ajuchitlán, Mexico ; Centro Nacional de Investigación Disciplinaria en Fisiología y Mejoramiento Animal, Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias - Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación Ajuchitlán, Mexico
| | - Curtis P Van Tassell
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, United State Department of Agriculture Beltsville, MD, USA
| | - Hugo H Montaldo
- Facultad de Estudios Superiores Cuautitlán, Universidad Nacional Autónoma de México Ajuchitlán, Mexico
| | - Heather J Huson
- Animal Genomics and Improvement Laboratory, Agricultural Research Service, United State Department of Agriculture Beltsville, MD, USA ; Department of Animal Science, Cornell University Ithaca, NY, USA
| |
Collapse
|