1
|
Qi L, Xiao L, Fu R, Nie Q, Zhang X, Luo W. Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens. Poult Sci 2024; 103:103863. [PMID: 38810566 PMCID: PMC11166977 DOI: 10.1016/j.psj.2024.103863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2024] [Revised: 04/25/2024] [Accepted: 05/13/2024] [Indexed: 05/31/2024] Open
Abstract
The introduction of exotic breeds and the cultivation of new lines by breeding companies have posed challenges to native chickens in South China, including loss of breed characteristics, decreased genetic diversity, and declining purity. Understanding the population genetic structure and genetic diversity of native chickens in South China is crucial for further advancements in breeding efforts. In this study, we analyzed the population genetic structure and genetic diversity of 321 individuals from 10 different breeds in South China. By comparing commercial chickens with native ones, we identified selection signatures occurring between local chickens and commercial breeds. The analysis of population genetic structure revealed that the native chicken populations in South China exhibited a considerable level of genetic diversity. Moreover, the commercial lines of Xiaobai chicken and Huangma chicken displayed even higher levels of genetic diversity, which distinguished them from other native varieties at the clustering level. However, certain individuals within these commercial varieties showed a discernible genetic relationship with the native populations. Notably, both commercial varieties also retained a significant degree of genetic similarity to their respective native counterparts. In order to investigate the genomic changes occurring during the commercialization of native chickens, we employed 4 methods (Fst, ROD, XPCLR, and XPEHH) to identify potential candidate regions displaying selective signatures in Southern Chinese native chicken population. A total of 168 (identified by Fst and ROD) and 86 (identified by XPCLR and XPEHH) overlapping genes were discovered. Functional annotation analysis revealed that these genes may be associated with reproduction and growth (SAMSN1, HYLS1, ROBO3, FGF14, PRSS23), musculoskeletal development (DNER, MYBPC1, DGKB, ORC1, KLF10), disease resistance and environmental adaptability (PUS3, CRB2, CALD1, USP15, SGCD, LTBP1), as well as egg production (ADGRB3, ACSF3). Overall, native chickens in South China harbor numerous selective sweep regions compared to commercial chickens, enriching valuable genomic resources for future genetic research and breeding conservation.
Collapse
Affiliation(s)
- Lin Qi
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou Guangzhou 510642, China
| | - Liangchao Xiao
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou Guangzhou 510642, China
| | - Rong Fu
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou Guangzhou 510642, China
| | - Qinghua Nie
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou Guangzhou 510642, China
| | - Xiquan Zhang
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou Guangzhou 510642, China
| | - Wen Luo
- State Key Laboratory of Livestock and Poultry Breeding, & Lingnan Guangdong Laboratory of Agriculture, South China Agricultural University, Guangzhou 510642, China; Guangdong Provincial Key Lab of Agro-Animal Genomics and Molecular Breeding, and Key Lab of Chicken Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affair, South China Agricultural University, 510642, China; Department of Animal Genetics, Breeding and Reproduction, College of Animal Science, South China Agricultural University, Guangzhou Guangzhou 510642, China.
| |
Collapse
|
2
|
Lukic B, Curik I, Drzaic I, Galić V, Shihabi M, Vostry L, Cubric-Curik V. Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds. J Anim Sci Biotechnol 2023; 14:142. [PMID: 37932811 PMCID: PMC10626677 DOI: 10.1186/s40104-023-00936-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 09/04/2023] [Indexed: 11/08/2023] Open
Abstract
BACKGROUND The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type. RESULTS We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis. CONCLUSIONS Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.
Collapse
Affiliation(s)
- Boris Lukic
- Faculty of Agrobiotechnical Sciences Osijek, J.J, Strossmayer University of Osijek, Vladimira Preloga 1, 31000, Osijek, Croatia.
| | - Ino Curik
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia.
| | - Ivana Drzaic
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia
| | - Vlatko Galić
- Department of Maize Breeding and Genetics, Agricultural Institute Osijek, Južno predgrađe 17, 31000, Osijek, Croatia
| | - Mario Shihabi
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia
| | - Luboš Vostry
- Czech University of Life Sciences Prague, Kamýcká 129, 165 00, Praque, Czech Republic
| | - Vlatka Cubric-Curik
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia
| |
Collapse
|
3
|
Tanaka T, Hayakawa T, Teshima KM. Power of neutrality tests for detecting natural selection. G3 (BETHESDA, MD.) 2023; 13:jkad161. [PMID: 37481468 PMCID: PMC10542275 DOI: 10.1093/g3journal/jkad161] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/09/2023] [Revised: 06/09/2023] [Accepted: 07/19/2023] [Indexed: 07/24/2023]
Abstract
Detection of natural selection is one of the main interests in population genetics. Thus, many tests have been developed for detecting natural selection using genomic data. Although it is recognized that the utility of tests depends on several evolutionary factors, such as the timing of selection, strength of selection, frequency of selected alleles, demographic events, and initial frequency of selected allele when selection started acting (softness of selection), the relationships between such evolutionary factors and the power of tests are not yet entirely clear. In this study, we investigated the power of 4 tests: Tajiama's D, Fay and Wu's H, relative extended haplotype homozygosity (rEHH), and integrated haplotype score (iHS), under ranges of evolutionary parameters and demographic models to quantitatively expand the understanding of approaches for detecting selection. The results show that each test detects selection within a limited parameter range, and there are still wide ranges of parameters for which none of these tests work effectively. In addition, the parameter space in which each test shows the highest power overlaps the empirical results of previous research. These results indicate that our present perspective of adaptation is limited to only a part of actual adaptation.
Collapse
Affiliation(s)
- Tomotaka Tanaka
- Graduate School of System Life Science, Kyushu University, Fukuoka 819-0395, Japan
| | - Toshiyuki Hayakawa
- Graduate School of System Life Science, Kyushu University, Fukuoka 819-0395, Japan
- Faculty of Arts and Science, Kyushu University, Fukuoka 819-0395, Japan
| | - Kosuke M Teshima
- Department of Biology, Faculty of Science, Kyushu University, Fukuoka 819-0395, Japan
| |
Collapse
|
4
|
Dixit SP, Bhatia AK, Ganguly I, Singh S, Dash S, Sharma A, Anandkumar N, Dang AK, Jayakumar S. Genome analyses revealed genetic admixture and selection signatures in Bos indicus. Sci Rep 2021; 11:21924. [PMID: 34753978 PMCID: PMC8578574 DOI: 10.1038/s41598-021-01144-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2021] [Accepted: 10/21/2021] [Indexed: 11/20/2022] Open
Abstract
The genomic diversity and relationship among seven diverse cattle breeds viz. Sahiwal, Tharparkar, Gir, Vechur, Ongole, Kangayam and Hariana were investigated in 132 random samples based on high density SNP array comprising > 777 K SNPs. A total of 1993 SNPs (0.25% of the total) having greater power (FST ≥ 0.20) to differentiate these cattle populations were identified, and utilized to partition genome of each animal into a predefined number of clusters. The structure of these cattle indicated shared ancestry of dairy breeds viz. Gir, Tharparkar and Sahiwal. Most of the animals (> 76%) of different populations under study except Vechur clustered into their own group of animals called breed. Vechur population retained highest rate of admixture, consistent with its crossing with other breeds. Ongole, Kangayam and Hariana shared comparatively less of their genome (≤ 15%) with other breeds. The study indicated that all seven breeds evolved from their independent ancestry but there was intermixing of these breeds in the recent past. The selection signatures identified between draft (Kangayam) and dairy breeds included several genes like FAM19A2, RAB31P, BEST3, DGKA, AHCY, PIGU and PFKP which are involved in immune response, metabolic pathway, transportation of glucose and sugars, signaling pathways, cellular processes, cell division and glycolysis regulation, respectively. Moreover, these genomic regions also harbour QTLs affecting milk performance traits. The signatures were also identified even between the dairy breeds. In comparison to large-sized cattle, there were significant differences in the number of QTLs affecting production (body weight, growth rate etc.) and morphological traits (height) in short-statured Vechur breed. The presence of HMGA2 gene in the selection signature on chromosome 5 may explain the variations in stature between these cattle.
Collapse
Affiliation(s)
- S P Dixit
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India.
| | - A K Bhatia
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Indrajit Ganguly
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Sanjeev Singh
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Soumya Dash
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - Anurodh Sharma
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - N Anandkumar
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| | - A K Dang
- ICAR - National Dairy Research Institute, Karnal, Haryana, 132001, India
| | - S Jayakumar
- ICAR - National Bureau of Animal Genetic Resources, Karnal, Haryana, 132001, India
| |
Collapse
|
5
|
Rostamzadeh Mahdabi E, Esmailizadeh A, Ayatollahi Mehrgardi A, Asadi Fozi M. A genome-wide scan to identify signatures of selection in two Iranian indigenous chicken ecotypes. Genet Sel Evol 2021; 53:72. [PMID: 34503452 PMCID: PMC8428137 DOI: 10.1186/s12711-021-00664-9] [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: 01/11/2021] [Accepted: 08/25/2021] [Indexed: 11/10/2022] Open
Abstract
Background Various regions of the chicken genome have been under natural and artificial selection for thousands of years. The substantial diversity that exits among chickens from different geographic regions provides an excellent opportunity to investigate the genomic regions under selection which, in turn, will increase our knowledge about the mechanisms that underlie chicken diversity and adaptation. Several statistics have been developed to detect genomic regions that are under selection. In this study, we applied approaches based on differences in allele or haplotype frequencies (FST and hapFLK, respectively) between populations, differences in long stretches of consecutive homozygous sequences (ROH), and differences in allele frequencies within populations (composite likelihood ratio (CLR)) to identify inter- and intra-populations traces of selection in two Iranian indigenous chicken ecotypes, the Lari fighting chicken and the Khazak or creeper (short-leg) chicken. Results Using whole-genome resequencing data of 32 individuals from the two chicken ecotypes, approximately 11.9 million single nucleotide polymorphisms (SNPs) were detected and used in genomic analyses after quality processing. Examination of the distribution of ROH in the two populations indicated short to long ROH, ranging from 0.3 to 5.4 Mb. We found 90 genes that were detected by at least two of the four applied methods. Gene annotation of the detected putative regions under selection revealed candidate genes associated with growth (DCN, MEOX2 and CACNB1), reproduction (ESR1 and CALCR), disease resistance (S1PR1, ALPK1 and MHC-B), behavior pattern (AGMO, GNAO1 and PSEN1), and morphological traits (IHH and NHEJ1). Conclusions Our findings show that these two phenotypically different indigenous chicken populations have been under selection for reproduction, immune, behavioral, and morphology traits. The results illustrate that selection can play an important role in shaping signatures of differentiation across the genomic landscape of two chicken populations. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00664-9.
Collapse
Affiliation(s)
- Elaheh Rostamzadeh Mahdabi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran
| | - Ali Esmailizadeh
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran
| | - Ahmad Ayatollahi Mehrgardi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran
| | - Masood Asadi Fozi
- Department of Animal Science, Faculty of Agriculture, Shahid Bahonar University of Kerman, 22 Bahman Blvd, Kerman, Iran.
| |
Collapse
|
6
|
Casto-Rebollo C, Argente MJ, García ML, Blasco A, Ibáñez-Escriche N. Selection for environmental variance of litter size in rabbits involves genes in pathways controlling animal resilience. Genet Sel Evol 2021; 53:59. [PMID: 34256696 PMCID: PMC8276493 DOI: 10.1186/s12711-021-00653-y] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2020] [Accepted: 07/06/2021] [Indexed: 11/16/2022] Open
Abstract
Background Environmental variance (VE) is partially under genetic control, which means that the VE of individuals that share the same environment can differ because they have different genotypes. Previously, a divergent selection experiment for VE of litter size (LS) during 13 generations in rabbit yielded a successful response and revealed differences in resilience between the divergent lines. The aim of the current study was to identify signatures of selection in these divergent lines to better understand the molecular mechanisms and pathways that control VE of LS and animal resilience. Three methods (FST, ROH and varLD) were used to identify signatures of selection in a set of 473 genotypes from these rabbit lines (377) and a base population (96). A whole-genome sequencing (WGS) analysis was performed on 54 animals to detect genes with functional mutations. Results By combining signatures of selection and WGS data, we detected 373 genes with functional mutations in their transcription units, among which 111 had functions related to the immune system, stress response, reproduction and embryo development, and/or carbohydrate and lipid metabolism. The genes TTC23L, FBXL20, GHDC, ENSOCUG00000031631, SLC18A1, CD300LG, MC2R, and ENSOCUG00000006264 were particularly relevant, since each one carried a functional mutation that was fixed in one of the rabbit lines and absent in the other line. In the 3ʹUTR region of the MC2R and ENSOCUG00000006264 genes, we detected a novel insertion/deletion (INDEL) variant. Conclusions Our findings provide further evidence in favour of VE as a measure of animal resilience. Signatures of selection were identified for VE of LS in genes that have a functional mutation in their transcription units and are mostly implicated in the immune response and stress response pathways. However, the real implications of these genes for VE and animal resilience will need to be assessed through functional analyses. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00653-y.
Collapse
Affiliation(s)
- Cristina Casto-Rebollo
- Institute for Animal Science and Technology, Universitat Politècnica de València, Valencia, Spain
| | - María José Argente
- Departamento de Tecnología Agroalimentaria, Universidad Miguel Hernández de Elche, Orihuela, Spain
| | - María Luz García
- Departamento de Tecnología Agroalimentaria, Universidad Miguel Hernández de Elche, Orihuela, Spain
| | - Agustín Blasco
- Institute for Animal Science and Technology, Universitat Politècnica de València, Valencia, Spain
| | - Noelia Ibáñez-Escriche
- Institute for Animal Science and Technology, Universitat Politècnica de València, Valencia, Spain.
| |
Collapse
|
7
|
Garud NR, Messer PW, Petrov DA. Detection of hard and soft selective sweeps from Drosophila melanogaster population genomic data. PLoS Genet 2021; 17:e1009373. [PMID: 33635910 PMCID: PMC7946363 DOI: 10.1371/journal.pgen.1009373] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2020] [Revised: 03/10/2021] [Accepted: 01/17/2021] [Indexed: 12/12/2022] Open
Abstract
Whether hard sweeps or soft sweeps dominate adaptation has been a matter of much debate. Recently, we developed haplotype homozygosity statistics that (i) can detect both hard and soft sweeps with similar power and (ii) can classify the detected sweeps as hard or soft. The application of our method to population genomic data from a natural population of Drosophila melanogaster (DGRP) allowed us to rediscover three known cases of adaptation at the loci Ace, Cyp6g1, and CHKov1 known to be driven by soft sweeps, and detected additional candidate loci for recent and strong sweeps. Surprisingly, all of the top 50 candidates showed patterns much more consistent with soft rather than hard sweeps. Recently, Harris et al. 2018 criticized this work, suggesting that all the candidate loci detected by our haplotype statistics, including the positive controls, are unlikely to be sweeps at all and that instead these haplotype patterns can be more easily explained by complex neutral demographic models. They also claim that these neutral non-sweeps are likely to be hard instead of soft sweeps. Here, we reanalyze the DGRP data using a range of complex admixture demographic models and reconfirm our original published results suggesting that the majority of recent and strong sweeps in D. melanogaster are first likely to be true sweeps, and second, that they do appear to be soft. Furthermore, we discuss ways to take this work forward given that most demographic models employed in such analyses are necessarily too simple to capture the full demographic complexity, while more realistic models are unlikely to be inferred correctly because they require a large number of free parameters.
Collapse
Affiliation(s)
- Nandita R. Garud
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, California, United States of America
- Department of Human Genetics, University of California, Los Angeles, California, United States of America
| | - Philipp W. Messer
- Department of Computational Biology, Cornell University, Ithaca, New York, United States of America
| | - Dmitri A. Petrov
- Department of Biology, Stanford University, Stanford, California, United States of America
| |
Collapse
|
8
|
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
|
9
|
Peripolli E, Reimer C, Ha NT, Geibel J, Machado MA, Panetto JCDC, do Egito AA, Baldi F, Simianer H, da Silva MVGB. Genome-wide detection of signatures of selection in indicine and Brazilian locally adapted taurine cattle breeds using whole-genome re-sequencing data. BMC Genomics 2020; 21:624. [PMID: 32917133 PMCID: PMC7488563 DOI: 10.1186/s12864-020-07035-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2020] [Accepted: 08/27/2020] [Indexed: 12/30/2022] Open
Abstract
BACKGROUND The cattle introduced by European conquerors during the Brazilian colonization period were exposed to a process of natural selection in different types of biomes throughout the country, leading to the development of locally adapted cattle breeds. In this study, whole-genome re-sequencing data from indicine and Brazilian locally adapted taurine cattle breeds were used to detect genomic regions under selective pressure. Within-population and cross-population statistics were combined separately in a single score using the de-correlated composite of multiple signals (DCMS) method. Putative sweep regions were revealed by assessing the top 1% of the empirical distribution generated by the DCMS statistics. RESULTS A total of 33,328,447 biallelic SNPs with an average read depth of 12.4X passed the hard filtering process and were used to access putative sweep regions. Admixture has occurred in some locally adapted taurine populations due to the introgression of exotic breeds. The genomic inbreeding coefficient based on runs of homozygosity (ROH) concurred with the populations' historical background. Signatures of selection retrieved from the DCMS statistics provided a comprehensive set of putative candidate genes and revealed QTLs disclosing cattle production traits and adaptation to the challenging environments. Additionally, several candidate regions overlapped with previous regions under selection described in the literature for other cattle breeds. CONCLUSION The current study reported putative sweep regions that can provide important insights to better understand the selective forces shaping the genome of the indicine and Brazilian locally adapted taurine cattle breeds. Such regions likely harbor traces of natural selection pressures by which these populations have been exposed and may elucidate footprints for adaptation to the challenging climatic conditions.
Collapse
Affiliation(s)
- Elisa Peripolli
- São Paulo State University (Unesp), School of Agricultural and Veterinarian Sciences, Jaboticabal, 14884-900, Brazil
| | - Christian Reimer
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
| | - Ngoc-Thuy Ha
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
| | - Johannes Geibel
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
| | - Marco Antonio Machado
- National Council for Scientific and Technological Development (CNPq), Lago Sul, 71605-001, Brazil
- Embrapa Dairy Cattle, Juiz de Fora, 36038-330, Brazil
| | | | | | - Fernando Baldi
- São Paulo State University (Unesp), School of Agricultural and Veterinarian Sciences, Jaboticabal, 14884-900, Brazil
| | - Henner Simianer
- Animal Breeding and Genetics Group, Department of Animal Sciences, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
- Center for Integrated Breeding Research, University of Goettingen, Albrecht-Thaer-Weg 3, 37075, Goettingen, Germany
| | | |
Collapse
|
10
|
Cádiz MI, López ME, Díaz-Domínguez D, Cáceres G, Yoshida GM, Gomez-Uchida D, Yáñez JM. Whole genome re-sequencing reveals recent signatures of selection in three strains of farmed Nile tilapia (Oreochromis niloticus). Sci Rep 2020; 10:11514. [PMID: 32661317 PMCID: PMC7359307 DOI: 10.1038/s41598-020-68064-5] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2019] [Accepted: 06/16/2020] [Indexed: 01/24/2023] Open
Abstract
Nile tilapia belongs to the second most cultivated group of fish in the world, mainly because of its favorable characteristics for production. Genetic improvement programs and domestication process of Nile tilapia may have modified the genome through selective pressure, leaving signals that can be detected at the molecular level. In this work, signatures of selection were identified using genome-wide SNP data, by two haplotype-based (iHS and Rsb) and one FST based method. Whole-genome re-sequencing of 326 individuals from three strains (A, B and C) of farmed tilapia maintained in Brazil and Costa Rica was carried out using Illumina HiSeq 2500 technology. After applying conventional SNP-calling and quality-control filters, ~ 1.3 M high-quality SNPs were inferred and used as input for the iHS, Rsb and FST based methods. We detected several candidate genes putatively subjected to selection in each strain. A considerable number of these genes are associated with growth (e.g. NCAPG, KLF3, TBC1D1, TTN), early development (e.g. FGFR3, PFKFB3), and immunity traits (e.g. NLRC3, PIGR, MAP1S). These candidate genes represent putative genomic landmarks that could be associated to traits of biological and commercial interest in farmed Nile tilapia.
Collapse
Affiliation(s)
- María I Cádiz
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Avenida Santa Rosa 11735, 8820808, La Pintana, Santiago, Chile.,Programa de Doctorado en Ciencias Silvoagropecuarias y Veterinarias, Campus Sur, Universidad de Chile, Santa Rosa 11315, 8820808, La Pintana, Santiago, Chile
| | - María E López
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Avenida Santa Rosa 11735, 8820808, La Pintana, Santiago, Chile.,Department of Animal Breeding and Genetics, Swedish University of Agriculturall Sciences, Uppsala, Sweden
| | - Diego Díaz-Domínguez
- Departamento de Ciencias de la Computación, Universidad de Chile, Santiago, Chile
| | - Giovanna Cáceres
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Avenida Santa Rosa 11735, 8820808, La Pintana, Santiago, Chile.,Programa de Doctorado en Ciencias Silvoagropecuarias y Veterinarias, Campus Sur, Universidad de Chile, Santa Rosa 11315, 8820808, La Pintana, Santiago, Chile
| | - Grazyella M Yoshida
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Avenida Santa Rosa 11735, 8820808, La Pintana, Santiago, Chile
| | - Daniel Gomez-Uchida
- Facultad de Ciencias Naturales y Oceanográficas, Universidad de Concepción, Concepción, Chile.,Núcleo Milenio INVASAL, Concepción, Chile
| | - José M Yáñez
- Facultad de Ciencias Veterinarias y Pecuarias, Universidad de Chile, Avenida Santa Rosa 11735, 8820808, La Pintana, Santiago, Chile. .,Núcleo Milenio INVASAL, Concepción, Chile.
| |
Collapse
|
11
|
Abstract
An experiment of divergent selection for intramuscular fat was carried out at Universitat Politècnica de València. The high response of selection in intramuscular fat content, after nine generations of selection, and a multidimensional scaling analysis showed a high degree of genomic differentiation between the two divergent populations. Therefore, local genomic differences could link genomic regions, encompassing selective sweeps, to the trait used as selection criterion. In this sense, the aim of this study was to identify genomic regions related to intramuscular fat through three methods for detection of selection signatures and to generate a list of candidate genes. The methods implemented in this study were Wright's fixation index, cross population composite likelihood ratio and cross population - extended haplotype homozygosity. Genomic data came from the 9th generation of the two populations divergently selected, 237 from Low line and 240 from High line. A high single nucleotide polymorphism (SNP) density array, Affymetrix Axiom OrcunSNP Array (around 200k SNPs), was used for genotyping samples. Several genomic regions distributed along rabbit chromosomes (OCU) were identified as signatures of selection (SNPs having a value above cut-off of 1%) within each method. In contrast, 8 genomic regions, harbouring 80 SNPs (OCU1, OCU3, OCU6, OCU7, OCU16 and OCU17), were identified by at least 2 methods and none by the 3 methods. In general, our results suggest that intramuscular fat selection influenced multiple genomic regions which can be a consequence of either only selection effect or the combined effect of selection and genetic drift. In addition, 73 genes were retrieved from the 8 selection signatures. After functional and enrichment analyses, the main genes into the selection signatures linked to energy, fatty acids, carbohydrates and lipid metabolic processes were ACER2, PLIN2, DENND4C, RPS6, RRAGA (OCU1), ST8SIA6, VIM (OCU16), RORA, GANC and PLA2G4B (OCU17). This genomic scan is the first study using rabbits from a divergent selection experiment. Our results pointed out a large polygenic component of the intramuscular fat content. Besides, promising positional candidate genes would be analysed in further studies in order to bear out their contributions to this trait and their feasible implications for rabbit breeding programmes.
Collapse
|
12
|
Comparative selection signature analyses identify genomic footprints in Reggiana cattle, the traditional breed of the Parmigiano-Reggiano cheese production system. Animal 2020; 14:921-932. [PMID: 31928542 DOI: 10.1017/s1751731119003318] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Reggiana is an autochthonous cattle breed reared mainly in the province of Reggio Emilia, located in the North of Italy. Reggiana cattle (originally a triple-purpose population largely diffused in the North of Italy) are characterised by a typical solid red coat colour. About 2500 cows of this breed are currently registered to its herd book. Reggiana is now considered a dual-purpose breed even if it is almost completely dedicated to the production of a mono-breed branded Protected Designation of Origin Parmigiano-Reggiano cheese, which is the main driver of the sustainable conservation of this local genetic resource. In this study, we provided the first overview of genomic footprints that characterise Reggiana and define the diversity of this local cattle breed. A total of 168 Reggiana sires (all bulls born over 35 years for which semen was available) and other 3321 sires from 3 cosmopolitan breeds (Brown, Holstein and Simmental) were genotyped with the Illumina BovineSNP50 panel. ADMIXTURE analysis suggested that Reggiana breed might have been influenced, at least in part, by the other three breeds included in this study. Selection signatures in the Reggiana genome were identified using three statistical approaches based on allele frequency differences among populations or on properties of haplotypes segregating in the populations (fixation index (FST); integrated haplotype score; cross-population extended haplotype homozygosity). We identified several regions under peculiar selection in the Reggiana breed, particularly on bovine chromosome (BTA) 6 in the KIT gene region, that is known to be involved in coat colour pattern distribution, and within the region of the LAP3, NCAPG and LCORL genes, that are associated with stature, conformation and carcass traits. Another already known region that includes the PLAG1 gene (BTA14), associated with conformation traits, showed a selection signature in the Reggiana cattle. On BTA18, a signal of selection included the MC1R gene that causes the red coat colour in cattle. Other selection sweeps were in regions, with high density of quantitative trait loci for milk production traits (on BTA20) and in several other large regions that might have contributed to shape and define the Reggiana genome (on BTA17 and BTA29). All these results, overall, indicate that the Reggiana genome might still contain several signs of its multipurpose and non-specialised utilisation, as already described for other local cattle populations, in addition to footprints derived by its ancestral origin and by its adaptation to the specialised Parmigiano-Reggiano cheese production system.
Collapse
|
13
|
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]
|
14
|
Mouresan EF, González-Rodríguez A, Cañas-Álvarez JJ, Munilla S, Altarriba J, Díaz C, Baró JA, Molina A, Lopez-Buesa P, Piedrafita J, Varona L. Mapping Recombination Rate on the Autosomal Chromosomes Based on the Persistency of Linkage Disequilibrium Phase Among Autochthonous Beef Cattle Populations in Spain. Front Genet 2019; 10:1170. [PMID: 31824571 PMCID: PMC6880760 DOI: 10.3389/fgene.2019.01170] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2018] [Accepted: 10/23/2019] [Indexed: 01/14/2023] Open
Abstract
In organisms with sexual reproduction, genetic diversity, and genome evolution are governed by meiotic recombination caused by crossing-over, which is known to vary within the genome. In this study, we propose a simple method to estimate the recombination rate that makes use of the persistency of linkage disequilibrium (LD) phase among closely related populations. The biological material comprised 171 triplets (sire/dam/offspring) from seven populations of autochthonous beef cattle in Spain (Asturiana de los Valles, Avileña-Negra Ibérica, Bruna dels Pirineus, Morucha, Pirenaica, Retinta, and Rubia Gallega), which were genotyped for 777,962 SNPs with the BovineHD BeadChip. After standard quality filtering, we reconstructed the haplotype phases in the parental individuals and calculated the LD by the correlation -r- between each pair of markers that had a genetic distance < 1 Mb. Subsequently, these correlations were used to calculate the persistency of LD phase between each pair of populations along the autosomal genome. Therefore, the distribution of the recombination rate along the genome can be inferred since the effect of the number of generations of divergence should be equivalent throughout the genome. In our study, the recombination rate was highest in the largest chromosomes and at the distal portion of the chromosomes. In addition, the persistency of LD phase was highly heterogeneous throughout the genome, with a ratio of 25.4 times between the estimates of the recombination rates from the genomic regions that had the highest (BTA18-7.1 Mb) and the lowest (BTA12-42.4 Mb) estimates. Finally, an overrepresentation enrichment analysis (ORA) showed differences in the enriched gene ontology (GO) terms between the genes located in the genomic regions with estimates of the recombination rate over (or below) the 95th (or 5th) percentile throughout the autosomal genome.
Collapse
Affiliation(s)
- Elena Flavia Mouresan
- Departamento de Anatomía, Embriología y Genética Animal, Universidad de Zaragoza, Zaragoza, Spain
| | | | | | - Sebastián Munilla
- Departamento de Anatomía, Embriología y Genética Animal, Universidad de Zaragoza, Zaragoza, Spain.,Departamento de Producción Animal, Facultad de Agronomía, Universidad de Buenos Aires, CONICET, Buenos Aires, Argentina
| | - Juan Altarriba
- Departamento de Anatomía, Embriología y Genética Animal, Universidad de Zaragoza, Zaragoza, Spain
| | - Clara Díaz
- Departamento de Mejora Genética Animal, Instituto Nacional de Investigación y Tecnología Agraria y Alimentaria (INIA), Madrid, Spain
| | - Jesús A Baró
- Instituto Agroalimentario de Aragón (IA2), Zaragoza, Spain
| | - Antonio Molina
- Departamento de Ciencias Agroforestales, Universidad de Valladolid, Valladolid, Spain
| | - Pascual Lopez-Buesa
- Departamento de Anatomía, Embriología y Genética Animal, Universidad de Zaragoza, Zaragoza, Spain
| | - Jesús Piedrafita
- Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Luis Varona
- Departamento de Anatomía, Embriología y Genética Animal, Universidad de Zaragoza, Zaragoza, Spain
| |
Collapse
|
15
|
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
|
16
|
Weigand H, Leese F. Detecting signatures of positive selection in non-model species using genomic data. Zool J Linn Soc 2018. [DOI: 10.1093/zoolinnean/zly007] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Affiliation(s)
- Hannah Weigand
- Aquatic Ecosystem Research, University of Duisburg-Essen, Universitätsstraße, Essen, Germany
| | - Florian Leese
- Aquatic Ecosystem Research, University of Duisburg-Essen, Universitätsstraße, Essen, Germany
- Centre for Water and Environmental Research (ZWU), University of Duisburg-Essen, Universitätsstraße, Essen, Germany
| |
Collapse
|
17
|
Ruiz-Larrañaga O, Langa J, Rendo F, Manzano C, Iriondo M, Estonba A. Genomic selection signatures in sheep from the Western Pyrenees. Genet Sel Evol 2018; 50:9. [PMID: 29566643 PMCID: PMC5865298 DOI: 10.1186/s12711-018-0378-x] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2017] [Accepted: 02/12/2018] [Indexed: 02/07/2023] Open
Abstract
BACKGROUND The current large spectrum of sheep phenotypic diversity results from the combined product of sheep selection for different production traits such as wool, milk and meat, and its natural adaptation to new environments. In this study, we scanned the genome of 25 Sasi Ardi and 75 Latxa sheep from the Western Pyrenees for three types of regions under selection: (1) regions underlying local adaptation of Sasi Ardi semi-feral sheep, (2) regions related to a long traditional dairy selection pressure in Latxa sheep, and (3) regions experiencing the specific effect of the modern genetic improvement program established for the Latxa breed during the last three decades. RESULTS Thirty-two selected candidate regions including 147 annotated genes were detected by using three statistical parameters: pooled heterozygosity H, Tajima's D, and Wright's fixation index Fst. For Sasi Ardi sheep, chromosomes Ovis aries (OAR)4, 6, and 22 showed the strongest signals and harbored several candidate genes related to energy metabolism and morphology (BBS9, ELOVL3 and LDB1), immunity (NFKB2), and reproduction (H2AFZ). The major genomic difference between Sasi Ardi and Latxa sheep was on OAR6, which is known to affect milk production, with highly selected regions around the ABCG2, SPP1, LAP3, NCAPG, LCORL, and MEPE genes in Latxa sheep. The effect of the modern genetic improvement program on Latxa sheep was also evident on OAR15, on which several olfactory genes are located. We also detected several genes involved in reproduction such as ESR1 and ZNF366 that were affected by this selection program. CONCLUSIONS Natural and artificial selection have shaped the genome of both Sasi Ardi and Latxa sheep. Our results suggest that Sasi Ardi traits related to energy metabolism, morphological, reproductive, and immunological features have been under positive selection to adapt this semi-feral sheep to its particular environment. The highly selected Latxa sheep for dairy production showed clear signatures of selection in genomic regions related to milk production. Furthermore, our data indicate that the selection criteria applied in the modern genetic improvement program affect immunity and reproduction traits.
Collapse
Affiliation(s)
- Otsanda Ruiz-Larrañaga
- Genetics, Physical Anthropology and Animal Physiology Department, University of the Basque Country (UPV/EHU), Leioa, Spain.
| | - Jorge Langa
- Genetics, Physical Anthropology and Animal Physiology Department, University of the Basque Country (UPV/EHU), Leioa, Spain
| | - Fernando Rendo
- Genetics, Sequencing and Genotyping Unit, Advanced Research Facilities (SGIker), University of the Basque Country (UPV/EHU), Leioa, Spain
| | - Carmen Manzano
- Genetics, Physical Anthropology and Animal Physiology Department, University of the Basque Country (UPV/EHU), Leioa, Spain
| | - Mikel Iriondo
- Genetics, Physical Anthropology and Animal Physiology Department, University of the Basque Country (UPV/EHU), Leioa, Spain
| | - Andone Estonba
- Genetics, Physical Anthropology and Animal Physiology Department, University of the Basque Country (UPV/EHU), Leioa, Spain
| |
Collapse
|
18
|
Gutiérrez-Gil B, Esteban-Blanco C, Wiener P, Chitneedi PK, Suarez-Vega A, Arranz JJ. High-resolution analysis of selection sweeps identified between fine-wool Merino and coarse-wool Churra sheep breeds. Genet Sel Evol 2017; 49:81. [PMID: 29115919 PMCID: PMC5674817 DOI: 10.1186/s12711-017-0354-x] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2016] [Accepted: 10/19/2017] [Indexed: 12/20/2022] Open
Abstract
BACKGROUND With the aim of identifying selection signals in three Merino sheep lines that are highly specialized for fine wool production (Australian Industry Merino, Australian Merino and Australian Poll Merino) and considering that these lines have been subjected to selection not only for wool traits but also for growth and carcass traits and parasite resistance, we contrasted the OvineSNP50 BeadChip (50 K-chip) pooled genotypes of these Merino lines with the genotypes of a coarse-wool breed, phylogenetically related breed, Spanish Churra dairy sheep. Genome re-sequencing datasets of the two breeds were analyzed to further explore the genetic variation of the regions initially identified as putative selection signals. RESULTS Based on the 50 K-chip genotypes, we used the overlapping selection signals (SS) identified by four selection sweep mapping analyses (that detect genetic differentiation, reduced heterozygosity and patterns of haplotype diversity) to define 18 convergence candidate regions (CCR), five associated with positive selection in Australian Merino and the remainder indicating positive selection in Churra. Subsequent analysis of whole-genome sequences from 15 Churra and 13 Merino samples identified 142,400 genetic variants (139,745 bi-allelic SNPs and 2655 indels) within the 18 defined CCR. Annotation of 1291 variants that were significantly associated with breed identity between Churra and Merino samples identified 257 intragenic variants that caused 296 functional annotation variants, 275 of which were located across 31 coding genes. Among these, four synonymous and four missense variants (NPR2_His847Arg, NCAPG_Ser585Phe, LCORL_Asp1214Glu and LCORL_Ile1441Leu) were included. CONCLUSIONS Here, we report the mapping and genetic variation of 18 selection signatures that were identified between Australian Merino and Spanish Churra sheep breeds, which were validated by an additional contrast between Spanish Merino and Churra genotypes. Analysis of whole-genome sequencing datasets allowed us to identify divergent variants that may be viewed as candidates involved in the phenotypic differences for wool, growth and meat production/quality traits between the breeds analyzed. The four missense variants located in the NPR2, NCAPG and LCORL genes may be related to selection sweep regions previously identified and various QTL reported in sheep in relation to growth traits and carcass composition.
Collapse
Affiliation(s)
- Beatriz Gutiérrez-Gil
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| | - Cristina Esteban-Blanco
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
- Fundación Centro Supercomputación de Castilla y León, Campus de Vegazana, León, 24071 Spain
| | - Pamela Wiener
- Roslin Institute and R(D)SVS, University of Edinburgh, Easter Bush Campus, Midlothian, EH25 9RG UK
| | - Praveen Krishna Chitneedi
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| | - Aroa Suarez-Vega
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| | - Juan-Jose Arranz
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad de León, Campus de Vegazana s/n, León, 24071 Spain
| |
Collapse
|
19
|
Li X, Yang S, Dong K, Tang Z, Li K, Fan B, Wang Z, Liu B. Identification of positive selection signatures in pigs by comparing linkage disequilibrium variances. Anim Genet 2017; 48:600-605. [DOI: 10.1111/age.12574] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/11/2017] [Indexed: 11/26/2022]
Affiliation(s)
- X. Li
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education; Key Laboratory of Pig Genetics and Breeding of Ministry of Agriculture; Huazhong Agricultural University; Wuhan Hubei 430070 China
- The Cooperative Innovation Center for Sustainable Pig Production; Wuhan Hubei 430070 China
- Department of Agricultural, Food and Nutritional Science; University of Alberta; Edmonton AB T6G 2P5 Canada
| | - S. Yang
- College of Animal Science and Technology; Zhejiang A&F University; Lin'an Zhejiang 311300 China
| | - K. Dong
- The Key Laboratory for Domestic Animal Genetic Resources and Breeding of Ministry of Agriculture of China; Institute of Animal Science; Chinese Academy of Agricultural Sciences; Beijing 100193 China
| | - Z. Tang
- The Key Laboratory for Domestic Animal Genetic Resources and Breeding of Ministry of Agriculture of China; Institute of Animal Science; Chinese Academy of Agricultural Sciences; Beijing 100193 China
| | - K. Li
- The Key Laboratory for Domestic Animal Genetic Resources and Breeding of Ministry of Agriculture of China; Institute of Animal Science; Chinese Academy of Agricultural Sciences; Beijing 100193 China
| | - B. Fan
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education; Key Laboratory of Pig Genetics and Breeding of Ministry of Agriculture; Huazhong Agricultural University; Wuhan Hubei 430070 China
| | - Z. Wang
- Department of Agricultural, Food and Nutritional Science; University of Alberta; Edmonton AB T6G 2P5 Canada
| | - B. Liu
- Key Laboratory of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education; Key Laboratory of Pig Genetics and Breeding of Ministry of Agriculture; Huazhong Agricultural University; Wuhan Hubei 430070 China
- The Cooperative Innovation Center for Sustainable Pig Production; Wuhan Hubei 430070 China
| |
Collapse
|
20
|
Mouresan E, González-Rodríguez A, Cañas-Álvarez J, Díaz C, Altarriba J, Baro J, Piedrafita J, Molina A, Toro M, Varona L. On the haplotype diversity along the genome in Spanish beef cattle populations. Livest Sci 2017. [DOI: 10.1016/j.livsci.2017.04.015] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|