1
|
Li C, Wang X, Li H, Ahmed Z, Luo Y, Qin M, Yang Q, Long Z, Lei C, Yi K. Whole-genome resequencing reveals diversity and selective signals in the Wuxue goat. Anim Genet 2024. [PMID: 38806279 DOI: 10.1111/age.13437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2023] [Revised: 04/14/2024] [Accepted: 04/17/2024] [Indexed: 05/30/2024]
Abstract
Animal genetic resources are crucial for ensuring global food security. However, in recent years, a noticeable decline in the genetic diversity of livestock has occurred worldwide. This decline is pronounced in developing countries, where the management of these resources is insufficient. In the current study, we performed whole genome sequencing for 20 Wuxue (WX) and five Guizhou White (GW) goats. Additionally, we utilized the published genomes of 131 samples representing five different goat breeds from various regions in China. We investigated and compared the genetic diversity and selection signatures of WX goats. Whole genome sequencing analysis of the WX and GW populations yielded 120 425 063 SNPs, which resided primarily in intergenic and intron regions. Population genetic structure revealed that WX exhibited genetic resemblance to GW, Chengdu Brown, and Jintang Black and significant differentiation from the other goat breeds. In addition, three methods (nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity) showed moderate genetic diversity in WX goats. We used nucleotide diversity and composite likelihood ratio methods to identify within-breed signatures of positive selection in WX goats. A total of 369 genes were identified using both detection methods, including genes related to reproduction (GRID2, ZNF276, TCF25, and SPIRE2), growth (HMGA2 and GJA3), and immunity (IRF3 and SRSF3). Overall, this study explored the adaptability of WX goats, shedding light on their genetic richness and potential to thrive in challenges posed by climatic changes and diseases. Further investigations are warranted to harness these insights to enhance more efficient and sustainable goat breeding initiatives.
Collapse
Affiliation(s)
- Chuanqing Li
- Hunan Institute of Animal and Veterinary Science, Changsha, China
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Xianglin Wang
- Animal Husbandry and Aquatic Products Affairs Center of Xiangxi Autonomous Prefecture, Jishou, China
| | - Haobang Li
- Hunan Institute of Animal and Veterinary Science, Changsha, China
| | - Zulfiqar Ahmed
- Faculty of Veterinary and Animal Sciences, University of Poonch Rawalakot, Rawalakot, Pakistan
| | - Yang Luo
- Hunan Institute of Animal and Veterinary Science, Changsha, China
| | - Mao Qin
- Animal Husbandry and Aquatic Products Affairs Center of Xiangxi Autonomous Prefecture, Jishou, China
| | - Qiong Yang
- Animal Husbandry and Aquatic Products Affairs Center of Xiangxi Autonomous Prefecture, Jishou, China
| | - Zhangcheng Long
- Animal Husbandry and Aquatic Products Affairs Center of Xiangxi Autonomous Prefecture, Jishou, China
| | - Chuzhao Lei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Kangle Yi
- Hunan Institute of Animal and Veterinary Science, Changsha, China
| |
Collapse
|
2
|
Xie K, Ning C, Yang A, Zhang Q, Wang D, Fan X. Resequencing Analyses Revealed Genetic Diversity and Selection Signatures during Rabbit Breeding and Improvement. Genes (Basel) 2024; 15:433. [PMID: 38674368 PMCID: PMC11049387 DOI: 10.3390/genes15040433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 03/23/2024] [Accepted: 03/26/2024] [Indexed: 04/28/2024] Open
Abstract
Domestication has shaped the diverse characteristics of rabbits, including coat color, fur structure, body size, and various physiological traits. Utilizing whole-genome resequencing (DNBSEQ-T7), we analyzed the genetic diversity, population structure, and genomic selection across 180 rabbits from 17 distinct breeds to uncover the genetic basis of these traits. We conducted whole-genome sequencing on 17 rabbit breeds, identifying 17,430,184 high-quality SNPs and analyzing genomic diversity, patterns of genomic variation, population structure, and selection signatures related to coat color, coat structure, long hair, body size, reproductive capacity, and disease resistance. Through PCA and NJ tree analyses, distinct clusters emerged among Chinese indigenous rabbits, suggesting varied origins and domestication histories. Selective sweep testing pinpointed regions and genes linked to domestication and key morphological and economic traits, including those affecting coat color (TYR, ASIP), structure (LIPH), body size (INSIG2, GLI3), fertility (EDNRA, SRD5A2), heat stress adaptation (PLCB1), and immune response (SEC31A, CD86, LAP3). Our study identified key genomic signatures of selection related to traits such as coat color, fur structure, body size, and fertility; these findings highlight the genetic basis underlying phenotypic diversification in rabbits and have implications for breeding programs aiming to improve productive, reproductive, and adaptive traits. The detected genomic signatures of selection also provide insights into rabbit domestication and can aid conservation efforts for indigenous breeds.
Collapse
Affiliation(s)
- Kerui Xie
- College of Animal Science and Veterinary Medicine, Shandong Agricultural University, Tai’an 271018, China;
| | - Chao Ning
- Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Shandong Agricultural University, 61 Daizong Street, Tai’an 271018, China; (C.N.); (Q.Z.)
| | - Aiguo Yang
- Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Shandong Agricultural University, 61 Daizong Street, Tai’an 271018, China; (C.N.); (Q.Z.)
| | - Qin Zhang
- Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Shandong Agricultural University, 61 Daizong Street, Tai’an 271018, China; (C.N.); (Q.Z.)
| | - Dan Wang
- Key Laboratory of Efficient Utilization of Non-Grain Feed Resources (Co-Construction by Ministry and Province), Ministry of Agriculture and Rural Affairs, College of Animal Science and Veterinary Medicine, Shandong Agricultural University, Tai’an 271018, China
| | - Xinzhong Fan
- College of Animal Science and Veterinary Medicine, Shandong Agricultural University, Tai’an 271018, China;
| |
Collapse
|
3
|
Sarviaho K, Uimari P, Martikainen K. Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population. J Dairy Sci 2024:S0022-0302(24)00548-4. [PMID: 38490540 DOI: 10.3168/jds.2024-24105] [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: 08/21/2023] [Accepted: 02/15/2024] [Indexed: 03/17/2024]
Abstract
The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits, e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.
Collapse
Affiliation(s)
- Katri Sarviaho
- Department of Agricultural Sciences, University of Helsinki, Helsinki, Finland.
| | - Pekka Uimari
- Department of Agricultural Sciences, University of Helsinki, Helsinki, Finland
| | - Katja Martikainen
- Department of Agricultural Sciences, University of Helsinki, Helsinki, Finland
| |
Collapse
|
4
|
Liu J, Xiao Y, Ren P, Zhang S, Liu Y, Zhu M. Integrating genomics and transcriptomics to identify candidate genes for high egg production in Wulong geese (Anser cygnoides orientalis). BMC Genomics 2023; 24:481. [PMID: 37620752 PMCID: PMC10464066 DOI: 10.1186/s12864-023-09603-y] [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: 05/12/2023] [Accepted: 08/18/2023] [Indexed: 08/26/2023] Open
Abstract
BACKGROUND Wulong geese (Anser cygnoides orientalis) are known for their excellent egg-laying performance. However, they show considerable population differences in egg-laying behavior. This study combined genome-wide selection signal analysis with transcriptome analysis (RNA-seq) to identify the genes related to high egg production in Wulong geese. RESULTS A total of 132 selected genomic regions were screened using genome-wide selection signal analysis, and 130 genes related to high egg production were annotated in these regions. These selected genes were enriched in pathways related to egg production, including oocyte meiosis, the estrogen signaling pathway, the oxytocin signaling pathway, and progesterone-mediated oocyte maturation. Furthermore, a total of 890 differentially expressed genes (DEGs), including 340 up-regulated and 550 down-regulated genes, were identified by RNA-seq. Two genes - GCG and FAP - were common to the list of selected genes and DEGs. A non-synonymous single nucleotide polymorphism was identified in an exon of FAP. CONCLUSIONS Based on genome-wide selection signal analysis and transcriptome data, GCG and FAP were identified as candidate genes associated with high egg production in Wulong geese. These findings could promote the breeding of Wulong geese with high egg production abilities and provide a theoretical basis for exploring the mechanisms of reproductive regulation in poultry.
Collapse
Affiliation(s)
- Jingjing Liu
- College of Agronomy and Agricultural Engineering, Liaocheng University, Liaocheng, 252000, China
| | - Yu Xiao
- College of Agronomy and Agricultural Engineering, Liaocheng University, Liaocheng, 252000, China
| | - Pengwei Ren
- College of Agronomy and Agricultural Engineering, Liaocheng University, Liaocheng, 252000, China
| | - Shuer Zhang
- Shandong Animal Husbandry General Station, Jinan, 250010, China
| | - Yang Liu
- College of Agronomy and Agricultural Engineering, Liaocheng University, Liaocheng, 252000, China
| | - Mingxia Zhu
- College of Agronomy and Agricultural Engineering, Liaocheng University, Liaocheng, 252000, China.
| |
Collapse
|
5
|
Hewett AM, Stoffel MA, Peters L, Johnston SE, Pemberton JM. Selection, recombination and population history effects on runs of homozygosity (ROH) in wild red deer (Cervus elaphus). Heredity (Edinb) 2023; 130:242-250. [PMID: 36801920 PMCID: PMC10076382 DOI: 10.1038/s41437-023-00602-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 02/07/2023] [Accepted: 02/07/2023] [Indexed: 02/19/2023] Open
Abstract
The distribution of runs of homozygosity (ROH) may be shaped by a number of interacting processes such as selection, recombination and population history, but little is known about the importance of these mechanisms in shaping ROH in wild populations. We combined an empirical dataset of >3000 red deer genotyped at >35,000 genome-wide autosomal SNPs and evolutionary simulations to investigate the influence of each of these factors on ROH. We assessed ROH in a focal and comparison population to investigate the effect of population history. We investigated the role of recombination using both a physical map and a genetic linkage map to search for ROH. We found differences in ROH distribution between both populations and map types indicating that population history and local recombination rate have an effect on ROH. Finally, we ran forward genetic simulations with varying population histories, recombination rates and levels of selection, allowing us to further interpret our empirical data. These simulations showed that population history has a greater effect on ROH distribution than either recombination or selection. We further show that selection can cause genomic regions where ROH is common, only when the effective population size (Ne) is large or selection is particularly strong. In populations having undergone a population bottleneck, genetic drift can outweigh the effect of selection. Overall, we conclude that in this population, genetic drift resulting from a historical population bottleneck is most likely to have resulted in the observed ROH distribution, with selection possibly playing a minor role.
Collapse
Affiliation(s)
- Anna M Hewett
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK.
| | - Martin A Stoffel
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
| | - Lucy Peters
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
- GenPhySE, Université de Toulouse, INRAE, ENVT, 31326, Castanet-Tolosan, France
| | - Susan E Johnston
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
| | - Josephine M Pemberton
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, EH9 3FL, UK
| |
Collapse
|
6
|
Chen J, Wang Y, Qi X, Cheng H, Chen N, Ahmed Z, Chen Q, Lei C, Yang X. Genome-wide analysis emancipates genomic diversity and signature of selection in Altay white-headed cattle of Xinjiang, China. Front Genet 2023; 14:1144249. [PMID: 37065480 PMCID: PMC10098193 DOI: 10.3389/fgene.2023.1144249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Accepted: 03/17/2023] [Indexed: 04/01/2023] Open
Abstract
Altay white-headed cattle have not received enough attention for several reasons. Due to irrational breeding and selection practices, the number of pure Altay white-headed cattle has decreased significantly and the breed is now on the eve of extinction. The genomic characterization will be a crucial step towards understanding the genetic basis of productivity and adaptability to survival under native Chinese agropastoral systems; nevertheless, no attempt has been made in Altay white-headed cattle. In the current study, we compared the genomes of 20 Altay white-headed cattle to the genomes of 144 individuals in representative breeds. Population genetic diversity revealed that the nucleotide diversity of Altay white-headed cattle was less than that of indicine breeds and comparable to that of Chinese taurus cattle. Using population structure analysis, we also found that Altay white-headed cattle carried the ancestry of the European and East Asian cattle lineage. In addition, we used three different methods (FST, θπ ratio and XP-EHH) to investigate the adaptability and white-headed phenotype of Altay white-headed cattle and compared it with Bohai black cattle. We found EPB41L5, SCG5 and KIT genes on the list of the top one percent genes, these genes might have an association with environmental adaptability and the white-headed phenotype for this breed. Our research reveals the distinctive genomic features of Altay white-headed cattle at the genome-wide level.
Collapse
Affiliation(s)
- Jialei Chen
- Life Science College, Luoyang Normal University, Luoyang, China
- College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Yushu Wang
- College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Xingshan Qi
- Biyang Xianan Cattle Technology and Development Company Ltd., Biyang, China
| | - Haijian Cheng
- Shandong Key Lab of Animal Disease Control and Breeding, Institute of Animal Science and Veterinary Medicine, Shandong Academy of Agricultural Sciences, Jinan, China
| | - Ningbo Chen
- College of Animal Science and Technology, Northwest A&F University, Yangling, China
| | - Zulfiqar Ahmed
- Key Laboratory of Swine Genetics and Breeding of Ministry of Agriculture and Rural Affairs, and Key Lab of Agricultural Animal Genetics, Breeding and Reproduction of Ministry of Education, Huazhong Agricultural University, Wuhan, China
- Faculty of Veterinary and Animal Sciences, University of Poonch Rawalakot, Shabestar, Pakistan
| | - Qiuming Chen
- College of Animal Science, Xinjiang Agricultural University, Urumqi, Xinjiang, China
| | - Chuzhao Lei
- College of Animal Science and Technology, Northwest A&F University, Yangling, China
- *Correspondence: Chuzhao Lei, ; Xueyi Yang,
| | - Xueyi Yang
- Life Science College, Luoyang Normal University, Luoyang, China
- *Correspondence: Chuzhao Lei, ; Xueyi Yang,
| |
Collapse
|
7
|
Wang H, Wang Q, Tan X, Wang J, Zhang J, Zheng M, Zhao G, Wen J. Estimation of genetic variability and identification of regions under selection based on runs of homozygosity in Beijing-You Chickens. Poult Sci 2022; 102:102342. [PMID: 36470032 PMCID: PMC9719870 DOI: 10.1016/j.psj.2022.102342] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2021] [Revised: 11/08/2022] [Accepted: 11/09/2022] [Indexed: 11/16/2022] Open
Abstract
The genetic composition of populations is the result of a long-term process of selection and adaptation to specific environments and ecosystems. Runs of homozygosity (ROHs) are homozygous segments of the genome where the 2 haplotypes inherited from the parents are identical. The detection of ROH can be used to describe the genetic variability and quantify the level of inbreeding in an individual. Here, we investigated the occurrence and distribution of ROHs in 40 Beijing-You Chickens from the random breeding population (BJY_C) and 40 Beijing-You Chickens from the intramuscular fat (IMF) selection population (BJY_S). Principal component analysis (PCA) and maximum likelihood (ML) analyses showed that BJY_C was completely separated from the BJY_S. The nucleotide diversity of BJY_C was higher than that of BJY_S, and the decay rate of LD of BJY_C was faster. The ROHs were identified for a total of 7,101 in BJY_C and 9,273 in BJY_S, respectively. The ROH-based inbreeding estimate (FROH) of BJY_C was 0.079, which was significantly lower than that of BJY_S (FROH = 0.114). The results were the same as the estimates of the inbreeding coefficients calculated based on homozygosity (FHOM), the correlation between uniting gametes (FUNI), and the genomic relationship matrix (FGRM). Additionally, the distribution and number of ROH islands in chromosomes of BJY_C and BJY_S were significantly different. The ROH islands of BJY_S that included genes associated with lipid metabolism and fat deposition, such as CIDEA and S1PR1, were absent in BJY_C. However, GPR161 was detected in both populations, which is a candidate gene for the formation of the unique five-finger trait in Beijing-You chickens. Our findings contributed to the understanding of the genetic diversity of random or artificially selected populations, and allowed the accurate monitoring of population inbreeding using genomic information, as well as the detection of genomic regions that affect traits under selection.
Collapse
Affiliation(s)
- Hailong Wang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Qiao Wang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Xiaodong Tan
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Jie Wang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Jin Zhang
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Maiqing Zheng
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Guiping Zhao
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China
| | - Jie Wen
- Chinese Academy of Agricultural Science, State Key Laboratory of Animal Nutrition, Beijing 100193, China.
| |
Collapse
|
8
|
Wang G, Wang F, Pei H, Li M, Bai F, Lei C, Dang R. Genome-wide analysis reveals selection signatures for body size and drought adaptation in Liangzhou donkey. Genomics 2022; 114:110476. [PMID: 36057425 DOI: 10.1016/j.ygeno.2022.110476] [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/28/2021] [Revised: 07/21/2022] [Accepted: 08/12/2022] [Indexed: 01/14/2023]
Abstract
Liangzhou donkey is a domestic animal breed distributed on the edge of the Tengger Desert in Gansu Province of China. It has small body size and strong adaptability to dry environments. Here, we sequenced 10 Liangzhou donkey genomes and compared them to the 55 genomes of 8 representative donkey breeds worldwide. The population structure analysis revealed that Liangzhou donkey harboured the ancestry with the Asian domestic donkeys (0.863) and European domestic donkeys (0.137). Three methods (nucleotide diversity, linkage disequilibrium decay and runs of homozygosity) implied the genetic diversity in Liangzhou donkey. In addition, we analyzed the genetic basis of the small body size and drought adaptation of Liangzhou donkey by using Fst, θπ-ratio, XP-EHH, CLR and θπ methods. We found that the NCAPG-LCORL on chromosome 3 may be a candidate region for small body size trait of Liangzhou donkey. The CYP4A11 gene located on chromosome 5 showed strong sign of selection sweep. CYP4A11 can convert arachidonic acid into 19(S)-HETE, which can promote water reabsorption in renal tubule and enhance the ability of Liangzhou donkey to adapt to dry environment. These results contribute to a better understanding of the underlying population structure of Liangzhou donkeys and provides a valuable resource for future research on donkey breeding in response to climate change.
Collapse
Affiliation(s)
- Gang Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang 712100, Shaanxi Province, China.
| | - Fuwen Wang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang 712100, Shaanxi Province, China.
| | - Haoyu Pei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang 712100, Shaanxi Province, China.
| | - Mei Li
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang 712100, Shaanxi Province, China
| | - Fuxia Bai
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang 712100, Shaanxi Province, China.
| | - Chuzhao Lei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang 712100, Shaanxi Province, China.
| | - Ruihua Dang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, Yangling, Xianyang 712100, Shaanxi Province, China.
| |
Collapse
|
9
|
Assessment of Genetic Diversity, Runs of Homozygosity, and Signatures of Selection in Tropical Milking Criollo Cattle Using Pedigree and Genomic Data. Genes (Basel) 2022; 13:genes13101896. [PMID: 36292782 PMCID: PMC9602073 DOI: 10.3390/genes13101896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Revised: 10/15/2022] [Accepted: 10/17/2022] [Indexed: 11/10/2022] Open
Abstract
The objective of this study was to evaluate the genetic diversity of the Tropical Milking Criollo cattle (TMC) breed in Mexico through parameters derived from pedigree and genomic information assessment. The pedigree file consisted of 3780 animals. Seventy-nine bovines were genotyped with the medium-density single nucleotide polymorphism chip and considered a reference population for pedigree analysis. The effective population size and the probability of gene origin used to assess the evolution of genetic diversity were calculated with pedigree information. Inbreeding coefficients were evaluated based on pedigree (FPed), the genomic relationship matrix (FGRM), and runs of homozygosity (FROH) of different length classes. The average inbreeding was 2.82 ± 2.66%, −0.7 ± 3.8%, and 10.9 ± 3.0% for FPED, FGRM, and FROH, respectively. Correlation between FPED and FROH was significant only for runs of homozygosity > 4 Mb, indicating the FPED of a population with an average equivalent complete generation of five only recovers the most recent inbreeding. The parameters of the probability of gene origin indicated the existence of genetic bottlenecks and the loss of genetic diversity in the history of the TMC cattle population; however, pedigree and genomic information revealed the existence of current sufficient genetic diversity to design a sustainable breeding program.
Collapse
|
10
|
Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China. PLoS One 2022; 17:e0271718. [PMID: 36006904 PMCID: PMC9409551 DOI: 10.1371/journal.pone.0271718] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Accepted: 07/06/2022] [Indexed: 11/19/2022] Open
Abstract
Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of FROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.
Collapse
|
11
|
Berghöfer J, Khaveh N, Mundlos S, Metzger J. Simultaneous testing of rule- and model-based approaches for runs of homozygosity detection opens up a window into genomic footprints of selection in pigs. BMC Genomics 2022; 23:564. [PMID: 35933356 PMCID: PMC9357325 DOI: 10.1186/s12864-022-08801-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2022] [Accepted: 07/29/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes. RESULTS Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation. CONCLUSIONS The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes.
Collapse
Affiliation(s)
- Jan Berghöfer
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Department of Biology, Chemistry and Pharmacy, Institute of Chemistry and Biochemistry, Freie Universität Berlin, Berlin, Germany
| | - Nadia Khaveh
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany
| | - Stefan Mundlos
- Research Group Development & Disease, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Institute of Medical Genetics and Human Genetics, Charité-Universitätsmedizin Berlin, Berlin, Germany.,Charité-Universitätsmedizin Berlin, BCRT, Berlin Institute of Health Center for Regenerative Therapies, Berlin, Germany
| | - Julia Metzger
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany. .,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany.
| |
Collapse
|
12
|
Dadousis C, Ablondi M, Cipolat-Gotet C, van Kaam JT, Marusi M, Cassandro M, Sabbioni A, Summer A. Genomic inbreeding coefficients using imputed genotypes: Assessing different estimators in Holstein-Friesian dairy cows. J Dairy Sci 2022; 105:5926-5945. [DOI: 10.3168/jds.2021-21125] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2021] [Accepted: 03/08/2022] [Indexed: 11/19/2022]
|
13
|
Criscione A, Mastrangelo S, D'Alessandro E, Tumino S, Di Gerlando R, Zumbo A, Marletta D, Bordonaro S. Genome-wide survey on three local horse populations with a focus on runs of homozygosity pattern. J Anim Breed Genet 2022; 139:540-555. [PMID: 35445758 PMCID: PMC9541879 DOI: 10.1111/jbg.12680] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Revised: 02/04/2022] [Accepted: 04/04/2022] [Indexed: 11/26/2022]
Abstract
Purosangue Orientale Siciliano, Sanfratellano and Siciliano represent the Sicilian equine genetic resource. This study aimed to investigate the genetic diversity, population structure and the pattern of autozygosity of Sicilian horse populations using genome‐wide single‐nucleotide polymorphism (SNP) data generated with the Illumina Equine SNP70 array. The genotyping data of 17 European and Middle East populations were also included in the study. The patterns of genetic differentiation, model‐based clustering and Neighbour‐Net showed the expected positioning of Sicilian populations within the wide analysed framework and the close connections between the Purosangue Orientale Siciliano and the Arab as well as between Sanfratellano, Siciliano and Maremmano. The highest expected heterozygosity (He) and contemporary effective population size (cNe) were reported in Siciliano (He = 0.323, cNe = 397), and the lowest were reported in Purosangue Orientale Siciliano (He = 0.277, cNe = 10). The analysis of the runs of homozygosity and the relative derived inbreeding revealed high internal homogeneity in Purosangue Orientale Siciliano and Arab horses, intermediate values in Maremmano and Sanfratellano and high heterogeneity in the Siciliano population. The genome‐wide SNP analysis showed the selective pressure on Purosangue Orientale Siciliano towards traits related to endurance performance. Our results underline the importance of planning adequate conservation and exploitation programmes to reduce the level of inbreeding and, therefore, the loss of genetic diversity.
Collapse
Affiliation(s)
- Andrea Criscione
- Dipartimento di Agricoltura, Alimentazione e Ambiente, Università di Catania, Catania, Italy
| | - Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università di Palermo, Palermo, Italy
| | | | - Serena Tumino
- Dipartimento di Agricoltura, Alimentazione e Ambiente, Università di Catania, Catania, Italy
| | - Rosalia Di Gerlando
- Dipartimento Scienze Agrarie, Alimentari e Forestali, Università di Palermo, Palermo, Italy
| | - Alessandro Zumbo
- Dipartimento di Scienze Veterinarie, Università di Messina, Messina, Italy
| | - Donata Marletta
- Dipartimento di Agricoltura, Alimentazione e Ambiente, Università di Catania, Catania, Italy
| | - Salvatore Bordonaro
- Dipartimento di Agricoltura, Alimentazione e Ambiente, Università di Catania, Catania, Italy
| |
Collapse
|
14
|
Mulim HA, Brito LF, Pinto LFB, Ferraz JBS, Grigoletto L, Silva MR, Pedrosa VB. Characterization of runs of homozygosity, heterozygosity-enriched regions, and population structure in cattle populations selected for different breeding goals. BMC Genomics 2022; 23:209. [PMID: 35291953 PMCID: PMC8925140 DOI: 10.1186/s12864-022-08384-0] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 02/10/2022] [Indexed: 01/12/2023] Open
Abstract
Background A decline in the level of genetic diversity in livestock can result in reduced response to selection, greater incidence of genetic defects, and inbreeding depression. In this context, various metrics have been proposed to assess the level of genetic diversity in selected populations. Therefore, the main goals of this study were to: 1) investigate the population structure of 16 cattle populations from 15 different pure breeds or composite populations, which have been selected for different breeds goals; and, 2) identify and compare runs of homozygosity (ROH) and heterozygosity-enriched regions (HER) based on different single nucleotide polymorphism (SNP) panels and whole-genome sequence data (WGS), followed by functional genomic analyses. Results A total of 24,187 ROH were found across all cattle populations, with 55% classified in the 2-4 Mb size group. Fourteen homozygosity islands were found in five populations, where four ROH islands located on BTA1, BTA5, BTA16, and BTA19 overlapped between the Brahman (BRM) and Gyr (GIR) breeds. A functional analysis of the genes found in these islands revealed candidate genes known to play a role in the melanogenesis, prolactin signaling, and calcium signaling pathways. The correlations between inbreeding metrics ranged from 0.02 to 0.95, where the methods based on homozygous genotypes (FHOM), uniting of gametes (FUNI), and genotype additive variance (FGRM) showed strong correlations among them. All methods yielded low to moderate correlations with the inbreeding coefficients based on runs of homozygosity (FROH). For the HER, 3576 runs and 26 islands, distributed across all autosomal chromosomes, were found in regions containing genes mainly related to the immune system, indicating potential balancing selection. Although the analyses with WGS did not enable detection of the same island patterns, it unraveled novel regions not captured when using SNP panel data. Conclusions The cattle populations that showed the largest amount of ROH and HER were Senepol (SEN) and Montana (MON), respectively. Overlapping ROH islands were identified between GIR and BRM breeds, indicating a possible historical connection between the populations. The distribution and pattern of ROH and HER are population specific, indicating that different breeds have experienced divergent selection processes or different genetic processes. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08384-0.
Collapse
Affiliation(s)
| | - Luiz F Brito
- Department of Animal Science, Purdue University, West Lafayette, Indiana, USA
| | | | - José Bento Sterman Ferraz
- Department of Animal Sciences, College of Animal Sciences and Food Engineering, University of São Paulo, Pirassununga, São Paulo, Brazil
| | - Lais Grigoletto
- Department of Animal Science, Purdue University, West Lafayette, Indiana, USA.,Department of Animal Sciences, College of Animal Sciences and Food Engineering, University of São Paulo, Pirassununga, São Paulo, Brazil
| | | | - Victor Breno Pedrosa
- Department of Animal Science, Federal University of Bahia, Salvador, Bahia, Brazil. .,Department of Animal Science, State University of Ponta Grossa, Av. General Carlos Cavalcanti, 4748 - Uvaranas, Ponta Grossa, PR, 84030-900, Brazil.
| |
Collapse
|
15
|
Fountain ED, Zhou LC, Karklus A, Liu QX, Meyers J, Fontanilla IKC, Rafael EF, Yu JY, Zhang Q, Zhu XL, Pei EL, Yuan YH, Banes GL. Cross-Species Application of Illumina iScan Microarrays for Cost-Effective, High-Throughput SNP Discovery. Front Ecol Evol 2021. [DOI: 10.3389/fevo.2021.629252] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Microarrays can be a cost-effective alternative to high-throughput sequencing for discovering novel single-nucleotide polymorphisms (SNPs). Illumina’s iScan platform dominates the market, but their commercial microarray products are designed for model organisms. Further, the platform outputs data in a proprietary format. This cannot be easily converted to human-readable genotypes or be merged with pre-existing data. To address this, we present and validate a novel pipeline to facilitate data analysis from cross-species application of Illumina microarrays. This facilitates the generation of a compatible VCF from iScan data and the merging of this with a second VCF comprising genotypes derived from other samples and sources. Our pipeline includes a custom script, iScanVCFMerge (presented as a Python package), which we validate using iScan data from three great ape genera. We conclude that cross-species application of microarrays can be a rapid, cost-effective approach for SNP discovery in non-model organisms. Our pipeline surmounts the common challenges of integrating iScan genotypes with pre-existing data.
Collapse
|
16
|
Rodríguez-Ramilo S, Reverter A, Legarra A. Islands of runs of homozygosity indicate selection signatures in Ovis aries 6 (OAR6) of French dairy sheep. JDS COMMUNICATIONS 2021; 2:132-136. [PMID: 36339500 PMCID: PMC9623631 DOI: 10.3168/jdsc.2020-0011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 02/05/2021] [Indexed: 11/21/2022]
Abstract
The presence of runs of homozygosity is not randomly distributed across the genome. Islands of runs of homozygosity may be the result of selection pressure. Concordance existed between islands of runs of homozygosity and selection signatures on OAR6. Candidate genes NCAPG and LCORL on OAR6 have agricultural and adaptive importance.
Runs of homozygosity (ROH) are contiguous homozygous segments of the genome where the haplotypes inherited from each parent are identical. The occurrence of ROH is not randomly distributed across the genome, and ROH islands across many animals may be the result of selective pressure. The objective of this study was to demonstrate that the presence of ROH islands may be indicative of selection signatures in French dairy sheep breeds and subpopulations. The data set available included animals (artificial insemination males) from various breeds and subpopulations: Basco-Béarnaise breed (321 individuals), Manech Tête Noire breed (329 individuals), Manech Tête Rousse breed (1,906 individuals), Lacaune Confederation subpopulation (3,030 individuals), and Lacaune Ovitest subpopulation (3,114 individuals). Animals were genotyped with the Illumina OvineSNP50 BeadChip. After applying filtering criteria, the genomic data included 38,287 autosomal SNP distributed across 26 chromosomes and 8,700 individuals. One island of ROH was detected on OAR6 in the same genomic position across animals (between 30 and 40 Mb). Global Wright's differentiation coefficients for 2 SNP within this ROH island were high (0.67–0.68). The linkage disequilibrium between both SNP was also elevated (0.98). The divergence in allele frequencies in those SNP grouped Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds in one cluster and Lacaune Confederation and Lacaune Ovitest subpopulations in another cluster. The closest candidate genes are NCAPG and LCORL, which have been reported to be under positive selection and suggested to control weight and height in sheep. The preliminary identification of ROH suggests the presence of selection. However, for the identification of potential candidate genes, ROH detection should be combined with other approaches to improve mapping accuracy.
Collapse
Affiliation(s)
- S.T. Rodríguez-Ramilo
- GenPhySE, Université de Toulouse, INRAE, ENVT, F-31326, Castanet Tolosan, France
- Corresponding author
| | - A. Reverter
- CSIRO Agriculture & Food, Brisbane, QLD 4067, Australia
| | - A. Legarra
- GenPhySE, Université de Toulouse, INRAE, ENVT, F-31326, Castanet Tolosan, France
| |
Collapse
|
17
|
Fang Y, Hao X, Xu Z, Sun H, Zhao Q, Cao R, Zhang Z, Ma P, Sun Y, Qi Z, Wei Q, Wang Q, Pan Y. Genome-Wide Detection of Runs of Homozygosity in Laiwu Pigs Revealed by Sequencing Data. Front Genet 2021; 12:629966. [PMID: 33995477 PMCID: PMC8116706 DOI: 10.3389/fgene.2021.629966] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2020] [Accepted: 03/08/2021] [Indexed: 11/13/2022] Open
Abstract
Laiwu pigs, distinguished by their high intramuscular fat of 7-9%, is an indigenous pig breed of China, and recent studies also found that Laiwu pigs showed high resistance to Porcine circovirus type 2. However, with the introduction of commercial varieties, the population of Laiwu pigs has declined, and some lineages have disappeared, which could result in inbreeding. Runs of homozygosity (ROH) can be used as a good measure of individual inbreeding status and is also normally used to detect selection signatures so as to map the candidate genes associated with economically important traits. In this study, we used data from Genotyping by Genome Reducing and Sequencing to investigate the number, length, coverage, and distribution patterns of ROH in 93 Chinese Laiwu pigs and identified genomic regions with a high ROH frequency. The average inbreeding coefficient calculated by pedigree was 0.021, whereas that estimated by all detected ROH segments was 0.133. Covering 13.4% of the whole genome, a total of 7,508 ROH segments longer than 1 Mb were detected, whose average length was 3.76 Mb, and short segments (1-5 Mb) dominated. For individuals, the coverage was in the range between 0.56 and 36.86%. For chromosomes, SSC6 had the largest number (n = 688), and the number of ROH in SSC12 was the lowest (n = 215). Thirteen ROH islands were detected in our study, and 86 genes were found within those regions. Some of these genes were correlated with economically important traits, such as meat quality (ECI1, LRP12, NDUFA4L2, GIL1, and LYZ), immunity capacity (IL23A, STAT2, STAT6, TBK1, IFNG, and ITH2), production (DCSTAMP, RDH16, and GDF11), and reproduction (ODF1 and CDK2). A total of six significant Gene Ontology terms and nine significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were identified, most of which were correlated with disease resistance and biosynthesis processes, and one KEGG pathway was related to lipid metabolism. In addition, we aligned all of the ROH islands to the pig quantitative trait loci (QTL) database and finally found eight QTL related to the intramuscular fat trait. These results may help us understand the characteristics of Laiwu pigs and provide insight for future breeding strategies.
Collapse
Affiliation(s)
- Yifei Fang
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Xinyu Hao
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Zhong Xu
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Hao Sun
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Qingbo Zhao
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Rui Cao
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | - Zhe Zhang
- Department of Animal Science, College of Animal Sciences, Zhejiang University, Hangzhou, China
| | - Peipei Ma
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai, China
| | | | | | | | - Qishan Wang
- Department of Animal Science, College of Animal Sciences, Zhejiang University, Hangzhou, China
| | - Yuchun Pan
- Department of Animal Science, College of Animal Sciences, Zhejiang University, Hangzhou, China
| |
Collapse
|
18
|
Macciotta NPP, Colli L, Cesarani A, Ajmone-Marsan P, Low WY, Tearle R, Williams JL. The distribution of runs of homozygosity in the genome of river and swamp buffaloes reveals a history of adaptation, migration and crossbred events. Genet Sel Evol 2021; 53:20. [PMID: 33639853 PMCID: PMC7912491 DOI: 10.1186/s12711-021-00616-3] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2020] [Accepted: 02/17/2021] [Indexed: 01/03/2023] Open
Abstract
Background Water buffalo is one of the most important livestock species in the world. Two types of water buffalo exist: river buffalo (Bubalus bubalis bubalis) and swamp buffalo (Bubalus bubalis carabanensis). The buffalo genome has been recently sequenced, and thus a new 90 K single nucleotide polymorphism (SNP) bead chip has been developed. In this study, we investigated the genomic population structure and the level of inbreeding of 185 river and 153 swamp buffaloes using runs of homozygosity (ROH). Analyses were carried out jointly and separately for the two buffalo types. Results The SNP bead chip detected in swamp about one-third of the SNPs identified in the river type. In total, 18,116 ROH were detected in the combined data set (17,784 SNPs), and 16,251 of these were unique. ROH were present in both buffalo types mostly detected (~ 59%) in swamp buffalo. The number of ROH per animal was larger and genomic inbreeding was higher in swamp than river buffalo. In the separated datasets (46,891 and 17,690 SNPs for river and swamp type, respectively), 19,760 and 10,581 ROH were found in river and swamp, respectively. The genes that map to the ROH islands are associated with the adaptation to the environment, fitness traits and reproduction. Conclusions Analysis of ROH features in the genome of the two water buffalo types allowed their genomic characterization and highlighted differences between buffalo types and between breeds. A large ROH island on chromosome 2 was shared between river and swamp buffaloes and contained genes that are involved in environmental adaptation and reproduction. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00616-3.
Collapse
Affiliation(s)
| | - Licia Colli
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti-DIANA, Università Cattolica del Sacro Cuore, Piacenza, Italia.,Centro di Ricerca sulla Biodiversità e sul DNA Antico-BioDNA, Università Cattolica del Sacro Cuore, Piacenza, Italia
| | - Alberto Cesarani
- Dipartimento di Agraria, Università degli Studi di Sassari, Sassari, Italia. .,Department of Animal and Dairy Science, University of Georgia, Athens, GA, USA.
| | - Paolo Ajmone-Marsan
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti-DIANA, Università Cattolica del Sacro Cuore, Piacenza, Italia.,Centro di Ricerca Nutrigenomica e Proteomica-PRONUTRIGEN, Università Cattolica del Sacro Cuore, Piacenza, Italia
| | - Wai Y Low
- The Davies Research Centre, School of Animal and Veterinary Sciences, University of Adelaide, Roseworthy, SA, 5371, Australia
| | - Rick Tearle
- The Davies Research Centre, School of Animal and Veterinary Sciences, University of Adelaide, Roseworthy, SA, 5371, Australia
| | - John L Williams
- Dipartimento di Scienze Animali, della Nutrizione e degli Alimenti-DIANA, Università Cattolica del Sacro Cuore, Piacenza, Italia.,The Davies Research Centre, School of Animal and Veterinary Sciences, University of Adelaide, Roseworthy, SA, 5371, Australia
| |
Collapse
|
19
|
Laseca N, Anaya G, Peña Z, Pirosanto Y, Molina A, Demyda Peyrás S. Impaired Reproductive Function in Equines: From Genetics to Genomics. Animals (Basel) 2021; 11:393. [PMID: 33546520 PMCID: PMC7913728 DOI: 10.3390/ani11020393] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2021] [Revised: 01/29/2021] [Accepted: 01/31/2021] [Indexed: 12/15/2022] Open
Abstract
Fertility is one of the key factors in the economic and productive success of the equine industry. Despite this, studies on the genetic causes affecting reproductive performance are scarce, especially in mares, where the genetic architecture of the reproductive traits is extremely complex. Today, with the increasing availability of new genomic methodologies for this species, we are presented with an interesting opportunity to understand the genetic basis of equine reproductive disorders. These include, among others, novel techniques for detecting chromosomal abnormalities, whose association with infertility in horses was established over 50 years ago; new sequencing technologies permitting an accurate detection of point mutations influencing fertility, as well as the study of inbreeding and molecular homozygosity, which has been widely suggested as one of the main causes of low reproductive performance in horses. Finally, over the last few years, reproductive performance has also been associated with copy number variants and candidate genes detected by genome-wide association studies on fertility traits. However, such studies are still scarce, probably because they depend on the existence of large and accurate phenotypic datasets of reproductive and/or fertility traits, which are still difficult to obtain in equines.
Collapse
Affiliation(s)
- Nora Laseca
- Departamento de genética, Universidad de Córdoba, Campus de Rabanales Ctra, Madrid-Cádiz, km 396, 14071 Córdoba, Spain; (N.L.); (G.A.); (Z.P.); (A.M.)
| | - Gabriel Anaya
- Departamento de genética, Universidad de Córdoba, Campus de Rabanales Ctra, Madrid-Cádiz, km 396, 14071 Córdoba, Spain; (N.L.); (G.A.); (Z.P.); (A.M.)
| | - Zahira Peña
- Departamento de genética, Universidad de Córdoba, Campus de Rabanales Ctra, Madrid-Cádiz, km 396, 14071 Córdoba, Spain; (N.L.); (G.A.); (Z.P.); (A.M.)
| | - Yamila Pirosanto
- Departamento de Producción Animal, Facultad de Ciencias Veterinarias, Universidad Nacional de La Plata, La Plata 1900, Argentina;
- Consejo Superior de Investigaciones Científicas y Tecnológicas (CONICET), CCT-La Plata, La Plata 1900, Argentina
| | - Antonio Molina
- Departamento de genética, Universidad de Córdoba, Campus de Rabanales Ctra, Madrid-Cádiz, km 396, 14071 Córdoba, Spain; (N.L.); (G.A.); (Z.P.); (A.M.)
| | - Sebastián Demyda Peyrás
- Departamento de Producción Animal, Facultad de Ciencias Veterinarias, Universidad Nacional de La Plata, La Plata 1900, Argentina;
- Consejo Superior de Investigaciones Científicas y Tecnológicas (CONICET), CCT-La Plata, La Plata 1900, Argentina
| |
Collapse
|
20
|
Xia X, Zhang S, Zhang H, Zhang Z, Chen N, Li Z, Sun H, Liu X, Lyu S, Wang X, Li Z, Yang P, Xu J, Ding X, Shi Q, Wang E, Ru B, Xu Z, Lei C, Chen H, Huang Y. Assessing genomic diversity and signatures of selection in Jiaxian Red cattle using whole-genome sequencing data. BMC Genomics 2021; 22:43. [PMID: 33421990 PMCID: PMC7796570 DOI: 10.1186/s12864-020-07340-0] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Accepted: 12/22/2020] [Indexed: 01/20/2023] Open
Abstract
Background Native cattle breeds are an important source of genetic variation because they might carry alleles that enable them to adapt to local environment and tough feeding conditions. Jiaxian Red, a Chinese native cattle breed, is reported to have originated from crossbreeding between taurine and indicine cattle; their history as a draft and meat animal dates back at least 30 years. Using whole-genome sequencing (WGS) data of 30 animals from the core breeding farm, we investigated the genetic diversity, population structure and genomic regions under selection of Jiaxian Red cattle. Furthermore, we used 131 published genomes of world-wide cattle to characterize the genomic variation of Jiaxian Red cattle. Results The population structure analysis revealed that Jiaxian Red cattle harboured the ancestry with East Asian taurine (0.493), Chinese indicine (0.379), European taurine (0.095) and Indian indicine (0.033). Three methods (nucleotide diversity, linkage disequilibrium decay and runs of homozygosity) implied the relatively high genomic diversity in Jiaxian Red cattle. We used θπ, CLR, FST and XP-EHH methods to look for the candidate signatures of positive selection in Jiaxian Red cattle. A total number of 171 (θπ and CLR) and 17 (FST and XP-EHH) shared genes were identified using different detection strategies. Functional annotation analysis revealed that these genes are potentially responsible for growth and feed efficiency (CCSER1), meat quality traits (ROCK2, PPP1R12A, CYB5R4, EYA3, PHACTR1), fertility (RFX4, SRD5A2) and immune system response (SLAMF1, CD84 and SLAMF6). Conclusion We provide a comprehensive overview of sequence variations in Jiaxian Red cattle genomes. Selection signatures were detected in genomic regions that are possibly related to economically important traits in Jiaxian Red cattle. We observed a high level of genomic diversity and low inbreeding in Jiaxian Red cattle. These results provide a basis for further resource protection and breeding improvement of this breed. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-020-07340-0.
Collapse
Affiliation(s)
- Xiaoting Xia
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Shunjin Zhang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Huaju Zhang
- Pingdingshan animal husbandry technology promotion station, Pingdingshan, 467000, Henan, China
| | - Zijing Zhang
- Institute of Animal Husbandry and Veterinary Science, Henan Academy of Agricultural Sciences, Zhengzhou, 450002, Henan, China
| | - Ningbo Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Zhigang Li
- Pingdingshan animal husbandry technology promotion station, Pingdingshan, 467000, Henan, China
| | - Hongxia Sun
- Pingdingshan animal husbandry technology promotion station, Pingdingshan, 467000, Henan, China
| | - Xian Liu
- Henan Provincial Animal Husbandry General Station, Zhengzhou, 450008, Henan, China
| | - Shijie Lyu
- Institute of Animal Husbandry and Veterinary Science, Henan Academy of Agricultural Sciences, Zhengzhou, 450002, Henan, China
| | - Xianwei Wang
- Henan Provincial Animal Husbandry General Station, Zhengzhou, 450008, Henan, China
| | - Zhiming Li
- Henan Provincial Animal Husbandry General Station, Zhengzhou, 450008, Henan, China
| | - Peng Yang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Jiawei Xu
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Xiaoting Ding
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Qiaoting Shi
- Institute of Animal Husbandry and Veterinary Science, Henan Academy of Agricultural Sciences, Zhengzhou, 450002, Henan, China
| | - Eryao Wang
- Institute of Animal Husbandry and Veterinary Science, Henan Academy of Agricultural Sciences, Zhengzhou, 450002, Henan, China
| | - Baorui Ru
- Henan Provincial Animal Husbandry General Station, Zhengzhou, 450008, Henan, China
| | - Zejun Xu
- Henan Provincial Animal Husbandry General Station, Zhengzhou, 450008, Henan, China
| | - Chuzhao Lei
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Hong Chen
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China
| | - Yongzhen Huang
- Key Laboratory of Animal Genetics, Breeding and Reproduction of Shaanxi Province, College of Animal Science and Technology, Northwest A&F University, No. 22 Xinong Road, Yangling, 712100, Shaanxi, China.
| |
Collapse
|
21
|
Dzomba EF, Chimonyo M, Pierneef R, Muchadeyi FC. Runs of homozygosity analysis of South African sheep breeds from various production systems investigated using OvineSNP50k data. BMC Genomics 2021; 22:7. [PMID: 33407115 PMCID: PMC7788743 DOI: 10.1186/s12864-020-07314-2] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Accepted: 12/07/2020] [Indexed: 12/28/2022] Open
Abstract
Background Population history, production system and within-breed selection pressure impacts the genome architecture resulting in reduced genetic diversity and increased frequency of runs of homozygosity islands. This study tested the hypothesis that production systems geared towards specific traits of importance or natural or artificial selection pressures influenced the occurrence and distribution of runs of homozygosity (ROH) in the South African sheep population. The Illumina OvineSNP50 BeadChip was used to genotype 400 sheep belonging to 13 breeds from South Africa representing mutton, pelt and mutton and wool dual-purpose breeds, including indigenous non-descript breeds that are reared by smallholder farmers. To get more insight into the autozygosity and distribution of ROH islands of South African breeds relative to global populations, 623 genotypes of sheep from worldwide populations were included in the analysis. Runs of homozygosity were computed at cut-offs of 1–6 Mb, 6–12 Mb, 12–24 Mb, 24–48 Mb and > 48 Mb, using the R package detectRUNS. The Golden Helix SVS program was used to investigate the ROH islands. Results A total of 121,399 ROH with mean number of ROH per animal per breed ranging from 800 (African White Dorper) to 15,097 (Australian Poll Dorset) were obtained. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were in the short (1–6 Mb) category (88.2%). Most animals had no ROH > 48 Mb. Of the South African breeds, the Nguni and the Blackhead Persian displayed high ROH based inbreeding (FROH) of 0.31 ± 0.05 and 0.31 ± 0.04, respectively. Highest incidence of common runs per SNP across breeds was observed on chromosome 10 with over 250 incidences of common ROHs. Mean proportion of SNPs per breed per ROH island ranged from 0.02 ± 0.15 (island ROH224 on chromosome 23) to 0.13 ± 0.29 (island ROH175 on chromosome 15). Seventeen (17) of the islands had SNPs observed in single populations (unique ROH islands). The MacArthur Merino (MCM) population had five unique ROH islands followed by Blackhead Persian and Nguni with three each whilst the South African Mutton Merino, SA Merino, White Vital Swakara, Karakul, Dorset Horn and Chinese Merino each had one unique ROH island. Genes within ROH islands were associated with predominantly metabolic and immune response traits and predomestic selection for traits such as presence or absence of horns. Conclusions Overall, the frequency and patterns of distribution of ROH observed in this study corresponds to the breed history and implied selection pressures exposed to the sheep populations under study. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-020-07314-2.
Collapse
Affiliation(s)
- E F Dzomba
- Discipline of Genetics, School of Life Sciences, University of KwaZulu-Natal, Private Bag X01, Scottsville, 3209, South Africa.
| | - M Chimonyo
- Discipline of Animal & Poultry Science; School of Agricultural, Earth & Environmental Sciences, University of KwaZulu-Natal, Private Bag X01, Scottsville, 3209, South Africa
| | - R Pierneef
- Agricultural Research Council, Biotechnology Platform, Private Bag X5, Onderstepoort, 0110, South Africa
| | - F C Muchadeyi
- Agricultural Research Council, Biotechnology Platform, Private Bag X5, Onderstepoort, 0110, South Africa
| |
Collapse
|
22
|
Gorssen W, Meyermans R, Janssens S, Buys N. A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species. Genet Sel Evol 2021; 53:2. [PMID: 33397285 PMCID: PMC7784028 DOI: 10.1186/s12711-020-00599-7] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2020] [Accepted: 12/11/2020] [Indexed: 01/18/2023] Open
Abstract
BACKGROUND Runs of homozygosity (ROH) have become the state-of-the-art method for analysis of inbreeding in animal populations. Moreover, ROH are suited to detect signatures of selection via ROH islands and are used in other applications, such as genomic prediction and genome-wide association studies (GWAS). Currently, a vast amount of single nucleotide polymorphism (SNP) data is available online, but most of these data have never been used for ROH analysis. Therefore, we performed a ROH analysis on large medium-density SNP datasets in eight animal species (cat, cattle, dog, goat, horse, pig, sheep and water buffalo; 442 different populations) and make these results publicly available. RESULTS The results include an overview of ROH islands per population and a comparison of the incidence of these ROH islands among populations from the same species, which can assist researchers when studying other (livestock) populations or when looking for similar signatures of selection. We were able to confirm many known ROH islands, for example signatures of selection for the myostatin (MSTN) gene in sheep and horses. However, our results also included multiple other ROH islands, which are common to many populations and not identified to date (e.g. on chromosomes D4 and E2 in cats and on chromosome 6 in sheep). CONCLUSIONS We are confident that our repository of ROH islands is a valuable reference for future studies. The discovered ROH island regions represent a unique starting point for new studies or can be used as a reference for future studies. Furthermore, we encourage authors to add their population-specific ROH findings to our repository.
Collapse
Affiliation(s)
- Wim Gorssen
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium
| | - Roel Meyermans
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium
| | - Steven Janssens
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium
| | - Nadine Buys
- Livestock Genetics, Department of Biosystems, KU Leuven, Kasteelpark Arenberg 30, Box 2472, 3001, Leuven, Belgium.
| |
Collapse
|
23
|
Abstract
Runs of homozygosity (ROH) are chromosomal stretches that in a diploid genome appear in a homozygous state and display identical alleles at multiple contiguous loci. This study aimed to systematically compare the genomic distribution of the ROH islands among five populations of wild vs commercial chickens of both layer and broiler type. To this end, we analyzed whole genome sequences of 115 birds including white layer (WL, n = 25), brown layer (BL, n = 25), broiler line A (BRA, n = 20), broiler line B (BRB, n = 20) and Red Junglefowl (RJF, n = 25). The ROH segments varied in size markedly among populations, ranging from 0.3 to 21.83 Mb reflecting their past genealogy. White layers contained the largest portion of the genome in homozygous state with an average ROH length of 432.1 Mb (±18.7) per bird, despite carrying it in short segments (0.3-1 Mb). Population-wise inbreeding measures based on Wright's (Fis) and genomic (FROH) metrics revealed highly inbred genome of layer lines relative to the broilers and Red Junglefowl. We further revealed the ROH islands, among commercial lines overlapped with QTL related to limb development (GREM1, MEOX2), body weight (Meis2a.1, uc_338), eggshell color (GLCCI1, ICA1, UMAD1), antibody response to Newcastle virus (ROBO2), and feather pecking. Comparison of ROH landscape in sequencing resolution demonstrated that a sizable portion of genome of commercial lines segregates in homozygote state, reflecting many generations of assortative mating and intensive selection in their recent history. In contrary, wild birds carry shorter ROH segments, likely suggestive of older evolutionary events.
Collapse
|
24
|
Schiavo G, Bovo S, Bertolini F, Dall'Olio S, Nanni Costa L, Tinarelli S, Gallo M, Fontanesi L. Runs of homozygosity islands in Italian cosmopolitan and autochthonous pig breeds identify selection signatures in the porcine genome. Livest Sci 2020. [DOI: 10.1016/j.livsci.2020.104219] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
|
25
|
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
|
26
|
Himmelbauer J, Mészáros G, Sölkner J. Detection of Autosomal Hemizygous Regions in the Fleckvieh Population Based on SNP-chip Data and Parent Offspring Pairs. ACTA UNIVERSITATIS AGRICULTURAE ET SILVICULTURAE MENDELIANAE BRUNENSIS 2019. [DOI: 10.11118/actaun201967061447] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
|
27
|
Luigi-Sierra MG, Cardoso TF, Martínez A, Pons A, Bermejo LA, Jordana J, Delgado JV, Adán S, Ugarte E, Arranz JJ, Casellas J, Amills M. Low genome-wide homozygosity in 11 Spanish ovine breeds. Anim Genet 2019; 50:501-511. [PMID: 31393638 DOI: 10.1111/age.12832] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 06/04/2019] [Indexed: 12/16/2022]
Abstract
The population of Spanish sheep has decreased from 24 to 15 million heads in the last 75 years due to multiple social and economic factors. Such a demographic reduction might have caused an increase in homozygosity and inbreeding, thus limiting the viability of local breeds with excellent adaptations to harsh ecosystems. The main goal of our study was to investigate the homozygosity patterns of 11 Spanish ovine breeds and to elucidate the relationship of these Spanish breeds with reference populations from Europe, Africa and the Near East. By using Ovine SNP50 BeadChip data retrieved from previous publications, we have found that the majority of studied Spanish ovine breeds have close genetic relatedness with other European populations; the one exception is the Canaria de Pelo breed, which is similar to North African breeds. Our analysis has also demonstrated that, with few exceptions, the genomes of Spanish sheep harbor fewer than 50 runs of homozygosity (ROH) with a total length of less than 350 Mb. Moreover, the frequencies of very long ROH (>30 Mb) are very low, and the inbreeding coefficients (FROH ) are generally small (FROH < 0.10), ranging from 0.008 (Rasa Aragonesa) to 0.086 (Canaria de Pelo). The low levels of homozygosity observed in the 11 Spanish sheep under analysis might be due to their extensive management and the high number of small to medium farms.
Collapse
Affiliation(s)
- M G Luigi-Sierra
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| | - T F Cardoso
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain.,CAPES Foundation, Ministry of Education of Brazil, Brasilia D.F, 70.040-020, Brazil
| | - A Martínez
- Departamento de Genética, Universidad de Córdoba, Córdoba, 14071, Spain
| | - A Pons
- Unitat de Races Autòctones, Servei de Millora Agrària i Pesquera (SEMILLA), Son Ferriol, 07198, Spain
| | - L A Bermejo
- Departamento de Ingeniería, Producción y Economía Agrarias, Universidad de La Laguna, La Laguna, Tenerife, 38071, Spain
| | - J Jordana
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| | - J V Delgado
- Departamento de Genética, Universidad de Córdoba, Córdoba, 14071, Spain
| | - S Adán
- Pazo de Fontefiz, Federación de Razas Autóctonas de Galicia (BOAGA), 32152, Coles, Ourense, Spain
| | - E Ugarte
- Neiker-Tecnalia, Campus Agroalimentario de Arkaute, apdo 46 E-01080, Vitoria-Gazteiz (Araba), Spain
| | - J J Arranz
- Departamento de Producción Animal, Universidad de León, León, 24071, Spain
| | - J Casellas
- Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| | - M Amills
- Department of Animal Genetics, Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain.,Departament de Ciència Animal i dels Aliments, Facultat de Veterinària, Universitat Autònoma de Barcelona, Bellaterra, 08193, Spain
| |
Collapse
|
28
|
Nandolo W, Mészáros G, Banda LJ, Gondwe TN, Lamuno D, Mulindwa HA, Nakimbugwe HN, Wurzinger M, Utsunomiya YT, Woodward-Greene MJ, Liu M, Liu G, Van Tassell CP, Curik I, Rosen BD, Sölkner J. Timing and Extent of Inbreeding in African Goats. Front Genet 2019; 10:537. [PMID: 31214253 PMCID: PMC6558083 DOI: 10.3389/fgene.2019.00537] [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: 10/31/2018] [Accepted: 05/17/2019] [Indexed: 11/13/2022] Open
Abstract
Genetic characterization of African goats is one of the current priorities in the improvement of goats in the continent. This study contributes to the characterization effort by determining the levels and number of generations to common ancestors ("age") associated with inbreeding in African goat breeds and identifies regions that contain copy number variation mistyped as being homozygous. Illumina 50k single nucleotide polymorphism genotype data for 608 goats from 31 breeds were used to compute the level and age of inbreeding at both local (marker) and global levels (FG) using a model-based approach based on a hidden Markov model. Runs of homozygosity (ROH) segments detected using the Viterbi algorithm led to ROH-based inbreeding coefficients for all ROH (FROH) and for ROH longer than 2 Mb (FROH > 2Mb). Some of the genomic regions identified as having ROH are likely to be hemizygous regions (copy number deletions) mistyped as homozygous regions. Although the proportion of these miscalled ROH is small and does not substantially affect estimates of levels of inbreeding for individual animals, the inbreeding metrics were adjusted by removing these regions from the ROH. All the inbreeding metrics varied widely across breeds, with overall means of 0.0408, 0.0370, and 0.0691 and medians of 0.0125, 0.0098, and 0.0366 for FROH, FROH > 2Mb, and FG, respectively. Several breeds (including Menabe and Sofia from Madagascar) had high proportions of recent inbreeding, while Small East African, Ethiopian, and most of the West African breeds (including West African Dwarf) had more ancient inbreeding.
Collapse
Affiliation(s)
- Wilson Nandolo
- Division of Livestock Sciences, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences, Vienna, Vienna, Austria.,Department of Animal Science, Faculty of Agriculture, Lilongwe University of Agriculture and Natural Resources, Lilongwe, Malawi.,Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agriculture Research Service, Beltsville, MD, United States
| | - Gábor Mészáros
- Division of Livestock Sciences, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences, Vienna, Vienna, Austria
| | - Liveness Jessica Banda
- Department of Animal Science, Faculty of Agriculture, Lilongwe University of Agriculture and Natural Resources, Lilongwe, Malawi
| | - Timothy N Gondwe
- Department of Animal Science, Faculty of Agriculture, Lilongwe University of Agriculture and Natural Resources, Lilongwe, Malawi
| | - Doreen Lamuno
- Division of Livestock Sciences, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences, Vienna, Vienna, Austria
| | | | | | - Maria Wurzinger
- Division of Livestock Sciences, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences, Vienna, Vienna, Austria
| | - Yuri T Utsunomiya
- School of Agricultural and Veterinarian Sciences, Jaboticabal, São Paulo State University (UNESP), São Paulo, Brazil
| | - M Jennifer Woodward-Greene
- Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agriculture Research Service, Beltsville, MD, United States
| | - Mei Liu
- Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agriculture Research Service, Beltsville, MD, United States
| | - George Liu
- Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agriculture Research Service, Beltsville, MD, United States
| | - Curtis P Van Tassell
- Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agriculture Research Service, Beltsville, MD, United States
| | - Ino Curik
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Zagreb, Croatia
| | - Benjamin D Rosen
- Animal Genomics and Improvement Laboratory, United States Department of Agriculture, Agriculture Research Service, Beltsville, MD, United States
| | - Johann Sölkner
- Division of Livestock Sciences, Department of Sustainable Agricultural Systems, University of Natural Resources and Life Sciences, Vienna, Vienna, Austria
| |
Collapse
|
29
|
Runs of Homozygosity and Population History of Three Horse Breeds With Small Population Size. J Equine Vet Sci 2018. [DOI: 10.1016/j.jevs.2018.09.004] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
|