Genetic diversity and population history of eight Italian beef cattle breeds using measures of autozygosity

Uncovering genomic diversity and signatures of selection in red Angus × Chinese red steppe crossbred cattle population

Crossbreeding is a cornerstone of modern livestock improvement, combining desirable traits to enhance productivity and environmental resilience. This study conducts the first comprehensive genomic analysis of Red Angus × Chinese Red Steppe (RACS) crossbred cattle, evaluating their genetic architecture, diversity, and selection signatures relative to founder breeds (Red Angus and Chinese Red Steppe) and global populations. A total of 119 cattle, comprising 104 RACS crossbreds and 15 Chinese Red Steppes cattle, were genotyped using the GGP Bovine 100k SNP array. Additionally, the public available genotypic data generated using the BovineSNP50 chip from 550 animals across eight beef breeds (Angus, Hereford, Limousin, Charolais, Mongolian, Shorthorn, Red Angus, and Simmental) and one dairy breed (Holstein) were incorporated into the analysis. We aimed to (1) define the population structure of RACS cattle, (2) quantify their genomic diversity and inbreeding levels, and (3) pinpoint regions under selection linked to adaptive and economic traits. We employed runs of homozygosity (ROH) and population differentiation (Fst) analyses to detect selection signals. The results revealed that the crossbred (RACS), Angus, and Red Angus breeds exhibited similar clustering patterns in principal component analysis (PCA), but the crossbred population showed the highest nucleotide diversity and lowest inbreeding coefficients compared to other breeds. Notably, candidate regions associated with immune response, cold adaptation, and carcass traits were identified within the RACS population. These findings enhance our understanding of the genetic makeup of crossbred beef cattle and highlight their potential for genetic improvement, informing future selection and breeding strategies aimed at optimizing beef production in challenging environments.

Investigation of selection signatures of dairy goats using whole-genome sequencing data

Dairy goats, a livestock species with a long history of milk production, are essential for the economic advancement of nations, particularly in regions experiencing growth. In this study, we gathered whole-genome resequencing data of 58 goats, including 34 dairy goats and 24 wild goats (Bezoar), to explore the selection signatures linked to milk production traits using ROH (Runs of homozygosity), CLR (composite likelihood ratio), Fst (Fixation index), XP-EHH (Ex-tended haplotype homozygosity across populations) and XP-CLR(Cross-population composite likelihood ratio test) methods. Analysis of five tests of selection signatures for dairy goats revealed a total of 210 genes, with 24 genes consistently identified in at least two approaches. These genes are associated with milk fat, milk protein, and fat yield. Gene enrichment analysis highlighted important GO and KEGG pathways related to milk production, such as the "acyl-CoA metabolic process", "glycerolipid biosynthetic process", "cellular response to fatty ac-id", "hormone metabolic process", "Galactose metabolism". Additionally, genes linked to repro-duction, immune response, and environmental adaptation were identified in dairy goats. The findings from our study offer profound understanding into the critical economic features of dairy goats and offer practical guidance for the improvement and development of crossbreeding initiatives across different dairy goat breeds.

Examination of homozygosity runs and selection signatures in native goat breeds of Henan, China

Understanding the genomic characteristics of livestock is crucial for improving breeding efficiency and conservation efforts. However, there is a relative lack of information on the genetic makeup of local goat breeds in Henan, China. In this study, we identified runs of homozygosity (ROH), genomic inbreeding coefficients (FROH), and selection signatures in four breeds including Funiu White (FNW), Huai (HG), Lushan Bullleg (LS), and Taihang black (THB). The genomic analysis utilized a dataset of 46,278 SNP markers and 102 animals. A total of 342, 567, 1285, and 180 ROH segments were detected in FNW, HG, LS, and THB, respectively, with an average of 15.55, 29.84, 32.95, and 8.18 segments per individual. The lengths of ROH segments varied from 69.36 Mb in THB to 417.06 Mb in LS, with the most common lengths being 2-4 Mb and 4-8 Mb. The highest number of longest ROH segments (> 16 Mb) were found in LS (328) and the highest average FROH value was observed in LS (0.173), followed by HG (0.128), while the lowest FROH values were in THB (0.029) and FNW (0.070). Furthermore, the analysis of ROH islands and Composite Likelihood Ratio (CLR) identified a total of 175 significant genes. Among these, 25 genes were found to overlap, detected by both methods. These genes were associated with a diverse range of traits including reproductive ability (GPRIN3), weight (CCSER1), immune response (HERC5 and TIGD2), embryo development (NAP1L5), environmental adaptation (KLHL3, TRHDE, and IFNGR1), and milk characteristics (FAM13A). Significant Gene Ontology (GO) terms related to embryo skeletal system morphogenesis, brain ventricle development, and growth were also identified. This study helps reveal the genetic architecture of Henan goat breeds and provides valuable insights for the effective conservation and breeding programs of local goat breeds in Henan.

Genomic prediction and genome-wide association study using combined genotypic data from different genotyping systems: application to apple fruit quality traits

With advances in next-generation sequencing technologies, various marker genotyping systems have been developed for genomics-based approaches such as genomic selection (GS) and genome-wide association study (GWAS). As new genotyping platforms are developed, data from different genotyping platforms must be combined. However, the potential use of combined data for GS and GWAS has not yet been clarified. In this study, the accuracy of genomic prediction (GP) and the detection power of GWAS increased for most fruit quality traits of apples when using combined data from different genotyping systems, Illumina Infinium single-nucleotide polymorphism array and genotyping by random amplicon sequencing-direct (GRAS-Di) systems. In addition, the GP model, which considered the inbreeding effect, further improved the accuracy of the seven fruit traits. Runs of homozygosity (ROH) islands overlapped with the significantly associated regions detected by the GWAS for several fruit traits. Breeders may have exploited these regions to select promising apples by breeders, increasing homozygosity. These results suggest that combining genotypic data from different genotyping platforms benefits the GS and GWAS of fruit quality traits in apples. Information on inbreeding could be beneficial for improving the accuracy of GS for fruit traits of apples; however, further analysis is required to elucidate the relationship between the fruit traits and inbreeding depression (e.g. decreased vigor).

Investigation of cosmopolitan and local Italian beef cattle breeds uncover common patterns of heterozygosity

The analysis of livestock heterozygosity is less common compared to the study of homozygous patterns. Heterozygous-Rich Regions (HRRs) may harbor significant loci for functional traits such as immune response, survival rate, and fertility. For this reason, this study was conducted to investigate and characterize the heterozygosity patterns of four beef cattle breeds, which included two cosmopolitan breeds (Limousine and Charolaise) and two local breeds (Sarda and Sardo Bruna). Our analysis identified regions with a high degree of heterozygosity using a consecutive runs approach, the Tajima D test, nucleotide diversity estimation, and Hardy Weinberg equilibrium test. These regions exhibited recurrent heterozygosity peaks and were consistently found on specific chromosomes across all breeds, specifically autosomes 15, 16, 20, and 23. The cosmopolitan and Sardo Bruna breeds also displayed peaks on autosomes 2 and 21, respectively. Thirty-five top runs shared by more than 25% of the populations were identified. These genomic fragments encompassed 18 genes, two of which are directly linked to male fertility, while four are associated with lactation. Two other genes play roles in survival and immune response. Our study also detected a region related to growth and carcass traits in Limousine breed. Our analysis of heterozygosity-rich regions revealed particular segments of the cattle genome linked to various functional traits. It appears that balancing selection is occurring in specific regions within the four examined breeds, and unexpectedly, they are common across cosmopolitan and local breeds. The genes identified hold potential for applications in breeding programs and conservation studies to investigate the phenotypes associated with these heterozygous genotypes. In addition, Tajima D test, Nucleotide diversity, and Hardy Weinberg equilibrium test confirmed the presence of heterozygous fragments found with Heterozygous-Rich Regions analysis.

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.

Genome-wide landscape of runs of homozygosity and differentiation across Egyptian goat breeds

Understanding the genomic features of livestock is essential for successful breeding programs and conservation. This information is scarce for local goat breeds in Egypt. In the current study, genomic regions with selection signatures were identified as well as runs of homozygosity (ROH), genomic inbreeding coefficients (FROH) and fixation index (FST) were detected in Egyptian Nubian, Damascus, Barki and Boer goat breeds. A total of 46,268 SNP markers and 337 animals were available for the genomic analyses. On average, 145.44, 42.02, 87.90 and 126.95 ROHs were detected per individual in the autosomal genome of the respective breeds. The mean accumulative ROH lengths ranged from 46.5 Mb in Damascus to 360 Mb in Egyptian Nubian. The short ROH segments (< 2 Mb) were most frequent in all breeds, while the longest ROH segments (> 16 Mb) were exclusively found in the Egyptian Nubian. The highest average FROH was observed in Egyptian Nubian (~ 0.12) followed by Boer (~ 0.11), while the lowest FROH was found in Damascus (~ 0.05) and Barki breed (~ 0.03). The estimated mean FST was 0.14 (Egyptian Nubian and Boer), 0.077 (Egyptian Nubian and Barki), 0.075 (Egyptian Nubian and Damascus), 0.071 (Barki and Boer), 0.064 (Damascus and Boer), and 0.015 (Damascus and Barki), for each pair of breeds. Interestingly, multiple SNPs that accounted for high FST values were observed on chromosome 6 in regions harboring ALPK1 and KCNIP4. Genomic regions overlapping both FST and ROH harbor genes related to immunity (IL4R, PHF23, GABARAP, GPS2, and CD68), reproduction (SPATA2L, TNFSF12, TMEM95, and RNF17), embryonic development (TCF25 and SOX15) and adaptation (MC1R, KDR, and KIT), suggesting potential genetic adaptations to local environmental conditions. Our results contribute to the understanding of the genetic architecture of different goat breeds and may provide valuable information for effective preservation and breeding programs of local goat breeds in Egypt.

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

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 (F_ROH) and homozygous loci-based inbreeding estimate (F_HOM).

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 F_ROH, F_HOM, and pedigree-based inbreeding estimate (F_PED) 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.

Assessment of Genetic Diversity, Runs of Homozygosity, and Signatures of Selection in Tropical Milking Criollo Cattle Using Pedigree and Genomic Data

The objective of this study was to evaluate the genetic diversity of the Tropical Milking Criollo cattle (TMC) breed in Mexico through parameters derived from pedigree and genomic information assessment. The pedigree file consisted of 3780 animals. Seventy-nine bovines were genotyped with the medium-density single nucleotide polymorphism chip and considered a reference population for pedigree analysis. The effective population size and the probability of gene origin used to assess the evolution of genetic diversity were calculated with pedigree information. Inbreeding coefficients were evaluated based on pedigree (FPed), the genomic relationship matrix (FGRM), and runs of homozygosity (FROH) of different length classes. The average inbreeding was 2.82 ± 2.66%, −0.7 ± 3.8%, and 10.9 ± 3.0% for FPED, FGRM, and FROH, respectively. Correlation between FPED and FROH was significant only for runs of homozygosity > 4 Mb, indicating the FPED of a population with an average equivalent complete generation of five only recovers the most recent inbreeding. The parameters of the probability of gene origin indicated the existence of genetic bottlenecks and the loss of genetic diversity in the history of the TMC cattle population; however, pedigree and genomic information revealed the existence of current sufficient genetic diversity to design a sustainable breeding program.

Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems

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 (&lt;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.

Heterozygosity and homozygosity regions affect reproductive success and the loss of reproduction: a case study with litter traits in pigs

Runs of heterozygosity (ROHet) and homozygosity (ROH) harbor useful information related to traits of interest. There is a lack of investigating the effect of ROHet and ROH on reproductive success and the loss of reproduction in mammals. Here, we detected and characterized the ROHet and ROH patterns in the genomes of Chinese indigenous pigs (i.e., Jinhua, Chun’an, Longyou Black, and Shengxian Spotted pigs), revealing the similar genetic characteristics of indigenous pigs. Later, we highlighted the underlying litter traits-related ROHet and ROH using association analysis with linear model in these four indigenous pig breeds. To pinpoint the promising candidate genes associated with litter traits, we further in-depth explore the selection patterns of other five pig breeds (i.e., Erhualian, Meishan, Minzhu, Rongchang, and Diqing pigs) with different levels of reproduction performance at the underlying litter traits-related ROHet and ROH using F_ST and genetic diversity ratio. Then, we identified a set of known and novel candidate genes associated with reproductive performance in pigs. For the novel candidate genes (i.e., CCDC91, SASH1, SAMD5, MACF1, MFSD2A, EPC2, and MBD5), we obtained public available datasets and performed multi-omics analyses integrating transcriptome-wide association studies and comparative single-cell RNA-seq analyses to uncover the roles of them in mammalian reproductive performance. The genes have not been widely reported to be fertility-related genes and can be complementally considered as prior biological information to modify genomic selections models that benefits pig genetic improvement of litter traits. Besides, our findings provide new insights into the function of ROHet and ROH in mammals.

Genome-Wide Estimates of Runs of Homozygosity, Heterozygosity, and Genetic Load in Two Chinese Indigenous Goat Breeds

Runs of homozygosity (ROH) and heterozygosity (ROHet) are windows into population demographic history and adaptive evolution. Numerous studies have shown that deleterious mutations are enriched in the ROH of humans, pigs, cattle, and chickens. However, the relationship of deleterious variants to ROH and the pattern of ROHet in goats have been largely understudied. Here, 240 Guangfeng and Ganxi goats from Jiangxi Province, China, were genotyped using the Illumina GoatSNP50 BeadChip and genome-wide ROH, ROHet, and genetic load analyses were performed in the context of 32 global goat breeds. The classes with the highest percentage of ROH and ROHet were 0.5–2 Mb and 0.5–1 Mb, respectively. The results of inbreeding coefficients (based on SNP and ROH) and ROHet measurements showed that Guangfeng goats had higher genetic variability than most Chinese goats, while Ganxi goats had a high degree of inbreeding, even exceeding that of commercial goat breeds. Next, the predicted damaging homozygotes were more enriched in long ROHs, especially in Guangfeng goats. Therefore, we suggest that information on damaging alleles should also be incorporated into the design of breeding and conservation programs. A list of genes related to fecundity, growth, and environmental adaptation were identified in the ROH hotspots of two Jiangxi goats. A sense-related ROH hotspot (chromosome 12: 50.55–50.81 Mb) was shared across global goat breeds and may have undergone selection prior to goat domestication. Furthermore, an identical ROHet hotspot (chromosome 1: 132.21–132.54 Mb) containing two genes associated with embryonic development (STAG1 and PCCB) was detected in domestic goat breeds worldwide. Tajima’s D and BetaScan2 statistics indicated that this region may be caused by long-term balancing selection. These findings not only provide guidance for the design of conservation strategies for Jiangxi goat breeds but also enrich our understanding of the adaptive evolution of goats.