Genome-wide analysis of homozygosity regions in european simmental bulls

Unraveling inbreeding patterns and selection signals in Alpine Grey cattle

Inbreeding plays a crucial role in livestock breeding, influencing genetic diversity and phenotypic traits. Genomic data have helped address limitations posed by incomplete pedigrees, providing deeper insights into breed genetic diversity. This study assesses inbreeding levels via pedigree and genomic approaches and analyzes old and recent inbreeding using runs of homozygosity (ROH), and selection signals in Alpine Grey cattle. Pedigree data from 165 575 individuals, analyzed with INBUPGF90 software, computed inbreeding coefficients. Genomic-based coefficients derived from PLINK v1.9. or DetectRUNS R package analyses of 1 180 individuals' genotypes. Common single nucleotide polymorphisms within ROH pinpointed genomic regions, aggregating into "ROH islands" indicative of selection pressure. Overlaps with USCS Genome Browser unveiled gene presence. Moderate correlations (0.20-0.54) existed between pedigree and genomic coefficients, with most genomic estimators having higher (>0.8) correlation values. Inbreeding averaged 0.04 in < 8 Mb ROH segments, and 0.03 in > 16 Mb segments; > 90% of ROHs were < 8 Mb, indicating ancient inbreeding prevalence. Recent inbreeding proved less detrimental than in cosmopolitan breeds. Two major ROH islands on chromosomes 6 and 7 harbored genes linked to immune response, disease resistance (PYURF, HERC3), and fertility (EIF4EBP3, SRA1). This study underscores the need for detailed inbreeding analyses to understand genetic characteristics and historical changes in local breeds like Alpine Grey cattle. Genomic insights, especially from ROH, facilitated overcoming pedigree limitations, illuminating breed genetic diversity. Our findings reveal ancient inbreeding's enduring genetic impact and ROH islands potential for selective sweeps, elucidating traits in Alpine Grey cattle.

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

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.

Genome-Wide Detection for Runs of Homozygosity in Baoshan Pigs Using Whole Genome Resequencing

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.

Detection of Runs of Homozygosity and Identification of Candidate Genes in the Whole Genome of Tunchang Pigs

Tunchang pigs are an indigenous pig population in China known for their high tolerance to roughage, delicious meat, and fecundity. However, the number of Tunchang pigs has been declining due to the influence of commercial breeds and African swine fever, which could potentially lead to inbreeding. To assess the inbreeding level and the genetic basis of important traits in Tunchang pigs, our research investigated the patterns in "runs of homozygosity" (ROHs) using whole genome resequencing data from 32 Tunchang pigs. The study aimed to determine the length, number, coverage, and distribution model of ROHs in Tunchang pigs, as well as genomic regions with high ROH frequencies. The results of the study revealed that a total of 20,499,374 single-nucleotide polymorphisms (SNPs) and 1953 ROH fragments were recognized in 32 individuals. The ROH fragments in Tunchang pigs were predominantly short, ranging from 0.5 to 1 megabases (Mb) in length. Furthermore, the coverage of ROHs varied across chromosomes, with chromosome 3 having the highest coverage and chromosome 11 having the lowest coverage. The genetic diversity of Tunchang pigs was found to be relatively high based on the values of HE (expected heterozygosity), HO (observed heterozygosity), pi (nucleotide diversity), Ne (effective population size), and MAF (minor allele frequency). The average inbreeding coefficients of Tunchang pigs, as determined by three different methods (FHOM, FGRM, and FROH), were 0.019, 0.0138, and 0.0304, respectively. These values indicate that the level of inbreeding in Tunchang pigs is currently low. Additionally, the study identified a total of 13 ROH islands on all chromosomes, which in total contained 38,913 SNPs and 120 genes. These ROH islands included genes associated with economically important traits, including meat quality (GYS1, PHLPP1, SLC27A5, and CRTC1), growth and development (ANKS1A, TAF11, SPDEF, LHB, and PACSIN1), and environmental adaptation (SLC26A7). The findings of this research offer valuable perspectives on the present status of Tunchang pig resources and offer a reference for breeding conservation plans and the efficient utilization of Tunchang pigs in the future. By understanding the inbreeding level and genetic basis of important traits in Tunchang pigs, conservation efforts can be targeted towards maintaining genetic diversity and promoting the sustainable development of this indigenous pig population.

Analysis of runs of homozygosity of cattle living in different climate zones

Aim of this study was to analyze the distribution and characteristics of runs of homozygosity in Bos taurus taurus and Bos taurus indicus breeds, as well as their crosses, farmed all around the world. With this aim in view, we used single-nucleotide polymorphisms (SNP) genotypes for 3,263 cattle belonging to 204 different breeds. After quality control, 23,311 SNPs were retained for the analysis. Animals were divided into seven different groups: 1) continental taurus, 2) temperate taurus, 3) temperate indicus, 4) temperate composite, 5) tropical taurus, 6) tropical indicus, and 7) tropical composite. The climatic zones were created according to the latitude of the breeds' country of origin: i) continental, latitude ≥ 45°; ii) temperate, 45°< Latitude >23.26°; iii) tropics, latitude ≤ 23.26°. Runs of homozygosity were computed as 15 SNPs spanning in at least 2 Mb; number of ROH per animal (nROH), average ROH length (meanMb), and ROH-based inbreeding coefficients (FROH) were also computed. Temperate indicus showed the largest nROH, whereas Temperate taurus the lowest value. Moreover, the largest meanMb was observed for Temperate taurus, whereas the lowest value for Tropics indicus. Temperate indicus breeds showed the largest FROH values. Genes mapped in the identified ROH were reported to be associated with the environmental adaptation, disease resistance, coat color determinism, and production traits. Results of the present study confirmed that runs of homozygosity could be used to identify genomic signatures due to both artificial and natural selection.

Relationship between inbreeding and milk production traits in two Italian dairy sheep breeds

The effects of inbreeding in livestock species breeds have been well documented and they have a negative impact on profitability. The objective of this study was to evaluate the levels of inbreeding in Sarda (SAR, n = 785) and Valle del Belice (VdB, n = 473) dairy sheep breeds and their impact on milk production traits. Two inbreeding coefficients (F) were estimated: using pedigree (F_PED ), or runs of homozygosity (ROH; F_ROH ) at different minimum ROH lengths and different ROH classes. After the quality control, 38,779 single nucleotide polymorphisms remained for further analyses. A mixed-linear model was used to evaluate the impact of inbreeding coefficients on production traits within each breed. VdB showed higher inbreeding coefficients compared to SAR, with both breeds showing lower estimates as the minimum ROH length increased. Significant inbreeding depression was found only for milk yield, with a loss of around 7 g/day (for SAR) and 9 g/day (VdB) for a 1% increase of F_ROH . The present study confirms how the use of genomic information can be used to manage intra-breed diversity and to calculate the effects of inbreeding on phenotypic traits.

Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China

Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of F_ROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.

Assessment of Heterozygosity and Genome-Wide Analysis of Heterozygosity Regions in Two Duroc Pig Populations

Heterozygosity can effectively reflect the diverse models of population structure and demographic history. However, the genomic distribution of heterozygotes and the correlation between regions of heterozygosity (runs of heterozygosity, ROHet) and phenotypes are largely understudied in livestock. The objective of this study was to identify ROHet in the Duroc pig genome, and investigate the relationships between ROHet and eight important economic traits. Here, we genotyped 3,770 American Duroc (S21) and 2,096 Canadian Duroc (S22) pigs using 50 K single nucleotide polymorphism array to analyze heterozygosity. A total of 145,010 and 84,396 ROHets were characterized for S21 and S22 populations, respectively. ROHet segments were mostly enriched in 1–2 Mb length classification (75.48% in S21 and 72.25% in S22). The average genome length covered by ROHet was 66.53 ± 12.20 Mb in S21 and 73.32 ± 13.77 Mb in S22 pigs. Additionally, we detected 20 and 13 ROHet islands in S21 and S22 pigs. Genes in these genomic regions were mainly involved in the biological processes of immunity and reproduction. Finally, the genome-wide ROHet-phenotypes association analysis revealed that 130 ROHets of S21 and 84 ROHets of S22 were significantly associated with eight economic traits. Among the candidate genes in the significant ROHet regions, 16 genes related to growth, metabolism, and meat quality were considered as candidate genes for important economic traits of pigs. This work preliminarily explores the effect of heterozygosity-rich regions in the pig genome on production performance and provides new insights for subsequent research on pig genetic improvement.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.

Genomics Confirm an Alarming Status of the Genetic Diversity of Belgian Red and Belgian White Red Cattle

Simple Summary

Genetic diversity of livestock is vitally important as it enables the adaptation of future populations to changing environments. Therefore, preserving a sufficiently large genetic diversity is key. However, for many local cattle populations, little is known about their genetic diversity such as inbreeding level, effective size etc. We studied the genetic diversity of two local Belgian red cattle populations (Belgian Red and Belgian White Red cattle) using state-of-the-art genomic techniques. These tools assessed diversity at the population and individual level, and allowed the positioning of these two breeds in an international context of 52 other (European) cattle breeds. Accordingly, we contribute to the general knowledge of European red cattle, and more specifically we help the breeders, breed organization and the government to manage the genetic diversity of both breeds.

Abstract

Genetic diversity is increasingly important for researchers and society. Small and local populations deserve more attention especially, as they may harbor important characteristics. Moreover, small populations are at greater risk and their genetic management is often more challenging. Likewise, European red cattle populations are threatened, as they are outcompeted by more specialized cattle breeds. In this study, we investigate the genetic diversity of two local Belgian red cattle breeds: Belgian Red and Belgian White Red cattle. A total of 270 animals were genotyped via medium density SNP arrays. Genetic diversity was assessed using runs of homozygosity screening, effective population size estimation and F_st analyses. Genomic inbreeding coefficients based on runs of homozygosity were estimated at 7.0% for Belgian Red and 6.1% for Belgian White Red cattle, and both populations had a low effective population size (68 and 86, respectively). PCA, F_st and admixture analyses revealed the relationship to 52 other international breeds, where they were closest related to some Belgian, French, Scandinavian and Dutch breeds. Moreover, F_st analyses revealed for Belgian Red cattle a signature of selection on BTA6, adjacent to the KIT gene. This study gains important knowledge on the genetic diversity of these two small local red cattle breeds, and will aid in their (genetic) management.

Reverse Genetic Screen for Deleterious Recessive Variants in the Local Simmental Cattle Population of Switzerland

We herein report the result of a large-scale reverse genetic screen in the Swiss Simmental population, a local dual-purpose cattle breed. We aimed to detect possible recessively inherited variants affecting protein-coding genes, as such deleterious variants can impair fertility and rearing success significantly. We used 115,000 phased SNP data of almost 10 thousand cattle with pedigree data. This revealed evidence for 11 genomic regions of 1.17 Mb on average, with haplotypes (SH1 to SH11) showing a significant depletion in homozygosity and an allele frequency between 3.2 and 10.6%. For the proposed haplotypes, it was unfortunately not possible to evaluate associations with fertility traits as no corresponding data were available. For each haplotype region, possible candidate genes were listed based on their known function in development and disease. Subsequent mining of single-nucleotide variants and short indels in the genomes of 23 sequenced haplotype carriers allowed us to identify three perfectly linked candidate causative protein-changing variants: a SH5-related DIS3:p.Ile678fs loss-of-function variant, a SH8-related CYP2B6:p.Ile313Asn missense variant, and a SH9-related NUBPL:p.Ser143Tyr missense variant. None of these variants occurred in homozygous state in any of more than 5200 sequenced cattle of various breeds. Selection against these alleles in order to reduce reproductive failure and animal loss is recommended.

Inbreeding Calculated with Runs of Homozygosity Suggests Chromosome-Specific Inbreeding Depression Regions in Line 1 Hereford

The goal of this study was to evaluate inbreeding in a closed beef cattle population and assess phenotype prediction accuracy using inbreeding information. Effects of inbreeding on average daily gain phenotype in the Line 1 Hereford cattle population were assessed in this study. Genomic data were used to calculate inbreeding based on runs of homozygosity (ROH), and pedigree information was used to calculate the probability of an allele being identical by descent. Prediction ability of phenotypes using inbreeding coefficients calculated based on pedigree information and runs of homozygosity over the whole genome was close to 0, even in the case of significant inbreeding coefficient effects. On the other hand, inbreeding calculated per individual chromosomes' ROH yielded higher accuracies of prediction. Additionally, including only ROH from chromosomes with higher predicting ability further increased prediction accuracy. Phenotype prediction accuracy, inbreeding depression, and the effects of chromosome-specific ROHs varied widely across the genome. The results of this study suggest that inbreeding should be evaluated per individual regions of the genome. Moreover, mating schemes to avoid inbreeding depression should focus more on specific ROH with negative effects. Finally, using ROH as added information may increase prediction of the genetic merit of animals in a genomic selection program.

Whole-Genome Sequencing and Characterization of Buffalo Genetic Resources: Recent Advances and Future Challenges

The buffalo was domesticated around 3000-6000 years ago and has substantial economic significance as a meat, dairy, and draught animal. The buffalo has remained underutilized in terms of the development of a well-annotated and assembled reference genome de novo. It is mandatory to explore the genetic architecture of a species to understand the biology that helps to manage its genetic variability, which is ultimately used for selective breeding and genomic selection. Morphological and molecular data have revealed that the swamp buffalo population has strong geographical genomic diversity with low gene flow but strong phenotypic consistency, while the river buffalo population has higher phenotypic diversity with a weak phylogeographic structure. The availability of recent high-quality reference genome and genotyping marker panels has invigorated many genome-based studies on evolutionary history, genetic diversity, functional elements, and performance traits. The increasing molecular knowledge syndicate with selective breeding should pave the way for genetic improvement in the climatic resilience, disease resistance, and production performance of water buffalo populations globally.

Genetic Background and Inbreeding Depression in Romosinuano Cattle Breed in Mexico

Simple Summary

The objective of this study was to evaluate the genetic background and inbreeding depression in the Mexican Romosinuano cattle using pedigree and genomic information. Inbreeding was estimated using pedigree (FPED) and genomic information based on the genomic relationship matrix (FGRM) and runs of homozygosity (FROH). Linkage disequilibrium (LD) was evaluated using the correlation between pairs of loci, and the effective population size (Ne) was calculated based on LD and pedigree information. The pedigree file consisted of 4875 animals; 71 had genotypes. LD decreased with the increase in distance between markers, and Ne estimated using genomic information decreased from 610 to 72 animals (from 109 to 1 generation ago), the Ne estimated using pedigree information was 86.44. The number of runs of homozygosity per animal ranged between 18 and 102 segments with an average of 55. The average inbreeding was 2.98 ± 2.81, 2.98 ± 4.01, and 7.28 ± 3.68% for FPED, FGRM, and FROH, respectively. A 1% increase in inbreeding decreased birth weight by 0.103 kg and weaning weight by 0.685 kg. A strategy such as optimum genetic contributions to maximize selection response and manage the long-term genetic variability and inbreeding could lead to sustainable breeding programs for the Mexican Romosinuano cattle breed.

Abstract

The ultimate goal of genetic selection is to improve genetic progress by increasing favorable alleles in the population. However, with selection, homozygosity, and potentially harmful recessive alleles can accumulate, deteriorating genetic variability and hampering continued genetic progress. Such potential adverse side effects of selection are of particular interest in populations with a small effective population size like the Romosinuano beef cattle in Mexico. The objective of this study was to evaluate the genetic background and inbreeding depression in Mexican Romosinuano cattle using pedigree and genomic information. Inbreeding was estimated using pedigree (FPED) and genomic information based on the genomic relationship matrix (FGRM) and runs of homozygosity (FROH) of different length classes. Linkage disequilibrium (LD) was evaluated using the correlation between pairs of loci, and the effective population size (Ne) was calculated based on LD and pedigree information. The pedigree file consisted of 4875 animals born between 1950 and 2019, of which 71 had genotypes. LD decreased with the increase in distance between markers, and Ne estimated using genomic information decreased from 610 to 72 animals (from 109 to 1 generation ago), the Ne estimated using pedigree information was 86.44. The reduction in effective population size implies the existence of genetic bottlenecks and the decline of genetic diversity due to the intensive use of few individuals as parents of the next generations. The number of runs of homozygosity per animal ranged between 18 and 102 segments with an average of 55. The shortest and longest segments were 1.0 and 36.0 Mb long, respectively, reflecting ancient and recent inbreeding. The average inbreeding was 2.98 ± 2.81, 2.98 ± 4.01, and 7.28 ± 3.68% for FPED, FGRM, and FROH, respectively. The correlation between FPED and FGRM was −0.25, and the correlations among FPED and FROH of different length classes were low (from 0.16 to 0.31). The correlations between FGRM and FROH of different length classes were moderate (from 0.44 to 0.58), indicating better agreement. A 1% increase in population inbreeding decreased birth weight by 0.103 kg and weaning weight by 0.685 kg. A strategy such as optimum genetic contributions to maximize selection response and manage the long-term genetic variability and inbreeding could lead to more sustainable breeding programs for the Mexican Romosinuano beef cattle breed.

Litter Size of Sheep (Ovis aries): Inbreeding Depression and Homozygous Regions

Ovine litter size (LS) is an important trait showing variability within breeds. It remains largely unknown whether inbreeding depression on LS exists based on genomic homozygous regions, and whether the homozygous regions resulted from inbreeding are significantly associated with LS in sheep. We here reanalyze a set of single nucleotide polymorphism (SNP) chip of six breeds to characterize the patterns of runs of homozygosity (ROH), to evaluate inbreeding levels and inbreeding depressions on LS, and to identify candidate homozygous regions responsible for LS. Consequently, unique ROH patterns were observed among six sheep populations. Inbreeding depression on LS was only found in Hu sheep, where a significant reduction of 0.016, 0.02, and 0.02 per 1% elevated inbreeding FROH4-8, FROH>8 and the total inbreeding measure was observed, respectively. Nine significantly homozygous regions were found for LS in Hu sheep, where some promising genes for LS possibly via regulation of the development of oocytes (NGF, AKT1, and SYCP1), fertilization (SPAG17, MORC1, TDRD9, ZFYVE21, ADGRB3, and CKB), embryo implantation (PPP1R13B, INF2, and VANGL1) and development (DPPA2, DPPA4, CDCA4, CSDE1, and ADSSL1), and reproductive health (NRG3, BAG5, CKB, and XRCC3) were identified. These results from the present study would provide insights into the genetic management and complementary understandings of LS in sheep.