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; 55:575-587. [PMID: 38806279 DOI: 10.1111/age.13437] [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: 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
|
Bem RD, Benfica LF, Silva DA, Carrara ER, Brito LF, Mulim HA, Borges MS, Cyrillo JNSG, Canesin RC, Bonilha SFM, Mercadante MEZ. Assessing different metrics of pedigree and genomic inbreeding and inbreeding effect on growth, fertility, and feed efficiency traits in a closed-herd Nellore cattle population. BMC Genomics 2024; 25:738. [PMID: 39080557 PMCID: PMC11290228 DOI: 10.1186/s12864-024-10641-3] [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: 04/05/2024] [Accepted: 07/19/2024] [Indexed: 08/02/2024] Open
Abstract
BACKGROUND The selection of individuals based on their predicted breeding values and mating of related individuals can increase the proportion of identical-by-descent alleles. In this context, the objectives of this study were to estimate inbreeding coefficients based on alternative metrics and data sources such as pedigree (FPED), hybrid genomic relationship matrix H (FH), and ROH of different length (FROH); and calculate Pearson correlations between the different metrics in a closed Nellore cattle population selected for body weight adjusted to 378 days of age (W378). In addition to total FROH (all classes) coefficients were also estimated based on the size class of the ROH segments: FROH1 (1-2 Mb), FROH2 (2-4 Mb), FROH3 (4-8 Mb), FROH4 (8-16 Mb), and FROH5 (> 16 Mb), and for each chromosome (FROH_CHR). Furthermore, we assessed the effect of each inbreeding metric on birth weight (BW), body weights adjusted to 210 (W210) and W378, scrotal circumference (SC), and residual feed intake (RFI). We also evaluated the chromosome-specific effects of inbreeding on growth traits. RESULTS The correlation between FPED and FROH was 0.60 while between FH and FROH and FH and FPED were 0.69 and 0.61, respectively. The annual rate of inbreeding was 0.16% for FPED, 0.02% for FH, and 0.16% for FROH. A 1% increase in FROH5 resulted in a reduction of up to -1.327 ± 0.495 kg in W210 and W378. Four inbreeding coefficients (FPED, FH, FROH2, and FROH5) had a significant effect on W378, with reductions of up to -3.810 ± 1.753 kg per 1% increase in FROH2. There was an unfavorable effect of FPED on RFI (0.01 ± 0.0002 kg dry matter/day) and of FROH on SC (-0.056 ± 0.022 cm). The FROH_CHR coefficients calculated for BTA3, BTA5, and BTA8 significantly affected the growth traits. CONCLUSIONS Inbreeding depression was observed for all traits evaluated. However, these effects were greater for the criterion used for selection of the animals (i.e., W378). The increase in the genomic inbreeding was associated with a higher inbreeding depression on the traits evaluated when compared to pedigree-based inbreeding. Genomic information should be used as a tool during mating to optimize control of inbreeding and, consequently, minimize inbreeding depression in Nellore cattle.
Collapse
Affiliation(s)
- Ricardo D Bem
- Institute of Animal Science, Sertãozinho, SP, Brazil.
- Department of Animal Science, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil.
| | - Lorena F Benfica
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907, USA
| | - Delvan A Silva
- Department of Animal Sciences, Federal University of Viçosa, Viçosa, MG, Brazil
| | - Eula R Carrara
- Department of Animal Sciences, Federal University of Viçosa, Viçosa, MG, Brazil
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907, USA
| | - Henrique A Mulim
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907, USA
| | - Marcelo S Borges
- Department of Pathology, Reproduction and One Health, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil
| | | | | | | | | |
Collapse
|
3
|
Xiao Z, Li Y, Yang L, Cui M, Wang Z, Sun W, Wang J, Chen S, Lai S, Jia X. Genome-Wide Association Studies of Growth Trait Heterosis in Crossbred Meat Rabbits. Animals (Basel) 2024; 14:2096. [PMID: 39061558 PMCID: PMC11273820 DOI: 10.3390/ani14142096] [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: 06/07/2024] [Revised: 07/07/2024] [Accepted: 07/12/2024] [Indexed: 07/28/2024] Open
Abstract
The application of heterosis can not only effectively improve the disease resistance and meat quality of livestock, but also significantly enhance the reproduction and growth of livestock and poultry. We conducted genome-wide association studies using data from F2 crossbred meat rabbits to screen out candidate genes with significant dominant effects associated with economic trait variation. High-throughput sequencing technology was used to obtain SNPs covering the whole genome to evaluate the homozygosity of the population genome, and analyze the number, length, frequency, and distribution of ROHs in the population. Candidate genes related to economic traits of meat rabbits were searched based on high-frequency ROH regions. After quality control filtering of genotype data, 380 F2 crossbred rabbits were identified with 78,579 SNPs and 42,018 ROHs on the autosomes. The fitting of the Logistic growth curve model showed that 49-day-old rabbits were a growth inflection point. Then, through genome-wide association studies, 10 SNP loci and seven growth trait candidate genes were found to be significantly related to body weight in meat rabbits at 84 days of age. In addition, we revealed the functional roles and locations of 20 candidate genes in the high-frequency ROH region associated with economic traits in meat rabbits. This study identified potential genes associated with growth and development in the high-frequency ROH region of meat rabbits. In this study, the identified candidate genes can be used as molecular markers for assisted selection in meat rabbits. At the same time, the inbreeding situation based on ROH assessment can provide reference for breeding and breeding preservation of meat rabbits.
Collapse
Affiliation(s)
- Zhanjun Xiao
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China; (Z.X.); (M.C.); (Z.W.); (W.S.); (S.C.); (S.L.)
| | - Yuchao Li
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 610000, China; (Y.L.); (L.Y.); (J.W.)
| | - Li Yang
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 610000, China; (Y.L.); (L.Y.); (J.W.)
| | - Mingyan Cui
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China; (Z.X.); (M.C.); (Z.W.); (W.S.); (S.C.); (S.L.)
| | - Zicheng Wang
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China; (Z.X.); (M.C.); (Z.W.); (W.S.); (S.C.); (S.L.)
| | - Wenqiang Sun
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China; (Z.X.); (M.C.); (Z.W.); (W.S.); (S.C.); (S.L.)
| | - Jie Wang
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 610000, China; (Y.L.); (L.Y.); (J.W.)
| | - Shiyi Chen
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China; (Z.X.); (M.C.); (Z.W.); (W.S.); (S.C.); (S.L.)
| | - Songjia Lai
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China; (Z.X.); (M.C.); (Z.W.); (W.S.); (S.C.); (S.L.)
| | - Xianbo Jia
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu 611130, China; (Z.X.); (M.C.); (Z.W.); (W.S.); (S.C.); (S.L.)
| |
Collapse
|
4
|
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
|
5
|
Cars BS, Kessler C, Hoffman EA, Côté SD, Koelsch D, Shafer ABA. Island demographics and trait associations in white-tailed deer. Heredity (Edinb) 2024; 133:1-10. [PMID: 38802598 PMCID: PMC11222433 DOI: 10.1038/s41437-024-00685-2] [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: 08/03/2023] [Revised: 04/22/2024] [Accepted: 04/23/2024] [Indexed: 05/29/2024] Open
Abstract
When a population is isolated and composed of few individuals, genetic drift is the paramount evolutionary force and results in the loss of genetic diversity. Inbreeding might also occur, resulting in genomic regions that are identical by descent, manifesting as runs of homozygosity (ROHs) and the expression of recessive traits. Likewise, the genes underlying traits of interest can be revealed by comparing fixed SNPs and divergent haplotypes between affected and unaffected individuals. Populations of white-tailed deer (Odocoileus virginianus) on islands of Saint Pierre and Miquelon (SPM, France) have high incidences of leucism and malocclusions, both considered genetic defects; on the Florida Keys islands (USA) deer exhibit smaller body sizes, a polygenic trait. Here we aimed to reconstruct island demography and identify the genes associated with these traits in a pseudo case-control design. The two island populations showed reduced levels of genomic diversity and a build-up of deleterious mutations compared to mainland deer; there was also significant genome-wide divergence in Key deer. Key deer showed higher inbreeding levels, but not longer ROHs, consistent with long-term isolation. We identified multiple trait-related genes in ROHs including LAMTOR2 which has links to pigmentation changes, and NPVF which is linked to craniofacial abnormalities. Our mixed approach of linking ROHs, fixed SNPs and haplotypes matched a high number (~50) of a-priori body size candidate genes in Key deer. This suite of biomarkers and candidate genes should prove useful for population monitoring, noting all three phenotypes show patterns consistent with a complex trait and non-Mendelian inheritance.
Collapse
Affiliation(s)
- Brooklyn S Cars
- Environmental and Life Sciences Graduate Program, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada
- Department of Forensics, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada
| | - Camille Kessler
- Environmental and Life Sciences Graduate Program, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada
| | - Eric A Hoffman
- Department of Biology, University of Central Florida, 4000, Central Florida Blvd, Orlando, FL, USA
| | - Steeve D Côté
- Département de Biologie and Centre d'Études Nordiques, Université Laval, Québec, QC, G1V 0A6, Canada
| | - Daniel Koelsch
- Fédération des chasseurs de Saint-Pierre et Miquelon, Saint-Pierre et Miquelon, France
- Direction des Territoires de l'Alimentation et de la Mer, service Biodiversité, Saint-Pierre et Miquelon, France
| | - Aaron B A Shafer
- Environmental and Life Sciences Graduate Program, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada.
- Department of Forensics, Trent University, 2140 East Bank Drive, Peterborough, ON, K9J 7B8, Canada.
| |
Collapse
|
6
|
Sievers J, Distl O. The Estimation of Genetic Parameters for Chronic Progressive Lymphedema and Body Traits in the Rhenish German Draught Horse. Animals (Basel) 2024; 14:1214. [PMID: 38672362 PMCID: PMC11047459 DOI: 10.3390/ani14081214] [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: 02/23/2024] [Revised: 04/14/2024] [Accepted: 04/16/2024] [Indexed: 04/28/2024] Open
Abstract
Chronic progressive lymphedema (CPL) is a prevalent and progressive disease in Rhenish German draught horses. The objective of our follow-up study was to evaluate the heritability of this disease in Rhenish German draught horses using pedigree-based and genomic relationship matrices. We employed linear and threshold animal models. Models included the random animal effect and effects of breeding association, coat colour, sex, and age within sex, and farm-related factors, on CPL scores. In addition, we estimated heritabilities in models assuming censoring for data when horses were below an age of 1-15 years. The heritabilities of CPL scores across all ages were 0.595 ± 0.131 and 0.482 ± 0.105 in the threshold and linear animal model with pedigree-based relationship matrices, respectively. The restriction of data to horses with a minimum age at examination or accounting for censored data in younger animals showed an increase in heritabilities of CPL scores up to 0.788 ± 0.168 (threshold model) and 0.752 ± 0.153 (linear model) at an age of 7-8 years. Analyses including genomic relationship matrices yielded very similar estimates, but with smaller standard errors than pedigree-based analyses. Heritabilities in threshold models for CPL prevalence (CPL-bin-score) and the number of affected limbs (CPL-bin-sum) were 0.176-0.189 ± 0.061-0.064 and 0.375-0.433 ± 0.164-0.170, respectively. We were able to show moderately to highly positive genetic correlations between the CPL score and cannon bone circumference (0.529-0.825), height at withers (0.338-0.555), and skinfold thickness (0.241-0.517). Using the dichotomous trait for the CPL score and the genomic relationship matrix resulted in corresponding estimates of 0.868, 0.793, and 0.784, respectively. This study showed the great importance of additive genetic variation influencing the expression of chronic progressive lymphedema in Rhenish German draught horses. Therefore, further research is warranted to implement breeding programmes in a small breeding population that exploit the potential of additive genetic differences among animals for reducing the prevalence and severity of lesions of this incurable disease.
Collapse
Affiliation(s)
| | - Ottmar Distl
- Institute of Animal Breeding and Genetics, University of Veterinary Medicine Hannover (Foundation), 30559 Hannover, Germany;
| |
Collapse
|
7
|
Ojeda-Marín C, Gutiérrez JP, Formoso-Rafferty N, Goyache F, Cervantes I. Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight. J Anim Breed Genet 2024; 141:193-206. [PMID: 37990938 DOI: 10.1111/jbg.12835] [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/24/2023] [Revised: 09/11/2023] [Accepted: 11/09/2023] [Indexed: 11/23/2023]
Abstract
Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.
Collapse
Affiliation(s)
- Candela Ojeda-Marín
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
| | - Juan Pablo Gutiérrez
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
| | | | - Félix Goyache
- Departamento de Producción Agraria, E.T.S. Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Madrid, Spain
| | - Isabel Cervantes
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
| |
Collapse
|
8
|
Kar D, Ganguly I, Singh S, Bhatia AK, Dixit SP. Genome-wide runs of homozygosity signatures in diverse Indian goat breeds. 3 Biotech 2024; 14:81. [PMID: 38375512 PMCID: PMC10874352 DOI: 10.1007/s13205-024-03921-y] [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/22/2023] [Accepted: 01/05/2024] [Indexed: 02/21/2024] Open
Abstract
The present study analyzed ROH and consensus ROH regions in 102 animals of eleven diverse Indian goat (Capra hircus) breeds using whole genome sequencing. A total of 51,705 ROH and 21,271 consensus regions were identified. The mean number of ROH per animal was highest in the meat breed, Jharkhand Black (2693) and lowest in the pashmina breed, Changthangi (60). The average length of ROH (ALROH) was maximum in Kanniadu (974.11 Kb) and minimum in Tellicherry (146.98 Kb). Long ROH is typically associated with more recent inbreeding, whereas short ROH is connected to more ancient inbreeding. The overall ROH-based genomic inbreeding (FROH) was highest for Jharkhand Black (0.602) followed by Kanniadu (0.120) and Sangamneri (0.108) among all breeds. FROH of Jharkhand Black was higher than Kanniadu up to 5 Mb ROH length category. However, in > 20 Mb ROH length category, Kanniadu (0.98) exhibited significantly higher FROH than Jharkhand Black (0.46). This implies that Kanniadu had higher levels of recent inbreeding than Jharkhand Black. Despite this, due to the presence of both recent and ancient inbreeding, Jharkhand Black demonstrated higher overall FROH compared to Kanniadu. ROH patterns revealed dual purpose (meat and dairy) and pashmina breeds as less consanguineous while recent inbreeding was apparent in meat breeds. Analysis of ROH consensus regions identified selection sweeps in key genes governing intramuscular fat deposition, meat tenderisation, lean meat production and carcass weight (CDK4, ALOX15, CASP9, PRDM16, DVL1) in meat breeds; milk fat percentage and mammary gland development (POLD1, NOTCH2, ARHGAP35) in dual purpose (meat and dairy) breeds; while cold adaptation and hair follicle development (APOBEC1, DNAJC3, F2RL1, FGF9) in pashmina breed. MAPK, RAS, BMP and Wnt signaling pathways associated with hair follicle morphogenesis in Changthangi were also identified. PCA analysis based on ROH consensus regions revealed that meat breeds are more diverse than other goat breeds/populations. Supplementary Information The online version contains supplementary material available at 10.1007/s13205-024-03921-y.
Collapse
Affiliation(s)
- Dibyasha Kar
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
- Division of Animal Genetics and Breeding, ICAR-National Dairy Research Institute, Karnal, Haryana 132001 India
| | - Indrajit Ganguly
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| | - Sanjeev Singh
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| | - Avnish Kumar Bhatia
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| | - S. P. Dixit
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal, Haryana 132001 India
| |
Collapse
|
9
|
Sarlo Davila KM, Boggiatto P, Olsen S, Lippolis JD, Crooker BA, Putz EJ. Effect of selection genotype on immune response to Brucella abortus RB51 in Holstein cattle. Anim Genet 2024; 55:47-54. [PMID: 37946616 DOI: 10.1111/age.13372] [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: 07/10/2023] [Revised: 09/11/2023] [Accepted: 10/16/2023] [Indexed: 11/12/2023]
Abstract
Genetic selection for milk production traits in US Holsteins has affected numerous genes associated with reproduction and immunity. This study compares the transcriptomic response of peripheral blood mononuclear cells to an in vitro Brucella abortus strain RB51 (RB51) bacterial challenge between contemporary Holsteins and Holsteins that have not been selected for milk production traits since the mid-1960s. Total RNA was extracted from peripheral blood mononuclear cells from four contemporary and four unselected lactating, primiparous cows following 24-h incubation with or without stimulation with RB51 bacteria. RNA was sequenced and reads analyzed using tools from galaxy.scinet.usda.gov. A total of 412 differentially expressed genes (false discovery rate p < 0.05, log fold change > |1|) were identified. The upregulated genes (genes with higher expression in contemporary than unselected cattle) were enriched for 19 terms/pathways, including alanine, aspartate, and glutamate metabolism, indicating a cellular stress response. Downregulated genes (genes with higher expression in unselected than contemporary cows) were enriched for 37 terms/pathways, representing diverse immune responses, including natural killer cell-mediated immunity, interferon-γ production, negative regulation of interleukin-10 production, and cytokine receptor activity indicating a broad immune response with an emphasis on immune defense. These results provide evidence that differences exist between the two genotypes in response to in vitro bacterial challenge. This suggests that contemporary cows, genetically selected for milk production, may have reduced immune function, including limitations in response to intracellular bacteria.
Collapse
Affiliation(s)
- Kaitlyn M Sarlo Davila
- Infectious Bacterial Diseases Research Unit, National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa, USA
- Ruminant Diseases and Immunology Research Unit, National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa, USA
| | - Paola Boggiatto
- Infectious Bacterial Diseases Research Unit, National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa, USA
| | - Steven Olsen
- Infectious Bacterial Diseases Research Unit, National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa, USA
| | - John D Lippolis
- Ruminant Diseases and Immunology Research Unit, National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa, USA
| | - Brian A Crooker
- Department of Animal Science, University of Minnesota, St. Paul, Minnesota, USA
| | - Ellie J Putz
- Infectious Bacterial Diseases Research Unit, National Animal Disease Center, Agricultural Research Service, United States Department of Agriculture, Ames, Iowa, USA
| |
Collapse
|
10
|
Illa SK, Mumtaz S, Nath S, Mukherjee S, Mukherjee A. Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle. J Appl Genet 2024; 65:167-180. [PMID: 38110827 DOI: 10.1007/s13353-023-00816-1] [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: 02/10/2023] [Revised: 11/29/2023] [Accepted: 12/05/2023] [Indexed: 12/20/2023]
Abstract
Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.
Collapse
Affiliation(s)
- Satish Kumar Illa
- Livestock Research Station, Garividi, Sri Venkateswara Veterinary University, Tirupati, Andhra Pradesh State, India
| | - Shabahat Mumtaz
- Animal Husbandry Department, Kolkata, West Bengal State, India
| | - Sapna Nath
- College of Veterinary Science, Garividi, Sri Venkateswara Veterinary University, Tirupati, Andhra Pradesh State, India
| | - Sabyasachi Mukherjee
- Animal Genetics & Breeding Division, Indian Council of Agricultural Research (ICAR)-National Dairy Research Institute (NDRI), Karnal, Haryana State, India.
| | - Anupama Mukherjee
- Animal Genetics & Breeding Division, Indian Council of Agricultural Research (ICAR)-National Dairy Research Institute (NDRI), Karnal, Haryana State, India.
| |
Collapse
|
11
|
Falchi L, Cesarani A, Criscione A, Hidalgo J, Garcia A, Mastrangelo S, Macciotta NPP. Effect of genotyping density on the detection of runs of homozygosity and heterozygosity in cattle. J Anim Sci 2024; 102:skae147. [PMID: 38798158 PMCID: PMC11197001 DOI: 10.1093/jas/skae147] [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/24/2024] [Accepted: 05/24/2024] [Indexed: 05/29/2024] Open
Abstract
Runs of homozygosity (ROHom) are contiguous stretches of homozygous regions of the genome. In contrast, runs of heterozygosity (ROHet) are heterozygosity-rich regions. The detection of these two types of genomic regions (ROHom and ROHet) is influenced by the parameters involved in their identification and the number of available single-nucleotide polymorphisms (SNPs). The present study aimed to test the effect of chip density in detecting ROHom and ROHet in the Italian Simmental cattle breed. A sample of 897 animals were genotyped at low density (50k SNP; 397 individuals), medium density (140k SNP; 348 individuals), or high density (800k SNP; 152 individuals). The number of ROHom and ROHet per animal (nROHom and nROHet, respectively) and their average length were calculated. ROHom or ROHet shared by more than one animal and the number of times a particular SNP was inside a run were also computed (SNPROHom and SNPROHet). As the chip density increased, the nROHom increased, whereas their average length decreased. In contrast, the nROHet decreased and the average length increased as the chip density increased. The most repeated ROHom harbored no genes, whereas in the most repeated ROHet four genes (SNRPN, SNURF, UBE3A, and ATP10A) previously associated with reproductive traits were found. Across the 3 datasets, 31 SNP, located on Bos taurus autosome (BTA) 6, and 37 SNP (located on BTA21) exceeded the 99th percentile in the distribution of the SNPROHom and SNPROHet, respectively. The genomic region on BTA6 mapped the SLIT2, PACRGL, and KCNIP4 genes, whereas 19 and 18 genes were mapped on BTA16 and BTA21, respectively. Interestingly, most of genes found through the ROHet analysis were previously reported to be related to health, reproduction, and fitness traits. The results of the present study confirm that the detection of ROHom is more reliable when the chip density increases, whereas the ROHet trend seems to be the opposite. Genes and quantitative trait loci (QTL) mapped in the highlighted regions confirm that ROHet can be due to balancing selection, thus related to fitness traits, health, and reproduction, whereas ROHom are mainly involved in production traits. The results of the present study strengthened the usefulness of these parameters in analyzing the genomes of livestock and their biological meaning.
Collapse
Affiliation(s)
- Laura Falchi
- Dipartimento di Agraria, Università degli Studi di Sassari, Sassari 07100, Italy
| | - Alberto Cesarani
- Dipartimento di Agraria, Università degli Studi di Sassari, Sassari 07100, Italy
- Department of Animal and Dairy Science, University of Georgia, Athens 30602, USA
| | - Andrea Criscione
- Dipartimento di Agricoltura, Alimentazione e Ambiente, Università degli Studi di Catania, Catania 95123, Italy
| | - Jorge Hidalgo
- Department of Animal and Dairy Science, University of Georgia, Athens 30602, USA
| | - Andre Garcia
- American Angus Association, Angus Genetics Inc., Saint Joseph, MO, USA
| | - Salvatore Mastrangelo
- Dipartimento di Scienze Agrarie, Alimentari, e Forestali, Università degli Studi di Palermo, Palermo 90128, Italy
| | | |
Collapse
|
12
|
Smaragdov MG. Identification of homozygosity-rich regions in the Holstein genome. Vavilovskii Zhurnal Genet Selektsii 2023; 27:471-479. [PMID: 37808215 PMCID: PMC10556852 DOI: 10.18699/vjgb-23-57] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2022] [Revised: 02/03/2023] [Accepted: 02/27/2023] [Indexed: 10/10/2023] Open
Abstract
In this study, 371 Holstein cows from six herds and 26 Holstein bulls, which were used in these herds, were genotyped by the Illumina BovineSNP50 array. For runs of homozygosity (ROH) identification, consecutive and sliding runs were performed by the detectRUNS and Plink software. The missing calls did not significantly affect the ROH data. The mean number of ROH identified by consecutive runs was 95.4 ± 2.7, and that by sliding runs was 86.0 ± 2.6 in cows, while this number for Holstein bulls was lower 58.9 ± 1.9. The length of the ROH segments varied from 1 Mb to over 16 Mb, with the largest number of ROH having a length of 1-2 Mb. Of the 29 chromosomes, BTA 14, BTA 16, and BTA 7 were the most covered by ROH. The mean coefficient of inbreeding across the herds was 0.111 ± 0.003 and 0.104 ± 0.004 based on consecutive and sliding runs, respectively, and 0.078 ± 0.005 for bulls based on consecutive runs. These values do not exceed those for Holstein cattle in North America. The results of this study confirmed the more accurate identification of ROH by consecutive runs, and also that the number of allowed heterozygous SNPs may have a significant effect on ROH data.
Collapse
Affiliation(s)
- M G Smaragdov
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Science Center for Animal Husbandry, St. Petersburg, Pushkin, Russia
| |
Collapse
|
13
|
Pacheco HA, Rossoni A, Cecchinato A, Peñagaricano F. Identification of runs of homozygosity associated with male fertility in Italian Brown Swiss cattle. Front Genet 2023; 14:1227310. [PMID: 37485336 PMCID: PMC10356982 DOI: 10.3389/fgene.2023.1227310] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2023] [Accepted: 06/23/2023] [Indexed: 07/25/2023] Open
Abstract
Intensive selection for improved productivity has been accompanied by an increase in inbreeding rates and a reduction in genetic diversity. The increase in inbreeding tends to impact performance, especially fitness-related traits such as male fertility. Inbreeding can be monitored using runs of homozygosity (ROH), defined as contiguous lengths of homozygous genotypes observed in an individual's chromosome. The goal of this study was to evaluate the presence of ROH in Italian Brown Swiss cattle and assess its association with bull fertility. First, we evaluated the association between ROH and male fertility using 1,102 Italian Brown Swiss bulls with sire conception rate records and 572 K SNPs spanning the entire genome. Second, we split the entire population into 100 high-fertility and 100 low-fertility bulls to investigate the potential enrichment of ROH segments in the low-fertility group. Finally, we mapped the significant ROH regions to the bovine genome to identify candidate genes associated with sperm biology and male fertility. Notably, there was a negative association between bull fertility and the amount of homozygosity. Four different ROH regions located in chromosomes 6, 10, 11, and 24 were significantly overrepresented in low-fertility bulls (Fisher's exact test, p-value <0.01). Remarkably, these four genomic regions harbor many genes such as WDR19, RPL9, LIAS, UBE2K, DPF3, 5S-rRNA, 7SK, U6, and WDR7 that are related to sperm biology and male fertility. Overall, our findings suggest that inbreeding and increased homozygosity have a negative impact on male fertility in Italian Brown Swiss cattle. The quantification of ROH can contribute to minimizing the inbreeding rate and avoid its negative effect on fitness-related traits, such as male fertility.
Collapse
Affiliation(s)
- Hendyel A. Pacheco
- Department of Animal and Dairy Sciences, University of Wisconsin-Madison, Madison, WI, United States
| | | | - Alessio Cecchinato
- Department of Agronomy, Food, Natural Resources, Animals and Environment, University of Padova, Padua, Italy
| | - Francisco Peñagaricano
- Department of Animal and Dairy Sciences, University of Wisconsin-Madison, Madison, WI, United States
| |
Collapse
|
14
|
Tian S, Tang W, Zhong Z, Wang Z, Xie X, Liu H, Chen F, Liu J, Han Y, Qin Y, Tan Z, Xiao Q. Identification of Runs of Homozygosity Islands and Functional Variants in Wenchang Chicken. Animals (Basel) 2023; 13:ani13101645. [PMID: 37238076 DOI: 10.3390/ani13101645] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/14/2023] [Revised: 05/11/2023] [Accepted: 05/13/2023] [Indexed: 05/28/2023] Open
Abstract
Wenchang chickens, a native breed in the Hainan province of China, are famous for their meat quality and adaptability to tropical conditions. For effective management and conservation, in the present study, we systematically investigated the characteristics of genetic variations and runs of homozygosity (ROH) along the genome using re-sequenced whole-genome sequencing data from 235 Wenchang chickens. A total of 16,511,769 single nucleotide polymorphisms (SNPs) and 53,506 ROH segments were identified in all individuals, and the ROH of Wenchang chicken were mainly composed of short segments (0-1 megabases (Mb)). On average, 5.664% of the genome was located in ROH segments across the Wenchang chicken samples. According to several parameters, the genetic diversity of the Wenchang chicken was relatively high. The average inbreeding coefficient of Wenchang chickens based on FHOM, FGRM, and FROH was 0.060 ± 0.014, 0.561 ± 0.020, and 0.0566 ± 0.01, respectively. A total of 19 ROH islands containing 393 genes were detected on 9 different autosomes. Some of these genes were putatively associated with growth performance (AMY1a), stress resistance (THEMIS2, PIK3C2B), meat traits (MBTPS1, DLK1, and EPS8L2), and fat deposition (LANCL2, PPARγ). These findings provide a better understanding of the degree of inbreeding in Wenchang chickens and the hereditary basis of the characteristics shaped under selection. These results are valuable for the future breeding, conservation, and utilization of Wenchang and other chicken breeds.
Collapse
Affiliation(s)
- Shuaishuai Tian
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Wendan Tang
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Ziqi Zhong
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Ziyi Wang
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Xinfeng Xie
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Hong Liu
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Fuwen Chen
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Jiaxin Liu
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Yuxin Han
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Yao Qin
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Zhen Tan
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| | - Qian Xiao
- Hainan Key Laboratory of Tropical Animal Reproduction & Breeding and Epidemic Disease Research, College of Animal Science and Technology, Hainan University, Haikou 570228, China
| |
Collapse
|
15
|
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: 6] [Impact Index Per Article: 6.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
|
16
|
Wang Q, Zhang J, Wang H, Wang Z, Li Q, Zhao G, Zheng M, Wen J. Estimates of genomic inbreeding and identification of candidate regions in Beijing-You chicken populations. Anim Genet 2023; 54:155-165. [PMID: 36541281 DOI: 10.1111/age.13286] [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: 04/12/2022] [Revised: 10/09/2022] [Accepted: 12/02/2022] [Indexed: 12/24/2022]
Abstract
Runs of homozygosity (ROHs) has become an effective method for analysing inbreeding in livestock populations. Moreover, ROHs is well-suited to detect signatures of selection via ROH islands. This study aimed to investigate the occurrence and distribution of ROHs, compare the genomic inbreeding coefficients and identify the genomic regions with high ROH frequencies in different Beijing-You chicken (BY) populations, including a random conservation population (BY_R), a pedigree conservation population (BY_P), and a commercial population obtained from the market (BY_C). Among them, BY_R in 2010 and 2019 were BY_R1 and BY_R2 respectively. A total of 27 916 ROHs were identified. The average number of ROHs per individual across the three BY populations ranged from 213 (BY_P) to 161 (BY_C), and the average length of ROHs ranged from 0.432 Mb (BY_R2) to 0.451 Mb (BY_P). The highest inbreeding coefficient calculated based on ROHs (FROH ) was 0.1 in BY_P, whereas the lowest FROH was 0.0743 in BY_C. In addition, the inbreeding coefficient of BY_R2 (FROH = 0.0798) was higher than that of BY_R1 (FROH = 0.0579). Furthermore, the highest proportion of long ROH fragments (>4 Mb) was observed in BY_P and BY_C. This study showed the top 10 ROH islands of each population, and these ROH islands harboured 53 genes, some of which were related to limb development, body size and immune response. These findings contribute to the understanding of genetic diversity and population demography, and might help improve breeding and conservation strategies for BY populations.
Collapse
Affiliation(s)
- Qiao Wang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Jin Zhang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Hailong Wang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Zixuan Wang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Qinghe Li
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Guiping Zhao
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Maiqing Zheng
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Jie Wen
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| |
Collapse
|
17
|
Visser C, Lashmar SF, Reding J, Berry DP, van Marle-Köster E. Pedigree and genome-based patterns of homozygosity in the South African Ayrshire, Holstein, and Jersey breeds. Front Genet 2023; 14:1136078. [PMID: 37007942 PMCID: PMC10063850 DOI: 10.3389/fgene.2023.1136078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2023] [Accepted: 03/02/2023] [Indexed: 03/19/2023] Open
Abstract
The erosion of genetic diversity limits long-term genetic gain and impedes the sustainability of livestock production. In the South African (SA) dairy industry, the major commercial dairy breeds have been applying estimated breeding values (EBVs) and/or have been participating in Multiple Across Country Evaluations (MACE). The transition to genomic estimated breeding values (GEBVs) in selection strategies requires monitoring of the genetic diversity and inbreeding of current genotyped animals, especially considering the comparatively small population sizes of global dairy breeds in SA. This study aimed to perform a homozygosity-based evaluation of the SA Ayrshire (AYR), Holstein (HST), and Jersey (JER) dairy cattle breeds. Three sources of information, namely 1) single nucleotide polymorphism (SNP) genotypes (3,199 animals genotyped for 35,572 SNPs) 2) pedigree records (7,885 AYR; 28,391 HST; 18,755 JER), and 3) identified runs of homozygosity (ROH) segments were used to quantify inbreeding related parameters. The lowest pedigree completeness was for the HST population reducing from a value of 0.990 to 0.186 for generation depths of one to six. Across all breeds, 46.7% of the detected ROH were between 4 megabase pairs (Mb) and 8 Mb in length. Two conserved homozygous haplotypes were identified in more than 70% of the JER population on Bos taurus autosome (BTA) 7. The JER breed displayed the highest level of inbreeding across all inbreeding coefficients. The mean (± standard deviation) pedigree-based inbreeding coefficient (FPED) ranged from 0.051 (±0.020) for AYR to 0.062 (±0.027) for JER, whereas SNP-based inbreeding coefficients (FSNP) ranged from 0.020 (HST) to 0.190 (JER) and ROH-based inbreeding coefficients, considering all ROH segment coverage (FROH), ranged from 0.053 (AYR) to 0.085 (JER). Within-breed Spearman correlations between pedigree-based and genome-based estimates ranged from weak (AYR: 0.132 between FPED and FROH calculated for ROH <4Mb in size) to moderate (HST: 0.584 between FPED and FSNP). Correlations strengthened between FPED and FROH as the ROH length category was considered lengthened, suggesting a dependency on breed-specific pedigree depth. The genomic homozygosity-based parameters studied proved useful in investigating the current inbreeding status of reference populations genotyped to implement genomic selection in the three most prominent South African dairy cattle breeds.
Collapse
Affiliation(s)
- Carina Visser
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
- *Correspondence: Carina Visser,
| | - Simon Frederick Lashmar
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
| | - Jason Reding
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
| | - Donagh P. Berry
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
- Animal and Grassland Research and Innovation Centre, Teagasc, Co. Cork, Ireland
| | - Esté van Marle-Köster
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
| |
Collapse
|
18
|
Sarviaho K, Uimari P, Martikainen K. Estimating inbreeding rate and effective population size in the Finnish Ayrshire population in the era of genomic selection. J Anim Breed Genet 2023; 140:343-353. [PMID: 36808142 DOI: 10.1111/jbg.12762] [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: 01/27/2022] [Accepted: 01/26/2023] [Indexed: 02/23/2023]
Abstract
Genomic selection has been applied in dairy cattle breeding over the last decade. Using genomic information may speed up genetic gain as breeding values can be predicted reasonably accurately directly after birth. However, genetic diversity may decrease if the inbreeding rate per generation increases and the effective population size decreases. Despite many positive qualities of the Finnish Ayrshire, for example, high average protein yield and fertility, over time the breed has lost its place as the most common dairy breed in Finland. Thus, maintaining the genetic variability of the breed is becoming more important. The aim of our research was to estimate the impact of genomic selection on inbreeding rate and effective population size using both pedigree and genomic data. The genomic data included 46,914 imputed single nucleotide polymorphism (SNP) variants from 75,038 individuals, and the pedigree data included 2,770,025 individuals. All animals in the data were born between 2000 and 2020. Genomic inbreeding coefficients were estimated as the proportion of SNPs in runs of homozygosity (ROH) out of the total number of SNPs. The inbreeding rate was estimated by regressing the mean genomic inbreeding coefficients on birth years. Effective population size was then estimated based on the inbreeding rate. Additionally, effective population size was estimated from the mean increase in individual inbreeding using pedigree data. Introduction of genomic selection was assumed to have taken place gradually; years 2012-2014 were treated as a transition period from the traditional phenotype-based breeding value estimation to genomic-based estimation. The median length of the identified homozygous segments was 5.5 Mbp, and a slight increase in the proportion of segments over 10 Mbp was observed after 2010. The inbreeding rate decreased from 2000 to 2011 and subsequently increased slightly. The pedigree- and genomic-based estimates of inbreeding rate were similar to each other. The estimates of effective population size based on the regression method were very sensitive to the number of years considered; thus, the estimates were not very reliable. The effective population size estimated from the mean increase in individual inbreeding reached its highest value of 160 in 2011 and decreased to 150 after that. In addition, the generation interval in the sire path has decreased from 5.5 years to 3.5 years after genomic selection was implemented. Based on our results, after the implementation of genomic selection, the proportion of long ROH stretches has increased, the generation interval in the sire path has decreased, the inbreeding rate has increased and the effective population size has decreased. However, the effective population size is still at a good level, allowing for an efficient selection scheme in the Finnish Ayrshire breed.
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
|
19
|
Genome-Wide Genetic Structure of Henan Indigenous Chicken Breeds. Animals (Basel) 2023; 13:ani13040753. [PMID: 36830540 PMCID: PMC9952073 DOI: 10.3390/ani13040753] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2023] [Revised: 02/14/2023] [Accepted: 02/16/2023] [Indexed: 02/22/2023] Open
Abstract
There are five indigenous chicken breeds in Henan Province, China. These breeds have their own unique phenotypic characteristics in terms of morphology, behavior, skin and feather color, and productive performance, but their genetic basis is not well understood. Therefore, we analyzed the genetic structure, genomic diversity, and migration history of Henan indigenous chicken populations and the selection signals and genes responsible for Henan gamecock unique phenotypes using whole genome resequencing. The results indicate that Henan native chickens clustered most closely with the chicken populations in neighboring provinces. Compared to other breeds, Henan gamecock's inbreeding and selection intensity were more stringent. TreeMix analysis revealed the gene flow from southern chicken breeds into the Zhengyang sanhuang chicken and from the Xichuan black-bone chicken into the Gushi chicken. Selective sweep analysis identified several genes and biological processes/pathways that were related to body size, head control, muscle development, reproduction, and aggression control. Additionally, we confirmed the association between genotypes of SNPs in the strong selective gene LCORL and body size and muscle development in the Gushi-Anka F2 resource population. These findings made it easier to understand the traits of the germplasm and the potential for using the Henan indigenous chicken.
Collapse
|
20
|
Gao C, Du W, Tian K, Wang K, Wang C, Sun G, Kang X, Li W. Analysis of Conservation Priorities and Runs of Homozygosity Patterns for Chinese Indigenous Chicken Breeds. Animals (Basel) 2023; 13:ani13040599. [PMID: 36830386 PMCID: PMC9951684 DOI: 10.3390/ani13040599] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2022] [Revised: 01/30/2023] [Accepted: 02/06/2023] [Indexed: 02/11/2023] Open
Abstract
To achieve sustainable development of the poultry industry, the effective conservation of genetic resources has become increasingly important. In the present study, we systematically elucidated the population structure, conservation priority, and runs of homozygosity (ROH) patterns of Chinese native chicken breeds. We used a high-density genotyping dataset of 157 native chickens from eight breeds. The population structure showed different degrees of population stratification among the breeds. Chahua chicken was the most differentiated breed from the other breeds (Nei = 0.0813), and the Wannan three-yellow chicken (WanTy) showed the lowest degree of differentiation (Nei = 0.0438). On the basis of contribution priority, Xiaoshan chicken had the highest contribution to the total gene diversity (1.41%) and the maximum gene diversity of the synthetic population (31.1%). WanTy chicken showed the highest contribution to the total allelic diversity (1.31%) and the maximum allelic diversity of the syntenic population (17.0%). A total of 5242 ROH fragments and 5 ROH island regions were detected. The longest ROH fragment was 41.51 Mb. A comparison of the overlapping genomic regions between the ROH islands and QTLs in the quantitative trait loci (QTL) database showed that the annotated candidate genes were involved in crucial economic traits such as immunity, carcass weight, drumstick and leg muscle development, egg quality and egg production, abdominal fat precipitation, body weight, and feed intake. In conclusion, our findings revealed that Chahua, Xiaoshan, and WanTy should be the priority conservation breeds, which will help optimize the conservation and breeding programs for Chinese indigenous chicken breeds.
Collapse
Affiliation(s)
- Chaoqun Gao
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
- The Shennong Laboratory, Zhengzhou 450002, China
| | - Wenping Du
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
- The Shennong Laboratory, Zhengzhou 450002, China
| | - Kaiyuan Tian
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
- The Shennong Laboratory, Zhengzhou 450002, China
| | - Kejun Wang
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
| | - Chunxiu Wang
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
| | - Guirong Sun
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
- The Shennong Laboratory, Zhengzhou 450002, China
| | - Xiangtao Kang
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
- The Shennong Laboratory, Zhengzhou 450002, China
- Correspondence: (X.K.); (W.L.)
| | - Wenting Li
- Henan Innovative Engineering Research Center of Poultry Germplasm Resource, College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450002, China
- The Shennong Laboratory, Zhengzhou 450002, China
- Correspondence: (X.K.); (W.L.)
| |
Collapse
|
21
|
Lavanchy E, Goudet J. Effect of reduced genomic representation on using runs of homozygosity for inbreeding characterization. Mol Ecol Resour 2023; 23:787-802. [PMID: 36626297 DOI: 10.1111/1755-0998.13755] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2022] [Revised: 12/22/2022] [Accepted: 01/05/2023] [Indexed: 01/11/2023]
Abstract
Genomic measures of inbreeding based on identical-by-descent (IBD) segments are increasingly used to measure inbreeding and mostly estimated on SNP arrays and whole-genome sequencing (WGS) data. However, some softwares recurrently used for their estimation assume that genomic positions which have not been genotyped are nonvariant. This might be true for WGS data, but not for reduced genomic representations and can lead to spurious IBD segments estimation. In this project, we simulated the outputs of WGS, two SNP arrays of different sizes and RAD-sequencing for three populations with different sizes and histories. We compare the results of IBD segments estimation with two softwares: runs of homozygosity (ROHs) estimated with PLINK and homozygous-by-descent (HBD) segments estimated with RZooRoH. We demonstrate that to obtain meaningful estimates of inbreeding, RZooRoH requires a SNPs density 11 times smaller compared to PLINK: ranks of inbreeding coefficients were conserved among individuals above 22 SNPs/Mb for PLINK and 2 SNPs/Mb for RZooRoH. We also show that in populations with simple demographic histories, distribution of ROHs and HBD segments are correctly estimated with both SNP arrays and WGS. PLINK correctly estimated distribution of ROHs with SNP densities above 22 SNPs/Mb, while RZooRoH correctly estimated distribution of HBD segments with SNPs densities above 11 SNPs/Mb. However, in a population with a more complex demographic history, RZooRoH resulted in better distribution of IBD segments estimation compared to PLINK even with WGS data. Consequently, we advise researchers to use either methods relying on excess homozygosity averaged across SNPs or model-based HBD segments calling methods for inbreeding estimations.
Collapse
Affiliation(s)
- Eléonore Lavanchy
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, University of Lausanne, Lausanne, Switzerland
| | - Jérôme Goudet
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, University of Lausanne, Lausanne, Switzerland
| |
Collapse
|
22
|
Neumann GB, Korkuć P, Arends D, Wolf MJ, May K, König S, Brockmann GA. Genomic diversity and relationship analyses of endangered German Black Pied cattle (DSN) to 68 other taurine breeds based on whole-genome sequencing. Front Genet 2023; 13:993959. [PMID: 36712857 PMCID: PMC9875303 DOI: 10.3389/fgene.2022.993959] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2022] [Accepted: 11/28/2022] [Indexed: 01/05/2023] Open
Abstract
German Black Pied cattle (Deutsches Schwarzbuntes Niederungsrind, DSN) are an endangered dual-purpose cattle breed originating from the North Sea region. The population comprises about 2,500 cattle and is considered one of the ancestral populations of the modern Holstein breed. The current study aimed at defining the breeds closest related to DSN cattle, characterizing their genomic diversity and inbreeding. In addition, the detection of selection signatures between DSN and Holstein was a goal. Relationship analyses using fixation index (FST), phylogenetic, and admixture analyses were performed between DSN and 68 other breeds from the 1000 Bull Genomes Project. Nucleotide diversity, observed heterozygosity, and expected heterozygosity were calculated as metrics for genomic diversity. Inbreeding was measured as excess of homozygosity (FHom) and genomic inbreeding (FRoH) through runs of homozygosity (RoHs). Region-wide FST and cross-population-extended haplotype homozygosity (XP-EHH) between DSN and Holstein were used to detect selection signatures between the two breeds, and RoH islands were used to detect selection signatures within DSN and Holstein. DSN showed a close genetic relationship with breeds from the Netherlands, Belgium, Northern Germany, and Scandinavia, such as Dutch Friesian Red, Dutch Improved Red, Belgian Red White Campine, Red White Dual Purpose, Modern Angler, Modern Danish Red, and Holstein. The nucleotide diversity in DSN (0.151%) was higher than in Holstein (0.147%) and other breeds, e.g., Norwegian Red (0.149%), Red White Dual Purpose (0.149%), Swedish Red (0.149%), Hereford (0.145%), Angus (0.143%), and Jersey (0.136%). The FHom and FRoH values in DSN were among the lowest. Regions with high FST between DSN and Holstein, significant XP-EHH regions, and RoH islands detected in both breeds harbor candidate genes that were previously reported for milk, meat, fertility, production, and health traits, including one QTL detected in DSN for endoparasite infection resistance. The selection signatures between DSN and Holstein provide evidence of regions responsible for the dual-purpose properties of DSN and the milk type of Holstein. Despite the small population size, DSN has a high level of diversity and low inbreeding. FST supports its relatedness to breeds from the same geographic origin and provides information on potential gene pools that could be used to maintain diversity in DSN.
Collapse
Affiliation(s)
- Guilherme B. Neumann
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Paula Korkuć
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Danny Arends
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt-Universität zu Berlin, Berlin, Germany,Department of Applied Sciences, Northumbria University, Newcastle Upon Tyne, United Kingdom
| | - Manuel J. Wolf
- Institute of Animal Breeding and Genetics, Justus-Liebig-Universität, Giessen, Germany
| | - Katharina May
- Institute of Animal Breeding and Genetics, Justus-Liebig-Universität, Giessen, Germany
| | - Sven König
- Institute of Animal Breeding and Genetics, Justus-Liebig-Universität, Giessen, Germany
| | - Gudrun A. Brockmann
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt-Universität zu Berlin, Berlin, Germany,*Correspondence: Gudrun A. Brockmann,
| |
Collapse
|
23
|
Legacies of domestication, Neolithic diffusion and trade between Indian subcontinent and Island Southeast Asia shape maternal genetic diversity of Andaman cattle. PLoS One 2022; 17:e0278681. [PMID: 36490290 PMCID: PMC9733863 DOI: 10.1371/journal.pone.0278681] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2022] [Accepted: 11/22/2022] [Indexed: 12/13/2022] Open
Abstract
Andaman cattle is a precious indigenous livestock species endemic to Andaman and Nicobar Islands, India. Till date, origin and genetic makeup of the breed which is warranted for breed conservation is not known. Moreover, the spread of zebu cattle from Indus valley to different parts of Island Southeast Asia (ISEA) is not properly understood. Here, we report the genetic diversity, population structure of Andaman cattle and their evolution in the context of epicentre of zebu domestication and ISEA. High genetic diversity in complete mitochondrial D-loop sequences indicated the ability of the breed to withstand impending climate change. Total 81 haplotypes were detected and all of them except three belonged to Bos indicus. The presence of taurine haplotypes in Andaman cattle indicate introgression by European-derived cattle. A poor phylogenetic signal of Andaman cattle with genetic affinities with cattle of Indian subcontinent and ISEA was observed. The poor phylogenetic structure may be due to multidirectional gene flow from Indian subcontinent and ISEA, with which Andaman shares a close cultural and trade relationship from Neolithic age. We hypothesize that Andaman cattle is the outcome of Neolithic diffusion from centre of zebu domestication along with multidirectional commercial exchange between Indian subcontinent and ISEA.
Collapse
|
24
|
Assessment and Distribution of Runs of Homozygosity in Horse Breeds Representing Different Utility Types. Animals (Basel) 2022; 12:ani12233293. [PMID: 36496815 PMCID: PMC9736150 DOI: 10.3390/ani12233293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2022] [Revised: 11/18/2022] [Accepted: 11/22/2022] [Indexed: 11/29/2022] Open
Abstract
The present study reports runs of homozygosity (ROH) distribution in the genomes of six horse breeds (571 horses in total) representing three horse types (primitive, light, and draft horses) based on the 65k Equine BeadChip assay. Of major interest was the length, quantity, and frequency of ROH characteristics, as well as differences between horse breeds and types. Noticeable differences in the number, length and distribution of ROH between breeds were observed, as well as in genomic inbreeding coefficients. We also identified regions of the genome characterized by high ROH coverage, known as ROH islands, which may be signals of recent selection events. Eight to fourteen ROH islands were identified per breed, which spanned multiple genes. Many were involved in important horse breed characteristics, including WFIKNN2, CACNA1G, STXBP4, NOG, FAM184B, QDPR, LCORL, and the zinc finger protein family. Regions of the genome with zero ROH occurrences were also of major interest in specific populations. Depending on the breed, we detected between 2 to 57 no-ROH regions and identified 27 genes in these regions that were common for five breeds. These genes were involved in, e.g., muscle contractility (CACNA1A) and muscle development (miR-23, miR-24, miR-27). To sum up, the obtained results can be furthered analyzed in the topic of identification of markers unique for specific horse breed characteristics.
Collapse
|
25
|
Runs of Homozygosity Analysis Reveals Genomic Diversity and Population Structure of an Indigenous Cattle Breed in Southwest China. Animals (Basel) 2022; 12:ani12233239. [PMID: 36496760 PMCID: PMC9737016 DOI: 10.3390/ani12233239] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2022] [Revised: 10/06/2022] [Accepted: 11/11/2022] [Indexed: 11/24/2022] Open
Abstract
In aiming to achieve sustainable development goals in the livestock industry, it is becoming increasingly necessary and important for the effective conservation of genetic resources. There are some indigenous cattle breeds in Sichuan, southwest China, for which, however, the genetic diversity and population structures still remain unknown because of the unavailability of systematic breeding programs and pedigree information. Xieka cattle are an indigenous breed locally distributed in southeastern Sichuan and have a long-term evolutionary adaptation to local environments and climates. In this study, we obtained 796,828 single nucleotide polymorphisms (SNPs) through sequencing the genomes of 30 Xieka cattle and used them for analyzing the genetic diversity and runs of homozygosity (ROH). The mean nucleotide diversity was 0.28 and 72% of SNPs were found to be in the heterozygous states. A total of 4377 ROH were detected with even distribution among all autosomes, and 74% of them were lower than 1 Mb in length. Meanwhile, only five ROH were found longer than 5 Mb. We further determined 19 significant genomic regions that were obviously enriched by ROH, in which 35 positional candidate genes were found. Some of these genes have been previously reported to be significantly associated with various production traits in cattle, such as meat quality, carcass performances, and diseases. In conclusion, the relatively high degree of genetic diversity of Xieka cattle was revealed using the genomic information, and the proposed candidate genes will help us optimize the breeding programs regarding this indigenous breed.
Collapse
|
26
|
Liu B, Tao W, Feng D, Wang Y, Heizatuola N, Ahemetbai T, Wu W. Revealing Genetic Diversity and Population Structure of Endangered Altay White-Headed Cattle Population Using 100 k SNP Markers. Animals (Basel) 2022; 12:3214. [PMID: 36428441 PMCID: PMC9686749 DOI: 10.3390/ani12223214] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 10/17/2022] [Accepted: 11/18/2022] [Indexed: 11/22/2022] Open
Abstract
Understanding the genetic basis of native cattle populations that have adapted to the local environment is of great significance for formulating appropriate strategies and programs for genetic improvement and protection. Therefore, it is necessary to understand the genetic diversity and population structure of Altay white-headed cattle so as to meet the current production needs under various environments, carry out continuous genetic improvement, and promote rapid adaptation to changing environments and breeding objectives. A total of 46 individual samples of endangered Xinjiang Altay white-headed cattle were collected in this study, including nine bulls and 37 cows. To collect genotype data, 100 k SNP markers were used, and then studies of genetic diversity, genetic structure, inbreeding degree, and family analysis were carried out. A total of 101,220 SNP loci were detected, and the genotype detection rate for individuals was ≥90%. There were 85,993 SNP loci that passed quality control, of which 93.5% were polymorphic. The average effective allele number was 0.036, the Polymorphism Information Content was 0.304 and the minimum allele frequency was 0.309, the average observed heterozygosity was 0.413, and the average expected heterozygosity was 0.403. The average genetic distance of Idengtical By State (IBS) was 0.3090, there were 461 ROH (genome-length homozygous fragments), 76.1% of which were between 1 and 5 MB in length, and the average inbreeding coefficient was 0.016. The 46 Altay white-headed cattle were divided into their families, and the individual numbers of each family were obviously different. To sum up, the Altay white-headed cattle conservation population had low heterozygosity, a high inbreeding degree, few families, and large differences in the number of individuals in each family, which can easily cause a loss of genetic diversity. In the follow-up seed conservation process, seed selection and matching should be carried out according to the divided families to ensure the long-term protection of Altay white-headed cattle genetic resources.
Collapse
Affiliation(s)
- Bo Liu
- Key Laboratory of Genetics Breeding and Reproduction of Xinjiang Wool Sheep & Cashmere Goat, Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi 830000, China
| | - Weikun Tao
- Key Laboratory of Genetics Breeding and Reproduction of Xinjiang Wool Sheep & Cashmere Goat, Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi 830000, China
| | - Donghe Feng
- Key Laboratory of Genetics Breeding and Reproduction of Xinjiang Wool Sheep & Cashmere Goat, Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi 830000, China
| | - Yue Wang
- Animal Husbandry Workstation in Altay Region, Altay 836000, China
| | | | - Tenes Ahemetbai
- Animal Husbandry Workstation in Altay Region, Altay 836000, China
| | - Weiwei Wu
- Key Laboratory of Genetics Breeding and Reproduction of Xinjiang Wool Sheep & Cashmere Goat, Institute of Animal Science, Xinjiang Academy of Animal Sciences, Urumqi 830000, China
| |
Collapse
|
27
|
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
|
28
|
Cortes-Hernández JG, Ruiz-López FJ, Vásquez-Peláez CG, García-Ruiz A. Runs of homocigosity and its association with productive traits in Mexican Holstein cattle. PLoS One 2022; 17:e0274743. [PMID: 36121861 PMCID: PMC9484644 DOI: 10.1371/journal.pone.0274743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 09/03/2022] [Indexed: 11/21/2022] Open
Abstract
The objective of this study was to describe the runs of homozygosity (ROH) detected in the Mexican Holstein population and to associate them with milk, fat and protein yields, and conformation final score. After imputation and genomic quality control, 4,227 genotyped animals with 100,806 SNPs markers each were used. ROH with a minimum length of 1 Mb and a minimum of 10 SNPs were included in the analysis. One heterozygous SNP marker and five missing genotypes per ROH were allowed. A total of 425,098 ROH were found in the studied population (71.83 ± 10.73 ROH per animal), with an average length and coverage of 4.80 ± 0.77 Mb, and 276.89 Mb, respectively. The average chromosome length covered by ROH was 10.40 ± 3.70 Mb. ROH between 1 and 2 Mb were the most frequent in the population (51.33%) while those between 14 and 16 Mb were the least frequent (1.20%). Long chromosomes showed a larger number of ROH. Chromosomes 10 and 20, had a greater percentage of their length covered by ROH because they presented a largest number of long ROH (>8 Mb). From the total ROH, 17 were detected in 1,847 animals and distributed among different chromosomes, and were associated with milk, fat and protein yield and percentage, and conformation final score. Of the ROH with effects on production traits, the majority were found with a length between 1 and 4 Mb. These results show evidence of genomic regions preserved by genetic selection and associated with the improvement of the productivity and functionality of dairy cattle.
Collapse
Affiliation(s)
- José G. Cortes-Hernández
- Programa de Maestría y Doctorado en Ciencias de la Producción y de la Salud Animal, Universidad Nacional Autónoma de México, Ciudad de México, México
| | - Felipe J. Ruiz-López
- Centro Nacional de Investigación Disciplinaria en Fisiología y Mejoramiento Animal, Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias, Ajuchitlán Colón Querétaro, México
| | - Carlos G. Vásquez-Peláez
- Departamento de Genética y Bioestadística de la Facultad de Medicina Veterinaria y Zootecnia, Universidad Nacional Autónoma de México, Ciudad de México, México
| | - Adriana García-Ruiz
- Centro Nacional de Investigación Disciplinaria en Fisiología y Mejoramiento Animal, Instituto Nacional de Investigaciones Forestales Agrícolas y Pecuarias, Ajuchitlán Colón Querétaro, México
- * E-mail:
| |
Collapse
|
29
|
Genetic Diversity and Selection Signatures in Jianchang Black Goats Revealed by Whole-Genome Sequencing Data. Animals (Basel) 2022; 12:ani12182365. [PMID: 36139225 PMCID: PMC9495118 DOI: 10.3390/ani12182365] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2022] [Revised: 08/28/2022] [Accepted: 09/07/2022] [Indexed: 11/17/2022] Open
Abstract
Understanding the genetic composition of indigenous goats is essential to promote the scientific conservation and sustainable utilization of these breeds. The Jianchang Black (JC) goat, a Chinese native breed, is solid black and exhibits crude feed tolerance, but is characterized by a low growth rate and small body size. Based on the whole-genome sequencing data for 30 JC, 41 Jintang Black (JT), and 40 Yunshang Black (YS) goats, and 21 Bezoar ibexes, here, we investigated the genetic composition of JC goats by conducting analyses of the population structure, runs of homozygosity (ROH), genomic inbreeding, and selection signature. Our results revealed that JT and YS showed a close genetic relationship with a non-negligible amount of gene flows but were genetically distant from JC, apart from Bezoars. An average of 2039 ROHs were present in the autosomal genome per individual. The ROH-based inbreeding estimates in JC goats generally showed moderate values ranging from 0.134 to 0.264, mainly due to rapid declines in the effective population size during recent generations. The annotated genes (e.g., IL2, IL7, and KIT) overlapping with ROH islands were significantly enriched in immune-related biological processes. Further, we found 61 genes (e.g., STIM1, MYO9A, and KHDRBS2) under positive selection in JC goats via three complementary approaches, which may underly genetic adaptations to local environmental conditions. Our findings provided references for the conservation and sustainable utilization of JC goats.
Collapse
|
30
|
Wu X, Zhou R, Wang Y, Zhang W, Zheng X, Zhao G, Zhang X, Yin Z, Ding Y. Genome‐wide scan for runs of homozygosity in Asian wild boars and Anqing six‐end‐white pigs. Anim Genet 2022; 53:867-871. [DOI: 10.1111/age.13250] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 04/26/2022] [Accepted: 08/02/2022] [Indexed: 01/28/2023]
Affiliation(s)
- Xudong Wu
- Anhui Provincial Laboratory of Local Animal Genetic Resource Conservation and Bio‐Breeding, College of Animal Science and Technology Anhui Agricultural University Hefei Anhui China
- Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences Hefei China
| | - Ren Zhou
- Anhui Provincial Laboratory of Local Animal Genetic Resource Conservation and Bio‐Breeding, College of Animal Science and Technology Anhui Agricultural University Hefei Anhui China
| | - Yuanlang Wang
- Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences Hefei China
| | - Wei Zhang
- Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences Hefei China
| | - Xianrui Zheng
- Anhui Provincial Laboratory of Local Animal Genetic Resource Conservation and Bio‐Breeding, College of Animal Science and Technology Anhui Agricultural University Hefei Anhui China
| | - Guiying Zhao
- College of Animal Science and Technology Yunnan Agricultural University Kunming China
| | - Xiaodong Zhang
- Anhui Provincial Laboratory of Local Animal Genetic Resource Conservation and Bio‐Breeding, College of Animal Science and Technology Anhui Agricultural University Hefei Anhui China
| | - Zongjun Yin
- Anhui Provincial Laboratory of Local Animal Genetic Resource Conservation and Bio‐Breeding, College of Animal Science and Technology Anhui Agricultural University Hefei Anhui China
| | - Yueyun Ding
- Anhui Provincial Laboratory of Local Animal Genetic Resource Conservation and Bio‐Breeding, College of Animal Science and Technology Anhui Agricultural University Hefei Anhui China
| |
Collapse
|
31
|
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
|
32
|
Shihabi M, Lukic B, Cubric-Curik V, Brajkovic V, Oršanić M, Ugarković D, Vostry L, Curik I. Identification of Selection Signals on the X-Chromosome in East Adriatic Sheep: A New Complementary Approach. Front Genet 2022; 13:887582. [PMID: 35615375 PMCID: PMC9126029 DOI: 10.3389/fgene.2022.887582] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2022] [Accepted: 03/21/2022] [Indexed: 12/13/2022] Open
Abstract
Sheep are one of the most important livestock species in Croatia, found mainly in the Mediterranean coastal and mountainous regions along the East Adriatic coast, well adapted to the environment and mostly kept extensively. Our main objective was therefore to map the positive selection of the X-chromosome (18,983 SNPs that passed quality control), since nothing is known about the adaptation genes on this chromosome for any of the breeds from the Balkan cluster. Analyses were performed on a sample of eight native Croatian breeds (101 females and 100 males) representing the East Adriatic metapopulation and on 10 mouflons (five females and males), all sampled in Croatia. Three classical within-population approaches (extreme Runs of Homozygosity islands, integrated Haplotype Score, and number of Segregating Sites by Length) were applied along with our new approach called Haplotype Richness Drop (HRiD), which uses only the information contained in male haplotypes. We have also shown that phylogenetic analyses, such as the Median-joining network, can provide additional information when performed with the selection signals identified by HRiD. Our new approach identifies positive selection signals by searching for genomic regions that exhibit a sudden decline in haplotype richness. In total, we identified 14 positive selection signals, 11 using the classical approach and three using the HRiD approach, all together containing 34 annotated genes. The most reliable selection signal was mapped by all four approaches in the same region, overlapping between 13.17 and 13.60 Mb, and assigned to the CA5B, ZRSR2, AP1S2, and GRPR genes. High repeatability (86%) of results was observed, as 12 identified selection signals were also confirmed in other studies with sheep. HRiD offers an interesting possibility to be used complementary to other approaches or when only males are genotyped, which is often the case in genomic breeding value estimations. These results highlight the importance of the X-chromosome in the adaptive architecture of domestic ruminants, while our novel HRiD approach opens new possibilities for research.
Collapse
Affiliation(s)
- Mario Shihabi
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Zagreb, Croatia
- *Correspondence: Mario Shihabi, ; Ino Curik,
| | - Boris Lukic
- Department for Animal Production and Biotechnology, Faculty of Agrobiotechnical Sciences Osijek, J.J. Strossmayer University of Osijek, Osijek, Croatia
| | - Vlatka Cubric-Curik
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Zagreb, Croatia
| | - Vladimir Brajkovic
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Zagreb, Croatia
| | - Milan Oršanić
- Department of Forest Ecology and Silviculture, Faculty of Forestry and Wood Technology, University of Zagreb, Zagreb, Croatia
| | - Damir Ugarković
- Department of Forest Ecology and Silviculture, Faculty of Forestry and Wood Technology, University of Zagreb, Zagreb, Croatia
| | - Luboš Vostry
- Department of Genetics and Breeding, Faculty Agrobiology, Food and Natural Resources, Czech University of Life Sciences, Prague, Czechia
| | - Ino Curik
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Zagreb, Croatia
- *Correspondence: Mario Shihabi, ; Ino Curik,
| |
Collapse
|
33
|
Nazari F, Seyedabadi HR, Noshary A, Emamjomeh-Kashan N, Banabazi MH. A genome-wide scan for signatures of selection in Kurdish horse breed. J Equine Vet Sci 2022; 113:103916. [PMID: 35218903 DOI: 10.1016/j.jevs.2022.103916] [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: 11/10/2021] [Revised: 02/20/2022] [Accepted: 02/21/2022] [Indexed: 11/17/2022]
Abstract
The genetic diversity and genomic regions being under putative natural selection in Kurdish horse population were studied. The samples from 72 horses were genotyped by using GGP Equine 70K SNP arrays. The Ne Slope (NeS) analyses revealed that a sharp decline in Ne has probably occurred around four generations ago, and high frequency of ROH with 2-4 Mbp in length suggested that the inbreeding has probably occurred around 20 generations ago. The effective population size (Ne) was 104 horses up to three generations ago and the average inbreeding (FROH) was 0.047(±0.045). Using de-correlated composite of multiple selection signals (DCMS) and runs of homozygosity (ROH) analyses the genomic regions being under putative selection were detected. By using DCMS, a total of 148 significant SNP (FDR < 0.05) were identified, 40% of which were located on ECA9, where the greatest peak was observed. This genomic region harbors several known QTL which are associated with withers height (body size). Also significant genomic regions (FDR<0.05), harboring QTL associated with insect bite hypersensitivity (IBH), hair density and coat texture, alternate gaits, guttural pouch tympany and temperament were identified. By using outputs of ROH analyses, two hotspot regions (i.e. 30% of individuals was considered as threshold), were identified on ECA7 (50.11-54.36 Mbp) and ECA11 (26.10-29.07 Mbp) harboring QTL associated with withers height, alternate gait and IBH. In summary, the genomic regions being under putative natural selection which harbors known QTL associated with body size and IBH, among others, were introduced. Nevertheless, additional functional and comparative studies are necessary to corroborate their effect on the observed genetic and phenotypic diversity of the Kurdish horses.
Collapse
Affiliation(s)
- Faezeh Nazari
- Department of Animal Science, Science and Research Branch, Islamic Azad University, Tehran, Iran
| | - Hamid-Reza Seyedabadi
- Animal Science Research Institute of Iran (ASRI), Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran.
| | - Alireza Noshary
- Department of Animal Science, Faculty of Agriculture, Karaj Branch, Islamic Azad University, Karaj, Iran
| | - Nasser Emamjomeh-Kashan
- Department of Animal Science, Science and Research Branch, Islamic Azad University, Tehran, Iran
| | - Mohamad-Hosein Banabazi
- Animal Science Research Institute of Iran (ASRI), Agricultural Research, Education and Extension Organization (AREEO), Karaj, Iran; Department of animal breeding and genetics (HGEN), Centre for Veterinary Medicine and Animal Science (VHC), Swedish University of Agricultural Sciences (SLU), Uppsala 75007, Sweden
| |
Collapse
|
34
|
Palma-Vera SE, Reyer H, Langhammer M, Reinsch N, Derezanin L, Fickel J, Qanbari S, Weitzel JM, Franzenburg S, Hemmrich-Stanisak G, Schoen J. Genomic characterization of the world's longest selection experiment in mouse reveals the complexity of polygenic traits. BMC Biol 2022; 20:52. [PMID: 35189878 PMCID: PMC8862358 DOI: 10.1186/s12915-022-01248-9] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 02/07/2022] [Indexed: 12/11/2022] Open
Abstract
BACKGROUND Long-term selection experiments are a powerful tool to understand the genetic background of complex traits. The longest of such experiments has been conducted in the Research Institute for Farm Animal Biology (FBN), generating extreme mouse lines with increased fertility, body mass, protein mass and endurance. For >140 generations, these lines have been maintained alongside an unselected control line, representing a valuable resource for understanding the genetic basis of polygenic traits. However, their history and genomes have not been reported in a comprehensive manner yet. Therefore, the aim of this study is to provide a summary of the breeding history and phenotypic traits of these lines along with their genomic characteristics. We further attempt to decipher the effects of the observed line-specific patterns of genetic variation on each of the selected traits. RESULTS Over the course of >140 generations, selection on the control line has given rise to two extremely fertile lines (>20 pups per litter each), two giant growth lines (one lean, one obese) and one long-distance running line. Whole genome sequencing analysis on 25 animals per line revealed line-specific patterns of genetic variation among lines, as well as high levels of homozygosity within lines. This high degree of distinctiveness results from the combined effects of long-term continuous selection, genetic drift, population bottleneck and isolation. Detection of line-specific patterns of genetic differentiation and structural variation revealed multiple candidate genes behind the improvement of the selected traits. CONCLUSIONS The genomes of the Dummerstorf trait-selected mouse lines display distinct patterns of genomic variation harbouring multiple trait-relevant genes. Low levels of within-line genetic diversity indicate that many of the beneficial alleles have arrived to fixation alongside with neutral alleles. This study represents the first step in deciphering the influence of selection and neutral evolutionary forces on the genomes of these extreme mouse lines and depicts the genetic complexity underlying polygenic traits.
Collapse
Affiliation(s)
- Sergio E Palma-Vera
- Institute of Reproductive Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany.
| | - Henry Reyer
- Institute of Genome Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Martina Langhammer
- Institute of Genetics and Biometry, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Norbert Reinsch
- Institute of Genetics and Biometry, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Lorena Derezanin
- Institute of Reproductive Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
- Department of Evolutionary Genetics, Research Institute for Zoo and Wildlife Research (IZW), Berlin, Germany
| | - Joerns Fickel
- Department of Evolutionary Genetics, Research Institute for Zoo and Wildlife Research (IZW), Berlin, Germany
- University of Potsdam, Institute for Biochemistry and Biology, Potsdam, Germany
| | - Saber Qanbari
- Institute of Genetics and Biometry, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Joachim M Weitzel
- Institute of Reproductive Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | | | | | - Jennifer Schoen
- Institute of Reproductive Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
- Department of Reproduction Biology, Research Institute for Zoo and Wildlife Research (IZW), Berlin, Germany
| |
Collapse
|
35
|
Ablondi M, Sabbioni A, Stocco G, Cipolat-Gotet C, Dadousis C, van Kaam JT, Finocchiaro R, Summer A. Genetic Diversity in the Italian Holstein Dairy Cattle Based on Pedigree and SNP Data Prior and After Genomic Selection. Front Vet Sci 2022; 8:773985. [PMID: 35097040 PMCID: PMC8792952 DOI: 10.3389/fvets.2021.773985] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Accepted: 11/30/2021] [Indexed: 01/09/2023] Open
Abstract
Genetic diversity has become an urgent matter not only in small local breeds but also in more specialized ones. While the use of genomic data in livestock breeding programs increased genetic gain, there is increasing evidence that this benefit may be counterbalanced by the potential loss of genetic variability. Thus, in this study, we aimed to investigate the genetic diversity in the Italian Holstein dairy cattle using pedigree and genomic data from cows born between 2002 and 2020. We estimated variation in inbreeding, effective population size, and generation interval and compared those aspects prior to and after the introduction of genomic selection in the breed. The dataset contained 84,443 single-nucleotide polymorphisms (SNPs), and 74,485 cows were analyzed. Pedigree depth based on complete generation equivalent was equal to 10.67. A run of homozygosity (ROH) analysis was adopted to estimate SNP-based inbreeding (FROH). The average pedigree inbreeding was 0.07, while the average FROH was more than double, being equal to 0.17. The pattern of the effective population size based on pedigree and SNP data was similar although different in scale, with a constant decrease within the last five generations. The overall inbreeding rate (ΔF) per year was equal to +0.27% and +0.44% for Fped and FROH throughout the studied period, which corresponded to about +1.35% and +2.2% per generation, respectively. A significant increase in the ΔF was found since the introduction of genomic selection in the breed. This study in the Italian Holstein dairy cattle showed the importance of controlling the loss of genetic diversity to ensure the long-term sustainability of this breed, as well as to guarantee future market demands.
Collapse
Affiliation(s)
- Michela Ablondi
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Alberto Sabbioni
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Giorgia Stocco
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Claudio Cipolat-Gotet
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
- *Correspondence: Claudio Cipolat-Gotet
| | - Christos Dadousis
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| | - Jan-Thijs van Kaam
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana, Cremona, Italy
| | - Raffaella Finocchiaro
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana, Cremona, Italy
| | - Andrea Summer
- Dipartimento di Scienze Medico-Veterinarie, University of Parma, Parma, Italy
| |
Collapse
|
36
|
Neumann GB, Korkuć P, Arends D, Wolf MJ, May K, Reißmann M, Elzaki S, König S, Brockmann GA. Design and performance of a bovine 200 k SNP chip developed for endangered German Black Pied cattle (DSN). BMC Genomics 2021; 22:905. [PMID: 34922441 PMCID: PMC8684242 DOI: 10.1186/s12864-021-08237-2] [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: 07/12/2021] [Accepted: 12/03/2021] [Indexed: 08/26/2023] Open
Abstract
BACKGROUND German Black Pied cattle (DSN) are an endangered dual-purpose breed which was largely replaced by Holstein cattle due to their lower milk yield. DSN cattle are kept as a genetic reserve with a current herd size of around 2500 animals. The ability to track sequence variants specific to DSN could help to support the conservation of DSN's genetic diversity and to provide avenues for genetic improvement. RESULTS Whole-genome sequencing data of 304 DSN cattle were used to design a customized DSN200k SNP chip harboring 182,154 variants (173,569 SNPs and 8585 indels) based on ten selection categories. We included variants of interest to DSN such as DSN unique variants and variants from previous association studies in DSN, but also variants of general interest such as variants with predicted consequences of high, moderate, or low impact on the transcripts and SNPs from the Illumina BovineSNP50 BeadChip. Further, the selection of variants based on haplotype blocks ensured that the whole-genome was uniformly covered with an average variant distance of 14.4 kb on autosomes. Using 300 DSN and 162 animals from other cattle breeds including Holstein, endangered local cattle populations, and also a Bos indicus breed, performance of the SNP chip was evaluated. Altogether, 171,978 (94.31%) of the variants were successfully called in at least one of the analyzed breeds. In DSN, the number of successfully called variants was 166,563 (91.44%) while 156,684 (86.02%) were segregating at a minor allele frequency > 1%. The concordance rate between technical replicates was 99.83 ± 0.19%. CONCLUSION The DSN200k SNP chip was proved useful for DSN and other Bos taurus as well as one Bos indicus breed. It is suitable for genetic diversity management and marker-assisted selection of DSN animals. Moreover, variants that were segregating in other breeds can be used for the design of breed-specific customized SNP chips. This will be of great value in the application of conservation programs for endangered local populations in the future.
Collapse
Affiliation(s)
- Guilherme B Neumann
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt Universität zu Berlin, Berlin, Germany
| | - Paula Korkuć
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt Universität zu Berlin, Berlin, Germany
| | - Danny Arends
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt Universität zu Berlin, Berlin, Germany
| | - Manuel J Wolf
- Institute of Animal Breeding and Genetics, Justus-Liebig-Universität, Gießen, Germany
| | - Katharina May
- Institute of Animal Breeding and Genetics, Justus-Liebig-Universität, Gießen, Germany
| | - Monika Reißmann
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt Universität zu Berlin, Berlin, Germany
| | - Salma Elzaki
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt Universität zu Berlin, Berlin, Germany.,Department of Genetics and Animal Breeding, Faculty of Animal Production, University of Khartoum, Khartoum North, Sudan
| | - Sven König
- Institute of Animal Breeding and Genetics, Justus-Liebig-Universität, Gießen, Germany
| | - Gudrun A Brockmann
- Animal Breeding Biology and Molecular Genetics, Albrecht Daniel Thaer-Institute for Agricultural and Horticultural Sciences, Humboldt Universität zu Berlin, Berlin, Germany.
| |
Collapse
|
37
|
Sanglard LP, Huang Y, Gray KA, Linhares DCL, Dekkers JCM, Niederwerder MC, Fernando RL, Serão NVL. Further host-genomic characterization of total antibody response to PRRSV vaccination and its relationship with reproductive performance in commercial sows: genome-wide haplotype and zygosity analyses. Genet Sel Evol 2021; 53:91. [PMID: 34875996 PMCID: PMC8650375 DOI: 10.1186/s12711-021-00676-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2021] [Accepted: 10/19/2021] [Indexed: 11/16/2022] Open
Abstract
Background The possibility of using antibody response (S/P ratio) to PRRSV vaccination measured in crossbred commercial gilts as a genetic indicator for reproductive performance in vaccinated crossbred sows has motivated further studies of the genomic basis of this trait. In this study, we investigated the association of haplotypes and runs of homozygosity (ROH) and heterozygosity (ROHet) with S/P ratio and their impact on reproductive performance. Results There was no association (P-value ≥ 0.18) of S/P ratio with the percentage of ROH or ROHet, or with the percentage of heterozygosity across the whole genome or in the major histocompatibility complex (MHC) region. However, specific ROH and ROHet regions were significantly associated (P-value ≤ 0.01) with S/P ratio on chromosomes 1, 4, 5, 7, 10, 11, 13, and 17 but not (P-value ≥ 0.10) with reproductive performance. With the haplotype-based genome-wide association study (GWAS), additional genomic regions associated with S/P ratio were identified on chromosomes 4, 7, and 9. These regions harbor immune-related genes, such as SLA-DOB, TAP2, TAPBP, TMIGD3, and ADORA. Four haplotypes at the identified region on chromosome 7 were also associated with multiple reproductive traits. A haplotype significantly associated with S/P ratio that is located in the MHC region may be in stronger linkage disequilibrium (LD) with the quantitative trait loci (QTL) than the previously identified single nucleotide polymorphism (SNP) (H3GA0020505) given the larger estimate of genetic variance explained by the haplotype than by the SNP. Conclusions Specific ROH and ROHet regions were significantly associated with S/P ratio. The haplotype-based GWAS identified novel QTL for S/P ratio on chromosomes 4, 7, and 9 and confirmed the presence of at least one QTL in the MHC region. The chromosome 7 region was also associated with reproductive performance. These results narrow the search for causal genes in this region and suggest SLA-DOB and TAP2 as potential candidate genes associated with S/P ratio on chromosome 7. These results provide additional opportunities for marker-assisted selection and genomic selection for S/P ratio as genetic indicator for litter size in commercial pig populations. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00676-5.
Collapse
Affiliation(s)
- Leticia P Sanglard
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Yijian Huang
- Smithfield Premium Genetic, Rose Hill, NC, 28458, USA
| | - Kent A Gray
- Smithfield Premium Genetic, Rose Hill, NC, 28458, USA
| | - Daniel C L Linhares
- Department of Veterinary Diagnostic & Production Animal Medicine, Iowa State University, Ames, IA, 50011, USA
| | - Jack C M Dekkers
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Megan C Niederwerder
- Department of Diagnostic Medicine/Pathobiology, Kansas State University, Manhattan, KS, 66506, USA
| | - Rohan L Fernando
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA
| | - Nick V L Serão
- Department of Animal Science, Iowa State University, Ames, IA, 50011, USA.
| |
Collapse
|
38
|
Makanjuola BO, Maltecca C, Miglior F, Marras G, Abdalla EA, Schenkel FS, Baes CF. Identification of unique ROH regions with unfavorable effects on production and fertility traits in Canadian Holsteins. Genet Sel Evol 2021; 53:68. [PMID: 34461820 PMCID: PMC8406729 DOI: 10.1186/s12711-021-00660-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2020] [Accepted: 08/19/2021] [Indexed: 12/02/2022] Open
Abstract
Background The advent of genomic information and the reduction in the cost of genotyping have led to the use of genomic information to estimate genomic inbreeding as an alternative to pedigree inbreeding. Using genomic measures, effects of genomic inbreeding on production and fertility traits have been observed. However, there have been limited studies on the specific genomic regions causing the observed negative association with the trait of interest. Our aim was to identify unique run of homozygosity (ROH) genotypes present within a given genomic window that display negative associations with production and fertility traits and to quantify the effects of these identified ROH genotypes. Methods In total, 50,575 genotypes based on a 50K single nucleotide polymorphism (SNP) array and 259,871 pedigree records were available. Of these 50,575 genotypes, 46,430 cows with phenotypic records for production and fertility traits and having a first calving date between 2008 and 2018 were available. Unique ROH genotypes identified using a sliding-window approach were fitted into an animal mixed model as fixed effects to determine their effect on production and fertility traits. Results In total, 133 and 34 unique ROH genotypes with unfavorable effects were identified for production and fertility traits, respectively, at a 1% genome-wise false discovery rate. Most of these ROH regions were located on bovine chromosomes 8, 13, 14 and 19 for both production and fertility traits. For production traits, the average of all the unfavorably identified unique ROH genotypes effects were estimated to decrease milk yield by 247.30 kg, fat yield by 11.46 kg and protein yield by 8.11 kg. Similarly, for fertility traits, an average 4.81-day extension in first service to conception, a 0.16 increase in number of services, and a − 0.07 incidence in 56-day non-return rate were observed. Furthermore, a ROH region located on bovine chromosome 19 was identified that, when homozygous, had a negative effect on production traits. Signatures of selection proximate to this region have implicated GH1 as a potential candidate gene, which encodes the growth hormone that binds the growth hormone receptor. This observed negative effect could be a consequence of unfavorable alleles in linkage disequilibrium with favorable alleles. Conclusions ROH genotypes with unfavorable effects on production and fertility traits were identified within and across multiple traits on most chromosomes. These identified ROH genotypes could be included in mate selection programs to minimize their frequency in future generations. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-021-00660-z.
Collapse
Affiliation(s)
- Bayode O Makanjuola
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.
| | - Christian Maltecca
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.,Department of Animal Science and Genetics Program, North Carolina State University, Raleigh, NC, 27607, USA
| | - Filippo Miglior
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | | | - Emhimad A Abdalla
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Flavio S Schenkel
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - Christine F Baes
- Centre for Genomic Improvement of Livestock, Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada.,Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| |
Collapse
|
39
|
Tarekegn GM, Khayatzadeh N, Liu B, Osama S, Haile A, Rischkowsky B, Zhang W, Tesfaye K, Dessie T, Mwai OA, Djikeng A, Mwacharo JM. Ethiopian indigenous goats offer insights into past and recent demographic dynamics and local adaptation in sub-Saharan African goats. Evol Appl 2021; 14:1716-1731. [PMID: 34295359 PMCID: PMC8287980 DOI: 10.1111/eva.13118] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Revised: 07/15/2020] [Accepted: 08/13/2020] [Indexed: 12/13/2022] Open
Abstract
Knowledge on how adaptive evolution and human socio-cultural and economic interests shaped livestock genomes particularly in sub-Saharan Africa remains limited. Ethiopia is in a geographic region that has been critical in the history of African agriculture with ancient and diverse human ethnicity and bio-climatic conditions. Using 52K genome-wide data analysed in 646 individuals from 13 Ethiopian indigenous goat populations, we observed high levels of genetic variation. Although runs of homozygosity (ROH) were ubiquitous genome-wide, there were clear differences in patterns of ROH length and abundance and in effective population sizes illustrating differences in genome homozygosity, evolutionary history, and management. Phylogenetic analysis incorporating patterns of genetic differentiation and gene flow with ancestry modelling highlighted past and recent intermixing and possible two deep ancient genetic ancestries that could have been brought by humans with the first introduction of goats in Africa. We observed four strong selection signatures that were specific to Arsi-Bale and Nubian goats. These signatures overlapped genomic regions with genes associated with morphological, adaptation, reproduction and production traits due possibly to selection under environmental constraints and/or human preferences. The regions also overlapped uncharacterized genes, calling for a comprehensive annotation of the goat genome. Our results provide insights into mechanisms leading to genome variation and differentiation in sub-Saharan Africa indigenous goats.
Collapse
Affiliation(s)
- Getinet M. Tarekegn
- Department of Animal Production and TechnologySchool of Animal Sciences and Veterinary MedicineBahir Dar UniversityBahir DarEthiopia
- Department of Animal Breeding and GeneticsSwedish University of Agricultural Sciences (SLU)UppsalaSweden
| | - Negar Khayatzadeh
- Department of Sustainable Agricultural SystemsDivision of Livestock SciencesUniversity of Natural Resources and Life SciencesViennaAustria
| | - Bin Liu
- Inner Mongolia Agricultural UniversityHohhotChina
| | - Sarah Osama
- The University of QueenslandSaint LuciaQLDAustralia
| | - Aynalem Haile
- Small Ruminant GenomicsInternational Centre for Agricultural Research in the Dry Areas (ICARDA)Addis AbabaEthiopia
| | - Barbara Rischkowsky
- Small Ruminant GenomicsInternational Centre for Agricultural Research in the Dry Areas (ICARDA)Addis AbabaEthiopia
| | | | - Kassahun Tesfaye
- Department of Microbial, Cellular and Molecular BiologyAddis Ababa UniversityAddis AbabaEthiopia
| | - Tadelle Dessie
- International Livestock Research Institute (ILRI)Addis AbabaEthiopia
| | - Okeyo A. Mwai
- International Livestock Research Institute (ILRI)NairobiKenya
| | - Appolinaire Djikeng
- Animal and Veterinary Sciences Group, SRUC and Centre for Tropical Livestock Genetics and Health (CTLGH)The Roslin InstituteEaster BushMidlothianUK
| | - Joram M. Mwacharo
- Small Ruminant GenomicsInternational Centre for Agricultural Research in the Dry Areas (ICARDA)Addis AbabaEthiopia
- Animal and Veterinary Sciences Group, SRUC and Centre for Tropical Livestock Genetics and Health (CTLGH)The Roslin InstituteEaster BushMidlothianUK
| |
Collapse
|
40
|
Deniskova TE, Dotsev AV, Selionova MI, Reyer H, Sölkner J, Fornara MS, Aybazov AMM, Wimmers K, Brem G, Zinovieva NA. SNP-Based Genotyping Provides Insight Into the West Asian Origin of Russian Local Goats. Front Genet 2021; 12:708740. [PMID: 34276802 PMCID: PMC8282346 DOI: 10.3389/fgene.2021.708740] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2021] [Accepted: 06/07/2021] [Indexed: 11/13/2022] Open
Abstract
Specific local environmental and sociocultural conditions have led to the creation of various goat populations in Russia. National goat diversity includes breeds that have been selected for down and mohair production traits as well as versatile local breeds for which pastoralism is the main management system. Effective preservation and breeding programs for local goat breeds are missing due to the lack of DNA-based data. In this work, we analyzed the genetic diversity and population structure of Russian local goats, including Altai Mountain, Altai White Downy, Dagestan Downy, Dagestan Local, Karachaev, Orenburg, and Soviet Mohair goats, which were genotyped with the Illumina Goat SNP50 BeadChip. In addition, we addressed genetic relationships between local and global goat populations obtained from the AdaptMap project. Russian goats showed a high level of genetic diversity. Although a decrease in historical effective population sizes was revealed, the recent effective population sizes estimated for three generations ago were larger than 100 in all studied populations. The mean runs of homozygosity (ROH) lengths ranged from 79.42 to 183.94 Mb, and the average ROH number varied from 18 to 41. Short ROH segments (<2 Mb) were predominant in all breeds, while the longest ROH class (>16 Mb) was the least frequent. Principal component analysis, Neighbor-Net graph, and Admixture clustering revealed several patterns in Russian local goats. First, a separation of the Karachaev breed from other populations was observed. Moreover, genetic connections between the Orenburg and Altai Mountain breeds were suggested and the Dagestan breeds were found to be admixed with the Soviet Mohair breed. Neighbor-Net analysis and clustering of local and global breeds demonstrated the close genetic relations between Russian local and Turkish breeds that probably resulted from past admixture events through postdomestication routes. Our findings contribute to the understanding of the genetic relationships of goats originating in West Asia and Eurasia and may be used to design breeding programs for local goats to ensure their effective conservation and proper management.
Collapse
Affiliation(s)
| | - Arsen V Dotsev
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Marina I Selionova
- Russian State Agrarian University - Moscow Timiryazev Agricultural Academy, Moscow, Russia
| | - Henry Reyer
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology, Dummerstorf, Germany
| | - Johann Sölkner
- Division of Livestock Sciences, University of Natural Resources and Life Sciences, Vienna, Vienna, Austria
| | | | - Ali-Magomed M Aybazov
- All-Russian Research Institute of Sheep and Goat Breeding - Branch of the Federal State Budgetary Scientific Institution, North Caucasian Agrarian Center, Stavropol, Russia
| | - Klaus Wimmers
- Institute of Genome Biology, Leibniz Institute for Farm Animal Biology, Dummerstorf, Germany
| | - Gottfried Brem
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia.,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Vienna, Vienna, Austria
| | | |
Collapse
|
41
|
Ablondi M, Malacarne M, Cipolat-Gotet C, van Kaam JT, Sabbioni A, Summer A. Genome-wide scan reveals genetic divergence in Italian Holstein cows bred within PDO cheese production chains. Sci Rep 2021; 11:12601. [PMID: 34131265 PMCID: PMC8206360 DOI: 10.1038/s41598-021-92168-1] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2020] [Accepted: 06/07/2021] [Indexed: 02/05/2023] Open
Abstract
Dairy cattle breeds have been exposed to intense artificial selection for milk production traits over the last fifty years. In Italy, where over 80% of milk is processed into cheese, selection has also focused on cheese-making traits. Due to a deep-rooted tradition in cheese-making, currently fifty Italian cheeses are marked with the Protected Designation of Origin (PDO) label as they proved traditional land of origin and procedures for milk transformation. This study aimed to explore from a genetic point of view if the presence of such diverse productive contexts in Italy have shaped in a different manner the genome of animals originally belonging to a same breed. We analyzed high density genotype data from 1000 Italian Holstein cows born between 2014 and 2018. Those animals were either farmed in one of four Italian PDO consortia or used for drinkable milk production only. Runs of Homozygosity, Bayesian Information Criterion and Discriminant Analysis of Principal Components were used to evaluate potential signs of genetic divergence within the breed. We showed that the analyzed Italian Holstein cows have genomic inbreeding level above 5% in all subgroups, reflecting the presence of ongoing artificial selection in the breed. Our study provided a comprehensive representation of the genetic structure of the Italian Holstein breed, highlighting the presence of potential genetic subgroups due to divergent dairy farming systems. This study can be used to further investigate genetic variants underlying adaptation traits in these subgroups, which in turn might be used to design more specialized breeding programs.
Collapse
Affiliation(s)
- Michela Ablondi
- grid.10383.390000 0004 1758 0937Dipartimento di Scienze Medico-Veterinarie, University of Parma, Via del Taglio 10, 43126 Parma, Italy
| | - Massimo Malacarne
- grid.10383.390000 0004 1758 0937Dipartimento di Scienze Medico-Veterinarie, University of Parma, Via del Taglio 10, 43126 Parma, Italy
| | - Claudio Cipolat-Gotet
- grid.10383.390000 0004 1758 0937Dipartimento di Scienze Medico-Veterinarie, University of Parma, Via del Taglio 10, 43126 Parma, Italy
| | - Jan-Thijs van Kaam
- Associazione Nazionale Allevatori della Razza Frisona e Jersey Italiana, Cremona, Italy
| | - Alberto Sabbioni
- grid.10383.390000 0004 1758 0937Dipartimento di Scienze Medico-Veterinarie, University of Parma, Via del Taglio 10, 43126 Parma, Italy
| | - Andrea Summer
- grid.10383.390000 0004 1758 0937Dipartimento di Scienze Medico-Veterinarie, University of Parma, Via del Taglio 10, 43126 Parma, Italy
| |
Collapse
|
42
|
Nascimento AV, Romero ARS, Nawaz MY, Cardoso DF, Santos DJA, Gondro C, Tonhati H. An updated Axiom buffalo genotyping array map and mapping of cattle quantitative trait loci to the new water buffalo reference genome assembly. Anim Genet 2021; 52:505-508. [PMID: 34106478 DOI: 10.1111/age.13103] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/26/2021] [Indexed: 11/30/2022]
Abstract
The objectives of this study were to provide the buffalo research community with an updated SNP map for the Axiom Buffalo Genotyping (ABG) array with genomic positions for SNP currently unmapped and to map all cattle QTL from the CattleQTLdb onto the buffalo reference assembly. To update the ABG array map, all SNP probe sequences from the ABG array were re-aligned against the UOA_WB_1 assembly. With the new map, the number of mapped markers increased by approximately 10% and went from 106 778 to 116 708, which reduced the average marker spacing by approximately 2 kb. A comparison of results between signatures of autozygosity study using the ABG and the new map showed that, when the additional markers were used there was an increase in the autozygosity peaks and additional peaks in BBU5 and BBU11 could be identified. After sequence alignment and quality control, 64 650 (UMD3.1) and 76 530 (ARS_UCD1.2) cattle QTL were mapped onto the buffalo genome. The mapping of the bovine QTL database onto the buffalo genome should be useful for genome-wide association studies in buffalo and, given the high homology between the two species, the positions of cattle QTL on the buffalo genome can serve as a stepping stone towards a water buffalo QTL database.
Collapse
Affiliation(s)
- A V Nascimento
- Department of Animal Science, São Paulo State University, Jaboticabal, 14884900, Brazil
| | - A R S Romero
- Department of Animal Science, São Paulo State University, Jaboticabal, 14884900, Brazil
| | - M Y Nawaz
- Genetics and Genome Sciences Program, Michigan State University, East Lansing, MI, 48823, USA
| | - D F Cardoso
- Department of Animal Biosciences, University of Guelph, Guelph, ON, N1G 2W1, Canada
| | - D J A Santos
- Department of Animal Science, University of Maryland, College Park, MD, 20740, USA
| | - C Gondro
- Department of Animal Science, Michigan State University, East Lansing, MI, 48823, USA
| | - H Tonhati
- Department of Animal Science, São Paulo State University, Jaboticabal, 14884900, Brazil
| |
Collapse
|
43
|
Peng Q, Ehlers CL. Long tracks of homozygosity predict the severity of alcohol use disorders in an American Indian population. Mol Psychiatry 2021; 26:2200-2211. [PMID: 33398086 PMCID: PMC8254832 DOI: 10.1038/s41380-020-00989-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/13/2020] [Revised: 11/30/2020] [Accepted: 12/07/2020] [Indexed: 11/20/2022]
Abstract
Runs of homozygosity (ROH) arise when an individual inherits two copies of the same haplotype segment. While ROH are ubiquitous across human populations, Native populations-with shared parental ancestry arising from isolation and endogamy-can carry a substantial enrichment for ROH. We have been investigating genetic and environmental risk factors for alcohol use disorders (AUD) in a group of American Indians (AI) who have higher rates of AUD than the general U. S. population. Here we explore whether ROH might be associated with incidence and severity of AUD in this admixed AI population (n = 742) that live on geographically contiguous reservations, using low-coverage whole genome sequences. We have found that the genomic regions in the ROH that were identified in this population had significantly elevated American Indian heritage compared with the rest of the genome. Increased ROH abundance and ROH burden are likely risk factors for AUD severity in this AI population, especially in those diagnosed with severe and moderate AUD. The association between ROH and AUD was mostly driven by ROH of moderate lengths between 1 and 2 Mb. An ROH island on chromosome 1p32.3 and a rare ROH pool on chromosome 3p12.3 were found to be significantly associated with AUD severity. They contain genes involved in lipid metabolism, oxidative stress and inflammatory responses; and OSBPL9 was found to reside on the consensus part of the ROH island. These data demonstrate that ROH are associated with risk for AUD severity in this AI population.
Collapse
Affiliation(s)
- Qian Peng
- Department of Neuroscience, The Scripps Research Institute, La Jolla, CA, 92037, USA.
| | - Cindy L Ehlers
- Department of Neuroscience, The Scripps Research Institute, La Jolla, CA, 92037, USA.
| |
Collapse
|
44
|
Genome-Wide Patterns of Homozygosity Reveal the Conservation Status in Five Italian Goat Populations. Animals (Basel) 2021; 11:ani11061510. [PMID: 34071004 PMCID: PMC8224610 DOI: 10.3390/ani11061510] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Revised: 05/11/2021] [Accepted: 05/20/2021] [Indexed: 02/07/2023] Open
Abstract
Simple Summary In the local populations, the increase in inbreeding is a relevant problem for the reduction in production, reproduction, and adaptive traits. The application of genomic technologies has facilitated the assessment of inbreeding in these populations. The current study aims to investigate the patterns of homozygosity in five Italian local goat populations. The results showed the different selection histories and breeding schemes of these goat populations. The analysis also indicated the importance of this information to avoid future loss of diversity and to produce information for designing optimal breeding and conservation programs. Abstract The application of genomic technologies has facilitated the assessment of genomic inbreeding based on single nucleotide polymorphisms (SNPs). In this study, we computed several runs of homozygosity (ROH) parameters to investigate the patterns of homozygosity using Illumina Goat SNP50 in five Italian local populations: Argentata dell’Etna (N = 48), Derivata di Siria (N = 32), Girgentana (N = 59), Maltese (N = 16) and Messinese (N = 22). The ROH results showed well-defined differences among the populations. A total of 3687 ROH segments >2 Mb were detected in the whole sample. The Argentata dell’Etna and Messinese were the populations with the lowest mean number of ROH and inbreeding coefficient values, which reflect admixture and gene flow. In the Girgentana, we identified an ROH pattern related with recent inbreeding that can endanger the viability of the breed due to reduced population size. The genomes of Derivata di Siria and Maltese breeds showed the presence of long ROH (>16 Mb) that could seriously impact the overall biological fitness of these breeds. Moreover, the results confirmed that ROH parameters are in agreement with the known demography of these populations and highlighted the different selection histories and breeding schemes of these goat populations. In the analysis of ROH islands, we detected harbored genes involved with important traits, such as for milk yield, reproduction, and immune response, and are consistent with the phenotypic traits of the studied goat populations. Finally, the results of this study can be used for implementing conservation programs for these local populations in order to avoid further loss of genetic diversity and to preserve the production and fitness traits. In view of this, the availability of genomic data is a fundamental resource.
Collapse
|
45
|
Wu X, Zhou R, Zhang W, Cao B, Xia J, Caiyun W, Zhang X, Chu M, Yin Z, Ding Y. Genome-wide scan for runs of homozygosity identifies candidate genes in Wannan Black pigs. Anim Biosci 2021; 34:1895-1902. [PMID: 33705632 PMCID: PMC8563231 DOI: 10.5713/ab.20.0679] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2020] [Accepted: 02/07/2021] [Indexed: 11/27/2022] Open
Abstract
Objective Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that can reveal inbreeding levels, selection pressure, and mating schemes. In this study, ROHs were evaluated in Wannan Black pigs to assess the inbreeding levels and the genome regions with high ROH frequency. Methods In a previous study, we obtained 501.52 GB of raw data from resequencing (10×) of the genome and identified 21,316,754 single-nucleotide variants in 20 Wannan Black pig samples. We investigated the number, length, and frequency of ROH using resequencing data to characterize the homozygosity in Wannan Black pigs and identified genomic regions with high ROH frequencies. Results In this work, 1,813 ROHs (837 ROHs in 100 to 500 kb, 449 ROHs in 500 to 1,000 kb, 527 ROHs in >1,000 kb) were identified in all samples, and the average genomic inbreeding coefficient (FROH) in Wannan Black pigs was 0.5234. Sixty-one regions on chromosomes 2, 3, 7, 8, 13, 15, and 16 harbored ROH islands. In total, 105 genes were identified in 42 ROH islands, among which some genes were related to production traits. Conclusion This is the first study to identify ROH across the genome of Wannan Black pigs, the Chinese native breed of the Anhui province. Overall, Wannan Black pigs have high levels of inbreeding due to the influence of ancient and recent inbreeding due to the genome. These findings are a reliable resource for future studies and contribute to save and use the germplasm resources of Wannan Black pigs.
Collapse
Affiliation(s)
- Xudong Wu
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| | - Ren Zhou
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| | - Wei Zhang
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, Anhui 230031, P.R. China
| | - Bangji Cao
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| | - Jing Xia
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| | - Wang Caiyun
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| | - Xiaodong Zhang
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| | - Mingxing Chu
- Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing,100193, P. R. China
| | - Zongjun Yin
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| | - Yueyun Ding
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, P. R. China.,Anhui province key laboratory of local livestock and poultry genetic resource conservation and bio-breeding, Anhui Agricultural University, Hefei, 230036, P.R. China
| |
Collapse
|
46
|
Ganteil A, Rodriguez-Ramilo ST, Ligonesche B, Larzul C. Characterization of Autozygosity in Pigs in Three-Way Crossbreeding. Front Genet 2021; 11:584556. [PMID: 33584790 PMCID: PMC7876413 DOI: 10.3389/fgene.2020.584556] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Accepted: 12/21/2020] [Indexed: 11/16/2022] Open
Abstract
Crossbreeding in livestock can be used to increase genetic diversity. The resulting increase in variability is related to the heterozygosity of the crossbred animal. The evolution of diversity during crossbreeding can be assessed using genomic data. The objective of this study was to describe patterns of runs of homozygosity (ROH) in animals resulting from three-way crossbreeding, from parental pure lines, and in their crossbred offspring. The crossbreeding scheme consisted of a first crossbreeding Pietrain boars and Large White sows, after which the offspring of the Pietrain × Large White were crossed with Duroc boars. The offspring of the second crossbreeding are called G0, the offspring of G0 boars and G0 sows are called G1. All the animals were genotyped using the Illumina SNP60 porcine chip. After filtering, analyses were performed with 2,336 animals and 48,579 autosomal single nucleotide polymorphism (SNP). The mean ROH-based inbreeding coefficients were shown to be 0.27 ± 0.05, 0.23 ± 0.04, and 0.26 ± 0.04 for Duroc, Large White, and Pietrain, respectively. ROH were detected in the Pietrain × Large White crossbred but the homozygous segments were fewer and smaller than in their parents. Similar results were obtained in the G0 crossbred. However, in the G1 crossbreds the number and the size of ROH were higher than in G0 parents. Similar ROH hotspots were detected on SSC1, SSC4, SSC7, SSC9, SSC13, SSC14, and SSC15 in both G0 and G1 animals. Long ROH (>16 Mb) were observed in G1 animals, suggesting regions with low recombination rates. The conservation of these homozygous segments in the three crossbred populations means that some haplotypes were shared between parental breeds. Gene annotation in ROH hotspots in G0 animals identified genes related to production traits including carcass composition and reproduction. These findings advance our understanding of how to manage genetic diversity in crossbred populations.
Collapse
Affiliation(s)
- Audrey Ganteil
- GenPhySE, Université de Toulouse, INRAE, ENVT, Castanet-Tolosan, France.,SAS NUCLEUS, Le Rheu, France
| | | | | | - Catherine Larzul
- GenPhySE, Université de Toulouse, INRAE, ENVT, Castanet-Tolosan, France
| |
Collapse
|
47
|
Pavan S, Delvento C, Mazzeo R, Ricciardi F, Losciale P, Gaeta L, D'Agostino N, Taranto F, Sánchez-Pérez R, Ricciardi L, Lotti C. Almond diversity and homozygosity define structure, kinship, inbreeding, and linkage disequilibrium in cultivated germplasm, and reveal genomic associations with nut and seed weight. HORTICULTURE RESEARCH 2021; 8:15. [PMID: 33423037 PMCID: PMC7797004 DOI: 10.1038/s41438-020-00447-1] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Revised: 11/04/2020] [Accepted: 11/13/2020] [Indexed: 05/04/2023]
Abstract
Almond [Prunus dulcis Miller (D.A. Webb)] is the main tree nut species worldwide. Here, genotyping-by-sequencing (GBS) was applied to 149 almond cultivars from the ex situ collections of the Italian Council for Agricultural Research (CREA) and the Spanish National Research Council (CSIC), leading to the detection of 93,119 single-nucleotide polymorphisms (SNPs). The study of population structure outlined four distinct genetic groups and highlighted diversification between the Mediterranean and Californian gene pools. Data on SNP diversity and runs of homozygosity (ROHs) allowed the definition of kinship, inbreeding, and linkage disequilibrium (LD) decay in almond cultivated germplasm. Four-year phenotypic observations, gathered on 98 cultivars of the CREA collection, were used to perform a genome-wide association study (GWAS) and, for the first time in a crop species, homozygosity mapping (HM), resulting in the identification of genomic associations with nut, shell, and seed weight. Both GWAS and HM suggested that loci controlling nut and seed weight are mostly independent. Overall, this study provides insights on the almond cultivation history and delivers information of major interest for almond genetics and breeding. In a broader perspective, our results encourage the use of ROHs in crop science to estimate inbreeding, choose parental combinations minimizing the risk of inbreeding depression, and identify genomic footprints of selection for specific traits.
Collapse
Affiliation(s)
- Stefano Pavan
- Department of Soil, Plant and Food Science, University of Bari Aldo Moro, Via Amendola 165/A, Bari, 70126, Italy.
- Institute of Biomedical Technologies, National Research Council (CNR), Via Amendola 122/D, Bari, 70126, Italy.
| | - Chiara Delvento
- Department of Soil, Plant and Food Science, University of Bari Aldo Moro, Via Amendola 165/A, Bari, 70126, Italy
| | - Rosa Mazzeo
- Department of the Sciences of Agriculture, Food and Environment, University of Foggia, Via Napoli 25, Foggia, 71100, Italy
| | - Francesca Ricciardi
- Department of the Sciences of Agriculture, Food and Environment, University of Foggia, Via Napoli 25, Foggia, 71100, Italy
| | - Pasquale Losciale
- Department of Soil, Plant and Food Science, University of Bari Aldo Moro, Via Amendola 165/A, Bari, 70126, Italy
| | - Liliana Gaeta
- Council for Agricultural Research and Economics-Research Centre for Agriculture and Environment (CREA-AA), Bari, 70125, Italy
| | - Nunzio D'Agostino
- Department of Agricultural Sciences, University of Naples Federico II, Portici, 80055, Italy
| | - Francesca Taranto
- Institute of Biosciences and Bioresources, National Research Council of Italy, Portici, 80055, Italy
| | | | - Luigi Ricciardi
- Department of Soil, Plant and Food Science, University of Bari Aldo Moro, Via Amendola 165/A, Bari, 70126, Italy
| | - Concetta Lotti
- Department of the Sciences of Agriculture, Food and Environment, University of Foggia, Via Napoli 25, Foggia, 71100, Italy.
| |
Collapse
|
48
|
Kumar S, Deng CH, Hunt M, Kirk C, Wiedow C, Rowan D, Wu J, Brewer L. Homozygosity Mapping Reveals Population History and Trait Architecture in Self-Incompatible Pear ( Pyrus spp.). FRONTIERS IN PLANT SCIENCE 2021; 11:590846. [PMID: 33469460 PMCID: PMC7813798 DOI: 10.3389/fpls.2020.590846] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/03/2020] [Accepted: 12/04/2020] [Indexed: 06/12/2023]
Abstract
Runs of homozygosity (ROH) have been widely used to study population history and trait architecture in humans and livestock species, but their application in self-incompatible plants has not been reported. The distributions of ROH in 199 accessions representing Asian pears (45), European pears (109), and interspecific hybrids (45) were investigated using genotyping-by-sequencing in this study. Fruit phenotypes including fruit weight, firmness, Brix, titratable acidity, and flavor volatiles were measured for genotype-phenotype analyses. The average number of ROH and the average total genomic length of ROH were 6 and 11 Mb, respectively, in Asian accessions, and 13 and 30 Mb, respectively, in European accessions. Significant associations between genomic inbreeding coefficients (FROH) and phenotypes were observed for 23 out of 32 traits analyzed. An overlap between ROH islands and significant markers from genome-wide association analyses was observed. Previously published quantitative trait loci for fruit traits and disease resistances also overlapped with some of the ROH islands. A prominent ROH island at the bottom of linkage group 17 overlapped with a recombination-supressed genomic region harboring the self-incompatibility locus. The observed ROH patterns suggested that systematic breeding of European pears would have started earlier than of Asian pears. Our research suggest that FROH would serve as a novel tool for managing inbreeding in gene-banks of self-incompatible plant species. ROH mapping provides a complementary strategy to unravel the genetic architecture of complex traits, and to evaluate differential selection in outbred plants. This seminal work would provide foundation for the ROH research in self-incompatible plants.
Collapse
Affiliation(s)
- Satish Kumar
- Hawke’s Bay Research Centre, The New Zealand Institute for Plant and Food Research Limited, Havelock North, New Zealand
| | - Cecilia Hong Deng
- Mount Albert Research Centre, The New Zealand Institute for Plant and Food Research Limited, Auckland, New Zealand
| | - Martin Hunt
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Chris Kirk
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Claudia Wiedow
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Daryl Rowan
- Palmerston North Research Centre, The New Zealand Institute for Plant and Food Research Limited, Palmerston North, New Zealand
| | - Jun Wu
- Centre of Pear Engineering Technology Research, Nanjing Agricultural University, Nanjing, China
| | - Lester Brewer
- Motueka Research Centre, The New Zealand Institute for Plant and Food Research Limited, Motueka, New Zealand
| |
Collapse
|
49
|
Genome-Wide Analysis Revealed Homozygosity and Demographic History of Five Chinese Sheep Breeds Adapted to Different Environments. Genes (Basel) 2020; 11:genes11121480. [PMID: 33317115 PMCID: PMC7764688 DOI: 10.3390/genes11121480] [Citation(s) in RCA: 10] [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/30/2020] [Revised: 12/06/2020] [Accepted: 12/07/2020] [Indexed: 01/18/2023] Open
Abstract
Homozygosity of long sequence genotypes are a result of parents transmitting identical haplotypes, which can be used to estimate their auto-zygosity. Therefore, we used high-density SNP Chip data to characterize the auto-zygosity of each breed according to the occurrence and distribution of runs of homozygosity (ROH). Subsequently, we identified the genomic regions with high runs of homozygosity frequencies within individuals of each breed. We selected 96 sheep samples from five local Chinese sheep breeds belonging to different geographical locations. We identified 3046 ROHs within the study breed individuals, among which the longer segments (>1–5 Mb) were dominant. On average, ROH segments covered about 12% of the genomes; the coverage rate of OAR20 was the lowest and that of OAR2 was the highest. The distribution analysis of runs of homozygosity showed that the detected ROH mainly distributed between >26 and 28 Mb. The Hetian and Hu sheep showed the lowest ROH distribution. The estimation of homozygosity level reflects the history of modern and ancient inbreeding, which may affect the genomes of Chinese indigenous sheep breeds and indicate that some animals have experienced recent self-pollination events (Yabuyi, Karakul and Wadi). In these sheep breeds, the genomic regions were assumed to be under selection signatures frequently in line with long ROH. These regions included candidate genes associated with disease resistance traits (5S_rRNA), the innate and adaptive immune response (HERC2 and CYFIP1), digestion and metabolism (CENPJ), growth (SPP1), body size and developments (GJB2 and GJA3). This study highlighted new insights into the ROH patterns and provides a basis for future breeding and conservation strategies of Chinese sheep breeds.
Collapse
|
50
|
Islam S, Reddy UK, Natarajan P, Abburi VL, Bajwa AA, Imran M, Zahoor MY, Abdullah M, Bukhari AM, Iqbal S, Ashraf K, Nadeem A, Rehman H, Rashid I, Shehzad W. Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight. PLoS One 2020; 15:e0242500. [PMID: 33232358 PMCID: PMC7685427 DOI: 10.1371/journal.pone.0242500] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2020] [Accepted: 11/03/2020] [Indexed: 11/20/2022] Open
Abstract
The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.
Collapse
Affiliation(s)
- Saher Islam
- Institute of Biochemistry and Biotechnology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Umesh K. Reddy
- Department of Biology, West Virginia State University, Institute, West Virginia, United States of America
| | - Purushothaman Natarajan
- Department of Biology, West Virginia State University, Institute, West Virginia, United States of America
| | - Venkata Lakshmi Abburi
- Department of Biology, West Virginia State University, Institute, West Virginia, United States of America
| | - Amna Arshad Bajwa
- Institute of Biochemistry and Biotechnology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Muhammad Imran
- Institute of Biochemistry and Biotechnology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Muhammad Yasir Zahoor
- Institute of Biochemistry and Biotechnology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Muhammad Abdullah
- Department of Livestock Production, University of Veterinary and Animal Sciences, Pattoki, Pakistan
| | - Aamir Mehmood Bukhari
- Semen Production Unit, Qadirabad, District Sahiwal, Pakistan
- Livestock and Dairy Development Department, Government of the Punjab, Lahore, Pakistan
| | - Sajid Iqbal
- Semen Production Unit, Qadirabad, District Sahiwal, Pakistan
- Livestock and Dairy Development Department, Government of the Punjab, Lahore, Pakistan
| | - Kamran Ashraf
- Department of Parasitology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Asif Nadeem
- Institute of Biochemistry and Biotechnology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Habibur Rehman
- Department of Physiology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Imran Rashid
- Department of Parasitology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| | - Wasim Shehzad
- Institute of Biochemistry and Biotechnology, University of Veterinary and Animal Sciences, Lahore, Pakistan
| |
Collapse
|