Runs of homozygosity analysis and genomic inbreeding estimation in Sumba Ongole cattle (Bos indicus) using a BovineSNP50K BeadChip.

Background and Aim: Runs of homozygosity (ROH) is a biocomputational technique for identifying homozygous regions in the genomics of livestock. This study aimed to determine the ROH in Sumba Ongole (SO) bulls (n = 48) using the BovineSNP50K BeadChip. Materials and Methods: GenomeStudio 2.0 software was used to generate the BovineSNP50K BeadChip output. The ROH and ROH-based inbreeding coefficients (FROH) were determined using the detect RUNS R v4.1.0 package. Using the following filtering criteria, PLINK v1.90 software was used to perform genotype quality control: (1) Individuals and single-nucleotide polymorphism (SNPs) had call rates >0.95; (2) more than 0.05 was the minor allele frequency; (3) the list contained only SNPs linked to autosomes; and (4) SNPs that strongly deviated (p < 1e-6) from Hardy-Weinberg equilibrium were removed. Subsequently, 25,252 autosomal SNP markers were included in the ROH and FROH analyses. Results: In general, the number and length of ROH segments in pool animals were 149.77 ± 16.02 Mb and 486.13 ± 156.11 Mb, respectively. Furthermore, the ROH segments in the animals under study can be discriminated into two classes of 1-4 Mb (83.33%) and 4-8 Mb (16.67%). Subsequently, Bos taurus autosomes (BTA) 1, BTA6, and BTA14 had significant homozygous segments comprising 13 genes. Despite this, the average FROH in pool animals was 0.20 ± 0.06. Conclusion: These findings indicate that a recent inbreeding event in SO cattle occurred many generations ago. Furthermore, the candidate genes identified from the ROH analysis indicate phenotypic attributes associated with environmental adaptation and economic traits.

Genome-Wide Analysis of Milk Production Traits and Selection Signatures in Serbian Holstein-Friesian Cattle.

To improve the genomic evaluation of milk-related traits in Holstein-Friesian (HF) cattle it is essential to identify the associated candidate genes. Novel SNP-based analyses, such as the genetic mapping of inherited diseases, GWAS, and genomic selection, have led to a new era of research. The aim of this study was to analyze the association of each individual SNP in Serbian HF cattle with milk production traits and inbreeding levels. The SNP 60 K chip Axiom Bovine BovMDv3 was deployed for the genotyping of 334 HF cows. The obtained genomic results, together with the collected phenotypic data, were used for a GWAS. Moreover, the identification of ROH segments was performed and served for inbreeding coefficient evaluation and ROH island detection. Using a GWAS, a polymorphism, rs110619097 (located in the intron of the CTNNA3 gene), was detected to be significantly (p < 0.01) associated with the milk protein concentration in the first lactation (adjusted to 305 days). The average genomic inbreeding value (FROH) was 0.079. ROH islands were discovered in proximity to genes associated with milk production traits and genomic regions under selection pressure for other economically important traits of dairy cattle. The findings of this pilot study provide useful information for a better understanding of the genetic architecture of milk production traits in Serbian HF dairy cows and can be used to improve lactation performances in Serbian HF cattle breeding programs.

Genome-wide runs of homozygosity signatures in diverse Indian goat breeds.

The present study analyzed ROH and consensus ROH regions in 102 animals of eleven diverse Indian goat (Capra hircus) breeds using whole genome sequencing. A total of 51,705 ROH and 21,271 consensus regions were identified. The mean number of ROH per animal was highest in the meat breed, Jharkhand Black (2693) and lowest in the pashmina breed, Changthangi (60). The average length of ROH (ALROH) was maximum in Kanniadu (974.11 Kb) and minimum in Tellicherry (146.98 Kb). Long ROH is typically associated with more recent inbreeding, whereas short ROH is connected to more ancient inbreeding. The overall ROH-based genomic inbreeding (FROH) was highest for Jharkhand Black (0.602) followed by Kanniadu (0.120) and Sangamneri (0.108) among all breeds. FROH of Jharkhand Black was higher than Kanniadu  up to 5 Mb ROH length category. However, in > 20 Mb ROH length category, Kanniadu (0.98)  exhibited significantly higher FROH than Jharkhand Black (0.46). This  implies that Kanniadu had higher levels of recent inbreeding than Jharkhand Black. Despite this, due to the presence of both recent and ancient inbreeding, Jharkhand Black  demonstrated higher overall FROH compared to Kanniadu. ROH patterns revealed dual purpose (meat and dairy) and pashmina breeds as less consanguineous while recent inbreeding was apparent in meat breeds. Analysis of ROH consensus regions identified selection sweeps in key genes governing intramuscular fat deposition, meat tenderisation, lean meat production and carcass weight (CDK4, ALOX15, CASP9, PRDM16, DVL1) in meat breeds; milk fat percentage and mammary gland development (POLD1, NOTCH2, ARHGAP35) in dual purpose (meat and dairy) breeds; while cold adaptation and hair follicle development (APOBEC1, DNAJC3, F2RL1, FGF9) in pashmina breed. MAPK, RAS, BMP and Wnt signaling pathways associated with hair follicle morphogenesis in Changthangi were also identified. PCA analysis based on ROH consensus regions revealed that meat breeds are more diverse than other goat breeds/populations. Supplementary Information: The online version contains supplementary material available at 10.1007/s13205-024-03921-y.

The role of inbreeding depression on productive performance in the Italian Holstein breed.

Inbreeding depression has become an urgent issue in cosmopolitan breeds where the massive genetic progress achieved in the latest generations is counterbalanced by a dramatic loss of genetic diversity causing increased health issues. Thus, the aim of this study was to estimate inbreeding depression on productive traits in Holstein dairy cattle. More precisely, we aimed to i) determine the level of inbreeding in 27,735 Italian Holstein dairy cows using pedigree and genotype data, ii) quantify the effect of inbreeding on 305-d in milk yield (MY; kg), fat yield (FY; kg), and protein yield (PY; kg) based on different statistical approaches, iii) determine if recent inbreeding has a more harmful impact than ancestral ones, and iv) quantify chromosomal homozygosity effect on productive traits. Quality control was performed on the autosomal chromosomes resulting in a final dataset of 84,443 single nucleotide polymorphisms. Four statistical models were used to evaluate the presence of inbreeding depression, which included linear regression analysis and division of FPED and FROH into percentile classes. Moreover, FROH was partitioned into i) length classes to assess the role of recent and ancestral inbreeding and ii) chromosome-specific contributions (FROH-CHR). Results evidenced that inbreeding negatively impacted the productive performance of Italian Holstein Friesian cows. However, differences between the estimated FPED and FROH coefficients resulted in different estimates of inbreeding depression. For instance, a 1% increase in FPED and FROH was associated with a decrease in MY of about 44 and 61 kg (P < 0.01). Further, when considering the extreme inbreeding percentile classes moving from the 5th lowest to the 95th highest, there was a reduction of -263 kg and -561 kg per lactation for FPED and FROH. Increased inbreeding, estimated by FPED and FROH, had also a negative effect on PY and FY, either fit as a regressor or percentile classes. When evaluating the impact of inbreeding based on runs of homozygosity (ROH) length classes, longer ROH (over 8 Mb) had a negative effect in all traits, indicating that recent inbreeding might be more harmful than the ancestral one. Finally, results within chromosome homozygosity highlighted specific chromosomes with a more deleterious effect on productive traits.

Inbreeding depression is a reduction in performance or health due to the mating of closely related individuals. The overall aim of this study was to investigate the level of inbreeding in the Italian Holstein dairy cow breed and quantify its negative effect on productive performances. The level of inbreeding was estimated by pedigree (FPED) and genomic data by looking at stretches of homozygosity (FROH). Both methods revealed a reduction in milk yield, fat yield, and protein yield when inbreeding increased. Moreover, the study demonstrated that FROH was able to capture more inbreeding depression compared to FPED. In addition, the more recent inbreeding had a stronger negative impact on productive performances compared to ancestral ones. Then, since the amount of runs of homozygosity can vary across the chromosomes of an individual, the effect of each chromosomal homozygosity region on productive traits was also evaluated. The chromosome-level results might be included in breeding programs to limit the accumulation of homozygosity in particular regions that appear to have a more detrimental effect on productive traits. Overall, this study highlights the importance of avoiding inbreeding in animal breeding programs to keep productive animals in the long term.

Genomic selection analysis of morphological and adaptation traits in Chinese indigenous dog breeds.

The significant morphological differences and abundant germplasm resources of Chinese indigenous dog breeds can be attributed to the diverse geographical environment, including plateaus, mountains, and a long history of raising dogs. The combination of both natural and artificial selection during the past several thousand years has led to hundreds of dog breeds with distinct morphological traits and environmental adaptations. China is one of the earliest countries to domesticate dogs and there are more than 50 ancient indigenous dog breeds. In this study, the run of homozygosity (ROH) and proportion of the autosomal genome covered by ROHs (FROH) were calculated for 10 dog breeds that are the most representative Chinese indigenous dogs based on 170K SNP microarray. The results of FROH showed that the Chuandong hound dogs (HCSSC) have the highest level of inbreeding among the tested breeds. The inbreeding in HCSSC occurred more recently than the Liangshan dogs (SCLSQ) dogs because of more numbers of long ROHs in HCSSC dogs, and the former also have higher inbreeding degree. In addition, there are significant differences in the inbreeding degree among different subpopulations of the same breed, such as the Thin dogs from Shaanxi and Shandong province. To explore genome-wide selection signatures among different breeds, including coat color, ear shape, and altitude adaptability, we performed genome selection analyses of FST and cross population extended haplotype homozygosity (XP-EHH). For the coat color, the FST analysis between Xiasi dogs (XSGZ) and HCSSC dogs was performed and identified multiple genes involved in coat color, hair follicle, and bone development, including MC1R, KITLG, SOX5, RSPO2, and TBX15. For the plateau adaptability, we performed FST and XP-EHH analyses between dogs from Tibet (Tibetan Mastiffs and Nyingchi dogs) and plain regions (Guangxi Biwei dogs GXBWQ and Guandong Sharpei dogs). The results showed the EPAS1 gene in dogs from Tibet undergo strong selection. Multiple genes identified for selection signals based on different usage of dogs. Furthermore, the results of ear shape analyses showed that MSRB3 was likely to be the main gene causing the drop ear of domestic dogs. Our study provides new insights into further understanding of Chinese indigenous dogs.

Identification of ROH Islands Conserved through Generations in Pigs Belonging to the Nero Lucano Breed.

The recovery of Nero Lucano (NL) pigs in the Basilicata region (Southern Italy) started in 2001 with the collaboration of several public authorities in order to preserve native breeds that can play a significant economic role both due to their remarkable ability to adapt to difficult environments and the value of typical products from their area of origin. In this study, by using the Illumina Porcine SNP60 BeadChip, we compared the genetic structures of NL pigs reared in a single farm in two different periods separated by a time interval corresponding to at least three generations. The results showed an increase in the percentage of polymorphic loci, a decrease in the inbreeding coefficient calculated according to ROH genome coverage (FROH), a reduction in the number of ROH longer than 16 Mb and an increase in ROH with a length between 2 and 4 Mb, highlighting a picture of improved genetic variability. In addition, ROH island analysis in the two groups allowed us to identify five conserved regions, located on chromosomes 1, 4, 8, 14 and 15, containing genes involved in biological processes affecting immune response, reproduction and production traits. Only the conserved ROH island on chromosome 14 contains markers which, according to the literature, are associated with QTLs affecting thoracic vertebra number, teat number, gestation length, age at puberty and mean platelet volume.

Declining autozygosity over time: An exploration in over 1 million individuals from three diverse cohorts.

Previous studies have hypothesized that autozygosity is decreasing over generational time. However, these studies were limited to relatively small samples (n < 11,000) lacking in diversity, which may limit the generalizability of their findings. We present data that partially support this hypothesis from three large cohorts of diverse ancestries, two from the US (All of Us, n = 82,474; the Million Veteran Program, n = 622,497) and one from the UK (UK Biobank, n = 380,899). Our results from a mixed-effect meta-analysis demonstrate an overall trend of decreasing autozygosity over generational time (meta-analyzed slope = -0.029, SE = 0.009, p = 6.03e-4). On the basis of our estimates, we would predict FROH to decline 0.29% for every 20-year increase in birth year. We determined that a model including an ancestry-by-country interaction term fit the data best, indicating that ancestry differences in this trend differ by country. We found further evidence to suggest a difference between the US and UK cohorts by meta-analyzing within country, observing a significant negative estimate in the US cohorts (meta-analyzed slope = -0.058, SE = 0.015, p = 1.50e-4) but a non-significant estimate in the UK (meta-analyzed slope = -0.001, SE = 0.008, p = 0.945). The association between autozygosity and birth year was substantially attenuated when accounting for educational attainment and income (meta-analyzed slope = -0.011, SE = 0.008, p = 0.167), suggesting they may partially account for decreasing autozygosity over time. Overall, we demonstrate decreasing autozygosity over time in a large, modern sample and speculate that this trend can be attributed to increases in urbanization and panmixia and differences in sociodemographic processes lead to country-specific differences in the rate of decline.

Using Pedigree and Genomic Data toward Better Management of Inbreeding in Italian Dairy Sheep and Goat Breeds.

The inbreeding coefficient is an important parameter for livestock management. Small ruminant breeders and associations mainly rely on pedigree information, but genomic tools are gaining relevance, overcoming possible pedigree inconsistencies. This study investigates the relationship between pedigree-based and genomic inbreeding in two goat and four sheep dairy breeds. Pedigree and genomic data (medium-density SNPchip) were obtained for 3107 goats and 2511 sheep. We estimated pedigree depth (number of fully traced generations, FullGen) and inbreeding (FPED), as well as two genomic inbreeding indexes, using runs of the homozygosity (FROH) and genomic relationship matrix (FGRM). The correlation between the inbreeding coefficients was assessed. A linear regression model (LRM) was fitted for estimating FPED from FROH. After quality control on genomic data, we retained 5085 animals. Mean inbreeding values were low, with higher FROH than FPED and FGRM. Breed differences can partially depend on different managements. The correlation between FPED and FROH was the highest and directly related to pedigree depth. The best LRM was chosen for FullGen ≥4 and ≥6 for goats and sheep, respectively; after excluding animals with extreme residuals, a new refined regression equation was calculated. Since massive genotyping is not affordable to small ruminant breeders, it is important to understand the distinction and relationship between differently calculated inbreeding coefficients, also in view of the introduction of genomic enhanced breeding values. Our study highlights the importance of accurate pedigree information and, especially if not obtainable, of calculating genomic-based inbreeding coefficients. A better estimation of animals' relatedness contributes to improve animal breeding and conservation.

Analysis of Autozygosity Using Whole-Genome Sequence Data of Full-Sib Families in Pikeperch (Sander lucioperca).

Pikeperch (Sander lucioperca) has emerged as a high value species to the aquaculture industry. However, its farming techniques are at an early stage and its production is often performed without a selective breeding program, potentially leading to high levels of inbreeding. In this study, we identified and characterized autozygosity based on genome-wide runs of homozygosity (ROH) on a sample of parental and offspring individuals, determined effective population size (N e ), and assessed relatedness among parental individuals. A mean of 2,235 ± 526 and 1,841 ± 363 ROH segments per individual, resulting in a mean inbreeding coefficient of 0.33 ± 0.06 and 0.25 ± 0.06 were estimated for the progeny and parents, respectively. N e was about 12 until four generations ago and at most 106 for 63 generations in the past, with varying genetic relatedness amongst the parents. This study shows the importance of genomic information when family relationships are unknown and the need of selective breeding programs for reproductive management decisions in the aquaculture industry.

Fine-scale estimation of inbreeding rates, runs of homozygosity and genome-wide heterozygosity levels in the Mangalarga Marchador horse breed.

With the availability of high-density SNP panels and the establishment of approaches for characterizing homozygosity and heterozygosity sites, it is possible to access fine-scale information regarding genomes, providing more than just comparisons of different inbreeding coefficients. This is the first study that seeks to access such information for the Mangalarga Marchador (MM) horse breed on a genomic scale. To this end, we aimed to assess inbreeding levels using different coefficients, as well as to characterize homozygous and heterozygous runs in the population. Using Axiom ® Equine Genotyping Array-670k SNP (Thermo Fisher), 192 horses were genotyped. Our results showed different estimates: inbreeding from genomic coefficients (FROH ) = 0.16; pedigree-based (FPED ) = 0.008; and a method based on excess homozygosity (FHOM ) = 0.010. The correlations between the inbreeding coefficients were low to moderate, and some comparisons showed negative correlations, being practically null. In total, 85,295 runs of homozygosity (ROH) and 10,016 runs of heterozygosity (ROHet) were characterized for the 31 horse autosomal chromosomes. The class with the highest percentage of ROH was 0-2 Mbps, with 92.78% of the observations. In the ROHet results, only the 0-2 class presented observations, with chromosome 11 highlighted in a region with high genetic variability. Three regions from the ROHet analyses showed genes with known functions: tripartite motif-containing 37 (TRIM37), protein phosphatase, Mg2+ /Mn2+ dependent 1E (PPM1E) and carbonic anhydrase 10 (CA10). Therefore, our findings suggest moderate inbreeding, possibly attributed to breed formation, annulling possible recent inbreeding. Furthermore, regions with high variability in the MM genome were identified (ROHet), associated with the recent selection and important events in the development and performance of MM horses over generations.

Genome-Wide Runs of Homozygosity Revealed Selection Signatures in Bos indicus.

Genome-wide runs of homozygosity (ROH) are suitable for understanding population history, calculating genomic inbreeding, deciphering genetic architecture of complex traits and diseases as well as identifying genes linked with agro-economic traits. Autozygosity and ROH islands, genomic regions with elevated ROH frequencies, were characterized in 112 animals of seven Indian native cattle breeds (B. indicus) using BovineHD BeadChip. In total, 4138 ROH were detected. The average number of ROH per animal was maximum in draft breed, Kangayam (63.62 ± 22.71) and minimum in dairy breed, Sahiwal (24.62 ± 11.03). The mean ROH length was maximum in Vechur (6.97 Mb) and minimum in Hariana (4.04 Mb). Kangayam revealed the highest ROH based inbreeding (FROH > 1Mb = 0.113 ± 0.059), whereas Hariana (FROH > 1Mb = 0.042 ± 0.031) and Sahiwal (FROH > 1Mb = 0.043 ± 0.048) showed the lowest. The high standard deviation observed in each breed highlights a considerable variability in autozygosity. Out of the total autozygous segments observed in each breed except Vechur, > 80% were of short length (< 8 Mb) and contributed almost 50% of the genome proportion under ROH. However, in Vechur cattle, long ROH contributed 75% of the genome proportion under ROH. ROH patterns revealed Hariana and Sahiwal breeds as less consanguineous, while recent inbreeding was apparent in Vechur. Maximum autozygosity observed in Kangayam is attributable to both recent and ancient inbreeding. The ROH islands were harbouring higher proportion of QTLs for production traits (20.68% vs. 14.64%; P≤ 0.05) but lower for reproductive traits (11.49% vs. 15.76%; P≤ 0.05) in dairy breeds compared to draft breed. In draft cattle, genes associated with resistant to diseases/higher immunity (LYZL1, SVIL, and GPX4) and stress tolerant (CCT4) were identified in ROH islands; while in dairy breeds, for milk production (PTGFR, CSN1S1, CSN2, CSN1S2, and CSN3). Significant difference in ROH islands among large and short statured breeds was observed at chromosome 3 and 5 involving genes like PTGFR and HMGA2 responsible for milk production and stature, respectively. PCA analysis on consensus ROH regions revealed distinct clustering of dairy, draft and short stature cattle breeds.