1
|
Macharia JK, Kim J, Kim M, Cho E, Munyaneza JP, Lee JH. Characterisation of runs of homozygosity and inbreeding coefficients in the red-brown Korean native chickens. Anim Biosci 2024; 37:1355-1366. [PMID: 38665087 PMCID: PMC11222857 DOI: 10.5713/ab.23.0514] [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: 12/10/2023] [Revised: 02/05/2024] [Accepted: 02/27/2024] [Indexed: 07/05/2024] Open
Abstract
OBJECTIVE The analysis of runs of homozygosity (ROH) has been applied to assess the level of inbreeding and identify selection signatures in various livestock species. The objectives of this study were to characterize the ROH pattern, estimate the rate of inbreeding, and identify signatures of selection in the red-brown Korean native chickens. METHODS The Illumina 60K single nucleotide polymorphism chip data of 651 chickens was used in the analysis. Runs of homozygosity were analysed using the PLINK v1.9 software. Inbreeding coefficients were estimated using the GCTA software and their correlations were examined. Genomic regions with high levels of ROH were explored to identify selection signatures. RESULTS A total of 32,176 ROH segments were detected in this study. The majority of the ROH segments were shorter than 4 Mb. The average ROH inbreeding coefficients (FROH) varied with the length of ROH segments. The means of inbreeding coefficients calculated from different methods were also variable. The correlations between different inbreeding coefficients were positive and highly variable (r = 0.18-1). Five ROH islands harbouring important quantitative trait loci were identified. CONCLUSION This study assessed the level of inbreeding and patterns of homozygosity in Red-brown native Korean chickens. The results of this study suggest that the level of recent inbreeding is low which indicates substantial progress in the conservation of red-brown Korean native chickens. Additionally, Candidate genomic regions associated with important production traits were detected in homozygous regions.
Collapse
Affiliation(s)
- John Kariuki Macharia
- Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134,
Korea
| | - Jaewon Kim
- Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134,
Korea
| | - Minjun Kim
- Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134,
Korea
| | - Eunjin Cho
- Department of Bio-AI Convergence, Chungnam National University, Daejeon 34134,
Korea
| | - Jean Pierre Munyaneza
- Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134,
Korea
| | - Jun Heon Lee
- Division of Animal and Dairy Science, Chungnam National University, Daejeon 34134,
Korea
- Department of Bio-AI Convergence, Chungnam National University, Daejeon 34134,
Korea
| |
Collapse
|
2
|
Ramljak J, Špehar M, Ceranac D, Držaić V, Pocrnić I, Barać D, Mioč B, Širić I, Barać Z, Ivanković A, Kasap A. Genomic Characterization of Local Croatian Sheep Breeds-Effective Population Size, Inbreeding & Signatures of Selection. Animals (Basel) 2024; 14:1928. [PMID: 38998043 PMCID: PMC11240672 DOI: 10.3390/ani14131928] [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: 06/05/2024] [Revised: 06/21/2024] [Accepted: 06/24/2024] [Indexed: 07/14/2024] Open
Abstract
The Istrian (IS) and the Pag sheep (PS) are local Croatian breeds which provide significant income for the regional economy and have a cultural and traditional importance for the inhabitants. The aim of this study was to estimate some important population specific genetic parameters in IS (N = 1293) and PS (N = 2637) based on genome wide SNPs. Estimates of linkage disequilibrium effective population size (Ne) evidenced more genetic variability in PS (Ne = 838) compared to IS (Ne = 197), regardless of historical time (both recent and ancient genetic variability). The discrepancy in the recent genetic variability between these breeds was additionally confirmed by the estimates of genomic inbreeding (FROH), which was estimated to be notably higher in IS (FROH>2 = 0.062) than in PS (FROH>2 = 0.029). The average FROH2-4, FROH4-8, FROH8-16, and FROH>16 were 0.26, 1.65, 2.14, and 3.72 for IS and 0.22, 0.61, 0.75, and 1.58 for PS, thus evidencing a high contribution of recent inbreeding in the overall inbreeding. One ROH island with > 30% of SNP incidence in ROHs was detected in IS (OAR6; 34,253,440-38,238,124 bp) while there was no ROH islands detected in PS. Seven genes (CCSER1, HERC3, LCORL, NAP1L5, PKD2, PYURF, and SPP1) involved in growth, feed intake, milk production, immune responses, and resistance were associated with the found autozygosity. The results of this study represent the first comprehensive insight into genomic variability of these two Croatian local sheep breeds and will serve as a baseline for setting up the most promising strategy of genomic Optimum Contribution Selection.
Collapse
Affiliation(s)
- Jelena Ramljak
- Faculty of Agriculture, University of Zagreb, 10000 Zagreb, Croatia; (V.D.); (B.M.); (I.Š.); (A.I.); (A.K.)
| | - Marija Špehar
- Croatian Agency for Agriculture and Food, 10000 Zagreb, Croatia; (M.Š.); (D.C.); (D.B.)
| | - Dora Ceranac
- Croatian Agency for Agriculture and Food, 10000 Zagreb, Croatia; (M.Š.); (D.C.); (D.B.)
| | - Valentino Držaić
- Faculty of Agriculture, University of Zagreb, 10000 Zagreb, Croatia; (V.D.); (B.M.); (I.Š.); (A.I.); (A.K.)
| | - Ivan Pocrnić
- The Roslin Institute, University of Edinburgh, Easter Bush Campus, Midlothian EH25 9RG, UK;
| | - Dolores Barać
- Croatian Agency for Agriculture and Food, 10000 Zagreb, Croatia; (M.Š.); (D.C.); (D.B.)
| | - Boro Mioč
- Faculty of Agriculture, University of Zagreb, 10000 Zagreb, Croatia; (V.D.); (B.M.); (I.Š.); (A.I.); (A.K.)
| | - Ivan Širić
- Faculty of Agriculture, University of Zagreb, 10000 Zagreb, Croatia; (V.D.); (B.M.); (I.Š.); (A.I.); (A.K.)
| | | | - Ante Ivanković
- Faculty of Agriculture, University of Zagreb, 10000 Zagreb, Croatia; (V.D.); (B.M.); (I.Š.); (A.I.); (A.K.)
| | - Ante Kasap
- Faculty of Agriculture, University of Zagreb, 10000 Zagreb, Croatia; (V.D.); (B.M.); (I.Š.); (A.I.); (A.K.)
| |
Collapse
|
3
|
Vostry L, Vostra-Vydrova H, Moravcikova N, Kasarda R, Margetin M, Rychtarova J, Drzaic I, Shihabi M, Cubric-Curik V, Sölkner J, Curik I. Genomic analysis of conservation status, population structure and admixture in local Czech and Slovak dairy goat breeds. J Dairy Sci 2024:S0022-0302(24)00937-8. [PMID: 38908686 DOI: 10.3168/jds.2023-24607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2023] [Accepted: 05/30/2024] [Indexed: 06/24/2024]
Abstract
While dairy goat production, characterized by traditional production on small farms, is an important source of income in the Czech Republic and Slovakia, locally adapted breeds have not been fully consolidated over the last 100 years due to large fluctuations in population size and inconsistent breeding programs that allowed for different crossbreeding strategies. Our main objective in this study was therefore to assess the conservation status of 4 Czech (Alpine Goat, White Shorthair, Brown Shorthair and Czech Landrace) and one Slovak (Slovak White Shorthair) local goat breeds, to analyze their population structure and admixture, and to estimate their relatedness to several neighboring breeds. Our analyses included 142 goats belonging to 5 local breeds genotyped with the Illumina 50K BeadChip and 618 previously genotyped animals representing 15 goat breeds from Austria and Switzerland (all analyses based on 46,862 autosomal SNPs and 760 animals). In general, the conservation status of the Czech and Slovak local goat breeds was satisfactory, with the exception of the Brown Shorthair goat, as the analyzed parameters (heterozygosity, haplotype richness, ROH-based inbreeding and effective population size) were mostly above the median of 20 breeds. However, for all 5 Czech and Slovakian breeds, an examination of historical effective population size indicated a substantial decline about 8 to 22 generations ago. In addition, our study revealed that the Czech and Slovakian breeds are not fully consolidated; for instance, White Shorthair and Brown Shorthair were not clearly distinguishable. Considerable admixture, especially in Czech Landrace (effective number of parental clusters equal to 4.2), and low but numerous migration rates from other Austrian and Swiss breeds were found. These results provide valuable insights for future breeding programs and genetic diversity management of local Czech and Slovak goat breeds.
Collapse
Affiliation(s)
- Lubos Vostry
- Czech University of Life Science Prague, Kamycka 129, 16500 Prague, Czech Republic.
| | - Hana Vostra-Vydrova
- Czech University of Life Science Prague, Kamycka 129, 16500 Prague, Czech Republic
| | - Nina Moravcikova
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Radovan Kasarda
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Milan Margetin
- Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovak Republic
| | - Jana Rychtarova
- Institute of Animal Science, Přátelství 815, 104 00 Prague, Czech Republic
| | - Ivana Drzaic
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Mario Shihabi
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Vlatka Cubric-Curik
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia
| | - Johan Sölkner
- University of Natural Resources & Life Sciences Vienna, Gregor-Mendel-Strasse 33, 1180 Vienna, Austria
| | - Ino Curik
- University of Zagreb, Faculty of Agriculture, Svetošimunska cesta 25, 10000 Zagreb, Croatia; Institute of Animal Sciences, Hungarian University of Agriculture and Life Sciences (MATE), Guba Sándor u. 40, 7400 Kaposvár, Hungary.
| |
Collapse
|
4
|
Giovannini S, Chessari G, Riggio S, Marletta D, Sardina MT, Mastrangelo S, Sarti FM. Insight into the current genomic diversity, conservation status and population structure of Tunisian Barbarine sheep breed. Front Genet 2024; 15:1379086. [PMID: 38881792 PMCID: PMC11176520 DOI: 10.3389/fgene.2024.1379086] [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/30/2024] [Accepted: 05/07/2024] [Indexed: 06/18/2024] Open
Abstract
Local livestock breeds play a crucial role in global biodiversity, connecting natural and human-influenced environments and contributing significantly to ecosystem services. While commercial breeds dominate industrial systems, local livestock breeds in developing countries, like Barbarine sheep in Tunisia, are vital for food security and community maintenance. The Tunisian Barbarine sheep, known for its adaptability and distinctive fat-tailed morphology, faces challenges due to historical crossbreeding. In this study, the Illumina Ovine SNP50K BeadChip array was used to perform a genome-wide characterization of Tunisian Barbarine sheep to investigate its genetic diversity, the genome structure, and the relationship within the context of Mediterranean breeds. The results show moderate genetic diversity and low inbreeding. Runs of Homozygosity analysis find genomic regions linked to important traits, including fat tail characteristics. Genomic relationship analysis shows proximity to Algerian thin-tailed breeds, suggesting crossbreeding impacts. Admixture analysis reveals unique genetic patterns, emphasizing the Tunisian Barbarine's identity within the Mediterranean context and its closeness to African breeds. Current results represent a starting point for the creation of monitoring and conservation plans. In summary, despite genetic dilution due to crossbreeding, the identification of genomic regions offers crucial insights for conservation. The study confirms the importance of preserving unique genetic characteristics of local breeds, particularly in the face of ongoing crossbreeding practices and environmental challenges. These findings contribute valuable insights for the sustainable management of this unique genetic reservoir, supporting local economies and preserving sheep species biodiversity.
Collapse
Affiliation(s)
- Samira Giovannini
- Dipartimento di Scienze Agrarie, Alimentari e Ambientali, University of Perugia, Perugia, Italy
| | - Giorgio Chessari
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Silvia Riggio
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Donata Marletta
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Maria Teresa Sardina
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Francesca Maria Sarti
- Dipartimento di Scienze Agrarie, Alimentari e Ambientali, University of Perugia, Perugia, Italy
| |
Collapse
|
5
|
Abdelmanova AA, Deniskova TE, Kharzinova VR, Chinarov RY, Boronetskaya OI, Sölkner J, Brem G, Ai H, Huang L, Trukhachev VI, Zinovieva NA. Tracing the Dynamical Genetic Diversity Changes of Russian Livni Pigs during the Last 50 Years with the Museum, Old, and Modern Samples. Animals (Basel) 2024; 14:1629. [PMID: 38891676 PMCID: PMC11171240 DOI: 10.3390/ani14111629] [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: 04/24/2024] [Revised: 05/27/2024] [Accepted: 05/27/2024] [Indexed: 06/21/2024] Open
Abstract
The pig industry is usually considered an intensive livestock industry, mainly supported by hybrid breeding between commercial pig breeds. However, people's pursuit of a more natural environment and higher meat quality has led to an increasing demand for eco-friendly and diverse pig feeding systems. Therefore, the importance of rearing and conserving local pig breeds is increasing. The Livni pig is a local breed with good adaptability to the environmental and fodder conditions in central Russia. In this study, we aimed to analyze the genetic diversity and population structure of Livni pigs using whole-genome single nucleotide polymorphism (SNP) data. We utilized the Porcine GGP HD BeadChip on genotype samples from old (n = 32, 2004) and modern (n = 32, 2019) populations of Livni pigs. For the museum samples of Livni pigs (n = 3), we extracted DNA from their teeth, performed genomic sequencing, and obtained SNP genotypes from the whole-genome sequences. SNP genotypes of Landrace (n = 32) and Large White (n = 32) pigs were included for comparative analysis. We observed that the allelic richness of Livni pigs was higher than those of Landrace and Large White pigs (AR = 1.775-1.798 vs. 1.703 and 1.668, respectively). The effective population size estimates (NE5 = 108 for Livni pigs, NE5 = 59 for Landrace and Large White pigs) confirmed their genetic diversity tendency. This was further supported by the length and number of runs of homozygosity, as well as the genomic inbreeding coefficient (almost twofold lower in Livni pigs compared to Landrace and Large White pigs). These findings suggest that the Livni pig population exhibits higher genetic diversity and experiences lower selection pressure compared to commercial pig populations. Furthermore, both principal component and network tree analyses demonstrated a clear differentiation between Livni pigs and transboundary commercial pigs. The TreeMix results indicated gene flow from Landrace ancestors to Livni pigs (2019) and from Large White ancestors to Livni pigs (2004), which was consistent with their respective historical breeding backgrounds. The comparative analysis of museum, old, and modern Livni pigs indicated that the modern Livni pig populations have preserved their historical genomic components, suggesting their potential suitability for future design selection programs.
Collapse
Affiliation(s)
- Alexandra A. Abdelmanova
- L.K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk Municipal District, Moscow Region, Podolsk 142132, Russia; (A.A.A.); (V.R.K.); (R.Y.C.)
| | - Tatiana E. Deniskova
- L.K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk Municipal District, Moscow Region, Podolsk 142132, Russia; (A.A.A.); (V.R.K.); (R.Y.C.)
| | - Veronika R. Kharzinova
- L.K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk Municipal District, Moscow Region, Podolsk 142132, Russia; (A.A.A.); (V.R.K.); (R.Y.C.)
| | - Roman Yu Chinarov
- L.K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk Municipal District, Moscow Region, Podolsk 142132, Russia; (A.A.A.); (V.R.K.); (R.Y.C.)
| | - Oksana I. Boronetskaya
- Museum of Livestock, Timiryazev Russian State Agrarian University—Moscow Agrarian Academy, 49, ul. Timiryazevskaya, Moscow 127550, Russia; (O.I.B.); (V.I.T.)
| | - Johann Sölkner
- Division of Livestock Sciences, University of Natural Resources and Life Sciences, 1180 Vienna, Austria;
| | - Gottfried Brem
- Institute of Animal Breeding and Genetics, University of Veterinary Medicine (VMU), Veterinärplatz, 1210 Vienna, Austria;
| | - Huashui Ai
- National Key Laboratory for Swine Genetic Improvement and Germplasm Innovation, Ministry of Science and Technology of China, Jiangxi Agricultural University, Nanchang 330045, China; (H.A.); (L.H.)
| | - Lusheng Huang
- National Key Laboratory for Swine Genetic Improvement and Germplasm Innovation, Ministry of Science and Technology of China, Jiangxi Agricultural University, Nanchang 330045, China; (H.A.); (L.H.)
| | - Vladimir I. Trukhachev
- Museum of Livestock, Timiryazev Russian State Agrarian University—Moscow Agrarian Academy, 49, ul. Timiryazevskaya, Moscow 127550, Russia; (O.I.B.); (V.I.T.)
| | - Natalia A. Zinovieva
- L.K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk Municipal District, Moscow Region, Podolsk 142132, Russia; (A.A.A.); (V.R.K.); (R.Y.C.)
| |
Collapse
|
6
|
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
|
7
|
Ristanic M, Zorc M, Glavinic U, Stevanovic J, Blagojevic J, Maletic M, Stanimirovic Z. Genome-Wide Analysis of Milk Production Traits and Selection Signatures in Serbian Holstein-Friesian Cattle. Animals (Basel) 2024; 14:669. [PMID: 38473054 DOI: 10.3390/ani14050669] [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: 01/05/2024] [Revised: 02/05/2024] [Accepted: 02/14/2024] [Indexed: 03/14/2024] Open
Abstract
To improve the genomic evaluation of milk-related traits in Holstein-Friesian (HF) cattle it is essential to identify the associated candidate genes. Novel SNP-based analyses, such as the genetic mapping of inherited diseases, GWAS, and genomic selection, have led to a new era of research. The aim of this study was to analyze the association of each individual SNP in Serbian HF cattle with milk production traits and inbreeding levels. The SNP 60 K chip Axiom Bovine BovMDv3 was deployed for the genotyping of 334 HF cows. The obtained genomic results, together with the collected phenotypic data, were used for a GWAS. Moreover, the identification of ROH segments was performed and served for inbreeding coefficient evaluation and ROH island detection. Using a GWAS, a polymorphism, rs110619097 (located in the intron of the CTNNA3 gene), was detected to be significantly (p < 0.01) associated with the milk protein concentration in the first lactation (adjusted to 305 days). The average genomic inbreeding value (FROH) was 0.079. ROH islands were discovered in proximity to genes associated with milk production traits and genomic regions under selection pressure for other economically important traits of dairy cattle. The findings of this pilot study provide useful information for a better understanding of the genetic architecture of milk production traits in Serbian HF dairy cows and can be used to improve lactation performances in Serbian HF cattle breeding programs.
Collapse
Affiliation(s)
- Marko Ristanic
- Department of Biology, Faculty of Veterinary Medicine, University of Belgrade, Bul. Oslobodjenja 18, 11000 Belgrade, Serbia
| | - Minja Zorc
- Department of Animal Science, Biotechnical Faculty, University of Ljubljana, Groblje 3, 1000 Ljubljana, Slovenia
| | - Uros Glavinic
- Department of Biology, Faculty of Veterinary Medicine, University of Belgrade, Bul. Oslobodjenja 18, 11000 Belgrade, Serbia
| | - Jevrosima Stevanovic
- Department of Biology, Faculty of Veterinary Medicine, University of Belgrade, Bul. Oslobodjenja 18, 11000 Belgrade, Serbia
| | - Jovan Blagojevic
- Department of Biology, Faculty of Veterinary Medicine, University of Belgrade, Bul. Oslobodjenja 18, 11000 Belgrade, Serbia
| | - Milan Maletic
- Department of Reproduction, Fertility and Artificial Insemination, Faculty of Veterinary Medicine, University of Belgrade, Bul. Oslobodjenja 18, 11000 Belgrade, Serbia
| | - Zoran Stanimirovic
- Department of Biology, Faculty of Veterinary Medicine, University of Belgrade, Bul. Oslobodjenja 18, 11000 Belgrade, Serbia
| |
Collapse
|
8
|
Dementieva NV, Shcherbakov YS, Ryabova AE, Vakhrameev AB, Makarova AV, Nikolaeva OA, Dysin AP, Azovtseva AI, Reinbah NR, Mitrofanova OV. Comparative peculiarities of genomic diversity in Gallus gallus domesticus chickens with decorative plumage: the muffs and beard phenotype. Vavilovskii Zhurnal Genet Selektsii 2024; 28:108-116. [PMID: 38465249 PMCID: PMC10917671 DOI: 10.18699/vjgb-24-13] [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: 03/28/2023] [Revised: 09/25/2023] [Accepted: 09/26/2023] [Indexed: 03/12/2024] Open
Abstract
Throughout history, humans have been attempting to develop the ornamental features of domestic animals in addition to their productive qualities. Many chicken breeds have developed tufts of elongated feathers that jut out from the sides and bottom of the beak, leading to the phenotype known as muffs and beard. It is an incomplete autosomal dominant phenotype determined by the Mb locus localised on chromosome GGA27. This project aimed to analyse the genetic diversity of chicken breeds using full genomic genotyping with the Chicken 60K BeadChip. A total of 53,313 Single Nucleotide Polymorphisms were analysed. DNA was obtained from breeds with the muffs and beard as a marker phenotype: Faverolles (n = 20), Ukrainian Muffed (n = 18), Orloff (n = 20), Novopavlov White (n = 20), and Novopavlov Coloured (n = 15). The Russian White (n = 20) was selected as an alternative breed without the muffs and beard phenotype. The chickens are owned by the Centre of Collective Use "Genetic Collection of Rare and Endangered Breeds of Chickens" (St. Petersburg region, Pushkin), and are also included in the Core Shared Research Facility (CSRF) and/or Large-Scale Research Facility (LSRF). Multidimensional scaling revealed that the Novopavlov White and the Novopavlov Coloured populations formed a separate group. The Ukrainian Muffed and the Orloff have also been combined into a separate group. Based on cluster analysis, with the cross-validation error and the most probable number of clusters K = 4 taken into account, the Orloff was singled out as a separate group. The Ukrainian Muffed exhibited a notable similarity with the Orloff under the same conditions. At K = 5, the populations of the Novopavlov White and the Novopavlov Coloured diverged. Only at K = 6, a distinct and separate cluster was formed by the Ukrainian Muffed. The Russian White had the greatest number of short (1-2 Mb) homozygous regions. If the HOXB8 gene is located between 3.402 and 3.404 Mb on chromosome GGA27, homozygous regions are rarely found in the chickens with the muffs and beard phenotype. Scanning the chicken genome with the Chicken 60K BeadChip provided enough information about the genetic diversity of the chicken breeds for the peculiarities of the development of the ornamental muffs and beard phenotypes in them to be understood. For example, Phoenix bantams, whose tail feathers grow throughout their lives, require greater consideration of husbandry conditions.
Collapse
Affiliation(s)
- N V Dementieva
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - Y S Shcherbakov
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - A E Ryabova
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - A B Vakhrameev
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - A V Makarova
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - O A Nikolaeva
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - A P Dysin
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - A I Azovtseva
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - N R Reinbah
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| | - O V Mitrofanova
- Russian Research Institute of Farm Animal Genetics and Breeding - Branch of the L.K. Ernst Federal Research Center for Animal Husbandry, Tyarlevo, St. Petersburg, Russia
| |
Collapse
|
9
|
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
|
10
|
Ojeda-Marín C, Cervantes I, Formoso-Rafferty N, Gutiérrez JP. Genomic inbreeding measures applied to a population of mice divergently selected for birth weight environmental variance. Front Genet 2023; 14:1303748. [PMID: 38155710 PMCID: PMC10752941 DOI: 10.3389/fgene.2023.1303748] [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/2023] [Accepted: 11/27/2023] [Indexed: 12/30/2023] Open
Abstract
This study aimed to compare different inbreeding measures estimated from pedigree and molecular data from two divergent mouse lines selected for environmental birth weight during 26 generations. Furthermore, the performance of different approaches and both molecular and pedigree data sources for estimating Ne were tested in this population. A total of 1,699 individuals were genotyped using a high-density genotyping array. Genomic relationship matrices were used to calculate molecular inbreeding: Nejati-Javaremi (F NEJ), Li and Horvitz (F L&H), Van Raden method 1 (F VR1) and method 2 (F VR2), and Yang (F YAN). Inbreeding based on runs of homozygosity (F ROH) and pedigree inbreeding (F PED) were also computed. F ROH, F NEJ, and F L&H were also adjusted for their average values in the first generation of selection and named F ROH0, F NEJ0, and F L&H0. ∆F was calculated from pedigrees as the individual inbreeding rate between the individual and his parents (∆F PEDt) and individual increases in inbreeding (∆F PEDi). Moreover, individual ∆F was calculated from the different molecular inbreeding coefficients (∆F NEJ0, ∆F L&H, ∆F L&H0, ∆F VR1, ∆F VR2, ∆F YAN, and ∆F ROH0). The Ne was obtained from different ∆F, such as Ne PEDt, Ne PEDi, Ne NEJ0, Ne L&H, Ne L&H0, Ne VR1, Ne VR2, Ne YAN, and Ne ROH0. Comparing with F PED , F ROH , F NEJ and F VR2 overestimated inbreeding while F NEJ0 , F L&H , F L&H0 , F VR1 and F YAN underestimated inbreeding. Correlations between inbreeding coefficients and ∆F were calculated. F ROH had the highest correlation with F PED (0.89); F YAN had correlations >0.95 with all the other molecular inbreeding coefficients. Ne PEDi was more reliable than Ne PEDt and presented similar behaviour to Ne L&H0 and Ne NEJ0. Stable trends in Ne were not observed until the 10th generation. In the 10th generation Ne PEDi was 42.20, Ne L&H0 was 45.04 and Ne NEJ0 was 45.05 and in the last generation these Ne were 35.65, 35.94 and 35.93, respectively F ROH presented the highest correlation with F PED, which addresses the identity by descent probability (IBD). The evolution of Ne L&H0 and Ne NEJ0 was the most similar to that of Ne PEDi. Data from several generations was necessary to reach a stable trend for Ne, both with pedigree and molecular data. This population was useful to test different approaches to computing inbreeding coefficients and Ne using molecular and pedigree data.
Collapse
Affiliation(s)
- Candela Ojeda-Marín
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
| | - Isabel Cervantes
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
| | - Nora Formoso-Rafferty
- Departamento de Producción Agraria, E.T.S. Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Madrid, Spain
| | - Juan Pablo Gutiérrez
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
| |
Collapse
|
11
|
Lukic B, Curik I, Drzaic I, Galić V, Shihabi M, Vostry L, Cubric-Curik V. Genomic signatures of selection, local adaptation and production type characterisation of East Adriatic sheep breeds. J Anim Sci Biotechnol 2023; 14:142. [PMID: 37932811 PMCID: PMC10626677 DOI: 10.1186/s40104-023-00936-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Accepted: 09/04/2023] [Indexed: 11/08/2023] Open
Abstract
BACKGROUND The importance of sheep breeding in the Mediterranean part of the eastern Adriatic has a long tradition since its arrival during the Neolithic migrations. Sheep production system is extensive and generally carried out in traditional systems without intensive systematic breeding programmes for high uniform trait production (carcass, wool and milk yield). Therefore, eight indigenous Croatian sheep breeds from eastern Adriatic treated here as metapopulation (EAS), are generally considered as multipurpose breeds (milk, meat and wool), not specialised for a particular type of production, but known for their robustness and resistance to certain environmental conditions. Our objective was to identify genomic regions and genes that exhibit patterns of positive selection signatures, decipher their biological and productive functionality, and provide a "genomic" characterization of EAS adaptation and determine its production type. RESULTS We identified positive selection signatures in EAS using several methods based on reduced local variation, linkage disequilibrium and site frequency spectrum (eROHi, iHS, nSL and CLR). Our analyses identified numerous genomic regions and genes (e.g., desmosomal cadherin and desmoglein gene families) associated with environmental adaptation and economically important traits. Most candidate genes were related to meat/production and health/immune response traits, while some of the candidate genes discovered were important for domestication and evolutionary processes (e.g., HOXa gene family and FSIP2). These results were also confirmed by GO and QTL enrichment analysis. CONCLUSIONS Our results contribute to a better understanding of the unique adaptive genetic architecture of EAS and define its productive type, ultimately providing a new opportunity for future breeding programmes. At the same time, the numerous genes identified will improve our understanding of ruminant (sheep) robustness and resistance in the harsh and specific Mediterranean environment.
Collapse
Affiliation(s)
- Boris Lukic
- Faculty of Agrobiotechnical Sciences Osijek, J.J, Strossmayer University of Osijek, Vladimira Preloga 1, 31000, Osijek, Croatia.
| | - Ino Curik
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia.
| | - Ivana Drzaic
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia
| | - Vlatko Galić
- Department of Maize Breeding and Genetics, Agricultural Institute Osijek, Južno predgrađe 17, 31000, Osijek, Croatia
| | - Mario Shihabi
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia
| | - Luboš Vostry
- Czech University of Life Sciences Prague, Kamýcká 129, 165 00, Praque, Czech Republic
| | - Vlatka Cubric-Curik
- Department of Animal Science, Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000, Zagreb, Croatia
| |
Collapse
|
12
|
Arias KD, Gutiérrez JP, Fernández I, Álvarez I, Goyache F. Approaching autozygosity in a small pedigree of Gochu Asturcelta pigs. Genet Sel Evol 2023; 55:74. [PMID: 37880572 PMCID: PMC10601182 DOI: 10.1186/s12711-023-00846-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2023] [Accepted: 10/03/2023] [Indexed: 10/27/2023] Open
Abstract
BACKGROUND In spite of the availability of single nucleotide polymorphism (SNP) array data, differentiation between observed homozygosity and that caused by mating between relatives (autozygosity) introduces major difficulties. Homozygosity estimators show large variation due to different causes, namely, Mendelian sampling, population structure, and differences among chromosomes. Therefore, the ascertainment of how inbreeding is reflected in the genome is still an issue. The aim of this research was to study the usefulness of genomic information for the assessment of genetic diversity in the highly endangered Gochu Asturcelta pig breed. Pedigree depth varied from 0 (founders) to 4 equivalent discrete generations (t). Four homozygosity parameters (runs of homozygosity, FROH; heterozygosity-rich regions, FHRR; Li and Horvitz's, FLH; and Yang and colleague's FYAN) were computed for each individual, adjusted for the variability in the base population (BP; six individuals) and further jackknifed over autosomes. Individual increases in homozygosity (depending on t) and increases in pairwise homozygosity (i.e., increase in the parents' mean) were computed for each individual in the pedigree, and effective population size (Ne) was computed for five subpopulations (cohorts). Genealogical parameters (individual inbreeding, individual increase in inbreeding, and Ne) were used for comparisons. RESULTS The mean F was 0.120 ± 0.074 and the mean BP-adjusted homozygosity ranged from 0.099 ± 0.081 (FLH) to 0.152 ± 0.075 (FYAN). After jackknifing, the mean values were slightly lower. The increase in pairwise homozygosity tended to be twofold higher than the corresponding individual increase in homozygosity values. When compared with genealogical estimates, estimates of Ne obtained using FYAN tended to have low root-mean-squared errors. However, Ne estimates based on increases in pairwise homozygosity using both FROH and FHRR estimates of genomic inbreeding had lower root-mean-squared errors. CONCLUSIONS Parameters characterizing homozygosity may not accurately depict losses of variability in small populations in which breeding policy prohibits matings between close relatives. After BP adjustment, the performance of FROH and FHRR was highly consistent. Assuming that an increase in homozygosity depends only on pedigree depth can lead to underestimating it in populations with shallow pedigrees. An increase in pairwise homozygosity computed from either FROH or FHRR is a promising approach for characterizing autozygosity.
Collapse
Affiliation(s)
- Katherine D Arias
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain
| | - Juan Pablo Gutiérrez
- Departamento de Producción Animal, Universidad Complutense de Madrid, Avda. Puerta de Hierro S/N, 28040, Madrid, Spain
| | - Iván Fernández
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain
| | - Isabel Álvarez
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain
| | - Félix Goyache
- Área de Genética y Reproducción Animal, SERIDA-Deva, Camino de Rioseco 1225, 33394, Gijón, Spain.
| |
Collapse
|
13
|
Sallam AM, Reyer H, Wimmers K, Bertolini F, Aboul-Naga A, Braz CU, Rabee AE. Genome-wide landscape of runs of homozygosity and differentiation across Egyptian goat breeds. BMC Genomics 2023; 24:573. [PMID: 37752425 PMCID: PMC10521497 DOI: 10.1186/s12864-023-09679-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 09/15/2023] [Indexed: 09/28/2023] Open
Abstract
Understanding the genomic features of livestock is essential for successful breeding programs and conservation. This information is scarce for local goat breeds in Egypt. In the current study, genomic regions with selection signatures were identified as well as runs of homozygosity (ROH), genomic inbreeding coefficients (FROH) and fixation index (FST) were detected in Egyptian Nubian, Damascus, Barki and Boer goat breeds. A total of 46,268 SNP markers and 337 animals were available for the genomic analyses. On average, 145.44, 42.02, 87.90 and 126.95 ROHs were detected per individual in the autosomal genome of the respective breeds. The mean accumulative ROH lengths ranged from 46.5 Mb in Damascus to 360 Mb in Egyptian Nubian. The short ROH segments (< 2 Mb) were most frequent in all breeds, while the longest ROH segments (> 16 Mb) were exclusively found in the Egyptian Nubian. The highest average FROH was observed in Egyptian Nubian (~ 0.12) followed by Boer (~ 0.11), while the lowest FROH was found in Damascus (~ 0.05) and Barki breed (~ 0.03). The estimated mean FST was 0.14 (Egyptian Nubian and Boer), 0.077 (Egyptian Nubian and Barki), 0.075 (Egyptian Nubian and Damascus), 0.071 (Barki and Boer), 0.064 (Damascus and Boer), and 0.015 (Damascus and Barki), for each pair of breeds. Interestingly, multiple SNPs that accounted for high FST values were observed on chromosome 6 in regions harboring ALPK1 and KCNIP4. Genomic regions overlapping both FST and ROH harbor genes related to immunity (IL4R, PHF23, GABARAP, GPS2, and CD68), reproduction (SPATA2L, TNFSF12, TMEM95, and RNF17), embryonic development (TCF25 and SOX15) and adaptation (MC1R, KDR, and KIT), suggesting potential genetic adaptations to local environmental conditions. Our results contribute to the understanding of the genetic architecture of different goat breeds and may provide valuable information for effective preservation and breeding programs of local goat breeds in Egypt.
Collapse
Affiliation(s)
- Ahmed M Sallam
- Animal and Poultry Breeding Department, Desert Research Center, Cairo, Egypt.
| | - Henry Reyer
- Research Institute for Farm Animal Biology (FBN), Wilhelm-Stahl-Allee 2, 18196, Dummerstorf, Germany
| | - Klaus Wimmers
- Research Institute for Farm Animal Biology (FBN), Wilhelm-Stahl-Allee 2, 18196, Dummerstorf, Germany
- Faculty of Agricultural and Environmental Sciences, University of Rostock, Justus-von-Liebig-Weg 6b, 18059, Rostock, Germany
| | - Francesca Bertolini
- Department of Agricultural and Food Sciences, University of Bologna, Bologna, Italy
| | - Adel Aboul-Naga
- Animal Production Research Institute, Agricultural Research Center, Dokki, Cairo, Egypt
| | - Camila U Braz
- Animal and Poultry Nutrition Department, Desert Research Center, Cairo, Egypt
| | - Alaa Emara Rabee
- Department of Animal Sciences, University of Illinois Urbana-Champaign, 1207 Gregory Dr, Urbana, 61801, USA
| |
Collapse
|
14
|
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
|
15
|
Bhardwaj A, Tandon G, Pal Y, Sharma NK, Nayan V, Soni S, Iquebal MA, Jaiswal S, Legha RA, Talluri TR, Bhattacharya TK, Kumar D, Rai A, Tripathi BN. Genome-Wide Single-Nucleotide Polymorphism-Based Genomic Diversity and Runs of Homozygosity for Selection Signatures in Equine Breeds. Genes (Basel) 2023; 14:1623. [PMID: 37628674 PMCID: PMC10454598 DOI: 10.3390/genes14081623] [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: 03/28/2023] [Revised: 07/31/2023] [Accepted: 08/01/2023] [Indexed: 08/27/2023] Open
Abstract
The horse, one of the most domesticated animals, has been used for several purposes, like transportation, hunting, in sport, or for agriculture-related works. Kathiawari, Marwari, Manipuri, Zanskari, Bhutia, Spiti, and Thoroughbred are the main breeds of horses, particularly due to their agroclimatic adaptation and role in any kind of strong physical activity, and these characteristics are majorly governed by genetic factors. The genetic diversity and phylogenetic relationship of these Indian equine breeds using microsatellite markers have been reported, but further studies exploring the SNP diversity and runs of homozygosity revealing the selection signature of breeds are still warranted. In our study, the identification of genes that play a vital role in muscle development is performed through SNP detection via the whole-genome sequencing approach. A total of 96 samples, categorized under seven breeds, and 620,721 SNPs were considered to ascertain the ROH patterns amongst all the seven breeds. Over 5444 ROH islands were mined, and the maximum number of ROHs was found to be present in Zanskari, while Thoroughbred was confined to the lowest number of ROHs. Gene enrichment of these ROH islands produced 6757 functional genes, with AGPAT1, CLEC4, and CFAP20 as important gene families. However, QTL annotation revealed that the maximum QTLs were associated with Wither's height trait ontology that falls under the growth trait in all seven breeds. An Equine SNP marker database (EqSNPDb) was developed to catalogue ROHs for all these equine breeds for the flexible and easy chromosome-wise retrieval of ROH along with the genotype details of all the SNPs. Such a study can reveal breed divergence in different climatic and ecological conditions.
Collapse
Affiliation(s)
- Anuradha Bhardwaj
- ICAR-National Research Centre on Equines, Sirsa Road, Hisar 125001, India; (Y.P.)
| | - Gitanjali Tandon
- Division of Agricultural Bioinformatics, ICAR-Indian Agricultural Statistics Research Institute, New Delhi 110012, India; (G.T.); (N.K.S.); (D.K.)
| | - Yash Pal
- ICAR-National Research Centre on Equines, Sirsa Road, Hisar 125001, India; (Y.P.)
| | - Nitesh Kumar Sharma
- Division of Agricultural Bioinformatics, ICAR-Indian Agricultural Statistics Research Institute, New Delhi 110012, India; (G.T.); (N.K.S.); (D.K.)
| | - Varij Nayan
- ICAR-Central Institute for Research on Buffaloes, Hisar 125001, India;
| | - Sonali Soni
- ICAR-National Research Centre on Equines, Sirsa Road, Hisar 125001, India; (Y.P.)
| | - Mir Asif Iquebal
- Division of Agricultural Bioinformatics, ICAR-Indian Agricultural Statistics Research Institute, New Delhi 110012, India; (G.T.); (N.K.S.); (D.K.)
| | - Sarika Jaiswal
- Division of Agricultural Bioinformatics, ICAR-Indian Agricultural Statistics Research Institute, New Delhi 110012, India; (G.T.); (N.K.S.); (D.K.)
| | - Ram Avatar Legha
- ICAR-National Research Centre on Equines, Sirsa Road, Hisar 125001, India; (Y.P.)
| | | | | | - Dinesh Kumar
- Division of Agricultural Bioinformatics, ICAR-Indian Agricultural Statistics Research Institute, New Delhi 110012, India; (G.T.); (N.K.S.); (D.K.)
| | - Anil Rai
- Division of Agricultural Bioinformatics, ICAR-Indian Agricultural Statistics Research Institute, New Delhi 110012, India; (G.T.); (N.K.S.); (D.K.)
| | - B. N. Tripathi
- ICAR-National Research Centre on Equines, Sirsa Road, Hisar 125001, India; (Y.P.)
- Indian Council of Agricultural Research, Krishi Bhawan, New Delhi 110001, India
| |
Collapse
|
16
|
Romanov MN, Abdelmanova AS, Fisinin VI, Gladyr EA, Volkova NA, Anshakov DV, Stanishevskaya OI, Vakhrameev AB, Dotsev AV, Griffin DK, Zinovieva NA. Whole Genome Screening Procures a Holistic Hold of the Russian Chicken Gene Pool Heritage and Demographic History. BIOLOGY 2023; 12:979. [PMID: 37508409 PMCID: PMC10376169 DOI: 10.3390/biology12070979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Revised: 07/01/2023] [Accepted: 07/08/2023] [Indexed: 07/30/2023]
Abstract
A study for genomic variation that may reflect putative selective signaling and be associated with economically important traits is instrumental for obtaining information about demographic and selection history in domestic animal species and populations. A rich variety of the Russian chicken gene pool breeds warrants a further detailed study. Specifically, their genomic features can derive implications from their genome architecture and selective footprints for their subsequent breeding and practical efficient exploitation. In the present work, whole genome genotyping of 19 chicken breeds (20 populations with up to 71 samples each) was performed using the Chicken 50 K BeadChip DNA chip. The studied breed sample included six native Russian breeds of chickens developed in the 17th-19th centuries, as well as eight Russian chicken breeds, including the Russian White (RW), created in the 20th century on the basis of improving local chickens using breeds of foreign selection. Five specialized foreign breeds of chickens, including the White Leghorn (WL), were used along with other breeds representing the Russian gene pool. The characteristics of the genetic diversity and phylogenetic relationships of the native breeds of chickens were represented in comparison with foreign breeds. It was established that the studied native breeds demonstrate their own genetic structure that distinguishes them from foreign breeds, and from each other. For example, we previously made an assumption on what could cause the differences between two RW populations, RW1 and RW2. From the data obtained here, it was verified that WL was additionally crossed to RW2, unlike RW1. Thus, inherently, RW1 is a purer population of this improved Russian breed. A significant contribution of the gene pool of native breeds to the global genetic diversity of chickens was shown. In general, based on the results of a multilateral survey of this sample of breeds, it can be concluded that phylogenetic relationships based on their genetic structure and variability robustly reflect the known, previously postulated and newly discovered patterns of evolution of native chickens. The results herein presented will aid selection and breeding work using this gene pool.
Collapse
Affiliation(s)
- Michael N Romanov
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia
- School of Biosciences, University of Kent, Canterbury CT2 7NJ, Kent, UK
| | - Alexandra S Abdelmanova
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia
| | - Vladimir I Fisinin
- Center "All-Russian Poultry Research and Technological Institute" of the Russian Academy of Sciences, Sergiev Posad 141311, Moscow Oblast, Russia
| | - Elena A Gladyr
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia
| | - Natalia A Volkova
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia
| | - Dmitry V Anshakov
- Breeding and Genetic Center "Zagorsk Experimental Breeding Farm"-Branch of the Federal Research Centre "All-Russian Poultry Research and Technological Institute" of the Russian Academy of Sciences, Sergiev Posad 141311, Moscow Oblast, Russia
| | - Olga I Stanishevskaya
- Russian Research Institute of Farm Animal Genetics and Breeding-Branch of the L. K. Ernst Federal Research Center for Animal Husbandry, Pushkin, Saint Petersburg 196601, Russia
| | - Anatoly B Vakhrameev
- Russian Research Institute of Farm Animal Genetics and Breeding-Branch of the L. K. Ernst Federal Research Center for Animal Husbandry, Pushkin, Saint Petersburg 196601, Russia
| | - Arsen V Dotsev
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia
| | - Darren K Griffin
- School of Biosciences, University of Kent, Canterbury CT2 7NJ, Kent, UK
| | - Natalia A Zinovieva
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia
| |
Collapse
|
17
|
Nishio M, Inoue K, Ogawa S, Ichinoseki K, Arakawa A, Fukuzawa Y, Okamura T, Kobayashi E, Taniguchi M, Oe M, Ishii K. Comparing pedigree and genomic inbreeding coefficients, and inbreeding depression of reproductive traits in Japanese Black cattle. BMC Genomics 2023; 24:376. [PMID: 37403068 DOI: 10.1186/s12864-023-09480-5] [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: 03/28/2023] [Accepted: 06/23/2023] [Indexed: 07/06/2023] Open
Abstract
BACKGROUND Pedigree-based inbreeding coefficients have been generally included in statistical models for genetic evaluation of Japanese Black cattle. The use of genomic data is expected to provide precise assessment of inbreeding level and depression. Recently, many measures have been used for genome-based inbreeding coefficients; however, with no consensus on which is the most appropriate. Therefore, we compared the pedigree- ([Formula: see text]) and multiple genome-based inbreeding coefficients, which were calculated from the genomic relationship matrix with observed allele frequencies ([Formula: see text]), correlation between uniting gametes ([Formula: see text]), the observed vs expected number of homozygous genotypes ([Formula: see text]), runs of homozygosity (ROH) segments ([Formula: see text]) and heterozygosity by descent segments ([Formula: see text]). We quantified inbreeding depression from estimating regression coefficients of inbreeding coefficients on three reproductive traits: age at first calving (AFC), calving difficulty (CD) and gestation length (GL) in Japanese Black cattle. RESULTS The highest correlations with [Formula: see text] were for [Formula: see text] (0.86) and [Formula: see text] (0.85) whereas [Formula: see text] and [Formula: see text] provided weak correlations with [Formula: see text], with range 0.33-0.55. Except for [Formula: see text] and [Formula: see text], there were strong correlations among genome-based inbreeding coefficients ([Formula: see text] 0.94). The estimates of regression coefficients of inbreeding depression for [Formula: see text] was 2.1 for AFC, 0.63 for CD and -1.21 for GL, respectively, but [Formula: see text] had no significant effects on all traits. Genome-based inbreeding coefficients provided larger effects on all reproductive traits than [Formula: see text]. In particular, for CD, all estimated regression coefficients for genome-based inbreeding coefficients were significant, and for GL, that for [Formula: see text] had a significant.. Although there were no significant effects when using overall genome-level inbreeding coefficients for AFC and GL, [Formula: see text] provided significant effects at chromosomal level in four chromosomes for AFC, three chromosomes for CD, and two chromosomes for GL. In addition, similar results were obtained for [Formula: see text]. CONCLUSIONS Genome-based inbreeding coefficients can capture more phenotypic variation than [Formula: see text]. In particular, [Formula: see text] and [Formula: see text] can be considered good estimators for quantifying inbreeding level and identifying inbreeding depression at the chromosome level. These findings might improve the quantification of inbreeding and breeding programs using genome-based inbreeding coefficients.
Collapse
Affiliation(s)
- Motohide Nishio
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan.
| | - Keiichi Inoue
- University of Miyazaki, Miyazaki, Miyazaki, 889-2192, Japan
- National Livestock Breeding Center, Nishigo, Fukushima, 961-8511, Japan
| | - Shinichiro Ogawa
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| | - Kasumi Ichinoseki
- National Livestock Breeding Center, Nishigo, Fukushima, 961-8511, Japan
| | - Aisaku Arakawa
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| | - Yo Fukuzawa
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| | - Toshihiro Okamura
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| | - Eiji Kobayashi
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| | - Masaaki Taniguchi
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| | - Mika Oe
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| | - Kazuo Ishii
- Institute of Livestock and Grassland Science, NARO, Tsukuba, Ibaraki, 3050901, Japan
| |
Collapse
|
18
|
Perfilyeva A, Bespalova K, Bespalov S, Begmanova M, Kuzovleva Y, Vishnyakova O, Nazarenko I, Abylkassymova G, Perfilyeva Y, Plakhov K, Djansugurova L, Bekmanov B. Homozygosity mapping in the Kazakh national dog breed Tazy. Sci Rep 2023; 13:10735. [PMID: 37400480 DOI: 10.1038/s41598-023-37990-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2023] [Accepted: 06/30/2023] [Indexed: 07/05/2023] Open
Abstract
The Tazy is a breed of sighthound common in Kazakhstan. The identification of runs of homozygosity (ROH) is an informative approach to assessing the history and possible patterns of directional selection pressure. To our knowledge, the present study is the first to provide an overview of the ROH pattern in the Tazy dogs from a genome-wide perspective. The ROH of the Tazy was found to be mainly composed of shorter segments (1-2 Mb), accounting for approximately 67% of the total ROH. The estimated ROH-based inbreeding coefficients (FROH) ranged from 0.028 to 0.058 with a mean of 0.057. Five genomic regions under positive selection were identified on chromosomes 18, 22, and 25. The regions on chromosomes 18 and 22 may be breed specific, while the region on chromosome 22 overlaps with regions of hunting traits in other hunting dog breeds. Among the 12 candidate genes located in these regions, the gene CAB39L may be a candidate that affects running speed and endurance of the Tazy dog. Eight genes could belong to an evolutionarily conserved complex as they were clustered in a large protein network with strong linkages. The results may enable effective interventions when incorporated into conservation planning and selection of the Tazy breed.
Collapse
Affiliation(s)
- Anastassiya Perfilyeva
- Department of Molecular Genetics, Institute of Genetics and Physiology, 050060, Almaty, Kazakhstan
| | - Kira Bespalova
- Department of Molecular Genetics, Institute of Genetics and Physiology, 050060, Almaty, Kazakhstan.
- Department of Biology and Biotechnology, Al-Farabi Kazakh National University, 050040, Almaty, Kazakhstan.
| | - Sergey Bespalov
- Department of Theriology, Institute of Zoology, 050060, Almaty, Kazakhstan
| | - Mamura Begmanova
- Department of Molecular Genetics, Institute of Genetics and Physiology, 050060, Almaty, Kazakhstan
| | - Yelena Kuzovleva
- Department of Molecular Genetics, Institute of Genetics and Physiology, 050060, Almaty, Kazakhstan
| | - Olga Vishnyakova
- Republican Federation of Public Associations of Hunters and Hunting Societies "Kansonar", 050008, Almaty, Kazakhstan
| | - Inna Nazarenko
- Republican Federation of Public Associations of Hunters and Hunting Societies "Kansonar", 020000, Astana, Kazakhstan
| | - Gulnar Abylkassymova
- Department of Molecular Genetics, Institute of Genetics and Physiology, 050060, Almaty, Kazakhstan
| | - Yuliya Perfilyeva
- Laboratory of Molecular Immunology and Immunobiotechnology, M.A. Aitkhozhin's Institute of Molecular Biology and Biochemistry, 050012, Almaty, Kazakhstan
| | - Konstantin Plakhov
- Laboratory of Biocenology and Hunting Management, Institute of Zoology, 050060, Almaty, Kazakhstan
| | - Leyla Djansugurova
- Department of Molecular Genetics, Institute of Genetics and Physiology, 050060, Almaty, Kazakhstan
| | - Bakhytzhan Bekmanov
- Department of Molecular Genetics, Institute of Genetics and Physiology, 050060, Almaty, Kazakhstan
- Department of Biology and Biotechnology, Al-Farabi Kazakh National University, 050040, Almaty, Kazakhstan
| |
Collapse
|
19
|
Pérez‐Pereira N, Quesada H, Caballero A. An empirical evaluation of the estimation of inbreeding depression from molecular markers under suboptimal conditions. Evol Appl 2023; 16:1302-1315. [PMID: 37492144 PMCID: PMC10363801 DOI: 10.1111/eva.13568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Revised: 05/30/2023] [Accepted: 05/30/2023] [Indexed: 07/27/2023] Open
Abstract
Inbreeding depression (ID), the reduction in fitness due to inbreeding, is typically measured by the regression of the phenotypic values of individuals for a particular trait on their corresponding inbreeding coefficients (F). While genealogical records can provide these coefficients, they may be unavailable or incomplete, making molecular markers a useful alternative. The power to detect ID and its accuracy depend on the variation of F values of individuals, the sample sizes available, and the accuracy in the estimation of individual fitness traits and F values. In this study, we used Drosophila melanogaster to evaluate the effectiveness of molecular markers in estimating ID under suboptimal conditions. We generated two sets of 100 pairs of unrelated individuals from a large panmictic population and mated them for two generations to produce non-inbred and unrelated individuals (F = 0) and inbred individuals (full-sib progeny; F = 0.25). Using these expected genealogical F values, we calculated inbreeding depression for two fitness-related traits, pupae productivity and competitive fitness. We then sequenced the males from 17 non-inbred pairs and 17 inbred pairs to obtain their genomic inbreeding coefficients and estimate ID for the two traits. The scenario assumed was rather restrictive in terms of estimation of ID because: (1) the individuals belonged to the same generation of a large panmictic population, leading to low variation in individual F coefficients; (2) the sample sizes were small; and (3) the traits measured depended on both males and females while only males were sequenced. Despite the challenging conditions of our study, we found that molecular markers provided estimates of ID that were comparable to those obtained from simple pedigree estimations with larger sample sizes. The results therefore suggest that genomic measures of inbreeding are useful to provide estimates of inbreeding depression even under very challenging scenarios.
Collapse
Affiliation(s)
- Noelia Pérez‐Pereira
- Centro de Investigación MariñaUniversidade de Vigo, Facultade de BioloxíaVigoSpain
| | - Humberto Quesada
- Centro de Investigación MariñaUniversidade de Vigo, Facultade de BioloxíaVigoSpain
| | - Armando Caballero
- Centro de Investigación MariñaUniversidade de Vigo, Facultade de BioloxíaVigoSpain
| |
Collapse
|
20
|
Wilmot H, Druet T, Hulsegge I, Gengler N, Calus M. Estimation of inbreeding, between-breed genomic relatedness and definition of sub-populations in red-pied cattle breeds. Animal 2023; 17:100793. [PMID: 37087997 DOI: 10.1016/j.animal.2023.100793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2022] [Revised: 03/16/2023] [Accepted: 03/21/2023] [Indexed: 03/31/2023] Open
Abstract
Currently, enhancing the collaboration between related breeds is of main importance to increase the competitivity and the sustainability of local breeds. One type of collaboration is the development of an across-breed reference population that will allow a better management of local breeds. For this purpose, the genomic relatedness between the local target breed and possible breeds to be included in the reference population should be estimated. In Europe, there are several local red-pied cattle breeds that would benefit from this kind of collaboration. However, how different red-pied cattle breeds from the Benelux are related to each other and can collaborate is still unclear. The objectives of this study were therefore: (1) to estimate the level of inbreeding of the East Belgian Red and White (EBRW), the Red-Pied of the Ösling (RPO) and Dutch red-pied cattle breeds; (2) to determine the genomic relatedness of several red-pied cattle breeds, with a special focus on two endangered breeds: the EBRW and the RPO, and (3) based on the second objective, to detect animals from other breeds that were genomically close enough to be considered as advantageous in the creation of an across-breed reference population of EBRW or RPO. The estimated inbreeding levels based on runs of homozygosity were relatively low for almost all the studied breeds and especially for the EBRW and RPO. This would imply that inbreeding is currently not an issue in these two endangered breeds and that their sustainability is not threatened by their level of inbreeding. The results from the principal component analysis, the phylogenetic tree and the clustering all highlighted that the EBRW and RPO breeds were included in the genomic continuum of the studied red-pied cattle breeds and can be therefore considered as genomically close to Dutch red-pied cattle breeds, highlighting the possibility of a collaboration between these breeds. Especially, EBRW animals were closely related to Deep Red and Improved Red animals while, to a lesser extent, the RPO animals were closely related to the Meuse-Rhine-Yssel breed. Based on these results, we could use distance measures, based either on the principal component analysis or clustering, to detect animals from Dutch breeds that were genomically closest to the EBRW or RPO breeds. This will finally allow the building of an across-breed reference population for EBRW or RPO for further genomic evaluations, considering these genomically closest animals from other breeds.
Collapse
|
21
|
Spatial genetic structure of European wild boar, with inferences on late-Pleistocene and Holocene demographic history. Heredity (Edinb) 2023; 130:135-144. [PMID: 36639700 PMCID: PMC9981775 DOI: 10.1038/s41437-022-00587-1] [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: 02/23/2021] [Revised: 12/13/2022] [Accepted: 12/14/2022] [Indexed: 01/15/2023] Open
Abstract
European wildlife has been subjected to intensifying levels of anthropogenic impact throughout the Holocene, yet the main genetic partitioning of many species is thought to still reflect the late-Pleistocene glacial refugia. We analyzed 26,342 nuclear SNPs of 464 wild boar (Sus scrofa) across the European continent to infer demographic history and reassess the genetic consequences of natural and anthropogenic forces. We found that population fragmentation, inbreeding and recent hybridization with domestic pigs have caused the spatial genetic structure to be heterogeneous at the local scale. Underlying local anthropogenic signatures, we found a deep genetic structure in the form of an arch-shaped cline extending from the Dinaric Alps, via Southeastern Europe and the Baltic states, to Western Europe and, finally, to the genetically diverged Iberian peninsula. These findings indicate that, despite considerable anthropogenic influence, the deeper, natural continental structure is still intact. Regarding the glacial refugia, our findings show a weaker signal than generally assumed, but are nevertheless suggestive of two main recolonization routes, with important roles for Southern France and the Balkans. Our results highlight the importance of applying genomic resources and framing genetic results within a species' demographic history and geographic distribution for a better understanding of the complex mixture of underlying processes.
Collapse
|
22
|
Romanov MN, Abdelmanova AS, Fisinin VI, Gladyr EA, Volkova NA, Koshkina OA, Rodionov AN, Vetokh AN, Gusev IV, Anshakov DV, Stanishevskaya OI, Dotsev AV, Griffin DK, Zinovieva NA. Selective footprints and genes relevant to cold adaptation and other phenotypic traits are unscrambled in the genomes of divergently selected chicken breeds. J Anim Sci Biotechnol 2023; 14:35. [PMID: 36829208 PMCID: PMC9951459 DOI: 10.1186/s40104-022-00813-0] [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: 06/28/2022] [Accepted: 11/27/2022] [Indexed: 02/26/2023] Open
Abstract
BACKGROUND The genomes of worldwide poultry breeds divergently selected for performance and other phenotypic traits may also be affected by, and formed due to, past and current admixture events. Adaptation to diverse environments, including acclimation to harsh climatic conditions, has also left selection footprints in breed genomes. RESULTS Using the Chicken 50K_CobbCons SNP chip, we genotyped four divergently selected breeds: two aboriginal, cold tolerant Ushanka and Orloff Mille Fleur, one egg-type Russian White subjected to artificial selection for cold tolerance, and one meat-type White Cornish. Signals of selective sweeps were determined in the studied breeds using three methods: (1) assessment of runs of homozygosity islands, (2) FST based population differential analysis, and (3) haplotype differentiation analysis. Genomic regions of true selection signatures were identified by two or more methods or in two or more breeds. In these regions, we detected 540 prioritized candidate genes supplemented them with those that occurred in one breed using one statistic and were suggested in other studies. Amongst them, SOX5, ME3, ZNF536, WWP1, RIPK2, OSGIN2, DECR1, TPO, PPARGC1A, BDNF, MSTN, and beta-keratin genes can be especially mentioned as candidates for cold adaptation. Epigenetic factors may be involved in regulating some of these important genes (e.g., TPO and BDNF). CONCLUSION Based on a genome-wide scan, our findings can help dissect the genetic architecture underlying various phenotypic traits in chicken breeds. These include genes representing the sine qua non for adaptation to harsh environments. Cold tolerance in acclimated chicken breeds may be developed following one of few specific gene expression mechanisms or more than one overlapping response known in cold-exposed individuals, and this warrants further investigation.
Collapse
Affiliation(s)
- Michael N. Romanov
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia ,grid.9759.20000 0001 2232 2818School of Biosciences, University of Kent, Canterbury, UK
| | - Alexandra S. Abdelmanova
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Vladimir I. Fisinin
- grid.4886.20000 0001 2192 9124Federal State Budget Scientific Institution Federal Research Centre “All-Russian Poultry Research and Technological Institute” of the Russian Academy of Sciences, Sergiev Posad, Moscow Region Russia
| | - Elena A. Gladyr
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Natalia A. Volkova
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Olga A. Koshkina
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Andrey N. Rodionov
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Anastasia N. Vetokh
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Igor V. Gusev
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Dmitry V. Anshakov
- grid.4886.20000 0001 2192 9124Breeding and Genetic Centre “Zagorsk Experimental Breeding Farm” – Branch of the Federal Research Centre “All-Russian Poultry Research and Technological Institute” of the Russian Academy of Sciences, Sergiev Posad, Moscow Region Russia
| | - Olga I. Stanishevskaya
- grid.473314.6Russian Research Institute of Farm Animal Genetics and Breeding – Branch of the L.K. Ernst Federal Research Centre for Animal Husbandry, St. Petersburg, Russia
| | - Arsen V. Dotsev
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| | - Darren K. Griffin
- grid.9759.20000 0001 2232 2818School of Biosciences, University of Kent, Canterbury, UK
| | - Natalia A. Zinovieva
- L.K. Ernst Federal Research Centre for Animal Husbandry, Dubrovitsy, Podolsk, Moscow Region Russia
| |
Collapse
|
23
|
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
|
24
|
The coefficients of inbreeding revealed by ROH study among inbred individuals belonging to each type of the first cousin marriage: A preliminary report from North India. Genes Genomics 2023; 45:813-825. [PMID: 36807878 DOI: 10.1007/s13258-023-01367-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 01/27/2023] [Indexed: 02/23/2023]
Abstract
BACKGROUND Genome-wide runs of homozygosity (ROH) are appropriate to estimate genomic inbreeding, determine population history, unravel the genetic architecture of complex traits and disorders. OBJECTIVE The study sought to investigate and compare the actual proportion of homozygosity or autozygosity in the genomes of progeny of four subtypes of first cousin mating in humans, using both pedigree and genomic measures for autosomes and sex chromosomes. METHODS For this purpose, Illumina Global Screening Array-24 v1.0 BeadChip followed by cyto-ROH analysis through Illumina Genome Studio was used to characterise the homozygosity in five participants from North Indian state (Uttar Pradesh). PLINK v.1.9 software was used to estimate the genomic inbreeding coefficients viz. ROH-based inbreeding estimate (FROH) and homozygous loci-based inbreeding estimate (FHOM). RESULTS A total of 133 ROH segments were detected with maximum number and genomic coverage in Matrilateral Parallel (MP) type and minimum in outbred individual. ROH pattern revealed that MP type has a higher degree of homozygosity than other subtypes. The comparison of FROH, FHOM, and pedigree-based inbreeding estimate (FPED) showed some difference in theoretical and realised proportion of homozygosity for sex-chromosomal loci but not for autosome for each type of consanguinity. CONCLUSIONS This is the very first study to compare and estimate the pattern of homozygosity among the kindreds of first cousin unions. However, a greater number of individuals from each type of marriage is required for statistical inference of no difference between theoretical and realized homozygosity among different degrees of inbreeding prevalent in humans worldwide.
Collapse
|
25
|
Ablondi M, Summer A, Stocco G, Finocchiaro R, van Kaam JT, Cassandro M, Dadousis C, Sabbioni A, Cipolat-Gotet C. The role of inbreeding depression on productive performance in the Italian Holstein breed. J Anim Sci 2023; 101:skad382. [PMID: 37983004 PMCID: PMC10693289 DOI: 10.1093/jas/skad382] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Accepted: 11/18/2023] [Indexed: 11/21/2023] Open
Abstract
Inbreeding depression has become an urgent issue in cosmopolitan breeds where the massive genetic progress achieved in the latest generations is counterbalanced by a dramatic loss of genetic diversity causing increased health issues. Thus, the aim of this study was to estimate inbreeding depression on productive traits in Holstein dairy cattle. More precisely, we aimed to i) determine the level of inbreeding in 27,735 Italian Holstein dairy cows using pedigree and genotype data, ii) quantify the effect of inbreeding on 305-d in milk yield (MY; kg), fat yield (FY; kg), and protein yield (PY; kg) based on different statistical approaches, iii) determine if recent inbreeding has a more harmful impact than ancestral ones, and iv) quantify chromosomal homozygosity effect on productive traits. Quality control was performed on the autosomal chromosomes resulting in a final dataset of 84,443 single nucleotide polymorphisms. Four statistical models were used to evaluate the presence of inbreeding depression, which included linear regression analysis and division of FPED and FROH into percentile classes. Moreover, FROH was partitioned into i) length classes to assess the role of recent and ancestral inbreeding and ii) chromosome-specific contributions (FROH-CHR). Results evidenced that inbreeding negatively impacted the productive performance of Italian Holstein Friesian cows. However, differences between the estimated FPED and FROH coefficients resulted in different estimates of inbreeding depression. For instance, a 1% increase in FPED and FROH was associated with a decrease in MY of about 44 and 61 kg (P < 0.01). Further, when considering the extreme inbreeding percentile classes moving from the 5th lowest to the 95th highest, there was a reduction of -263 kg and -561 kg per lactation for FPED and FROH. Increased inbreeding, estimated by FPED and FROH, had also a negative effect on PY and FY, either fit as a regressor or percentile classes. When evaluating the impact of inbreeding based on runs of homozygosity (ROH) length classes, longer ROH (over 8 Mb) had a negative effect in all traits, indicating that recent inbreeding might be more harmful than the ancestral one. Finally, results within chromosome homozygosity highlighted specific chromosomes with a more deleterious effect on productive traits.
Collapse
Affiliation(s)
- Michela Ablondi
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Andrea Summer
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Giorgia Stocco
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Raffaella Finocchiaro
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana (ANAFIBJ), 26100 Cremona, Italy
| | - Jan-Thijs van Kaam
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana (ANAFIBJ), 26100 Cremona, Italy
| | - Martino Cassandro
- Associazione Nazionale Allevatori della Razza Frisona Bruna e Jersey Italiana (ANAFIBJ), 26100 Cremona, Italy
- Department of Agronomy, Food, Natural resources, Animals and Environment, University of Padova, 35020 Legnaro, Italy
| | - Christos Dadousis
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | - Alberto Sabbioni
- Department of Veterinary Science, University of Parma, 43126 Parma, Italy
| | | |
Collapse
|
26
|
Chessari G, Criscione A, Tolone M, Bordonaro S, Rizzuto I, Riggio S, Macaluso V, Moscarelli A, Portolano B, Sardina MT, Mastrangelo S. High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context. Front Vet Sci 2023; 10:1127354. [PMID: 37205231 PMCID: PMC10185747 DOI: 10.3389/fvets.2023.1127354] [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/19/2022] [Accepted: 04/13/2023] [Indexed: 05/21/2023] Open
Abstract
Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise FST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and FST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.
Collapse
Affiliation(s)
- Giorgio Chessari
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Andrea Criscione
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Marco Tolone
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Salvatore Bordonaro
- Dipartimento Agricoltura, Alimentazione e Ambiente, University of Catania, Catania, Italy
| | - Ilaria Rizzuto
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Silvia Riggio
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Vito Macaluso
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Angelo Moscarelli
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Baldassare Portolano
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Maria Teresa Sardina
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Salvatore Mastrangelo
- Dipartimento Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
- *Correspondence: Salvatore Mastrangelo,
| |
Collapse
|
27
|
Ballan M, Schiavo G, Bovo S, Schiavitto M, Negrini R, Frabetti A, Fornasini D, Fontanesi L. Comparative analysis of genomic inbreeding parameters and runs of homozygosity islands in several fancy and meat rabbit breeds. Anim Genet 2022; 53:849-862. [PMID: 36073189 PMCID: PMC9826494 DOI: 10.1111/age.13264] [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: 01/24/2022] [Revised: 07/01/2022] [Accepted: 08/25/2022] [Indexed: 01/11/2023]
Abstract
Runs of homozygosity (ROH) are defined as long stretches of DNA homozygous at each polymorphic position. The proportion of genome covered by ROH and their length are indicators of the level and origin of inbreeding. In this study, we analysed SNP chip datasets (obtained using the Axiom OrcunSNP Array) of a total of 702 rabbits from 12 fancy breeds and four meat breeds to identify ROH with different approaches and calculate several genomic inbreeding parameters. The highest average number of ROH per animal was detected in Belgian Hare (~150) and the lowest in Italian Silver (~106). The average length of ROH ranged from 4.001 ± 0.556 Mb in Italian White to 6.268 ± 1.355 Mb in Ermine. The same two breeds had the lowest (427.9 ± 86.4 Mb, Italian White) and the highest (921.3 ± 179.8 Mb, Ermine) average values of the sum of all ROH segments. More fancy breeds had a higher level of genomic inbreeding (as defined by ROH) than meat breeds. Several ROH islands contain genes involved in body size, body length, pigmentation processes, carcass traits, growth, and reproduction traits (e.g.: AOX1, GPX5, IFRD1, ITGB8, NELL1, NR3C1, OCA2, TRIB1, TRIB2). Genomic inbreeding parameters can be useful to overcome the lack of information in the management of rabbit genetic resources. ROH provided information to understand, to some extent, the genetic history of rabbit breeds and to identify signatures of selection in the rabbit genome.
Collapse
Affiliation(s)
- Mohamad Ballan
- Division of Animal Sciences, Department of Agricultural and Food SciencesUniversity of BolognaBolognaItaly
| | - Giuseppina Schiavo
- Division of Animal Sciences, Department of Agricultural and Food SciencesUniversity of BolognaBolognaItaly
| | - Samuele Bovo
- Division of Animal Sciences, Department of Agricultural and Food SciencesUniversity of BolognaBolognaItaly
| | - Michele Schiavitto
- Associazione Nazionale Coniglicoltori Italiani (ANCI), Contrada Giancola SncVolturara AppulaItaly
| | | | | | | | - Luca Fontanesi
- Division of Animal Sciences, Department of Agricultural and Food SciencesUniversity of BolognaBolognaItaly
| |
Collapse
|
28
|
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
|
29
|
Assessing Genetic Diversity and Searching for Selection Signatures by Comparison between the Indigenous Livni and Duroc Breeds in Local Livestock of the Central Region of Russia. DIVERSITY 2022. [DOI: 10.3390/d14100859] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Indigenous pig breeds are mainly associated with the adaptive capacity that is necessary to respond adequately to climate change, food security, and livelihood needs, and natural resources conservation. Livni pigs are an indigenous fat-type breed farmed in a single farm in the Orel region and located in the Central European part of the Russian Federation. To determine the genomic regions and genes that are affected by artificial selection, we conducted the comparative study of two pig breeds with different breeding histories and breeding objectives, i.e., the native fat-type Livni and meat-type Duroc breeds using the Porcine GGP HD BeadChip, which contains ~80,000 SNPs. To check the Livni pigs for possible admixture, the Landrace and the Large White breeds were included into the study of genetic diversity as these breeds participated in the formation of the Livni pigs. We observed the highest level of genetic diversity in Livni pigs compared to commercial breeds (UHE = 0.409 vs. 0.319–0.359, p < 0.001; AR = 1.995 vs. 1.894–1.964, p < 0.001). A slight excess of heterozygotes was found in all of the breeds. We identified 291 candidate genes, which were localized within the regions under putative selection, including 22 and 228 genes, which were specific for Livni and Duroc breeds, respectively, and 41 genes common for both breeds. A detailed analysis of the molecular functions identified the genes, which were related to the formation of meat and fat traits, and adaptation to environmental stress, including extreme temperatures, which were different between breeds. Our research results are useful for conservation and sustainable breeding of Livni breed, which shows a high level of genetic diversity. This makes Livni one of the valuable national pig genetic resources.
Collapse
|
30
|
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: 6] [Impact Index Per Article: 3.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
|
31
|
Berghöfer J, Khaveh N, Mundlos S, Metzger J. Simultaneous testing of rule- and model-based approaches for runs of homozygosity detection opens up a window into genomic footprints of selection in pigs. BMC Genomics 2022; 23:564. [PMID: 35933356 PMCID: PMC9357325 DOI: 10.1186/s12864-022-08801-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2022] [Accepted: 07/29/2022] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes. RESULTS Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation. CONCLUSIONS The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes.
Collapse
Affiliation(s)
- Jan Berghöfer
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Department of Biology, Chemistry and Pharmacy, Institute of Chemistry and Biochemistry, Freie Universität Berlin, Berlin, Germany
| | - Nadia Khaveh
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany
| | - Stefan Mundlos
- Research Group Development & Disease, Max Planck Institute for Molecular Genetics, Berlin, Germany.,Institute of Medical Genetics and Human Genetics, Charité-Universitätsmedizin Berlin, Berlin, Germany.,Charité-Universitätsmedizin Berlin, BCRT, Berlin Institute of Health Center for Regenerative Therapies, Berlin, Germany
| | - Julia Metzger
- Research Group Veterinary Functional Genomics, Max Planck Institute for Molecular Genetics, Berlin, Germany. .,Institute of Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany.
| |
Collapse
|
32
|
Drzaic I, Curik I, Lukic B, Shihabi M, Li MH, Kantanen J, Mastrangelo S, Ciani E, Lenstra JA, Cubric-Curik V. High-Density Genomic Characterization of Native Croatian Sheep Breeds. Front Genet 2022; 13:940736. [PMID: 35910220 PMCID: PMC9337876 DOI: 10.3389/fgene.2022.940736] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2022] [Accepted: 06/16/2022] [Indexed: 11/13/2022] Open
Abstract
A recent comprehensive genomic analysis based on 50K SNP profiles has shown that the regional Balkan sheep populations have considerable genetic overlap but are distinctly different from surrounding breeds. All eight Croatian sheep breeds were represented by a small number of individuals per breed. Here, we genotyped 220 individuals representing the native Croatian sheep breeds (Istrian Sheep, Krk Island Sheep, Cres Island Sheep, Rab Island Sheep, Lika Pramenka, Pag Island Sheep, Dalmatian Pramenka, Dubrovnik Sheep) and mouflon using the Ovine Infinium® HD SNP BeadChip (606,006 SNPs). In addition, we included publicly available Balkan Pramenka and other Mediterranean sheep breeds. Our analyses revealed the complex population structure of Croatian sheep breeds and their origin and geographic barriers (island versus mainland). Migration patterns confirmed the historical establishment of breeds and the pathways of gene flow. Inbreeding coefficients (FROH>2 Mb) between sheep populations ranged from 0.025 to 0.070, with lower inbreeding coefficients observed in Dalmatian Pramenka and Pag Island Sheep and higher inbreeding in Dubrovnik sheep. The estimated effective population size ranged from 61 to 1039 for Krk Island Sheep and Dalmatian Pramenka, respectively. Higher inbreeding levels and lower effective population size indicate the need for improved conservation management to maintain genetic diversity in some breeds. Our results will contribute to breeding and conservation strategies of native Croatian sheep breeds.
Collapse
Affiliation(s)
- Ivana Drzaic
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
- *Correspondence: Ivana Drzaic, ; Vlatka Cubric-Curik,
| | - Ino Curik
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
| | - Boris Lukic
- Department of Animal Production and Biotechnology, Faculty of Agrobiotechnical Sciences Osijek, Chair for Domestic Animal Breeding and Genetics, J. J. Strossmayer University of Osijek, Osijek, Croatia
| | - Mario Shihabi
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
| | - Meng-Hua Li
- College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Juha Kantanen
- Production Systems, Natural Resources Institute Finland (Luke), Jokioinen, Finland
| | - Salvatore Mastrangelo
- Dipartimento di Scienze Agrarie, Alimentari e Forestali, University of Palermo, Palermo, Italy
| | - Elena Ciani
- Dipartimento di Bioscienze, Biotecnologie e Biofarmaceutica, Universita Degli Studi di Bari “Aldo Moro”, Bari, Italy
| | | | - Vlatka Cubric-Curik
- Department of Animal Science, University of Zagreb Faculty of Agriculture, Zagreb, Croatia
- *Correspondence: Ivana Drzaic, ; Vlatka Cubric-Curik,
| |
Collapse
|
33
|
Zhang Y, Zhuo Y, Ning C, Zhou L, Liu JF. Estimate of inbreeding depression on growth and reproductive traits in a Large White pig population. G3 (BETHESDA, MD.) 2022; 12:jkac118. [PMID: 35551391 PMCID: PMC9258530 DOI: 10.1093/g3journal/jkac118] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Accepted: 05/02/2022] [Indexed: 06/15/2023]
Abstract
With the broad application of genomic information, SNP-based measures of estimating inbreeding have been widely used in animal breeding, especially based on runs of homozygosity. Inbreeding depression is better estimated by SNP-based inbreeding coefficients than pedigree-based inbreeding in general. However, there are few comprehensive comparisons of multiple methods in pigs so far, to some extent limiting their application. In this study, to explore an appropriate strategy for estimating inbreeding depression on both growth traits and reproductive traits in a Large White pig population, we compared multiple methods for the inbreeding coefficient estimation based on both pedigree and genomic information. This pig population for analyzing the influence of inbreeding was from a pig breeding farm in the Inner Mongolia of China. There were 26,204 pigs with records of age at 100 kg (AGE) and back-fat thickness at 100 kg (BF), and 6,656 sows with reproductive records of the total number of piglets at birth (TNB), and the number of alive piglets at birth (NBA), and litter weight at birth. Inbreeding depression affected growth and reproductive traits. The results indicated that pedigree-based and SNP-based inbreeding coefficients had significant effects on AGE, TNB, and NBA, except for BF. However, only SNP-based inbreeding coefficients revealed a strong association with inbreeding depression on litter weight at birth. Runs of homozygosity-based methods showed a slight advantage over other methods in the correlation analysis of inbreeding coefficients and estimation of inbreeding depression. Furthermore, our results demonstrated that the model-based approach (RZooRoH) could avoid miscalculations of inbreeding and inbreeding depression caused by inappropriate parameters, which had a good performance on both AGE and reproductive traits. These findings might improve the extensive application of runs of homozygosity analysis in pig breeding and breed conservation.
Collapse
Affiliation(s)
- Yu Zhang
- National Engineering Laboratory for Animal Breeding; Key Laboratory of Animal Genetics, Breeding, and Reproduction, Ministry of Agriculture; College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Yue Zhuo
- National Engineering Laboratory for Animal Breeding; Key Laboratory of Animal Genetics, Breeding, and Reproduction, Ministry of Agriculture; College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Chao Ning
- Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, College of Animal Science and Technology, Shandong Agricultural University, Tai’an 271000, China
| | - Lei Zhou
- National Engineering Laboratory for Animal Breeding; Key Laboratory of Animal Genetics, Breeding, and Reproduction, Ministry of Agriculture; College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
| | - Jian-Feng Liu
- Corresponding author: College of Animal Science and Technology, China Agricultural University, 2 Yuanmingyuan West Road, Beijing 100193, China. ; Corresponding author: College of Animal Science and Technology, China Agricultural University, 2 Yuanmingyuan West Road, Beijing 100193, China.
| |
Collapse
|
34
|
Abdelmanova AS, Sermyagin AA, Dotsev AV, Rodionov AN, Stolpovsky YA, Zinovieva NA. Whole-Genomic Studies of the Population Structure of Russian Local Black-Pied Breeds. RUSS J GENET+ 2022. [DOI: 10.1134/s102279542207002x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
35
|
Chen Z, Zhang Z, Wang Z, Zhang Z, Wang Q, Pan Y. Heterozygosity and homozygosity regions affect reproductive success and the loss of reproduction: a case study with litter traits in pigs. Comput Struct Biotechnol J 2022; 20:4060-4071. [PMID: 35983229 PMCID: PMC9364102 DOI: 10.1016/j.csbj.2022.07.039] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2022] [Revised: 07/21/2022] [Accepted: 07/22/2022] [Indexed: 12/23/2022] Open
Abstract
Runs of heterozygosity (ROHet) and homozygosity (ROH) harbor useful information related to traits of interest. There is a lack of investigating the effect of ROHet and ROH on reproductive success and the loss of reproduction in mammals. Here, we detected and characterized the ROHet and ROH patterns in the genomes of Chinese indigenous pigs (i.e., Jinhua, Chun’an, Longyou Black, and Shengxian Spotted pigs), revealing the similar genetic characteristics of indigenous pigs. Later, we highlighted the underlying litter traits-related ROHet and ROH using association analysis with linear model in these four indigenous pig breeds. To pinpoint the promising candidate genes associated with litter traits, we further in-depth explore the selection patterns of other five pig breeds (i.e., Erhualian, Meishan, Minzhu, Rongchang, and Diqing pigs) with different levels of reproduction performance at the underlying litter traits-related ROHet and ROH using FST and genetic diversity ratio. Then, we identified a set of known and novel candidate genes associated with reproductive performance in pigs. For the novel candidate genes (i.e., CCDC91, SASH1, SAMD5, MACF1, MFSD2A, EPC2, and MBD5), we obtained public available datasets and performed multi-omics analyses integrating transcriptome-wide association studies and comparative single-cell RNA-seq analyses to uncover the roles of them in mammalian reproductive performance. The genes have not been widely reported to be fertility-related genes and can be complementally considered as prior biological information to modify genomic selections models that benefits pig genetic improvement of litter traits. Besides, our findings provide new insights into the function of ROHet and ROH in mammals.
Collapse
|
36
|
Mészárosová M, Mészáros G, Moravčíková N, Pavlík I, Margetín M, Kasarda R. Within- and between-Breed Selection Signatures in the Original and Improved Valachian Sheep. Animals (Basel) 2022; 12:ani12111346. [PMID: 35681809 PMCID: PMC9179888 DOI: 10.3390/ani12111346] [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/28/2022] [Revised: 05/13/2022] [Accepted: 05/23/2022] [Indexed: 11/16/2022] Open
Abstract
This study explored the genomic diversity and selection signatures in two Slovakian national breeds, the Original Valachian and the Improved Valachian sheep. As they are an important animal genetic resource within the country, but with decreasing population size, our aim is to identify potentially valuable genomic regions. A total of 97 sheep (18 male and 79 female) from the Original Valachian, and 69 sheep (25 male and 44 female) from the Improved Valachian populations were genotyped using the GeneSeek GGP Ovine 50 K chip. The inbreeding levels were assessed with runs of homozygosity (ROH). The selection signatures within breeds were identified based on the top 1% of most homozygous regions within the breed, the so-called ROH islands. The selection signatures between breeds were assessed based on variance in linkage disequilibrium. Overall, we have identified selection signatures with quantitative trait loci (QTL) and genes pointing towards all three production purposes of the Valachian sheep, milk, meat, and wool, including their quality characteristics. Another group with apparent large importance was the various traits related to health and resistance to parasites, which is well in line with the sturdy nature of this breed.
Collapse
Affiliation(s)
- Mária Mészárosová
- Faculty of Agrobiology and Food Resources, Institute of Nutrition and Genomics, Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovakia; (M.M.); (R.K.)
| | - Gábor Mészáros
- Department of Sustainable Agricultural Systems, Division of Livestock Sciences, University of Natural Resources and Life Sciences, Vienna, Gregor-Mendel-Straße 33, 1180 Vienna, Austria;
| | - Nina Moravčíková
- Faculty of Agrobiology and Food Resources, Institute of Nutrition and Genomics, Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovakia; (M.M.); (R.K.)
- Correspondence:
| | - Ivan Pavlík
- Research Institute of Animal Production—NPPC Slovakia, Hlohovecká 2, 95141 Nitra—Lužianky, Slovakia;
| | - Milan Margetín
- Faculty of Agrobiology and Food Resources, Institute of Animal Science, Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovakia;
| | - Radovan Kasarda
- Faculty of Agrobiology and Food Resources, Institute of Nutrition and Genomics, Slovak University of Agriculture in Nitra, Tr. A. Hlinku 2, 94976 Nitra, Slovakia; (M.M.); (R.K.)
| |
Collapse
|
37
|
Malatji DP. Breeding of African sheep reared under low-input/output smallholder production systems for trypanotolerance. Vet World 2022; 15:1031-1043. [PMID: 35698514 PMCID: PMC9178589 DOI: 10.14202/vetworld.2022.1031-1043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Accepted: 03/11/2022] [Indexed: 11/30/2022] Open
Abstract
Trypanosomiasis is a disease caused by unicellular protozoan parasites. Small ruminants succumb to trypanosomiasis in areas of high tsetse fly challenge, resulting in serious economic loss often to farmers in low-input smallholder systems. At present, trypanosomiasis is treated with trypanocidal drugs, but access to these can be limited, and increasing parasite resistance raises questions about their efficacy. The development of trypanotolerance in small ruminant flocks through targeted breeding strategies is considered a sustainable and economical option for controlling African trypanosomiasis. Recently, quantitative trait loci (QTLs) associated with trypanotolerance traits in sheep have been reported. The results of these studies form the basis for more studies to identify QTLs associated with trypanosomiasis resistance, particularly in African livestock species. For example, signatures of positive selection for trypanotolerance have been identified using genome-wide single-nucleotide polymorphism data. However, there are several challenges in performing genetic analyses using data from low-input smallholder systems, including a lack of recorded pedigree and production records and the need for large sample sizes when flock sizes are often fewer than 50 animals. Breeding strategies to improve trypanotolerance should also preserve existing genetic diversity as well as minimize excessive genetic introgression by trypanosusceptible breeds. This review discusses the possibilities of breeding for trypanosome tolerance/resistance in low-input/low-output small ruminant production systems. Potential challenges are outlined, and potential available genetic resources are described as a foundation for future work.
Collapse
Affiliation(s)
- Dikeledi P. Malatji
- Department of Agriculture and Animal Health, College of Agriculture and Environmental Sciences, University of South Africa, Johannesburg, Gauteng Province, South Africa
| |
Collapse
|
38
|
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
|
39
|
Liu SH, Ma XY, Hassan FU, Gao TY, Deng TX. Genome-wide analysis of runs of homozygosity in Italian Mediterranean buffalo. J Dairy Sci 2022; 105:4324-4334. [PMID: 35307184 DOI: 10.3168/jds.2021-21543] [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: 11/07/2021] [Accepted: 02/07/2022] [Indexed: 11/19/2022]
Abstract
Runs of homozygosity (ROH) are a powerful tool to explore patterns of genomic inbreeding in animal populations and detect signatures of selection. The present study used ROH analysis to evaluate the genome-wide patterns of homozygosity, inbreeding levels, and distribution of ROH islands using the SNP data sets from 899 Mediterranean buffaloes. A total of 42,433 ROH segments were identified, with an average of 47.20 segments per individual. The ROH comprising mostly shorter segments (1-4 Mb) accounted for approximately 72.29% of all ROH. In contrast, the larger ROH (>8 Mb) class accounted for only 7.97% of all ROH segments. Estimated inbreeding coefficients from ROH (FROH) ranged from 0.0201 to 0.0371. Pearson correlations between FROH and genomic relationship matrix increased with the increase of ROH length. We identified ROH hotspots in 12 genomic regions, located on chromosomes 1, 2, 3, 5, 17, and 19, harboring a total of 122 genes. Protein-protein interaction (PPI) analysis revealed the clustering of these genes into 7 PPI networks. Many genes located in these regions were associated with different production traits. In addition, 5 ROH islands overlapped with cattle quantitative trait loci that were mainly associated with milk traits. These findings revealed the genome-wide autozygosity patterns and inbreeding levels in Mediterranean buffalo. Our study identified many candidate genes related to production traits that could be used to assist in selective breeding for genetic improvement of buffalo.
Collapse
Affiliation(s)
- Shen-He Liu
- College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450046, China.
| | - Xiao-Ya Ma
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China
| | - Faiz-Ul Hassan
- Institute of Animal and Dairy Sciences, University of Agriculture, Faisalabad 38040, Pakistan
| | - Teng-Yun Gao
- College of Animal Science and Technology, Henan Agricultural University, Zhengzhou 450046, China
| | - Ting-Xian Deng
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning 530001, China.
| |
Collapse
|
40
|
Signer-Hasler H, Henkel J, Bangerter E, Bulut Z, Drögemüller C, Leeb T, Flury C. Runs of homozygosity in Swiss goats reveal genetic changes associated with domestication and modern selection. Genet Sel Evol 2022; 54:6. [PMID: 35073837 PMCID: PMC8785455 DOI: 10.1186/s12711-022-00695-w] [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: 05/26/2021] [Accepted: 01/06/2022] [Indexed: 11/16/2022] Open
Abstract
Background The domestication of goat (Capra hircus) started 11,000 years ago in the fertile crescent. Breed formation in the nineteenth century, establishment of herd books, and selection for specific traits resulted in 10 modern goat breeds in Switzerland. We analyzed whole-genome sequencing (WGS) data from 217 modern goats and nine wild Bezoar goats (Capra aegagrus). After quality control, 27,728,288 biallelic single nucleotide variants (SNVs) were used for the identification of runs of homozygosity (ROH) and the detection of ROH islands. Results Across the 226 caprine genomes from 11 populations, we detected 344 ROH islands that harbor 1220 annotated genes. We compared the ROH islands between the modern breeds and the Bezoar goats. As a proof of principle, we confirmed a signature of selection, which contains the ASIP gene that controls several breed-specific coat color patterns. In two other ROH islands, we identified two missense variants, STC1:p.Lys139Arg and TSHR:p.Ala239Thr, which might represent causative functional variants for domestication signatures. Conclusions We have shown that the information from ROH islands using WGS data is suitable for the analysis of signatures of selection and allowed the detection of protein coding variants that may have conferred beneficial phenotypes during goat domestication. We hypothesize that the TSHR:p.Ala239Thr variant may have played a role in changing the seasonality of reproduction in modern domesticated goats. The exact functional significance of the STC1:p.Lys139Arg variant remains unclear and requires further investigation. Nonetheless, STC1 might represent a new domestication gene affecting relevant traits such as body size and/or milk yield in goats. Supplementary Information The online version contains supplementary material available at 10.1186/s12711-022-00695-w.
Collapse
Affiliation(s)
- Heidi Signer-Hasler
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, 3052, Zollikofen, Switzerland.
| | - Jan Henkel
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| | - Erika Bangerter
- Swiss Goat Breeding Association SZZV, Schützenstrasse 10, 3052, Zollikofen, Switzerland
| | - Zafer Bulut
- Department of Biochemistry, Faculty of Veterinary Medicine, Selcuk University, Konya, Turkey
| | | | - Cord Drögemüller
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| | - Tosso Leeb
- Institute of Genetics, Vetsuisse Faculty, University of Bern, 3001, Bern, Switzerland
| | - Christine Flury
- School of Agricultural, Forest and Food Sciences, Bern University of Applied Sciences, 3052, Zollikofen, Switzerland
| |
Collapse
|
41
|
De los Ríos-Pérez L, Druet T, Goldammer T, Wittenburg D. Analysis of Autozygosity Using Whole-Genome Sequence Data of Full-Sib Families in Pikeperch (Sander lucioperca). Front Genet 2022; 12:786934. [PMID: 35111201 PMCID: PMC8801746 DOI: 10.3389/fgene.2021.786934] [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: 09/30/2021] [Accepted: 12/23/2021] [Indexed: 11/13/2022] Open
Abstract
Pikeperch (Sander lucioperca) has emerged as a high value species to the aquaculture industry. However, its farming techniques are at an early stage and its production is often performed without a selective breeding program, potentially leading to high levels of inbreeding. In this study, we identified and characterized autozygosity based on genome-wide runs of homozygosity (ROH) on a sample of parental and offspring individuals, determined effective population size (Ne), and assessed relatedness among parental individuals. A mean of 2,235 ± 526 and 1,841 ± 363 ROH segments per individual, resulting in a mean inbreeding coefficient of 0.33 ± 0.06 and 0.25 ± 0.06 were estimated for the progeny and parents, respectively. Ne was about 12 until four generations ago and at most 106 for 63 generations in the past, with varying genetic relatedness amongst the parents. This study shows the importance of genomic information when family relationships are unknown and the need of selective breeding programs for reproductive management decisions in the aquaculture industry.
Collapse
Affiliation(s)
- Lidia De los Ríos-Pérez
- Institute of Genetics and Biometry, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
| | - Tom Druet
- Unit of Animal Genomics, GIGA-R & Faculty of Veterinary Medicine, University of Liège, Liège, Belgium
| | - Tom Goldammer
- Institute of Genome Biology, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
- Molecular Biology and Fish Genetics, Faculty of Agriculture and Environmental Sciences, University of Rostock, Rostock, Germany
| | - Dörte Wittenburg
- Institute of Genetics and Biometry, Research Institute for Farm Animal Biology (FBN), Dummerstorf, Germany
- *Correspondence: Dörte Wittenburg,
| |
Collapse
|
42
|
Cesarani A, Gaspa G, Correddu F, Dimauro C, Macciotta NPP. Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity. J Anim Breed Genet 2022; 139:292-306. [PMID: 34984736 DOI: 10.1111/jbg.12666] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2021] [Revised: 12/15/2021] [Accepted: 12/22/2021] [Indexed: 02/04/2023]
Abstract
Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.
Collapse
Affiliation(s)
- Alberto Cesarani
- Department of Agricultural Sciences, University of Sassari, Sassari, Italy.,Department of Animal and Dairy Science, University of Georgia, Athens, Georgia, USA
| | - Giustino Gaspa
- Department of Agricultural, Forestry and Alimentary Sciences, University of Torino, Grugliasco, Italy
| | - Fabio Correddu
- Department of Agricultural Sciences, University of Sassari, Sassari, Italy
| | - Corrado Dimauro
- Department of Agricultural Sciences, University of Sassari, Sassari, Italy
| | | |
Collapse
|
43
|
Schiavo G, Bovo S, Ribani A, Moscatelli G, Bonacini M, Prandi M, Mancin E, Mantovani R, Dall'Olio S, Fontanesi L. Comparative analysis of inbreeding parameters and runs of homozygosity islands in 2 Italian autochthonous cattle breeds mainly raised in the Parmigiano-Reggiano cheese production region. J Dairy Sci 2021; 105:2408-2425. [PMID: 34955250 DOI: 10.3168/jds.2021-20915] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2021] [Accepted: 10/25/2021] [Indexed: 01/19/2023]
Abstract
Reggiana and Modenese are autochthonous cattle breeds, reared in the North of Italy, that can be mainly distinguished for their standard coat color (Reggiana is red, whereas Modenese is white with some pale gray shades). Almost all milk produced by these breeds is transformed into 2 mono-breed branded Parmigiano-Reggiano cheeses, from which farmers receive the economic incomes needed for the sustainable conservation of these animal genetic resources. After the setting up of their herd books in 1960s, these breeds experienced a strong reduction in the population size that was subsequently reverted starting in the 1990s (Reggiana) or more recently (Modenese) reaching at present a total of about 2,800 and 500 registered cows, respectively. Due to the small population size of these breeds, inbreeding is a very important cause of concern for their conservation programs. Inbreeding is traditionally estimated using pedigree data, which are summarized in an inbreeding coefficient calculated at the individual level (FPED). However, incompleteness of pedigree information and registration errors can affect the effectiveness of conservation strategies. High-throughput SNP genotyping platforms allow investigation of inbreeding using genome information that can overcome the limits of pedigree data. Several approaches have been proposed to estimate genomic inbreeding, with the use of runs of homozygosity (ROH) considered to be the more appropriate. In this study, several pedigree and genomic inbreeding parameters, calculated using the whole herd book populations or considering genotyping information (GeneSeek GGP Bovine 150K) from 1,684 Reggiana cattle and 323 Modenese cattle, were compared. Average inbreeding values per year were used to calculate effective population size. Reggiana breed had generally lower genomic inbreeding values than Modenese breed. The low correlation between pedigree-based and genomic-based parameters (ranging from 0.187 to 0.195 and 0.319 to 0.323 in the Reggiana and Modenese breeds, respectively) reflected the common problems of local populations in which pedigree records are not complete. The high proportion of short ROH over the total number of ROH indicates no major recent inbreeding events in both breeds. ROH islands spread over the genome of the 2 breeds (15 in Reggiana and 14 in Modenese) identified several signatures of selection. Some of these included genes affecting milk production traits, stature, body conformation traits (with a main ROH island in both breeds on BTA6 containing the ABCG2, NCAPG, and LCORL genes) and coat color (on BTA13 in Modenese containing the ASIP gene). In conclusion, this work provides an extensive comparative analysis of pedigree and genomic inbreeding parameters and relevant genomic information that will be useful in the conservation strategies of these 2 iconic local cattle breeds.
Collapse
Affiliation(s)
- Giuseppina Schiavo
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Samuele Bovo
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Anisa Ribani
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Giulia Moscatelli
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Massimo Bonacini
- Associazione Nazionale Allevatori Bovini di Razza Reggiana (ANABORARE), Via Masaccio 11, 42124 Reggio Emilia, Italy
| | - Marco Prandi
- Associazione Nazionale Allevatori Bovini di Razza Reggiana (ANABORARE), Via Masaccio 11, 42124 Reggio Emilia, Italy
| | - Enrico Mancin
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale dell'Università 16, 35020 Legnaro (PD), Italy
| | - Roberto Mantovani
- Department of Agronomy, Food, Natural Resources, Animals and Environment (DAFNAE), University of Padova, Viale dell'Università 16, 35020 Legnaro (PD), Italy
| | - Stefania Dall'Olio
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy
| | - Luca Fontanesi
- Department of Agricultural and Food Sciences, Division of Animal Sciences, University of Bologna, Viale Giuseppe Fanin 46, 40127 Bologna, Italy.
| |
Collapse
|
44
|
Krupa E, Moravčíková N, Krupová Z, Žáková E. Assessment of the Genetic Diversity of a Local Pig Breed Using Pedigree and SNP Data. Genes (Basel) 2021; 12:1972. [PMID: 34946921 PMCID: PMC8702119 DOI: 10.3390/genes12121972] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2021] [Revised: 12/06/2021] [Accepted: 12/08/2021] [Indexed: 12/04/2022] Open
Abstract
Herein, the genetic diversity of the local Přeštice Black-Pied pig breed was assessed by the simultaneous analysis of the pedigree and single nucleotide polymorphism (SNP) data. The information about sire line, dam, date of birth, sex, breeding line, and herd for 1971 individuals was considered in the pedigree analysis. The SNP analysis (n = 181) was performed using the Illumina PorcineSNP60 BeadChip kit. The quality of pedigree and SNPs and the inbreeding coefficients (F) and effective population size (Ne) were evaluated. The correlations between inbreeding based on the runs of homozygosity (FROH) and pedigree (FPED) were also calculated. The average FPED for all animals was 3.44%, while the FROH varied from 10.81% for a minimum size of 1 Mbp to 3.98% for a minimum size of 16 Mbp. The average minor allele frequency was 0.28 ± 0.11. The observed and expected within breed heterozygosities were 0.38 ± 0.13 and 0.37 ± 0.12, respectively. The Ne, obtained using both the data sources, reached values around 50 animals. Moderate correlation coefficients (0.49-0.54) were observed between FPED and FROH. It is necessary to make decisions that stabilize the inbreeding rate in the long-term using optimal contribution selection based on the available SNP data.
Collapse
Affiliation(s)
- Emil Krupa
- Institute of Animal Science, 104 00 Prague, Czech Republic; (Z.K.); (E.Ž.)
| | - Nina Moravčíková
- Faculty of Agrobiology and Food Resources, Institute of Nutrition and Genomics, Slovak University of Agriculture, 949 76 Nitra, Slovakia;
| | - Zuzana Krupová
- Institute of Animal Science, 104 00 Prague, Czech Republic; (Z.K.); (E.Ž.)
| | - Eliška Žáková
- Institute of Animal Science, 104 00 Prague, Czech Republic; (Z.K.); (E.Ž.)
| |
Collapse
|
45
|
Paul K, D'Ambrosio J, Phocas F. Temporal and region‐specific variations in genome‐wide inbreeding effects on female size and reproduction traits of rainbow trout. Evol Appl 2021; 15:645-662. [PMID: 35505890 PMCID: PMC9046919 DOI: 10.1111/eva.13308] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2020] [Revised: 07/11/2021] [Accepted: 09/14/2021] [Indexed: 02/01/2023] Open
Abstract
Recent studies have shown that current levels of inbreeding, estimated by runs of homozygosity (ROH), are moderate to high in farmed rainbow trout lines. Based on ROH metrics, the aims of our study were to (i) quantify inbreeding effects on female size (postspawning body weight, fork length) and reproduction traits (spawning date, coelomic fluid weight, spawn weight, egg number, average egg weight) in rainbow trout, and (ii) identify both the genomic regions and inbreeding events affecting performance. We analysed the performance of 1346 females under linear animal models including random additive and dominance genetics effects, with fixed covariates accounting for inbreeding effects at different temporal and genomic scales. A significant effect of genome‐wide inbreeding (F) was only observed for spawning date and egg weight, with performance variations of +12.3% and −3.8%, respectively, for 0.1 unit increase in F level. At different local genomic scales, we observed highly variable inbreeding effects on the seven traits under study, ranging from increasing to decreasing trait values. As widely reported in the literature, the main scenario observed during this study was a negative impact of recent inbreeding. However, other scenarios such as positive effects of recent inbreeding or negative impacts of old inbreeding were also observed. Although partial dominance appeared to be the main hypothesis explaining inbreeding depression for all the traits studied, the overdominance hypothesis might also play a significant role in inbreeding depression affecting fecundity (egg number and mass) traits in rainbow trout. These findings suggest that region‐specific inbreeding can strongly impact performance without necessarily observing genome‐wide inbreeding effects. They shed light on the genetic architecture of inbreeding depression and its evolution along the genome over time. The use of region‐specific metrics may enable breeders to more accurately manage the trade‐off between genetic merit and the undesirable side effects associated with inbreeding.
Collapse
Affiliation(s)
- Katy Paul
- Université Paris‐Saclay INRAE AgroParisTech GABIJouy‐en‐Josas France
| | - Jonathan D'Ambrosio
- Université Paris‐Saclay INRAE AgroParisTech GABIJouy‐en‐Josas France
- SYSAAFStation INRAE‐LPGPCampus de Beaulieu Rennes France
| | - Florence Phocas
- Université Paris‐Saclay INRAE AgroParisTech GABIJouy‐en‐Josas France
| |
Collapse
|
46
|
Liu D, Chen Z, Zhao W, Guo L, Sun H, Zhu K, Liu G, Shen X, Zhao X, Wang Q, Ma P, Pan Y. Genome-wide selection signatures detection in Shanghai Holstein cattle population identified genes related to adaption, health and reproduction traits. BMC Genomics 2021; 22:747. [PMID: 34654366 PMCID: PMC8520274 DOI: 10.1186/s12864-021-08042-x] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Accepted: 09/28/2021] [Indexed: 12/26/2022] Open
Abstract
Background Over several decades, a wide range of natural and artificial selection events in response to subtropical environments, intensive pasture and intensive feedlot systems have greatly changed the customary behaviour, appearance, and important economic traits of Shanghai Holstein cattle. In particular, the longevity of the Shanghai Holstein cattle population is generally short, approximately the 2nd to 3rd lactation. In this study, two complementary approaches, integrated haplotype score (iHS) and runs of homozygosity (ROH), were applied for the detection of selection signatures within the genome using genotyping by genome-reduced sequence data from 1092 cows. Results In total, 101 significant iHS genomic regions containing selection signatures encompassing a total of 256 candidate genes were detected. There were 27 significant |iHS| genomic regions with a mean |iHS| score > 2. The average number of ROH per individual was 42.15 ± 25.47, with an average size of 2.95 Mb. The length of 78 % of the detected ROH was within the range of 1–2 MB and 2–4 MB, and 99 % were shorter than 8 Mb. A total of 168 genes were detected in 18 ROH islands (top 1 %) across 16 autosomes, in which each SNP showed a percentage of occurrence > 30 %. There were 160 and 167 genes associated with the 52 candidate regions within health-related QTL intervals and 59 candidate regions within reproduction-related QTL intervals, respectively. Annotation of the regions harbouring clustered |iHS| signals and candidate regions for ROH revealed a panel of interesting candidate genes associated with adaptation and economic traits, such as IL22RA1, CALHM3, ITGA9, NDUFB3, RGS3, SOD2, SNRPA1, ST3GAL4, ALAD, EXOSC10, and MASP2. In a further step, a total of 1472 SNPs in 256 genes were matched with 352 cis-eQTLs in 21 tissues and 27 trans-eQTLs in 6 tissues. For SNPs located in candidate regions for ROH, a total of 108 cis-eQTLs in 13 tissues and 4 trans-eQTLs were found for 1092 SNPs. Eighty-one eGenes were significantly expressed in at least one tissue relevant to a trait (P value < 0.05) and matched the 256 genes detected by iHS. For the 168 significant genes detected by ROH, 47 gene-tissue pairs were significantly associated with at least one of the 37 traits. Conclusions We provide a comprehensive overview of selection signatures in Shanghai Holstein cattle genomes by combining iHS and ROH. Our study provides a list of genes associated with immunity, reproduction and adaptation. For functional annotation, the cGTEx resource was used to interpret SNP-trait associations. The results may facilitate the identification of genes relevant to important economic traits and can help us better understand the biological processes and mechanisms affected by strong ongoing natural or artificial selection in livestock populations. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-08042-x.
Collapse
Affiliation(s)
- Dengying Liu
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Zhenliang Chen
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Wei Zhao
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Longyu Guo
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Hao Sun
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China
| | - Kai Zhu
- Shanghai Dairy Cattle Breeding Centre Co., Ltd, 201901, Shanghai, P.R. China
| | - Guanglei Liu
- Shanghai Dairy Cattle Breeding Centre Co., Ltd, 201901, Shanghai, P.R. China
| | - Xiuping Shen
- Shanghai Agricultural Development Promotion Center, 200335, Shanghai, PR China
| | - Xiaoduo Zhao
- Shanghai Dairy Cattle Breeding Centre Co., Ltd, 201901, Shanghai, P.R. China
| | - Qishan Wang
- Department of Animal Breeding and Reproduction, College of Animal Science, Zhejiang University, 310058, Hangzhou, PR China
| | - Peipei Ma
- Department of Animal Science, School of Agriculture and Biology, Shanghai Jiao Tong University, 200240, Shanghai, PR China.
| | - Yuchun Pan
- Department of Animal Breeding and Reproduction, College of Animal Science, Zhejiang University, 310058, Hangzhou, PR China.
| |
Collapse
|
47
|
Santos WB, Schettini GP, Maiorano AM, Bussiman FO, Balieiro JCC, Ferraz GC, Pereira GL, Baldassini WA, Neto ORM, Oliveira HN, Curi RA. Genome-wide scans for signatures of selection in Mangalarga Marchador horses using high-throughput SNP genotyping. BMC Genomics 2021; 22:737. [PMID: 34645387 PMCID: PMC8515666 DOI: 10.1186/s12864-021-08053-8] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2020] [Accepted: 09/07/2021] [Indexed: 12/25/2022] Open
Abstract
BACKGROUND The detection of signatures of selection in genomic regions provides insights into the evolutionary process, enabling discoveries regarding complex phenotypic traits. In this research, we focused on identifying genomic regions affected by different selection pressures, mainly highlighting the recent positive selection, as well as understanding the candidate genes and functional pathways associated with the signatures of selection in the Mangalarga Marchador genome. Besides, we seek to direct the discussion about genes and traits of importance in this breed, especially traits related to the type and quality of gait, temperament, conformation, and locomotor system. RESULTS Three different methods were used to search for signals of selection: Tajima's D (TD), the integrated haplotype score (iHS), and runs of homozygosity (ROH). The samples were composed of males (n = 62) and females (n = 130) that were initially chosen considering well-defined phenotypes for gait: picada (n = 86) and batida (n = 106). All horses were genotyped using a 670 k Axiom® Equine Genotyping Array (Axiom MNEC670). In total, 27, 104 (chosen), and 38 candidate genes were observed within the signatures of selection identified in TD, iHS, and ROH analyses, respectively. The genes are acting in essential biological processes. The enrichment analysis highlighted the following functions: anterior/posterior pattern for the set of genes (GLI3, HOXC9, HOXC6, HOXC5, HOXC4, HOXC13, HOXC11, and HOXC10); limb morphogenesis, skeletal system, proximal/distal pattern formation, JUN kinase activity (CCL19 and MAP3K6); and muscle stretch response (MAPK14). Other candidate genes were associated with energy metabolism, bronchodilator response, NADH regeneration, reproduction, keratinization, and the immunological system. CONCLUSIONS Our findings revealed evidence of signatures of selection in the MM breed that encompass genes acting on athletic performance, limb development, and energy to muscle activity, with the particular involvement of the HOX family genes. The genome of MM is marked by recent positive selection. However, Tajima's D and iHS results point also to the presence of balancing selection in specific regions of the genome.
Collapse
Affiliation(s)
- Wellington B Santos
- Department of Animal Science, São Paulo State University (Unesp) - FCAV, Via de Acesso Professor Paulo Donato Castelane, NN, CEP: 14884-900, Jaboticabal, SP, Brazil.
| | - Gustavo P Schettini
- Department of Animal Science, São Paulo State University (Unesp) - FCAV, Via de Acesso Professor Paulo Donato Castelane, NN, CEP: 14884-900, Jaboticabal, SP, Brazil
| | - Amanda M Maiorano
- Department of Animal Science, São Paulo State University (Unesp) - FCAV, Via de Acesso Professor Paulo Donato Castelane, NN, CEP: 14884-900, Jaboticabal, SP, Brazil
| | - Fernando O Bussiman
- Department of Animal Science, University of São Paulo (USP) - FZEA, Pirassununga, Brazil
| | - Júlio C C Balieiro
- Department of Animal Science, University of São Paulo (USP) - FZEA, Pirassununga, Brazil
| | - Guilherme C Ferraz
- Department of Animal Science, São Paulo State University (Unesp) - FCAV, Via de Acesso Professor Paulo Donato Castelane, NN, CEP: 14884-900, Jaboticabal, SP, Brazil
| | - Guilherme L Pereira
- Department of Breeding and Animal Nutrition, São Paulo State University (Unesp) - FMVZ, Botucatu, Brazil
| | - Welder Angelo Baldassini
- Department of Breeding and Animal Nutrition, São Paulo State University (Unesp) - FMVZ, Botucatu, Brazil
| | - Otávio R M Neto
- Department of Breeding and Animal Nutrition, São Paulo State University (Unesp) - FMVZ, Botucatu, Brazil
| | - Henrique N Oliveira
- Department of Animal Science, São Paulo State University (Unesp) - FCAV, Via de Acesso Professor Paulo Donato Castelane, NN, CEP: 14884-900, Jaboticabal, SP, Brazil
| | - Rogério A Curi
- Department of Breeding and Animal Nutrition, São Paulo State University (Unesp) - FMVZ, Botucatu, Brazil
| |
Collapse
|
48
|
Genomic characterization of the nucleus for conservation of the Chino Santandereano breed using SNP markers. Trop Anim Health Prod 2021; 53:492. [PMID: 34596804 DOI: 10.1007/s11250-021-02936-4] [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/04/2021] [Accepted: 09/23/2021] [Indexed: 10/20/2022]
Abstract
The Chino Santandereano (CHS) creole cattle breed has undergone an adaptation process of more than 500 years to the Colombian mountainous tropics. Despite its advantages, the breed has been replaced by specialized cattle that put it at risk of extinction. The aim of this work was to estimate the genomic diversity and population structure of the nucleus of conservation of the CHS breed. Thirty-seven CHS and 20 Brahman animals were genotyped due to the possible introgression of genes with the GGP Bovine LD v3 chip. Quality control was performed, and linkage disequilibrium (LD), effective population size (Ne), ROH segments, homozygosity, and genomic inbreeding in the breed were estimated. Subsequently, 50 K genomic information of the Holstein (n = 30) and Romosinuano (n = 8) breeds were included to estimate the minor allele frequency (MAF) with common markers and constructing the graphs of the principal component analysis (PCA). Pairwise FSTs were estimated and a neighbor-joining tree was constructed using the IBS matrix. Admixture was used with k = 2 to 10 for the racial composition. LD (r2) was found up to a distance of 0.13 Mb, r2 > 0.3 at a distance of 340.3 kb, and Ne of 32 ± 1. ROH inbreeding was 5.36 ± 0.86%, with a higher contribution from recent inbreeding (4.55%). The PCA showed that the creole breeds were closer together, and the Brahman was more distant. The admixture analysis suggested k = 5 possible ancestral groups and shows that within the CHS breed, there seem to be two different groups with little Holstein and Brahman introgression. The genetic diversity parameters obtained in this work show minimal diversity in this breed and reinforce the need to protect this resource and the conservation banks.
Collapse
|
49
|
Hajihosseinlo A, Nejati-Javaremi A, Miraei-Ashtiani SR. Genetic structure analysis in several populations of cattle using SNP genotypes. Anim Biotechnol 2021; 34:288-300. [PMID: 34591729 DOI: 10.1080/10495398.2021.1960360] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
Parameters such as effective population size (Ne), runs of homozygosity (ROH), and inbreeding based on ROH (FROH) can give new insight into the level of genetic diversity for the population under selection. This research aimed to measure the extent of linkage disequilibrium (LD), effective population size (Ne), Haplotype Block Structure, and runs of homozygosity (ROHs) in several populations of cattle using SNP genotypes. In this study, that the average r2 decreased with the increasing distance of SNP pairs. A general decrease in Ne can be seen for all four populations, indicating a loss of genetic diversity. The Iranian Holstein had the lowest level of genomic inbreeding at an ROH of 1, 5, 10 Mb, while the French Holstein had the highest. The maximum number of ROH is seen at a distance of less than 1 Mb, and the lowest number of ROH is seen at a distance of 10 Mb. The number of ROH decreases with increasing distance due to the increased recombination rate. This is a concern as an increase in inbreeding leads to a reduction in the effective population size, which was also evident in the study populations.
Collapse
Affiliation(s)
- Abbas Hajihosseinlo
- Department of Animal Science, University of Tehran Aras International Campus, Jolfa, Iran
| | | | | |
Collapse
|
50
|
Runs of homozygosity analysis reveals consensus homozygous regions affecting production traits in Chinese Simmental beef cattle. BMC Genomics 2021; 22:678. [PMID: 34548021 PMCID: PMC8454143 DOI: 10.1186/s12864-021-07992-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2021] [Accepted: 09/07/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Genomic regions with a high frequency of runs of homozygosity (ROH) are related to important traits in farm animals. We carried out a comprehensive analysis of ROH and evaluated their association with production traits using the BovineHD (770 K) SNP array in Chinese Simmental beef cattle. RESULTS We detected a total of 116,953 homozygous segments with 2.47Gb across the genome in the studied population. The average number of ROH per individual was 99.03 and the average length was 117.29 Mb. Notably, we detected 42 regions with a frequency of more than 0.2. We obtained 17 candidate genes related to body size, meat quality, and reproductive traits. Furthermore, using Fisher's exact test, we found 101 regions were associated with production traits by comparing high groups with low groups in terms of production traits. Of those, we identified several significant regions for production traits (P < 0.05) by association analysis, within which candidate genes including ECT2, GABRA4, and GABRB1 have been previously reported for those traits in beef cattle. CONCLUSIONS Our study explored ROH patterns and their potential associations with production traits in beef cattle. These results may help to better understand the association between production traits and genome homozygosity and offer valuable insights into managing inbreeding by designing reasonable breeding programs in farm animals.
Collapse
|