1
|
Mahar K, Gurao A, Kumar A, Pratap Singh L, Chitkara M, Gowane GR, Ahlawat S, Niranjan SK, Pundir RK, Kataria RS, Dige MS. Genomic inbreeding analysis reveals resilience and genetic diversity in Indian yak populations. Gene 2024; 928:148787. [PMID: 39053660 DOI: 10.1016/j.gene.2024.148787] [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/15/2024] [Revised: 07/02/2024] [Accepted: 07/19/2024] [Indexed: 07/27/2024]
Abstract
The yak (Bos grunniens), renowned for its adaptability to extreme cold and hypoxic conditions, stands as a remarkable domestic animal crucial for sustaining livelihoods in harsh climates. We conducted a comprehensive analysis of the whole genome sequence data from three distinct Indian yak populations: Arunachali yak (n = 10), Himachali yak (n = 10), and Ladakhi yak (n = 10). The genomic data for Indian yaks were meticulously generated by our laboratory and compared with their Chinese counterpart, the Jinchuan yak (n = 8), for a more nuanced understanding. Our investigation revealed a total of 37,437 runs of homozygosity (ROH) segments in 34 animals representing four distinct yak populations. The Jinchuan yak population exhibited the highest proportion, constituting 80.8 % of total ROHs, predominantly as small segments (<0.1 Mb), accounting for 63 % of the overall ROHs. Further analysis uncovered a significantly higher degree of inbreeding in Chinese yaks compared to their Indian counterparts. The Indian yak populations, in contrast, demonstrated relatively lower and consistent levels of inbreeding. Moreover, we identified ROH hotspots that covered at least 60 % of individuals in our study, indicating their pivotal role in environmental adaptation. A total of five hotspot regions were detected, housing genes such as ENSBGRG00000015023 (WNT2), YIPF4, SPAST, TLN2, and DSG4. These genes are associated with traits including hair follicle initiation, nutrient stress response, microtubule assembly, development of cardiac muscle, hair follicle, and coat color. This observation strongly suggests that there is substantial selection acting on these genes, emphasizing their important role in environmental adaptation among yak populations.
Collapse
Affiliation(s)
- Karan Mahar
- Division of Animal Genetics and Breeding, ICAR-National Dairy Research Institute, Karnal 132001, Haryana, India
| | - Ankita Gurao
- Division of Animal Genetic Resources, ICAR-National Bureau of Animal Genetic Resources, Karnal 132001, Haryana, India
| | - Amod Kumar
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal 132001, Haryana, India
| | - Lalit Pratap Singh
- Division of Animal Biotechnology, ICAR-National Dairy Research Institute, Karnal 132001, Haryana, India
| | - Meenakshi Chitkara
- Division of Animal Biotechnology, ICAR-National Dairy Research Institute, Karnal 132001, Haryana, India
| | - Gopal R Gowane
- Division of Animal Genetics and Breeding, ICAR-National Dairy Research Institute, Karnal 132001, Haryana, India
| | - Sonika Ahlawat
- Division of Animal Biotechnology, ICAR-National Bureau of Animal Genetic Resources, Karnal 132001, Haryana, India
| | - S K Niranjan
- Division of Animal Genetics, ICAR-National Bureau of Animal Genetic Resources, Karnal 132001, Haryana, India
| | - R K Pundir
- Division of Animal Genetic Resources, ICAR-National Bureau of Animal Genetic Resources, Karnal 132001, Haryana, India
| | - R S Kataria
- Division of Animal Biotechnology, ICAR-National Bureau of Animal Genetic Resources, Karnal 132001, Haryana, India
| | - M S Dige
- Division of Animal Genetic Resources, ICAR-National Bureau of Animal Genetic Resources, Karnal 132001, Haryana, India.
| |
Collapse
|
2
|
Moreno-Grau S, Vernekar M, Lopez-Pineda A, Mas-Montserrat D, Barrabés M, Quinto-Cortés CD, Moatamed B, Lee MTM, Yu Z, Numakura K, Matsuda Y, Wall JD, Ioannidis AG, Katsanis N, Takano T, Bustamante CD. Polygenic risk score portability for common diseases across genetically diverse populations. Hum Genomics 2024; 18:93. [PMID: 39218908 PMCID: PMC11367857 DOI: 10.1186/s40246-024-00664-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2024] [Accepted: 08/19/2024] [Indexed: 09/04/2024] Open
Abstract
BACKGROUND Polygenic risk scores (PRS) derived from European individuals have reduced portability across global populations, limiting their clinical implementation at worldwide scale. Here, we investigate the performance of a wide range of PRS models across four ancestry groups (Africans, Europeans, East Asians, and South Asians) for 14 conditions of high-medical interest. METHODS To select the best-performing model per trait, we first compared PRS performances for publicly available scores, and constructed new models using different methods (LDpred2, PRS-CSx and SNPnet). We used 285 K European individuals from the UK Biobank (UKBB) for training and 18 K, including diverse ancestries, for testing. We then evaluated PRS portability for the best models in Europeans and compared their accuracies with respect to the best PRS per ancestry. Finally, we validated the selected PRS models using an independent set of 8,417 individuals from Biobank of the Americas-Genomelink (BbofA-GL); and performed a PRS-Phewas. RESULTS We confirmed a decay in PRS performances relative to Europeans when the evaluation was conducted using the best-PRS model for Europeans (51.3% for South Asians, 46.6% for East Asians and 39.4% for Africans). We observed an improvement in the PRS performances when specifically selecting ancestry specific PRS models (phenotype variance increase: 1.62 for Africans, 1.40 for South Asians and 0.96 for East Asians). Additionally, when we selected the optimal model conditional on ancestry for CAD, HDL-C and LDL-C, hypertension, hypothyroidism and T2D, PRS performance for studied populations was more comparable to what was observed in Europeans. Finally, we were able to independently validate tested models for Europeans, and conducted a PRS-Phewas, identifying cross-trait interplay between cardiometabolic conditions, and between immune-mediated components. CONCLUSION Our work comprehensively evaluated PRS accuracy across a wide range of phenotypes, reducing the uncertainty with respect to which PRS model to choose and in which ancestry group. This evaluation has let us identify specific conditions where implementing risk-prioritization strategies could have practical utility across diverse ancestral groups, contributing to democratizing the implementation of PRS.
Collapse
Affiliation(s)
- Sonia Moreno-Grau
- Galatea Bio, Inc, 14350 Commerce Way, Miami Lakes, FL, 33146, USA
- Department of Biomedical Data Science, Stanford University School of Medicine, 1265 Welch Road, Stanford, CA, 94305, USA
| | - Manvi Vernekar
- Genomelink, Inc, 2150 Shattuck Avenue, Berkeley, CA, 94704, USA
| | - Arturo Lopez-Pineda
- Galatea Bio, Inc, 14350 Commerce Way, Miami Lakes, FL, 33146, USA
- , Amphora Health. Batallon Independencia 80, Morelia, Michoacan, 58260, Mexico
- Escuela Nacional de Estudios Superiores, Unidad Morelia, Universidad Nacional Autonoma de México, Antigua Carretera a Pátzcuaro No. 8701, Col. Ex Hacienda de San José de la Huerta, Morelia, Michoacán, C.P. 58190, Mexico
| | | | - Míriam Barrabés
- Galatea Bio, Inc, 14350 Commerce Way, Miami Lakes, FL, 33146, USA
| | | | - Babak Moatamed
- Galatea Bio, Inc, 14350 Commerce Way, Miami Lakes, FL, 33146, USA
| | | | - Zhenning Yu
- Genomelink, Inc, 2150 Shattuck Avenue, Berkeley, CA, 94704, USA
| | | | - Yuta Matsuda
- Genomelink, Inc, 2150 Shattuck Avenue, Berkeley, CA, 94704, USA
| | - Jeffrey D Wall
- Galatea Bio, Inc, 14350 Commerce Way, Miami Lakes, FL, 33146, USA
| | - Alexander G Ioannidis
- Galatea Bio, Inc, 14350 Commerce Way, Miami Lakes, FL, 33146, USA
- Department of Biomedical Data Science, Stanford University School of Medicine, 1265 Welch Road, Stanford, CA, 94305, USA
- University of California Santa Cruz, 1156 High Street, Santa Cruz, CA, 95064, USA
| | | | - Tomohiro Takano
- Genomelink, Inc, 2150 Shattuck Avenue, Berkeley, CA, 94704, USA.
- Japan: Awakens Japan K.K. (Japanese subsidiary of Genomelink, Inc.), 2-11-3, Meguro, Meguro-ku, 1530063, Tokyo, Japan.
| | - Carlos D Bustamante
- Galatea Bio, Inc, 14350 Commerce Way, Miami Lakes, FL, 33146, USA.
- Department of Biomedical Data Science, Stanford University School of Medicine, 1265 Welch Road, Stanford, CA, 94305, USA.
| |
Collapse
|
3
|
Bem RD, Benfica LF, Silva DA, Carrara ER, Brito LF, Mulim HA, Borges MS, Cyrillo JNSG, Canesin RC, Bonilha SFM, Mercadante MEZ. Assessing different metrics of pedigree and genomic inbreeding and inbreeding effect on growth, fertility, and feed efficiency traits in a closed-herd Nellore cattle population. BMC Genomics 2024; 25:738. [PMID: 39080557 PMCID: PMC11290228 DOI: 10.1186/s12864-024-10641-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2024] [Accepted: 07/19/2024] [Indexed: 08/02/2024] Open
Abstract
BACKGROUND The selection of individuals based on their predicted breeding values and mating of related individuals can increase the proportion of identical-by-descent alleles. In this context, the objectives of this study were to estimate inbreeding coefficients based on alternative metrics and data sources such as pedigree (FPED), hybrid genomic relationship matrix H (FH), and ROH of different length (FROH); and calculate Pearson correlations between the different metrics in a closed Nellore cattle population selected for body weight adjusted to 378 days of age (W378). In addition to total FROH (all classes) coefficients were also estimated based on the size class of the ROH segments: FROH1 (1-2 Mb), FROH2 (2-4 Mb), FROH3 (4-8 Mb), FROH4 (8-16 Mb), and FROH5 (> 16 Mb), and for each chromosome (FROH_CHR). Furthermore, we assessed the effect of each inbreeding metric on birth weight (BW), body weights adjusted to 210 (W210) and W378, scrotal circumference (SC), and residual feed intake (RFI). We also evaluated the chromosome-specific effects of inbreeding on growth traits. RESULTS The correlation between FPED and FROH was 0.60 while between FH and FROH and FH and FPED were 0.69 and 0.61, respectively. The annual rate of inbreeding was 0.16% for FPED, 0.02% for FH, and 0.16% for FROH. A 1% increase in FROH5 resulted in a reduction of up to -1.327 ± 0.495 kg in W210 and W378. Four inbreeding coefficients (FPED, FH, FROH2, and FROH5) had a significant effect on W378, with reductions of up to -3.810 ± 1.753 kg per 1% increase in FROH2. There was an unfavorable effect of FPED on RFI (0.01 ± 0.0002 kg dry matter/day) and of FROH on SC (-0.056 ± 0.022 cm). The FROH_CHR coefficients calculated for BTA3, BTA5, and BTA8 significantly affected the growth traits. CONCLUSIONS Inbreeding depression was observed for all traits evaluated. However, these effects were greater for the criterion used for selection of the animals (i.e., W378). The increase in the genomic inbreeding was associated with a higher inbreeding depression on the traits evaluated when compared to pedigree-based inbreeding. Genomic information should be used as a tool during mating to optimize control of inbreeding and, consequently, minimize inbreeding depression in Nellore cattle.
Collapse
Affiliation(s)
- Ricardo D Bem
- Institute of Animal Science, Sertãozinho, SP, Brazil.
- Department of Animal Science, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil.
| | - Lorena F Benfica
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907, USA
| | - Delvan A Silva
- Department of Animal Sciences, Federal University of Viçosa, Viçosa, MG, Brazil
| | - Eula R Carrara
- Department of Animal Sciences, Federal University of Viçosa, Viçosa, MG, Brazil
| | - Luiz F Brito
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907, USA
| | - Henrique A Mulim
- Department of Animal Sciences, Purdue University, 270 S. Russell Street, West Lafayette, IN, 47907, USA
| | - Marcelo S Borges
- Department of Pathology, Reproduction and One Health, Faculty of Agricultural and Veterinary Sciences, Sao Paulo State University, Jaboticabal, SP, Brazil
| | | | | | | | | |
Collapse
|
4
|
Davoudi P, Do DN, Rathgeber B, Colombo S, Sargolzaei M, Plastow G, Wang Z, Miar Y. Identification of consensus homozygous regions and their associations with growth and feed efficiency traits in American mink. BMC Genom Data 2024; 25:68. [PMID: 38982354 PMCID: PMC11234557 DOI: 10.1186/s12863-024-01252-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Accepted: 07/01/2024] [Indexed: 07/11/2024] Open
Abstract
The recent chromosome-based genome assembly and the newly developed 70K single nucleotide polymorphism (SNP) array for American mink (Neogale vison) facilitate the identification of genetic variants underlying complex traits in this species. The objective of this study was to evaluate the association between consensus runs of homozygosity (ROH) with growth and feed efficiency traits in American mink. A subsample of two mink populations (n = 2,986) were genotyped using the Affymetrix Mink 70K SNP array. The identified ROH segments were included simultaneously, concatenated into consensus regions, and the ROH-based association studies were carried out with linear mixed models considering a genomic relationship matrix for 11 growth and feed efficiency traits implemented in ASReml-R version 4. In total, 298,313 ROH were identified across all individuals, with an average length and coverage of 4.16 Mb and 414.8 Mb, respectively. After merging ROH segments, 196 consensus ROH regions were detected and used for genome-wide ROH-based association analysis. Thirteen consensus ROH regions were significantly (P < 0.01) associated with growth and feed efficiency traits. Several candidate genes within the significant regions are known for their involvement in growth and body size development, including MEF2A, ADAMTS17, POU3F2, and TYRO3. In addition, we found ten consensus ROH regions, defined as ROH islands, with frequencies over 80% of the population. These islands harbored 12 annotated genes, some of which were related to immune system processes such as DTX3L, PARP9, PARP14, CD86, and HCLS1. This is the first study to explore the associations between homozygous regions with growth and feed efficiency traits in American mink. Our findings shed the light on the effects of homozygosity in the mink genome on growth and feed efficiency traits, that can be utilized in developing a sustainable breeding program for mink.
Collapse
Affiliation(s)
- Pourya Davoudi
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Duy Ngoc Do
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Bruce Rathgeber
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Stefanie Colombo
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada
| | - Mehdi Sargolzaei
- Department of Pathobiology, University of Guelph, Guelph, ON, Canada
- Select Sires Inc, Plain City, OH, USA
| | - Graham Plastow
- Department of Agricultural, Food and Nutritional Science, Livestock Gentec, University of Alberta, Edmonton, AB, Canada
| | - Zhiquan Wang
- Department of Agricultural, Food and Nutritional Science, Livestock Gentec, University of Alberta, Edmonton, AB, Canada
| | - Younes Miar
- Department of Animal Science and Aquaculture, Dalhousie University, Truro, NS, Canada.
| |
Collapse
|
5
|
Sarviaho K, Uimari P, Martikainen K. Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population. J Dairy Sci 2024; 107:4822-4832. [PMID: 38490540 DOI: 10.3168/jds.2024-24105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 02/15/2024] [Indexed: 03/17/2024]
Abstract
The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits (e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G) overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.
Collapse
Affiliation(s)
- Katri Sarviaho
- Department of Agricultural Sciences, University of Helsinki, Helsinki 00014, Finland.
| | - Pekka Uimari
- Department of Agricultural Sciences, University of Helsinki, Helsinki 00014, Finland
| | - Katja Martikainen
- Department of Agricultural Sciences, University of Helsinki, Helsinki 00014, Finland
| |
Collapse
|
6
|
Lawson JM, Shilton CA, Lindsay-McGee V, Psifidi A, Wathes DC, Raudsepp T, de Mestre AM. Does inbreeding contribute to pregnancy loss in Thoroughbred horses? Equine Vet J 2024; 56:711-718. [PMID: 38221707 DOI: 10.1111/evj.14057] [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: 03/13/2023] [Accepted: 12/29/2023] [Indexed: 01/16/2024]
Abstract
BACKGROUND Excessive inbreeding increases the probability of uncovering homozygous recessive genotypes and has been associated with an increased risk of retained placenta and lower semen quality. No genomic analysis has investigated the association between inbreeding levels and pregnancy loss. OBJECTIVES To compare genetic inbreeding coefficients (F) of naturally occurring Thoroughbred Early Pregnancy Loss (EPLs), Mid and Late term Pregnancy Loss (MLPL) and Controls. The F value was hypothesised to be higher in cases of pregnancy loss (EPLs and MLPLs) than Controls. STUDY DESIGN Observational case-control study. METHODS Allantochorion and fetal DNA from EPL (n = 37, gestation age 14-65 days), MLPL (n = 94, gestational age 70 days-24 h post parturition) and Controls (n = 58) were genotyped on the Axiom Equine 670K SNP Genotyping Array. Inbreeding coefficients using Runs of Homozygosity (FROH) were calculated using PLINK software. ROHs were split into size categories to investigate the recency of inbreeding. RESULTS MLPLs had significantly higher median number of ROH (188 interquartile range [IQR], 180.8-197.3), length of ROH (3.10, IQR 2.93-3.33), and total number of ROH (590.8, IQR 537.3-632.3), and FROH (0.26, IQR 0.24-0.28) when compared with the Controls and the EPLs (p < 0.05). There was no significant difference in any of the inbreeding indices between the EPLs and Controls. The MLPLs had a significantly higher proportion of long (>10 Mb) ROH (2.5%, IQR 1.6-3.6) than the Controls (1.7%, IQR 0.6-2.5), p = 0.001. No unique ROHs were found in the EPL or MLPL populations. MAIN LIMITATIONS SNP-array data does not allow analysis of every base in the sequence. CONCLUSIONS This first study of the effect of genomic inbreeding levels on pregnancy loss showed that inbreeding is a contributor to MLPL, but not EPL in the UK Thoroughbred population. Mating choices remain critical, because inbreeding may predispose to MLPL by increasing the risk of homozygosity for specific lethal allele(s).
Collapse
Affiliation(s)
- Jessica M Lawson
- Department of Pathobiology and Population Sciences, The Royal Veterinary College, University of London, Hatfield, UK
| | - Charlotte A Shilton
- Department of Comparative Biomedical Sciences, The Royal Veterinary College, University of London, London, UK
| | - Victoria Lindsay-McGee
- Department of Clinical Science and Services, The Royal Veterinary College, University of London, London, UK
| | - Androniki Psifidi
- Department of Clinical Science and Services, The Royal Veterinary College, University of London, London, UK
| | - D Claire Wathes
- Department of Pathobiology and Population Sciences, The Royal Veterinary College, University of London, Hatfield, UK
| | - Terje Raudsepp
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station, Texas, USA
| | - Amanda M de Mestre
- Department of Comparative Biomedical Sciences, The Royal Veterinary College, University of London, London, UK
| |
Collapse
|
7
|
Tan X, Liu L, Dong J, Huang M, Zhang J, Li Q, Wang H, Bai L, Cui M, Zhou Z, Wu D, Xiang Y, Li W, Wang D. Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens. BMC Genomics 2024; 25:485. [PMID: 38755540 PMCID: PMC11097469 DOI: 10.1186/s12864-024-10349-4] [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: 03/01/2024] [Accepted: 04/25/2024] [Indexed: 05/18/2024] Open
Abstract
BACKGROUND Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding. RESULTS To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively. CONCLUSION Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.
Collapse
Affiliation(s)
- Xiaodong Tan
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Lu Liu
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
- Jinhua Jinfan Feed Co., Ltd, Jinhua, Zhejiang, 321000, China
| | - Jie Dong
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Minjie Huang
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Jiawen Zhang
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - Qinghai Li
- Animal Husbandry Institute, Hangzhou Academy of Agricultural Sciences, Hangzhou, 310024, China
| | - Huanhuan Wang
- Animal Husbandry Institute, Hangzhou Academy of Agricultural Sciences, Hangzhou, 310024, China
| | - Lijuan Bai
- Zhejiang Animal Husbandry Technology Extension and Breeding Livestock and Poultry Monitoring Station, Hangzhou, 310020, China
| | - Ming Cui
- Zhejiang Animal Husbandry Technology Extension and Breeding Livestock and Poultry Monitoring Station, Hangzhou, 310020, China
| | - Zhenzhen Zhou
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China
| | - De Wu
- Postdoctoral Research Station, Jinhua Development Zone, Jinhua, Zhejiang, 321000, China
| | - Yun Xiang
- Jinhua Jinfan Feed Co., Ltd, Jinhua, Zhejiang, 321000, China.
| | - Weifen Li
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
| | - Deqian Wang
- Institute of Animal Husbandry and Veterinary Science, Zhejiang Academy of Agricultural Sciences, Hangzhou, 310021, China.
| |
Collapse
|
8
|
Romanov MN, Shakhin AV, Abdelmanova AS, Volkova NA, Efimov DN, Fisinin VI, Korshunova LG, Anshakov DV, Dotsev AV, Griffin DK, Zinovieva NA. Dissecting Selective Signatures and Candidate Genes in Grandparent Lines Subject to High Selection Pressure for Broiler Production and in a Local Russian Chicken Breed of Ushanka. Genes (Basel) 2024; 15:524. [PMID: 38674458 PMCID: PMC11050503 DOI: 10.3390/genes15040524] [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: 03/28/2024] [Revised: 04/16/2024] [Accepted: 04/20/2024] [Indexed: 04/28/2024] Open
Abstract
Breeding improvements and quantitative trait genetics are essential to the advancement of broiler production. The impact of artificial selection on genomic architecture and the genetic markers sought remains a key area of research. Here, we used whole-genome resequencing data to analyze the genomic architecture, diversity, and selective sweeps in Cornish White (CRW) and Plymouth Rock White (PRW) transboundary breeds selected for meat production and, comparatively, in an aboriginal Russian breed of Ushanka (USH). Reads were aligned to the reference genome bGalGal1.mat.broiler.GRCg7b and filtered to remove PCR duplicates and low-quality reads using BWA-MEM2 and bcftools software; 12,563,892 SNPs were produced for subsequent analyses. Compared to CRW and PRW, USH had a lower diversity and a higher genetic distinctiveness. Selective sweep regions and corresponding candidate genes were examined based on ZFST, hapFLK, and ROH assessment procedures. Twenty-seven prioritized chicken genes and the functional projection from human homologs suggest their importance for selection signals in the studied breeds. These genes have a functional relationship with such trait categories as body weight, muscles, fat metabolism and deposition, reproduction, etc., mainly aligned with the QTLs in the sweep regions. This information is pivotal for further executing genomic selection to enhance phenotypic traits.
Collapse
Affiliation(s)
- Michael N. Romanov
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia; (A.V.S.); (A.S.A.); (N.A.V.); (A.V.D.)
- School of Biosciences, University of Kent, Canterbury CT2 7NJ, UK;
| | - Alexey V. Shakhin
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia; (A.V.S.); (A.S.A.); (N.A.V.); (A.V.D.)
| | - Alexandra S. Abdelmanova
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia; (A.V.S.); (A.S.A.); (N.A.V.); (A.V.D.)
| | - Natalia A. Volkova
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia; (A.V.S.); (A.S.A.); (N.A.V.); (A.V.D.)
| | - Dmitry N. Efimov
- Federal State Budget Scientific Institution Federal Scientific Center “All-Russian Research and Technological Poultry Institute”, Sergiev Posad 141311, Moscow Oblast, Russia; (D.N.E.); (V.I.F.); (L.G.K.)
| | - Vladimir I. Fisinin
- Federal State Budget Scientific Institution Federal Scientific Center “All-Russian Research and Technological Poultry Institute”, Sergiev Posad 141311, Moscow Oblast, Russia; (D.N.E.); (V.I.F.); (L.G.K.)
| | - Liudmila G. Korshunova
- Federal State Budget Scientific Institution Federal Scientific Center “All-Russian Research and Technological Poultry Institute”, Sergiev Posad 141311, Moscow Oblast, Russia; (D.N.E.); (V.I.F.); (L.G.K.)
| | - Dmitry V. Anshakov
- Breeding and Genetic Center “Zagorsk Experimental Breeding Farm”—Branch of the Federal Research Center “All-Russian Poultry Research and Technological Institute”, Russian Academy of Sciences, Sergiev Posad 141311, Moscow Oblast, Russia;
| | - Arsen V. Dotsev
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia; (A.V.S.); (A.S.A.); (N.A.V.); (A.V.D.)
| | | | - Natalia A. Zinovieva
- L. K. Ernst Federal Research Center for Animal Husbandry, Dubrovitsy, Podolsk 142132, Moscow Oblast, Russia; (A.V.S.); (A.S.A.); (N.A.V.); (A.V.D.)
| |
Collapse
|
9
|
Liu X, Lin L, Sinding MHS, Bertola LD, Hanghøj K, Quinn L, Garcia-Erill G, Rasmussen MS, Schubert M, Pečnerová P, Balboa RF, Li Z, Heaton MP, Smith TPL, Pinto RR, Wang X, Kuja J, Brüniche-Olsen A, Meisner J, Santander CG, Ogutu JO, Masembe C, da Fonseca RR, Muwanika V, Siegismund HR, Albrechtsen A, Moltke I, Heller R. Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations. Nat Commun 2024; 15:2921. [PMID: 38609362 PMCID: PMC11014984 DOI: 10.1038/s41467-024-47015-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2023] [Accepted: 03/11/2024] [Indexed: 04/14/2024] Open
Abstract
The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.
Collapse
Affiliation(s)
- Xiaodong Liu
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Long Lin
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | | - Laura D Bertola
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Kristian Hanghøj
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Liam Quinn
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | | | | - Mikkel Schubert
- Novo Nordisk Foundation Center for Basic Metabolic Research, University of Copenhagen, Copenhagen, Denmark
| | | | - Renzo F Balboa
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Zilong Li
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Michael P Heaton
- USDA, ARS, U.S. Meat Animal Research Center (USMARC), Clay Center, NE, USA
| | - Timothy P L Smith
- USDA, ARS, U.S. Meat Animal Research Center (USMARC), Clay Center, NE, USA
| | - Rui Resende Pinto
- CIIMAR-Interdisciplinary Centre of Marine and Environmental Research-University of Porto, Porto, Portugal
- Section for Biodiversity, Globe Institute, University of Copenhagen, Copenhagen, Denmark
| | - Xi Wang
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Josiah Kuja
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | | - Jonas Meisner
- Novo Nordisk Foundation Center for Basic Metabolic Research, University of Copenhagen, Copenhagen, Denmark
- Copenhagen Research Centre for Mental Health, Copenhagen University Hospital, Copenhagen, Denmark
| | - Cindy G Santander
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | - Joseph O Ogutu
- Biostatistics Unit, Institute of Crop Science, University of Hohenheim, Stuttgart, Germany
| | - Charles Masembe
- Department of Zoology, Entomology and Fisheries Sciences, Makerere University, P. O. Box 7062, Kampala, Uganda
| | - Rute R da Fonseca
- CIIMAR-Interdisciplinary Centre of Marine and Environmental Research-University of Porto, Porto, Portugal
- Section for Biodiversity, Globe Institute, University of Copenhagen, Copenhagen, Denmark
| | - Vincent Muwanika
- Department of Environmental Management, Makerere University, PO Box 7062, Kampala, Uganda
| | - Hans R Siegismund
- Department of Biology, University of Copenhagen, Copenhagen, Denmark
| | | | - Ida Moltke
- Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| | - Rasmus Heller
- Department of Biology, University of Copenhagen, Copenhagen, Denmark.
| |
Collapse
|
10
|
Rajawat D, Ghildiyal K, Sonejita Nayak S, Sharma A, Parida S, Kumar S, Ghosh AK, Singh U, Sivalingam J, Bhushan B, Dutt T, Panigrahi M. Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle. Gene 2024; 901:148178. [PMID: 38242377 DOI: 10.1016/j.gene.2024.148178] [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: 09/26/2023] [Revised: 01/10/2024] [Accepted: 01/16/2024] [Indexed: 01/21/2024]
Abstract
The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from -0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075-35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.
Collapse
Affiliation(s)
- Divya Rajawat
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Kanika Ghildiyal
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Sonali Sonejita Nayak
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Anurodh Sharma
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Subhashree Parida
- Pharmacology & Toxicology Division, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Shive Kumar
- Department of Animal Genetics and Breeding, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Uttarakhand, India
| | - A K Ghosh
- Department of Animal Genetics and Breeding, Govind Ballabh Pant University of Agriculture and Technology, Pantnagar, Uttarakhand, India
| | - Umesh Singh
- ICAR Central Institute for Research on Cattle, Meerut, UP, India
| | | | - Bharat Bhushan
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Triveni Dutt
- Livestock Production and Management Section, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India
| | - Manjit Panigrahi
- Division of Animal Genetics, Indian Veterinary Research Institute, Izatnagar, Bareilly 243122, UP, India.
| |
Collapse
|
11
|
Zhang W, Lin K, Fu W, Xie J, Fan X, Zhang M, Luo H, Yin Y, Guo Q, Huang H, Chen T, Lin X, Yuan Y, Huang C, Du S. Insights for the Captive Management of South China Tigers Based on a Large-Scale Genetic Survey. Genes (Basel) 2024; 15:398. [PMID: 38674333 PMCID: PMC11049310 DOI: 10.3390/genes15040398] [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/17/2024] [Revised: 03/17/2024] [Accepted: 03/22/2024] [Indexed: 04/28/2024] Open
Abstract
There is an urgent need to find a way to improve the genetic diversity of captive South China tiger (SCT, Panthera tigris amoyensis), the most critically endangered taxon of living tigers, facing inbreeding depression. The genomes showed that 13 hybrid SCTs from Meihuashan were divided into two groups; one group included three individuals who had a closer relationship with pureblood SCTs than another group. The three individuals shared more that 40% of their genome with pureblood SCTs and might be potential individuals for genetic rescuing in SCTs. A large-scale genetic survey based on 319 pureblood SCTs showed that the mean microsatellite inbreeding coefficient of pureblood SCTs decreased significantly from 0.1789 to 0.0600 (p = 0.000009) and the ratio of heterozygous loci increased significantly from 38.5% to 43.2% (p = 0.02) after one individual of the Chongqing line joined the Suzhou line and began to breed in the mid-1980s, which is a reason why the current SCTs keep a moderate level of microsatellite heterozygosity and nucleotide diversity. However, it is important to establish a back-up population based on the three individuals through introducing one pureblood SCT into the back-up population every year. The back-up population should be an important reserve in case the pureblood SCTs are in danger in the future.
Collapse
Affiliation(s)
- Wenping Zhang
- Key Laboratory of Monitoring Biological Diversity in Minshan Mountain of National Park of Giant Pandas, College of Life Science & Biotechnology, Mianyang Normal University, Mianyang 621000, China; (W.Z.)
| | - Kaixiong Lin
- Fujian Meihuashan Institute of South China Tiger Breeding, Longyan 364201, China; (K.L.); (H.L.)
| | - Wenyuan Fu
- Longyan Geopark Protection and Development Center, Longyan 364201, China
| | - Junjin Xie
- Chengdu Research Base of Giant Panda Breeding, Chengdu 610081, China
| | - Xueyang Fan
- Chengdu Research Base of Giant Panda Breeding, Chengdu 610081, China
| | - Mingchun Zhang
- China Conservation and Research Center for the Giant Panda, Chengdu 611830, China;
| | - Hongxing Luo
- Fujian Meihuashan Institute of South China Tiger Breeding, Longyan 364201, China; (K.L.); (H.L.)
| | | | - Qiang Guo
- Key Laboratory of Monitoring Biological Diversity in Minshan Mountain of National Park of Giant Pandas, College of Life Science & Biotechnology, Mianyang Normal University, Mianyang 621000, China; (W.Z.)
| | - He Huang
- Chengdu Research Base of Giant Panda Breeding, Chengdu 610081, China
| | - Tengteng Chen
- Fujian Meihuashan Institute of South China Tiger Breeding, Longyan 364201, China; (K.L.); (H.L.)
| | - Xipan Lin
- Fujian Meihuashan Institute of South China Tiger Breeding, Longyan 364201, China; (K.L.); (H.L.)
| | | | - Cheng Huang
- Fujian Meihuashan Institute of South China Tiger Breeding, Longyan 364201, China; (K.L.); (H.L.)
| | - Shizhang Du
- Key Laboratory of Monitoring Biological Diversity in Minshan Mountain of National Park of Giant Pandas, College of Life Science & Biotechnology, Mianyang Normal University, Mianyang 621000, China; (W.Z.)
| |
Collapse
|
12
|
González-Cano R, González-Martínez A, Ramón M, González Serrano M, Moreno Millán M, Rubio de Juan A, Rodero Serrano E. Exploring the Effects of Robertsonian Translocation 1/29 (Rob (1;29)) on Genetic Diversity in Minor Breeds of Spanish Berrenda Cattle via Genome-Wide Analysis. Animals (Basel) 2024; 14:793. [PMID: 38473178 DOI: 10.3390/ani14050793] [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/24/2024] [Revised: 02/25/2024] [Accepted: 03/02/2024] [Indexed: 03/14/2024] Open
Abstract
Most of the previous studies on the genetic variability in Spanish "Berrenda" breeds have been carried out using DNA microsatellites. The present work aimed to estimate the genetic diversity, population structure, and potential genetic differences among individuals of both Berrenda breeds and groups based on the presence of the Robertsonian chromosomal translocation, rob (1;29). A total of 373 samples from animals belonging to the two breeds, including 169 cases diagnosed as rob (1;29)-positive, were genotyped using an SNP50K chip. The genetic diversity at the breed level did not show significant differences, but it was significantly lower in those subpopulations containing the rob (1;29). Runs of homozygosity identified a region of homozygosity on chromosome 6, where the KIT (KIT proto-oncogene, receptor tyrosine kinase) gene, which determines the typical spotted coat pattern in both breeds, is located. The four subpopulations considered showed minor genetic differences. The regions of the genome that most determined the differences between the breeds were observed on chromosomes 4, 6, 18, and 22. The presence of this Robertsonian translocation did not result in sub-structuring within each of the breeds considered. To improve the reproductive performance of Berrenda breeds, it would be necessary to implement strategies considering the involvement of potential breeding stock carrying rob (1;29).
Collapse
Affiliation(s)
- Rafael González-Cano
- Ministry of Agriculture, Fisheries and Food, Paseo Infanta Isabel 1, 28014 Madrid, Spain
- Regional Center of Animal Breeding and Reproduction (CERSYRA-IRIAF), Avenida del Vino 10, 13300 Ciudad Real, Spain
| | - Ana González-Martínez
- Department of Animal Production, Faculty of Veterinary Sciences, University of Cordoba, 14071 Córdoba, Spain
| | - Manuel Ramón
- Department of Animal Breeding and Genetics, National Institute for Agricultural and Food Research and Technology (INIA-CSIC), 28040 Madrid, Spain
| | - Miriam González Serrano
- Department of Animal Production, Faculty of Veterinary Sciences, University of Cordoba, 14071 Córdoba, Spain
| | - Miguel Moreno Millán
- Department of Genetic, Faculty of Veterinary Sciences, University of Cordoba, 14071 Córdoba, Spain
| | - Alejandro Rubio de Juan
- Regional Center of Animal Breeding and Reproduction (CERSYRA-IRIAF), Avenida del Vino 10, 13300 Ciudad Real, Spain
| | - Evangelina Rodero Serrano
- Department of Animal Production, Faculty of Veterinary Sciences, University of Cordoba, 14071 Córdoba, Spain
| |
Collapse
|
13
|
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
|
14
|
Li W, Wu X, Xiang D, Zhang W, Wu L, Meng X, Huo J, Yin Z, Fu G, Zhao G. Genome-Wide Detection for Runs of Homozygosity in Baoshan Pigs Using Whole Genome Resequencing. Genes (Basel) 2024; 15:233. [PMID: 38397222 PMCID: PMC10887577 DOI: 10.3390/genes15020233] [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: 01/16/2024] [Revised: 02/08/2024] [Accepted: 02/09/2024] [Indexed: 02/25/2024] Open
Abstract
Baoshan pigs (BS) are a local breed in Yunnan Province that may face inbreeding owing to its limited population size. To accurately evaluate the inbreeding level of the BS pig population, we used whole-genome resequencing to identify runs of homozygosity (ROH) regions in BS pigs, calculated the inbreeding coefficient based on pedigree and ROH, and screened candidate genes with important economic traits from ROH islands. A total of 22,633,391 SNPS were obtained from the whole genome of BS pigs, and 201 ROHs were detected from 532,450 SNPS after quality control. The number of medium-length ROH (1-5 Mb) was the highest (98.43%), the number of long ROH (>5 Mb) was the lowest (1.57%), and the inbreeding of BS pigs mainly occurred in distant generations. The inbreeding coefficient FROH, calculated based on ROH, was 0.018 ± 0.016, and the FPED, calculated based on the pedigree, was 0.027 ± 0.028, which were positively correlated. Forty ROH islands were identified, containing 507 genes and 891 QTLs. Several genes were associated with growth and development (IGFALS, PTN, DLX5, DKK1, WNT2), meat quality traits (MC3R, ACSM3, ECI1, CD36, ROCK1, CACNA2D1), and reproductive traits (NPW, TSHR, BMP7). This study provides a reference for the protection and utilization of BS pigs.
Collapse
Affiliation(s)
- Wenjun Li
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (W.L.); (L.W.); (X.M.); (J.H.); (G.F.)
| | - Xudong Wu
- Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei 230036, China; (X.W.); (W.Z.)
| | - Decai Xiang
- Institute of Pig and Animal Research, Yunnan Academy of Animal Husbandry and Veterinary Science, Kunming 650201, China;
| | - Wei Zhang
- Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei 230036, China; (X.W.); (W.Z.)
| | - Lingxiang Wu
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (W.L.); (L.W.); (X.M.); (J.H.); (G.F.)
| | - Xintong Meng
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (W.L.); (L.W.); (X.M.); (J.H.); (G.F.)
| | - Jinlong Huo
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (W.L.); (L.W.); (X.M.); (J.H.); (G.F.)
| | - Zongjun Yin
- College of Animal Science and Technology, Anhui Agricultural University, Hefei 230036, China;
| | - Guowen Fu
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (W.L.); (L.W.); (X.M.); (J.H.); (G.F.)
| | - Guiying Zhao
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming 650201, China; (W.L.); (L.W.); (X.M.); (J.H.); (G.F.)
| |
Collapse
|
15
|
Illa SK, Mumtaz S, Nath S, Mukherjee S, Mukherjee A. Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle. J Appl Genet 2024; 65:167-180. [PMID: 38110827 DOI: 10.1007/s13353-023-00816-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2023] [Revised: 11/29/2023] [Accepted: 12/05/2023] [Indexed: 12/20/2023]
Abstract
Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.
Collapse
Affiliation(s)
- Satish Kumar Illa
- Livestock Research Station, Garividi, Sri Venkateswara Veterinary University, Tirupati, Andhra Pradesh State, India
| | - Shabahat Mumtaz
- Animal Husbandry Department, Kolkata, West Bengal State, India
| | - Sapna Nath
- College of Veterinary Science, Garividi, Sri Venkateswara Veterinary University, Tirupati, Andhra Pradesh State, India
| | - Sabyasachi Mukherjee
- Animal Genetics & Breeding Division, Indian Council of Agricultural Research (ICAR)-National Dairy Research Institute (NDRI), Karnal, Haryana State, India.
| | - Anupama Mukherjee
- Animal Genetics & Breeding Division, Indian Council of Agricultural Research (ICAR)-National Dairy Research Institute (NDRI), Karnal, Haryana State, India.
| |
Collapse
|
16
|
Gislason H. SNP heterozygosity, relatedness and inbreeding of whole genomes from the isolated population of the Faroe Islands. BMC Genomics 2023; 24:707. [PMID: 37996805 PMCID: PMC10666429 DOI: 10.1186/s12864-023-09763-x] [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/06/2023] [Accepted: 10/23/2023] [Indexed: 11/25/2023] Open
Abstract
BACKGROUND The population of the Faroe Islands is an isolated population but very little is known about it from whole genome sequencing. The population of about 50000 people has a high incidence of rare diseases e.g., 1:300 for Primary Carnitine Deficiency. A screening programme was implemented, and eleven persons were also whole genome sequenced at x37 coverage for diagnostic purposes of those cases that were not affected by the known mutations. The purpose of our study is to utilize the high coverage data to explore the genomic variation and the ancestral history of the population. We study the SNP heterozygosity, the pairwise relatedness from kinship, the inbreeding from runs of homozygosity ROH, and we find the minor allele frequency distribution. We estimate the population ancestry and the timing of the founding event by using the whole genomes from eight consenting individuals. RESULTS We find the number of SNPs and the heterozygosity for the eight individual samples, and for merged samples, for which we also study the relatedness. We find close relatedness between the supposedly unrelated individuals. From ROH, we interpret the high relatedness as an ancient property of the isolated population. A bottleneck event is estimated starting between years [Formula: see text] with a maximum consanguineous population in year [Formula: see text] and similarly consanguineous between years [Formula: see text]. The ancestry analysis shows the population descends from founders of [Formula: see text] European and [Formula: see text] Admixed American ancestry. A distinct clustering near the central European and British populations of the 1000 Genome Project is likely the result of the population isolation and genetic drift. The minor allele frequency distribution suggests many rare variants. CONCLUSIONS The ancestry is mainly European while the inbreeding is higher compared to European populations and population isolates. The Faroese population has inbreeding more like ancient Europeans. We discovered a bottlenecked and consanguineous population event and estimated it starting in the 1st-4th century as compared to the oldest archaeological findings from the 4th-6th century.
Collapse
Affiliation(s)
- Hannes Gislason
- Faculty of Science and Technology, University of the Faroe Islands, Tórshavn, Faroe Islands.
| |
Collapse
|
17
|
Saif-Ur-Rehman M, Hassan FU, Reecy J, Deng T. Whole-genome SNP markers reveal runs of homozygosity in indigenous cattle breeds of Pakistan. Anim Biotechnol 2023; 34:1384-1396. [PMID: 35044288 DOI: 10.1080/10495398.2022.2026369] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/01/2022]
Abstract
The runs of homozygosity (ROH) were identified in 14 Pakistani cattle breeds (n = 105) by genotyping with the Illumina 50 K SNP BeadChip. These breeds were categorized into Dairy, Dual, and Draft breeds based on their utility and production performance. We identified a total of 10,936 ROHs which mainly consisted of a high number of shorter segments (1-4 Mb). Dairy group exhibited the highest level of inbreeding (FROH: 0.078 ± 0.028) while the lowest (FROH: 0.002 ± 0.008) was observed in Dual group. In 48 genomic regions identified with a high frequency of ROH, 207 genes were detected in the three breed groups. A substantially higher number of ROH islands detected in dairy breeds indicated the impact of the positive selection pressure over the years. Important candidate genes and QTL were detected in the ROH islands associated with economic traits like milk production, reproduction, meat, carcass, and health traits in dairy cattle.
Collapse
Affiliation(s)
| | - Faiz-Ul Hassan
- Institute of Animal and Dairy Sciences, University of Agriculture, Faisalabad, Pakistan
| | - James Reecy
- Department of Animal Science, Iowa State University, Ames, IA, USA
| | - Tingxian Deng
- Guangxi Provincial Key Laboratory of Buffalo Genetics, Breeding and Reproduction Technology, Buffalo Research Institute, Chinese Academy of Agricultural Sciences, Nanning, China
| |
Collapse
|
18
|
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
|
19
|
Lindsay-McGee V, Sanchez-Molano E, Banos G, Clark EL, Piercy RJ, Psifidi A. Genetic characterisation of the Connemara pony and the Warmblood horse using a within-breed clustering approach. Genet Sel Evol 2023; 55:60. [PMID: 37592264 PMCID: PMC10436415 DOI: 10.1186/s12711-023-00827-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2022] [Accepted: 07/07/2023] [Indexed: 08/19/2023] Open
Abstract
BACKGROUND The Connemara pony (CP) is an Irish breed that has experienced varied selection by breeders over the last fifty years, with objectives ranging from the traditional hardy pony to an agile athlete. We compared these ponies with well-studied Warmblood (WB) horses, which are also selectively bred for athletic performance but with a much larger census population. Using genome-wide single nucleotide polymorphism (SNP) and whole-genome sequencing data from 116 WB (94 UK WB and 22 European WB) and 36 CP (33 UK CP and 3 US CP), we studied the genomic diversity, inbreeding and population structure of these breeds. RESULTS The k-means clustering approach divided both the CP and WB populations into four genetic groups, among which the CP genetic group 1 (C1) associated with non-registered CP, C4 with US CP, WB genetic group 1 (W1) with Holsteiners, and W3 with Anglo European and British WB. Maximum and mean linkage disequilibrium (LD) varied significantly between the two breeds (mean from 0.077 to 0.130 for CP and from 0.016 to 0.370 for WB), but the rate of LD decay was generally slower in CP than WB. The LD block size distribution peaked at 225 kb for all genetic groups, with most of the LD blocks not exceeding 1 Mb. The top 0.5% harmonic mean pairwise fixation index (FST) values identified ontology terms related to cancer risk when the four CP genetic groups were compared. The four CP genetic groups were less inbred than the WB genetic groups, but C2, C3 and C4 had a lower proportion of shorter runs of homozygosity (ROH) (74 to 76% < 4 Mb) than the four WB genetic groups (80 to 85% < 4 Mb), indicating more recent inbreeding. The CP and WB genetic groups had a similar ratio of effective number of breeders (Neb) to effective population size (Ne). CONCLUSIONS Distinct genetic groups of individuals were revealed within each breed, and in WB these genetic groups reflected population substructure better than studbook or country of origin. Ontology terms associated with immune and inflammatory responses were identified from the signatures of selection between CP genetic groups, and while CP were less inbred than WB, the evidence pointed to a greater degree of recent inbreeding. The ratio of Neb to Ne was similar in CP and WB, indicating the influence of popular sires is similar in CP and WB.
Collapse
Affiliation(s)
- Victoria Lindsay-McGee
- Royal Veterinary College, London, UK
- Royal (Dick) School of Veterinary Studies, University of Edinburgh, Edinburgh, UK
| | | | | | - Emily L Clark
- The Roslin Institute, University of Edinburgh, Edinburgh, UK
| | | | | |
Collapse
|
20
|
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
|
21
|
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
|
22
|
Tolone M, Sardina MT, Criscione A, Lasagna E, Senczuk G, Rizzuto I, Riggio S, Moscarelli A, Macaluso V, Di Gerlando R, Cassandro M, Portolano B, Mastrangelo S. High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources. Poult Sci 2023; 102:102692. [PMID: 37120867 PMCID: PMC10172703 DOI: 10.1016/j.psj.2023.102692] [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/26/2023] [Revised: 03/21/2023] [Accepted: 03/29/2023] [Indexed: 05/02/2023] Open
Abstract
Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.
Collapse
Affiliation(s)
- Marco Tolone
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Maria Teresa Sardina
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Andrea Criscione
- Department of Agriculture, Food and the Environment, University of Catania, 95131 Catania, Italy
| | - Emiliano Lasagna
- Department of Agricultural, Food and Environmental Sciences, University of Perugia, 06121 Perugia, Italy
| | - Gabriele Senczuk
- Department of Agricultural, Environmental and Food Sciences, University of Molise, 86100 Campobasso, Italy
| | - Ilaria Rizzuto
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Silvia Riggio
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Angelo Moscarelli
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Vito Macaluso
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Rosalia Di Gerlando
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Martino Cassandro
- Department of Agronomy, Food, Natural Resources, Animals and Environment, University of Padova, 35020 Legnaro, Italy
| | - Baldassare Portolano
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy
| | - Salvatore Mastrangelo
- Department of Agricultural, Food and Forest Sciences, University of Palermo, 90128 Palermo, Italy.
| |
Collapse
|
23
|
Wang Q, Zhang J, Wang H, Wang Z, Li Q, Zhao G, Zheng M, Wen J. Estimates of genomic inbreeding and identification of candidate regions in Beijing-You chicken populations. Anim Genet 2023; 54:155-165. [PMID: 36541281 DOI: 10.1111/age.13286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Revised: 10/09/2022] [Accepted: 12/02/2022] [Indexed: 12/24/2022]
Abstract
Runs of homozygosity (ROHs) has become an effective method for analysing inbreeding in livestock populations. Moreover, ROHs is well-suited to detect signatures of selection via ROH islands. This study aimed to investigate the occurrence and distribution of ROHs, compare the genomic inbreeding coefficients and identify the genomic regions with high ROH frequencies in different Beijing-You chicken (BY) populations, including a random conservation population (BY_R), a pedigree conservation population (BY_P), and a commercial population obtained from the market (BY_C). Among them, BY_R in 2010 and 2019 were BY_R1 and BY_R2 respectively. A total of 27 916 ROHs were identified. The average number of ROHs per individual across the three BY populations ranged from 213 (BY_P) to 161 (BY_C), and the average length of ROHs ranged from 0.432 Mb (BY_R2) to 0.451 Mb (BY_P). The highest inbreeding coefficient calculated based on ROHs (FROH ) was 0.1 in BY_P, whereas the lowest FROH was 0.0743 in BY_C. In addition, the inbreeding coefficient of BY_R2 (FROH = 0.0798) was higher than that of BY_R1 (FROH = 0.0579). Furthermore, the highest proportion of long ROH fragments (>4 Mb) was observed in BY_P and BY_C. This study showed the top 10 ROH islands of each population, and these ROH islands harboured 53 genes, some of which were related to limb development, body size and immune response. These findings contribute to the understanding of genetic diversity and population demography, and might help improve breeding and conservation strategies for BY populations.
Collapse
Affiliation(s)
- Qiao Wang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Jin Zhang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Hailong Wang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Zixuan Wang
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Qinghe Li
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Guiping Zhao
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Maiqing Zheng
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| | - Jie Wen
- Institute of Animal Sciences of Chinese Academy of Agricultural Sciences, State Key Laboratory of Animal Nutrition, Beijing, China
| |
Collapse
|
24
|
Visser C, Lashmar SF, Reding J, Berry DP, van Marle-Köster E. Pedigree and genome-based patterns of homozygosity in the South African Ayrshire, Holstein, and Jersey breeds. Front Genet 2023; 14:1136078. [PMID: 37007942 PMCID: PMC10063850 DOI: 10.3389/fgene.2023.1136078] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/02/2023] [Accepted: 03/02/2023] [Indexed: 03/19/2023] Open
Abstract
The erosion of genetic diversity limits long-term genetic gain and impedes the sustainability of livestock production. In the South African (SA) dairy industry, the major commercial dairy breeds have been applying estimated breeding values (EBVs) and/or have been participating in Multiple Across Country Evaluations (MACE). The transition to genomic estimated breeding values (GEBVs) in selection strategies requires monitoring of the genetic diversity and inbreeding of current genotyped animals, especially considering the comparatively small population sizes of global dairy breeds in SA. This study aimed to perform a homozygosity-based evaluation of the SA Ayrshire (AYR), Holstein (HST), and Jersey (JER) dairy cattle breeds. Three sources of information, namely 1) single nucleotide polymorphism (SNP) genotypes (3,199 animals genotyped for 35,572 SNPs) 2) pedigree records (7,885 AYR; 28,391 HST; 18,755 JER), and 3) identified runs of homozygosity (ROH) segments were used to quantify inbreeding related parameters. The lowest pedigree completeness was for the HST population reducing from a value of 0.990 to 0.186 for generation depths of one to six. Across all breeds, 46.7% of the detected ROH were between 4 megabase pairs (Mb) and 8 Mb in length. Two conserved homozygous haplotypes were identified in more than 70% of the JER population on Bos taurus autosome (BTA) 7. The JER breed displayed the highest level of inbreeding across all inbreeding coefficients. The mean (± standard deviation) pedigree-based inbreeding coefficient (FPED) ranged from 0.051 (±0.020) for AYR to 0.062 (±0.027) for JER, whereas SNP-based inbreeding coefficients (FSNP) ranged from 0.020 (HST) to 0.190 (JER) and ROH-based inbreeding coefficients, considering all ROH segment coverage (FROH), ranged from 0.053 (AYR) to 0.085 (JER). Within-breed Spearman correlations between pedigree-based and genome-based estimates ranged from weak (AYR: 0.132 between FPED and FROH calculated for ROH <4Mb in size) to moderate (HST: 0.584 between FPED and FSNP). Correlations strengthened between FPED and FROH as the ROH length category was considered lengthened, suggesting a dependency on breed-specific pedigree depth. The genomic homozygosity-based parameters studied proved useful in investigating the current inbreeding status of reference populations genotyped to implement genomic selection in the three most prominent South African dairy cattle breeds.
Collapse
Affiliation(s)
- Carina Visser
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
- *Correspondence: Carina Visser,
| | - Simon Frederick Lashmar
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
| | - Jason Reding
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
| | - Donagh P. Berry
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
- Animal and Grassland Research and Innovation Centre, Teagasc, Co. Cork, Ireland
| | - Esté van Marle-Köster
- Department of Animal Science, Faculty of Natural and Agricultural Sciences, University of Pretoria, Pretoria, South Africa
| |
Collapse
|
25
|
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
|
26
|
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
|
27
|
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
|
28
|
Abumadini MS, Al Ghamdi KS, Alqahtani AH, Almedallah DK, Callans L, Jarad JA, Cyrus C, Koeleman BPC, Keating BJ, Pankratz N, Al-Ali AK. Genome-wide copy number variant screening of Saudi schizophrenia patients reveals larger deletions in cases versus controls. Front Mol Neurosci 2023; 16:1069375. [PMID: 36846569 PMCID: PMC9950097 DOI: 10.3389/fnmol.2023.1069375] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2022] [Accepted: 01/20/2023] [Indexed: 02/12/2023] Open
Abstract
Introduction Genome-wide association studies have discovered common polymorphisms in regions associated with schizophrenia. No genome-wide analyses have been performed in Saudi schizophrenia subjects. Methods Genome-wide genotyping data from 136 Saudi schizophrenia cases and 97 Saudi controls in addition to 4,625 American were examined for copy number variants (CNVs). A hidden Markov model approach was used to call CNVs. Results CNVs in schizophrenia cases were twice as large on average than CNVs in controls (p = 0.04). The analyses focused on extremely large >250 kilobases CNVs or homozygous deletions of any size. One extremely large deletion was noted in a single case (16.5 megabases on chromosome 10). Two cases had an 814 kb duplication of chromosome 7 spanning a cluster of genes, including circadian-related loci, and two other cases had 277 kb deletions of chromosome 9 encompassing an olfactory receptors gene family. CNVs were also seen in loci previously associated with schizophrenia, namely a 16p11 proximal duplication and two 22q11.2 deletions. Discussion Runs of homozygosity (ROHs) were analyzed across the genome to investigate correlation with schizophrenia risk. While rates and sizes of these ROHs were similar in cases and controls, we identified 10 regions where multiple cases had ROHs and controls did not.
Collapse
Affiliation(s)
- Mahdi S. Abumadini
- Department of Psychiatry, King Fahd Hospital of the University, Al-Khobar and College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
| | - Kholoud S. Al Ghamdi
- Department of Physiology, College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
| | - Abdullah H. Alqahtani
- Department of Psychiatry, King Fahd Hospital of the University, Al-Khobar and College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
| | - Dana K. Almedallah
- Department of Psychiatry, King Fahd Hospital of the University, Al-Khobar and College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
| | - Lauren Callans
- Department of Surgery, University of Pennsylvania School of Medicine, Philadelphia, PA, United States
| | - Jumanah A. Jarad
- Department of Psychiatry, King Fahd Hospital of the University, Al-Khobar and College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
| | - Cyril Cyrus
- Department of Clinical Biochemistry, College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia
| | - Bobby P. C. Koeleman
- Department of Genetics, Division Lab, University Medical Center Utrecht, Utrecht, Netherlands
| | - Brendan J. Keating
- Department of Surgery, University of Pennsylvania School of Medicine, Philadelphia, PA, United States
| | - Nathan Pankratz
- Department of Laboratory Medicine and Pathology, University of Minnesota, Minneapolis, MN, United States
| | - Amein K. Al-Ali
- Department of Clinical Biochemistry, College of Medicine, Imam Abdulrahman Bin Faisal University, Dammam, Saudi Arabia,*Correspondence: Amein K. Al-Ali, ✉
| |
Collapse
|
29
|
Pedigree reconstruction and population structure using SNP markers in Gir cattle. J Appl Genet 2023; 64:329-340. [PMID: 36645582 DOI: 10.1007/s13353-023-00747-x] [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: 07/20/2022] [Revised: 11/24/2022] [Accepted: 01/03/2023] [Indexed: 01/17/2023]
Abstract
Our objective was to establish a SNPs panel for pedigree reconstruction using microarrays of different densities and evaluate the genomic relationship coefficient of the inferred pedigree, in addition to analyzing the population structure based on genomic analyses in Gir cattle. For parentage analysis and genomic relationship, 16,205 genotyped Gir animals (14,458 females and 1747 males) and 1810 common markers to the four SNP microarrays were used. For population structure analyses, including linkage disequilibrium, effective population size, and runs of homozygosity (ROH), genotypes from 21,656 animals were imputed. Likelihood ratio (LR) approach was used to reconstruct the pedigree, deepening the pedigree and showing it is well established in terms of recent information. Coefficients for each relationship category of the inferred pedigree were adequate. Linkage disequilibrium showed rapid decay. We detected a decrease in the effective population size over the last 50 generations, with the average generation interval around 9.08 years. Higher ROH-based inbreeding coefficient in a class of short ROH segments, with moderate to high values, was also detected, suggesting bottlenecks in the Gir genome. Breeding strategies to minimize inbreeding and avoid massive use of few proven sires with high genetic value are suggested to maintain genetic variability in future generations. In addition, we recommend reducing the generation interval to maximize genetic progress and increase effective population size.
Collapse
|
30
|
Fieder M, Huber S. Facial attractiveness is only weakly linked to genome-wide heterozygosity. Front Psychol 2023; 14:1009962. [PMID: 37151335 PMCID: PMC10157054 DOI: 10.3389/fpsyg.2023.1009962] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2022] [Accepted: 03/30/2023] [Indexed: 05/09/2023] Open
Abstract
Introduction It has been frequently suggested that overall genomic heterozygosity and, particularly, heterozygosity of loci on the so-called major histocompatibility complex (MHC), which is responsible for the recognition of foreign substances/ pathogens and the recognition of self and non-self, is associated with better health and better resistance to infections and parasites. It has further been speculated that such a potentially beneficial heterozygosity can be detected through body odor and facial attractiveness. Methods On the basis of genome wide SNP data (713,014 SNPs) of participants from the Wisconsin Longitudinal Study, we therefore investigated whether homozygosity either on the MHC (measured as inbreeding coefficient) or genome-wide (measured as runs of homozygosity and inbreeding coefficient) is associated with rated facial attractiveness. Results Although we found that the genome-wide average length of homozygous segments and the genome-wide inbreeding coefficient are significantly negatively associated with some measures of facial attractiveness, if corrected for multiple testing, any significant association was no longer formally significant after correction. In addition, the variance in facial attractiveness explained by the genome wide homozygosity is very low (<0.15%). We did not find any significant association between the inbreeding coefficient on the MHC and facial attractiveness. Discussion We only find a weak association of genome- wide heterozygosity and facial attractiveness.
Collapse
Affiliation(s)
- Martin Fieder
- Department of Evolutionary Anthropology and Network of Human Evolution and Archeological Sciences (HEAS), University of Vienna, Vienna, Austria
- Religion and Transformation in Contemporary Society, University of Vienna, Vienna, Austria
- *Correspondence: Martin Fieder,
| | - Susanne Huber
- Department of Evolutionary Anthropology and Network of Human Evolution and Archeological Sciences (HEAS), University of Vienna, Vienna, Austria
| |
Collapse
|
31
|
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
|
32
|
Lashmar SF, Visser C, Okpeku M, Muchadeyi FC, Mapholi NO, van Marle-Köster E. A within- and across-country assessment of the genomic diversity and autozygosity of South African and eSwatini Nguni cattle. Trop Anim Health Prod 2022; 54:365. [DOI: 10.1007/s11250-022-03373-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Accepted: 10/24/2022] [Indexed: 11/06/2022]
|
33
|
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
|
34
|
Perdomo-González DI, Laseca N, Demyda-Peyrás S, Valera M, Cervantes I, Molina A. Fine-tuning genomic and pedigree inbreeding rates in equine population with a deep and reliable stud book: the case of the Pura Raza Española horse. J Anim Sci Biotechnol 2022; 13:127. [PMID: 36336696 PMCID: PMC9639299 DOI: 10.1186/s40104-022-00781-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2022] [Accepted: 09/13/2022] [Indexed: 11/09/2022] Open
Abstract
BACKGROUND Estimating inbreeding, which is omnipresent and inevitable in livestock populations, is a primary goal for management and animal breeding especially for those interested in mitigating the negative consequences of inbreeding. Inbreeding coefficients have been historically estimated by using pedigree information; however, over the last decade, genome-base inbreeding coefficients have come to the forefront in this field. The Pura Raza Española (PRE) horse is an autochthonous Spanish horse breed which has been recognised since 1912. The total PRE population (344,718 horses) was used to estimate Classical (F), Ballou's ancestral, Kalinowski's ancestral, Kalinowski's new and the ancestral history coefficient values. In addition, genotypic data from a selected population of 805 PRE individuals was used to determine the individual inbreeding coefficient using SNP-by-SNP-based techniques (methods of moments -FHOM-, the diagonal elements of the genomic -FG-, and hybrid matrixes -FH-) and ROH measures (FRZ). The analyse of both pedigree and genomic based inbreeding coefficients in a large and robust population such as the PRE horse, with proven parenteral information for the last 40 years and a high degree of completeness (over 90% for the last 70 years) will allow us to understand PRE genetic variability better and the correlations between the estimations will give the data greater reliability. RESULTS The mean values of the pedigree-based inbreeding coefficients ranged from 0.01 (F for the last 3 generations -F3-) to 0.44 (ancestral history coefficient) and the mean values of genomic-based inbreeding coefficients varied from 0.05 (FRZ for three generations, FH and FHOM) to 0.11 (FRZ for nine generations). Significant correlations were also found between pedigree and genomic inbreeding values, which ranged between 0.58 (F3 with FHOM) and 0.79 (F with FRZ). In addition, the correlations between FRZ estimated for the last 20 generations and the pedigree-based inbreeding highlight the fact that fewer generations of genomic data are required when comparing total inbreeding values, and the opposite when ancient values are calculated. CONCLUSIONS Ultimately, our results show that it is still useful to work with a deep and reliable pedigree in pedigree-based genetic studies with very large effective population sizes. Obtaining a satisfactory parameter will always be desirable, but the approximation obtained with a robust pedigree will allow us to work more efficiently and economically than with massive genotyping.
Collapse
Affiliation(s)
- Davinia Isabel Perdomo-González
- Departamento Agronomía, Escuela Técnica Superior de Ingeniería Agromómica, Universidad de Sevilla, Ctra Utrera Km 1, 41013, Sevilla, Spain.
| | - Nora Laseca
- Departamento de Genética, Universidad de Córdoba, Córdoba, Spain
| | - Sebastián Demyda-Peyrás
- Facultad de Ciencias Veterinarias, Universidad Nacional de La Plata, La Plata, Argentina.,Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), La Plata, Argentina
| | - Mercedes Valera
- Departamento Agronomía, Escuela Técnica Superior de Ingeniería Agromómica, Universidad de Sevilla, Ctra Utrera Km 1, 41013, Sevilla, Spain
| | - Isabel Cervantes
- Departamento de Producción Animal, Facultad de Veterinaria, Universidad Complutense de Madrid, Madrid, Spain
| | - Antonio Molina
- Departamento de Genética, Universidad de Córdoba, Córdoba, Spain
| |
Collapse
|
35
|
Zhang X, Jamwal K, Distl O. Tracking footprints of artificial and natural selection signatures in breeding and non-breeding cats. Sci Rep 2022; 12:18061. [PMID: 36302822 PMCID: PMC9613910 DOI: 10.1038/s41598-022-22155-7] [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: 03/17/2022] [Accepted: 10/10/2022] [Indexed: 01/24/2023] Open
Abstract
Stray non-breeding cats (stray) represent the largest heterogeneous cat population subject to natural selection, while populations of the Siamese (SIAM) and Oriental Shorthair (OSH) breeds developed through intensive artificial selection for aesthetic traits. Runs of homozygosity (ROH) and demographic measures are useful tools to discover chromosomal regions of recent selection and to characterize genetic diversity in domestic cat populations. To achieve this, we genotyped 150 stray and 26 household non-breeding cats (household) on the Illumina feline 63 K SNP BeadChip and compared them to SIAM and OSH. The 50% decay value of squared correlation coefficients (r2) in stray (0.23), household (0.25), OSH (0.24) and SIAM (0.25) corresponded to a mean marker distance of 1.12 Kb, 4.55 Kb, 62.50 Kb and 175.07 Kb, respectively. The effective population size (Ne) decreased in the current generation to 55 in stray, 11 in household, 9 in OSH and 7 in SIAM. In the recent generation, the increase in inbreeding per generation (ΔF) reached its maximum values of 0.0090, 0.0443, 0.0561 and 0.0710 in stray, household, OSH and SIAM, respectively. The genomic inbreeding coefficient (FROH) based on ROH was calculated for three length categories. The FROH was between 0.014 (FROH60) and 0.020 (FROH5) for stray, between 0.018 (FROH60) and 0.024 (FROH5) for household, between 0.048 (FROH60) and 0.069 (FROH5) for OSH and between 0.053 (FROH60) and 0.073 (FROH5) for SIAM. We identified nine unique selective regions for stray through genome-wide analyses for regions with reduced heterozygosity based on FST statistics. Genes in these regions have previously been associated with reproduction (BUB1B), motor/neurological behavior (GPHN, GABRB3), cold-induced thermogenesis (DIO2, TSHR), immune system development (TSHR), viral carcinogenesis (GTF2A1), host immune response against bacteria, viruses, chemoattractant and cancer cells (PLCB2, BAHD1, TIGAR), and lifespan and aging (BUB1B, FGF23). In addition, we identified twelve unique selective regions for OSH containing candidate genes for a wide range of coat colors and patterns (ADAMTS20, KITLG, TYR, TYRO3-a MITF regulator, GPNMB, FGF7, RAB38) as well as congenital heart defects (PDE4D, PKP2) and gastrointestinal disorders (NLGN1, ALDH1B1). Genes in stray that represent unique selective events indicate, at least in part, natural selection for environmental adaptation and resistance to infectious disease, and should be the subject of future research. Stray cats represent an important genetic resource and have the potential to become a research model for disease resistance and longevity, which is why we recommend preserving semen before neutering.
Collapse
Affiliation(s)
- Xuying Zhang
- grid.412970.90000 0001 0126 6191Institute for Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany
| | - Kokila Jamwal
- grid.412970.90000 0001 0126 6191Institute for Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany
| | - Ottmar Distl
- grid.412970.90000 0001 0126 6191Institute for Animal Breeding and Genetics, University of Veterinary Medicine Hannover, Hannover, Germany
| |
Collapse
|
36
|
Integrative QTL mapping and selection signatures in Groningen White Headed cattle inferred from whole-genome sequences. PLoS One 2022; 17:e0276309. [PMID: 36288367 PMCID: PMC9605288 DOI: 10.1371/journal.pone.0276309] [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: 08/04/2021] [Accepted: 10/04/2022] [Indexed: 11/04/2022] Open
Abstract
Here, we aimed to identify and characterize genomic regions that differ between Groningen White Headed (GWH) breed and other cattle, and in particular to identify candidate genes associated with coat color and/or eye-protective phenotypes. Firstly, whole genome sequences of 170 animals from eight breeds were used to evaluate the genetic structure of the GWH in relation to other cattle breeds by carrying out principal components and model-based clustering analyses. Secondly, the candidate genomic regions were identified by integrating the findings from: a) a genome-wide association study using GWH, other white headed breeds (Hereford and Simmental), and breeds with a non-white headed phenotype (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel, Dutch Belted, and Holstein Friesian); b) scans for specific signatures of selection in GWH cattle by comparison with four other Dutch traditional breeds (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel and Dutch Belted) and the commercial Holstein Friesian; and c) detection of candidate genes identified via these approaches. The alignment of the filtered reads to the reference genome (ARS-UCD1.2) resulted in a mean depth of coverage of 8.7X. After variant calling, the lowest number of breed-specific variants was detected in Holstein Friesian (148,213), and the largest in Deep Red (558,909). By integrating the results, we identified five genomic regions under selection on BTA4 (70.2-71.3 Mb), BTA5 (10.0-19.7 Mb), BTA20 (10.0-19.9 and 20.0-22.7 Mb), and BTA25 (0.5-9.2 Mb). These regions contain positional and functional candidate genes associated with retinal degeneration (e.g., CWC27 and CLUAP1), ultraviolet protection (e.g., ERCC8), and pigmentation (e.g. PDE4D) which are probably associated with the GWH specific pigmentation and/or eye-protective phenotypes, e.g. Ambilateral Circumocular Pigmentation (ACOP). Our results will assist in characterizing the molecular basis of GWH phenotypes and the biological implications of its adaptation.
Collapse
|
37
|
Wang X, Li G, Jiang Y, Tang J, Fan Y, Ren J. Genomic insights into the conservation and population genetics of two Chinese native goat breeds. J Anim Sci 2022; 100:skac274. [PMID: 35998083 PMCID: PMC9585554 DOI: 10.1093/jas/skac274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Accepted: 08/20/2022] [Indexed: 11/14/2022] Open
Abstract
Chinese goats are an important group of goats worldwide. However, there are few studies on the conservation priority, genetic relationship, and potential gene flow between Chinese and global goat breeds. Here, we genotyped 239 goats from conservation populations of the Chinese Guangfeng and Ganxi breeds using the GoatSNP50 BeadChip. The conservation priority, population structure, selection signatures and introgression of these goats were analyzed in the context of 36 global goat breeds. First, we showed that Guangfeng and Ganxi goats had the largest effective population sizes across the global breeds 13 generations ago. Nevertheless, Ganxi goats have recently experienced a high degree of inbreeding, resulting in their conservation priority based on total gene and allelic diversities being lower than that of most other Chinese breeds (including Guangfeng goats). Population structure and admixture analyses showed that an average of 18% of Guangfeng genomic components were introgressed from Boer goats approximately 18-yr ago. Next, we reconstructed the subfamily structure of the core populations of Guangfeng and Ganxi goats, and proposed reasonable conservation strategies for inbreeding management. Moreover, a list of candidate genes under selection for fertility, immunity, growth, and meat quality were detected in Guangfeng and Ganxi goats. Finally, we identified some genes related to body development and reproduction, which were introgressed from Boer goats and may be beneficial for improving performance and productivity of Guangfeng goats. In conclusion, this study not only provides new insights into the conservation and utilization of Guangfeng and Ganxi goats but also enriches our understanding of artificial introgression from exotic goats into Chinese local goats.
Collapse
Affiliation(s)
- Xiaopeng Wang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Guixin Li
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Yongchuang Jiang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| | - Jianhong Tang
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
- Laboratory Animal Engineering Research Center of Ganzhou, Gannan Medical University, Ganzhou 341000, China
| | - Yin Fan
- Department of Animal Science, Jiangxi Biotech Vocational College, Nanchang 330200, China
| | - Jun Ren
- Guangdong Laboratory for Lingnan Modern Agriculture, College of Animal Science, South China Agricultural University, Guangzhou 510642, China
| |
Collapse
|
38
|
Bosse M, van Loon S. Challenges in quantifying genome erosion for conservation. Front Genet 2022; 13:960958. [PMID: 36226192 PMCID: PMC9549127 DOI: 10.3389/fgene.2022.960958] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 08/09/2022] [Indexed: 11/18/2022] Open
Abstract
Massive defaunation and high extinction rates have become characteristic of the Anthropocene. Genetic effects of population decline can lead populations into an extinction vortex, where declining populations show lower genetic fitness, in turn leading to lower populations still. The lower genetic fitness in a declining population due to a shrinking gene pool is known as genetic erosion. Three different types of genetic erosion are highlighted in this review: overall homozygosity, genetic load and runs of homozygosity (ROH), which are indicative of inbreeding. The ability to quantify genetic erosion could be a very helpful tool for conservationists, as it can provide them with an objective, quantifiable measure to use in the assessment of species at risk of extinction. The link between conservation status and genetic erosion should become more apparent. Currently, no clear correlation can be observed between the current conservation status and genetic erosion. However, the high quantities of genetic erosion in wild populations, especially in those species dealing with habitat fragmentation and habitat decline, may be early signs of deteriorating populations. Whole genome sequencing data is the way forward to quantify genetic erosion. Extra screening steps for genetic load and hybridization can be included, since they could potentially have great impact on population fitness. This way, the information yielded from genetic sequence data can provide conservationists with an objective genetic method in the assessment of species at risk of extinction. However, the great complexity of genome erosion quantification asks for consensus and bridging science and its applications, which remains challenging.
Collapse
Affiliation(s)
- Mirte Bosse
- Amsterdam Institute for Life and Environment (A-LIFE), Section Ecology and Evolution, Vrije Universiteit Amsterdam, Amsterdam, Netherlands
- Animal Breeding and Genomics, Wageningen University and Research, Wageningen, Netherlands
| | - Sam van Loon
- Amsterdam Institute for Life and Environment (A-LIFE), Section Ecology and Evolution, Vrije Universiteit Amsterdam, Amsterdam, Netherlands
| |
Collapse
|
39
|
Lozada-Soto EA, Tiezzi F, Jiang J, Cole JB, VanRaden PM, Maltecca C. Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle. J Dairy Sci 2022; 105:8956-8971. [PMID: 36153159 DOI: 10.3168/jds.2022-22116] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 06/20/2022] [Indexed: 11/19/2022]
Abstract
Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (FPED), marker-based using a genomic relationship matrix (FGRM), and segment-based using ROH (FROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using FPED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using FGRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using FROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.
Collapse
Affiliation(s)
| | - Francesco Tiezzi
- Department of Agriculture, Food, Environment and Forestry (DAGRI), University of Florence, 50144 Florence, Italy
| | - Jicai Jiang
- Department of Animal Science, North Carolina State University, Raleigh 27607
| | | | - Paul M VanRaden
- Animal Genomics and Improvement Laboratory, Henry A. Wallace Beltsville Agricultural Research Service, USDA, Beltsville, MD 20705
| | - Christian Maltecca
- Department of Animal Science, North Carolina State University, Raleigh 27607
| |
Collapse
|
40
|
Meiring C, Schurz H, van Helden P, Hoal E, Tromp G, Kinnear C, Kleynhans L, Glanzmann B, van Schalkwyk L, Miller M, Möller M. African wild dogs (Lycaon pictus) from the Kruger National Park, South Africa are currently not inbred but have low genomic diversity. Sci Rep 2022; 12:14979. [PMID: 36056068 PMCID: PMC9440078 DOI: 10.1038/s41598-022-19025-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 08/23/2022] [Indexed: 12/24/2022] Open
Abstract
African wild dogs (Lycaon pictus) have undergone severe population reductions and are listed as endangered on the International Union for Conservation of Nature Red List. Small, isolated populations have the potential to suffer from threats to their genetic diversity that may impact species viability and future survival. This study provides the first set of population-wide genomic data to address conservation concerns for this endangered species. Whole genome sequencing data were generated for 71 free-ranging African wild dogs from the Kruger National Park (KNP), South Africa, and used to estimate important population genomic parameters. Genomic diversity metrics revealed that variation levels were low; however, this African wild dog population showed low levels of inbreeding. Very few first- and second-order relationships were observed in this cohort, with most relationships falling into the third-order or distant category. Patterns of homozygosity could have resulted from historical inbreeding or a loss in genome variation due to a population bottleneck. Although the results suggest that this stronghold African wild dog population maintains low levels of inbreeding, likely due to their cooperative breeding system, it may lead to a continuous population decline when a reduced number of suitable mates are available. Consequently, the low genomic variation may influence species viability over time. This study highlights the importance of assessing population genomic parameters to set conservation priorities. Future studies should include the investigation of the potential of this endangered species to adapt to environmental changes considering the low genomic diversity in this population.
Collapse
Affiliation(s)
- Christina Meiring
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa.
| | - Haiko Schurz
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
| | - Paul van Helden
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
| | - Eileen Hoal
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
| | - Gerard Tromp
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
- South African Tuberculosis Bioinformatics Initiative (SATBBI), Faculty of Medicine and Health Sciences, Stellenbosch University, Francie van Zijl Drive, PO Box 241, Cape Town, 7500, South Africa
| | - Craig Kinnear
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
- Genomics Centre, South African Medical Research Council, Francie van Zijl Drive, PO Box 19070, Cape Town, 7500, South Africa
| | - Léanie Kleynhans
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
| | - Brigitte Glanzmann
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
- Genomics Centre, South African Medical Research Council, Francie van Zijl Drive, PO Box 19070, Cape Town, 7500, South Africa
| | - Louis van Schalkwyk
- Department of Agriculture, Land Reform and Rural Development, PO Box 12, Skukuza, 1350, South Africa
- Department of Veterinary Tropical Diseases, Faculty of Veterinary Science, University of Pretoria, Private Bag X04, Soutpan Road, Pretoria, 0110, South Africa
- Department of Migration, Max Planck Institute of Animal Behavior, Am Obstberg 1, 78315, Radolfzell, Germany
| | - Michele Miller
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
| | - Marlo Möller
- DSI-NRF Centre of Excellence for Biomedical Tuberculosis Research, South African Medical Research Council Centre for Tuberculosis Research, Division of Molecular Biology and Human Genetics, Faculty of Medicine and Health Sciences, Stellenbosch University, PO Box 241, Francie van Zijl Drive, Cape Town, 7500, South Africa
- Centre for Bioinformatics and Computational Biology, Stellenbosch University, Private bag X1, Merriman Avenue, Stellenbosch, 7600, South Africa
| |
Collapse
|
41
|
|
42
|
Chokoe TC, Hadebe K, Muchadeyi FC, Nephawe KA, Dzomba EF, Mphahlele TD, Matelele TC, Mtileni BJ. Conservation status and historical relatedness of South African communal indigenous goat populations using a genome-wide single-nucleotide polymorphism marker. Front Genet 2022; 13:909472. [PMID: 36017496 PMCID: PMC9395594 DOI: 10.3389/fgene.2022.909472] [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: 03/31/2022] [Accepted: 06/27/2022] [Indexed: 11/18/2022] Open
Abstract
Indigenous goats form the majority of populations in smallholder, low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support of communal farming. Effective population size (Ne), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Across populations, the current Ne of Gauteng was the lowest at 371 animals, while the historical Ne across populations suggests that the ancestor Ne has decreased by 53.86%, 44.58%, 42.16%, and 41.16% in Free State (FS), North West (NW), Limpopo (LP), and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (FROH > 5 Mb) were highest in FS (0.08 ± 0.09) and lowest in the Eastern Cape (EC) (0.02 ± 0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR, and MAPK3 were identified (occurring in over 20% of the samples) in FS (n = 37), GP (n = 42), and NW (n = 2) populations only. The mean length of ROH across populations was 7.76 Mb and ranged from 1.61 Mb in KwaZulu-Natal (KZN) to 98.05 Mb (GP and NW). The distribution of ROH according to their size showed that the majority (n = 1949) of the detected ROH were > 5 Mb in length compared to the other categories. Assuming two hypothetical ancestral populations, the populations from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP share a common origin but have substantial admixture from the EC and NW populations. The findings revealed that the occurrence of high Ne and autozygosity varied largely across breeds in communal indigenous goat populations at recent and ancient events when a genome-wide single-nucleotide polymorphism (SNP) marker was used. The use of Illumina goat SNP50K BeadChip shows that there was a migration route of communal indigenous goat populations from the northern part (LP) of South Africa to the eastern areas of the KZN that confirmed their historical relatedness and coincides with the migration periods of the Bantu nation.
Collapse
Affiliation(s)
- T. C. Chokoe
- Farm Animal Genetic Resources, Department of Agriculture, Land Reform and Rural Development, Pretoria, South Africa
- School of Agriculture & Environmental Sciences, University of Limpopo, Polokwane, South Africa
- *Correspondence: T. C. Chokoe,
| | - K. Hadebe
- Biotechnology Platform, Agricultural Research Council, Pretoria, South Africa
| | - F. C. Muchadeyi
- Biotechnology Platform, Agricultural Research Council, Pretoria, South Africa
| | - K. A. Nephawe
- Department of Animal Sciences, Tshwane University of Technology, Pretoria, South Africa
| | - E. F. Dzomba
- Discipline of Genetics, School of Life Sciences, University of Kwazulu-Natal, Scottsville, South African
| | - T. D. Mphahlele
- Farm Animal Genetic Resources, Department of Agriculture, Land Reform and Rural Development, Pretoria, South Africa
| | - T. C. Matelele
- Farm Animal Genetic Resources, Department of Agriculture, Land Reform and Rural Development, Pretoria, South Africa
| | - B. J. Mtileni
- Department of Animal Sciences, Tshwane University of Technology, Pretoria, South Africa
| |
Collapse
|
43
|
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
|
44
|
Mulim HA, Brito LF, Batista Pinto LF, Moletta JL, Da Silva LR, Pedrosa VB. Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems. Front Genet 2022; 13:858970. [PMID: 35923708 PMCID: PMC9341487 DOI: 10.3389/fgene.2022.858970] [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: 01/20/2022] [Accepted: 06/10/2022] [Indexed: 11/13/2022] Open
Abstract
Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from −0.009 to 0.029. A low (−0.12–0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (<4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5–1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.
Collapse
Affiliation(s)
| | - Luiz F. Brito
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
| | | | | | | | - Victor Breno Pedrosa
- Department of Animal Science, Federal University of Bahia, Salvador, Brazil
- Department of Animal Sciences, Purdue University, West Lafayette, IN, United States
- Department of Animal Science, State University of Ponta Grossa, Ponta Grossa, Brazil
- *Correspondence: Victor Breno Pedrosa,
| |
Collapse
|
45
|
Wang P, Hou R, Wu Y, Zhang Z, Que P, Chen P. Genomic status of yellow-breasted bunting following recent rapid population decline. iScience 2022; 25:104501. [PMID: 35733787 PMCID: PMC9207672 DOI: 10.1016/j.isci.2022.104501] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Revised: 05/04/2022] [Accepted: 05/26/2022] [Indexed: 12/03/2022] Open
Abstract
Global biodiversity is facing serious threats. However, knowledge of the genomic consequences of recent rapid population declines of wild organisms is limited. Do populations experiencing recent rapid population decline have the same genomic status as wild populations that experience long-term declines? Yellow-breasted Bunting (Emberiza aureola) is a critically endangered species that has been experiencing a recent rapid population decline. To answer the question, we assembled and annotated the whole genome of Yellow-breasted Bunting. Furthermore, we found high genetic diversity, low linkage disequilibrium, and low proportion of long runs of homozygosity in Yellow-breasted Bunting, suggesting that the populations following recent rapid declines have different genomic statuses from the population that experienced long-term population decline.
Collapse
Affiliation(s)
- Pengcheng Wang
- Jiangsu Key Laboratory for Biodiversity and Biotechnology, College of Life Sciences, Nanjing Normal University, Nanjing 210023, P. R. China
| | - Rong Hou
- Chengdu Research Base of Giant Panda Breeding, Sichuan Key Laboratory of Conservation Biology for Endangered Wildlife, Chengdu, 610081, P. R. China
| | - Yang Wu
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, P. R. China
| | - Zhengwang Zhang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing 100875, P. R. China
| | - Pinjia Que
- Chengdu Research Base of Giant Panda Breeding, Sichuan Key Laboratory of Conservation Biology for Endangered Wildlife, Chengdu, 610081, P. R. China
- Sichuan Academy of Giant Panda, Chengdu 610086, P. R. China
| | - Peng Chen
- Chengdu Research Base of Giant Panda Breeding, Sichuan Key Laboratory of Conservation Biology for Endangered Wildlife, Chengdu, 610081, P. R. China
- Sichuan Academy of Giant Panda, Chengdu 610086, P. R. China
| |
Collapse
|
46
|
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
|
47
|
Toro-Ospina AM, Herrera Rios AC, Pimenta Schettini G, Vallejo Aristizabal VH, Bizarria dos Santos W, Zapata CA, Ortiz Morea EG. Identification of Runs of Homozygosity Islands and Genomic Estimated Inbreeding Values in Caqueteño Creole Cattle (Colombia). Genes (Basel) 2022; 13:genes13071232. [PMID: 35886015 PMCID: PMC9318017 DOI: 10.3390/genes13071232] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2022] [Revised: 06/29/2022] [Accepted: 07/01/2022] [Indexed: 02/04/2023] Open
Abstract
The Caqueteño Creole (CAQ) is a native breed of cattle from the Caquetá department (Colombia), adapted to tropical conditions, which is extremely important to production systems in those regions. However, CAQ is poorly studied. In this sense, population structure studies associated with runs of homozygosity (ROH) analysis would allow for a better understanding of CAQ. Through ROH analysis, it is possible to reveal genetic relationships between individuals, measure genome inbreeding levels, and identify regions associated with traits of economic interest. Samples from a CAQ population (n = 127) were genotyped with the Bovine HD BeadChip (777,000 SNPs) and analyzed with the PLINK 1.9 program to estimate FROH and ROH islands. We highlighted a decrease in inbreeding frequency for FROH 4−8 Mb, 8−16 Mb, and >16 Mb classes, indicating inbreeding control in recent matings. We also found genomic hotspot regions on chromosomes 3, 5, 6, 8, 16, 20, and 22, where chromosome 20 harbored four hotspots. Genes in those regions were associated with fertility and immunity traits, muscle development, and environmental resistance, which may be present in the CAQ breed due to natural selection. This indicates potential for production systems in tropical regions. However, further studies are necessary to elucidate the CAQ production objective.
Collapse
Affiliation(s)
- Alejandra M. Toro-Ospina
- Amazonian Research Center CIMAZ-MACAGUAL, Laboratory of Agrobiotechnology, University of the Amazon, Florencia 180002, Colombia; (A.C.H.R.); (V.H.V.A.); (C.A.Z.); (E.G.O.M.)
- Correspondence:
| | - Ana C. Herrera Rios
- Amazonian Research Center CIMAZ-MACAGUAL, Laboratory of Agrobiotechnology, University of the Amazon, Florencia 180002, Colombia; (A.C.H.R.); (V.H.V.A.); (C.A.Z.); (E.G.O.M.)
- Science and Humanities Faculty, Digital University Institute of Antioquia, IUDigital, Medellin, Antioquia 50010, Colombia
| | - Gustavo Pimenta Schettini
- Department of Animal and Poultry Science, Virginia Polytechnic Institute and State University, Blacksburg, VA 24061-0002, USA;
| | - Viviana H. Vallejo Aristizabal
- Amazonian Research Center CIMAZ-MACAGUAL, Laboratory of Agrobiotechnology, University of the Amazon, Florencia 180002, Colombia; (A.C.H.R.); (V.H.V.A.); (C.A.Z.); (E.G.O.M.)
| | - Wellington Bizarria dos Santos
- School of Agricultural and Veterinary Sciences (FCAV), São Paulo State University (UNESP), Jaboticabal, Sao Paulo 14884-900, Brazil;
| | - Cesar A. Zapata
- Amazonian Research Center CIMAZ-MACAGUAL, Laboratory of Agrobiotechnology, University of the Amazon, Florencia 180002, Colombia; (A.C.H.R.); (V.H.V.A.); (C.A.Z.); (E.G.O.M.)
| | - Edna Gicela Ortiz Morea
- Amazonian Research Center CIMAZ-MACAGUAL, Laboratory of Agrobiotechnology, University of the Amazon, Florencia 180002, Colombia; (A.C.H.R.); (V.H.V.A.); (C.A.Z.); (E.G.O.M.)
| |
Collapse
|
48
|
da Cruz PRS, Ananina G, Secolin R, Gil-da-Silva-Lopes VL, Lima CSP, de França PHC, Donatti A, Lourenço GJ, de Araujo TK, Simioni M, Lopes-Cendes I, Costa FF, de Melo MB. Demographic history differences between Hispanics and Brazilians imprint haplotype features. G3 GENES|GENOMES|GENETICS 2022; 12:6576632. [PMID: 35511163 PMCID: PMC9258545 DOI: 10.1093/g3journal/jkac111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Accepted: 04/27/2022] [Indexed: 11/24/2022]
Abstract
Admixture is known to greatly impact the genetic landscape of a population and, while genetic variation underlying human phenotypes has been shown to differ among populations, studies on admixed subjects are still scarce. Latin American populations are the result of complex demographic history, such as 2 or 3-way admixing events, bottlenecks and/or expansions, and adaptive events unique to the American continent. To explore the impact of these events on the genetic structure of Latino populations, we evaluated the following haplotype features: linkage disequilibrium, shared identity by descent segments, runs of homozygosity, and extended haplotype homozygosity (integrated haplotype score) in Latinos represented in the 1000 Genome Project along with array data from 171 Brazilians sampled in the South and Southeast regions of Brazil. We found that linkage disequilibrium decay relates to the amount of American and African ancestry. The extent of identity by descent sharing positively correlates with historical effective population sizes, which we found to be steady or growing, except for Puerto Ricans and Colombians. Long runs of homozygosity, a particular instance of autozygosity, was only enriched in Peruvians and Native Americans. We used simulations to account for random sampling and linkage disequilibrium to filter positive selection indexes and found 244 unique markers under selection, 26 of which are common to 2 or more populations. Some markers exhibiting positive selection signals had estimated time to the most recent common ancestor consistent with human adaptation to the American continent. In conclusion, Latino populations present highly divergent haplotype characteristics that impact genetic architecture and underlie complex phenotypes.
Collapse
Affiliation(s)
- Pedro Rodrigues Sousa da Cruz
- Laboratory of Human Genetics, Center for Molecular Biology and Genetic Engineering (CBMEG), University of Campinas—UNICAMP , Campinas, SP 13083-875, Brazil
| | - Galina Ananina
- Laboratory of Human Genetics, Center for Molecular Biology and Genetic Engineering (CBMEG), University of Campinas—UNICAMP , Campinas, SP 13083-875, Brazil
| | - Rodrigo Secolin
- Department of Medical Genetics and Genomic Medicine, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
- The Brazilian Institute of Neuroscience and Neurotechnology (BRAINN) , Campinas, SP 13083-887, Brazil
| | - Vera Lúcia Gil-da-Silva-Lopes
- Department of Medical Genetics and Genomic Medicine, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
| | - Carmen Silvia Passos Lima
- Clinical Oncology Service, Department of Internal Medicine, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
| | | | - Amanda Donatti
- Department of Medical Genetics and Genomic Medicine, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
- The Brazilian Institute of Neuroscience and Neurotechnology (BRAINN) , Campinas, SP 13083-887, Brazil
| | - Gustavo Jacob Lourenço
- Laboratory of Cancer Genetics, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
| | - Tânia Kawasaki de Araujo
- Department of Medical Genetics and Genomic Medicine, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
| | - Milena Simioni
- Department of Medical Genetics and Genomic Medicine, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
| | - Iscia Lopes-Cendes
- Department of Medical Genetics and Genomic Medicine, School of Medical Sciences, University of Campinas—UNICAMP , Campinas, SP 13083-887, Brazil
- The Brazilian Institute of Neuroscience and Neurotechnology (BRAINN) , Campinas, SP 13083-887, Brazil
| | - Fernando Ferreira Costa
- Hematology and Hemotherapy Center, University of Campinas—UNICAMP, Campinas, SP, 13083-878 , Brazil
| | - Mônica Barbosa de Melo
- Laboratory of Human Genetics, Center for Molecular Biology and Genetic Engineering (CBMEG), University of Campinas—UNICAMP , Campinas, SP 13083-875, Brazil
| |
Collapse
|
49
|
Jiang Y, Li X, Liu J, Zhang W, Zhou M, Wang J, Liu L, Su S, Zhao F, Chen H, Wang C. Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data. BMC Genomics 2022; 23:373. [PMID: 35581549 PMCID: PMC9115978 DOI: 10.1186/s12864-022-08583-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2021] [Accepted: 04/22/2022] [Indexed: 11/25/2022] Open
Abstract
Background Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (FROH) and homozygosity (FHOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein. Results The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10–20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. FROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, FROH and FHOM values were highly correlated (0.952–0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits. Conclusion Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08583-9.
Collapse
Affiliation(s)
- Yao Jiang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Xiaojin Li
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Jiali Liu
- Key Laboratory of Animal Genetics, Breeding and Reproduction (Poultry) of Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Wei Zhang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Mei Zhou
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Jieru Wang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Linqing Liu
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Shiguang Su
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China
| | - Fuping Zhao
- Key Laboratory of Animal Genetics, Breeding and Reproduction (Poultry) of Ministry of Agriculture, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
| | - Hongquan Chen
- College of Animal Science and Technology, Anhui Agricultural University, Hefei, 230036, China
| | - Chonglong Wang
- Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, 230031, China.
| |
Collapse
|
50
|
Berton MP, da Silva RP, Banchero G, Mourão GB, Ferraz JBS, Schenkel FS, Baldi F. Genomic integration to identify molecular biomarkers associated with indicator traits of gastrointestinal nematode resistance in sheep. J Anim Breed Genet 2022; 139:502-516. [PMID: 35535437 DOI: 10.1111/jbg.12682] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 04/21/2022] [Indexed: 12/19/2022]
Abstract
This study aimed to integrate GWAS and structural variants to propose possible molecular biomarkers related to gastrointestinal nematode resistance traits in Santa Inês sheep. The phenotypic records FAMACHA, haematocrit, white blood cell count, red blood cell count, haemoglobin, platelets and egg counts per gram of faeces were collected from 700 naturally infected animals, belonging to four Brazilian flocks. A total of 576 animals were genotyped using the Ovine SNP12k BeadChip and were imputed using a reference population with Ovine SNP50 BeadChip. The GWAS approaches were based on SNPs, haplotypes, CNVs and ROH. The overlapping between the significant genomic regions detected from all approaches was investigated, and the results were integrated using a network analysis. Genes related to the immune system were found, such as ABCB1, IL6, WNT5A and IRF5. Genomic regions containing candidate genes and metabolic pathways involved in immune responses, inflammatory processes and immune cells affecting parasite resistance traits were identified. The genomic regions, biological processes and candidate genes uncovered could lead to biomarkers for selecting more resilient sheep and improving herd welfare and productivity. The results obtained are the start point to identify molecular biomarkers related to indicator traits of gastrointestinal nematode resistance in Santa Inês sheep.
Collapse
Affiliation(s)
- Mariana Piatto Berton
- Departamento de Zootecnia, Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista, Jaboticabal, Brazil
| | - Rosiane Pereira da Silva
- Faculdade de Zootecnia e Engenharia de Alimentos, Universidade de São Paulo, Pirassununga, Brazil
| | - Georgget Banchero
- Instituto Nacional de Investigación Agropecuária (INIA), Colonia, Uruguay
| | - Gerson Barreto Mourão
- Departamento de Zootecnia, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo/ESALQ, Piracicaba, Brazil
| | | | | | - Fernando Baldi
- Departamento de Zootecnia, Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista, Jaboticabal, Brazil
| |
Collapse
|