Runs of homozygosity and population history in cattle

Insights from homozygous signatures of cervus nippon revealed genetic architecture for components of fitness

This study investigates the genomic landscape of Sika deer populations, emphasizing the detection and characterization of runs of homozygosity (ROH) and their contribution towards components of fitness. Using 85,001 high-confidence SNPs, the investigation into ROH distribution unveiled nuanced patterns of autozygosity across individuals especially in 2 out of the 8 farms, exhibiting elevated ROH levels and mean genome coverage under ROH segments. The prevalence of shorter ROH segments (0.5-4 Mb) suggests historical relatedness and potential selective pressures within these populations. Intriguingly, despite observed variations in ROH profiles, the overall genomic inbreeding coefficient (FROH) remained relatively low across all farms, indicating a discernible degree of genetic exchange and effective mitigation of inbreeding within the studied Sika deer populations. Consensus ROH (cROH) were found to harbor genes for important functions viz., EGFLAM gene which is involved in the vision function of the eye, SKP2 gene which regulates cell cycle, CAPSL involved in adipogenesis, SPEF2 which is essential for sperm flagellar assembly, DCLK3 involved in the heat stress. This first ever study on ROH in Sika deer, to shed light on the adaptive role of genes in these homozygous regions. The insights garnered from this study have broader implications in the management of genetic diversity in this vulnerable species.

Genomic inbreeding analysis reveals resilience and genetic diversity in Indian yak populations

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.

Analysis of genetic diversity and structure of endangered Dengchuan cattle population using a single-nucleotide polymorphism chip

This study aimed to evaluate the genetic diversity and structure within the Dengchuan cattle population and effectively protect and utilize their germplasm resources. Herein, the single-nucleotide polymorphisms (SNPs) of 100 Dengchuan cattle (46 bulls and 54 cows) were determined using the GGP Bovine 100K SNP Beadchip. The results showed that among the Dengchuan cattle, a total of 101,220 SNPs were detected, and there were 83,534 SNPs that passed quality control, of which 85.7% were polymorphic. The average genetic distance based on identity-by-state (IBS) within the conservation population of Dengchuan cattle was 0.26 ± 0.02. A total of 3,999 genome-length runs of homozygosity (ROHs) were detected in the Dengchuan cattle, with ROH lengths primarily concentrated in the range of 1-5 Mb, accounting for 87.02% of the total. The average inbreeding coefficient based on ROHs was 4.6%, within the conservation population of Dengchuan cattle, whereas it was 4.9% for bulls, and the Wright inbreeding coefficient (FIS) value was 2.4%, demonstrating a low level of inbreeding within the Dengchuan cattle population. Based on neighbor-joining tree analysis, the Dengchuan cattle could be divided into 16 families. In summary, the conservation population of Dengchuan cattle displays relatively abundant diversity and a moderate genetic relationship. Inbreeding was observed among a few individuals, but the overall inbreeding level of the population remained low. It is important to maintain this low level of inbreeding when introducing purebred bloodlines to expand the core group. This approach will ensure the long-term conservation of Dengchuan cattle germplasm resources and prevent loss of genetic diversity.

Assessment of inbreeding coefficients and inbreeding depression on complex traits from genomic and pedigree data in Nelore cattle

BACKGROUND

Nelore cattle play a key role in tropical production systems due to their resilience to harsh conditions, such as heat stress and seasonally poor nutrition. Monitoring their genetic diversity is essential to manage the negative impacts of inbreeding. Traditionally, inbreeding and inbreeding depression are assessed by pedigree-based coefficients (F), but recently, genetic markers have been preferred for their precision in capturing the inbreeding level and identifying animals at risk of reduced productive and reproductive performance. Hence, we compared the inbreeding and inbreeding depression for productive and reproductive performance traits in Nelore cattle using different inbreeding coefficient estimation methods from pedigree information (FPed), the genomic relationship matrix (FGRM), runs of homozygosity (FROH) of different lengths (> 1 Mb (genome), between 1 and 2 Mb - FROH 1-2; 2-4 Mb FROH 2-4 or > 8 Mb FROH >8) and excess homozygosity (FSNP).

RESULTS

The correlation between FPed and FROH was lower when the latter was based on shorter segments (r = 0.15 with FROH 1-2, r = 0.20 with FROH 2-4 and r = 0.28 with FROH 4-8). Meanwhile, the FPed had a moderate correlation with FSNP (r = 0.47) and high correlation with FROH >8 (r = 0.58) and FROH-genome (r = 0.60). The FROH-genome was highly correlated with inbreeding based on FROH>8 (r = 0.93) and FSNP (r = 0.88). The FGRM exhibited a high correlation with FROH-genome (r = 0.55) and FROH >8 (r = 0.51) and a lower correlation with other inbreeding estimators varying from 0.30 for FROH 2-4 to 0.37 for FROH 1-2. Increased levels of inbreeding had a negative impact on the productive and reproductive performance of Nelore cattle. The unfavorable inbreeding effect on productive and reproductive traits ranged from 0.12 to 0.51 for FPed, 0.19-0.59 for FGRM, 0.21-0.58 for FROH-genome, and 0.19-0.54 for FSNP per 1% of inbreeding scaled on the percentage of the mean. When scaling the linear regression coefficients on the standard deviation, the unfavorable inbreeding effect varied from 0.43 to 1.56% for FPed, 0.49-1.97% for FGRM, 0.34-2.2% for FROH-genome, and 0.50-1.62% for FSNP per 1% of inbreeding. The impact of the homozygous segments on reproductive and performance traits varied based on the chromosomes. This shows that specific homozygous chromosome segments can be signs of positive selection due to their beneficial effects on the traits.

CONCLUSIONS

The low correlation observed between FPed and genomic-based inbreeding estimates suggests that the presence of animals with one unknown parent (sire or dam) in the pedigree does not account for ancient inbreeding. The ROH hotspots surround genes related to reproduction, growth, meat quality, and adaptation to environmental stress. Inbreeding depression has adverse effects on productive and reproductive traits in Nelore cattle, particularly on age at puberty in young bulls and heifer calving at 30 months, as well as on scrotal circumference and body weight when scaled on the standard deviation of the trait.

[Runs of Homozygosity Decipher Genetic Diversity in Cattle Breed Dwelling in the Colder Regions of the World]

BACKGROUND

Our study focuses on Yakutian cattle, a Siberian native breed, examining its inbreeding and diversity through genome-wide analysis of runs of homozygosity (ROHs). Yakutian cattle are adapted to Siberia's harsh sub-arctic conditions, enduring temperatures below -70°C. However, the population genetics studies on this breed are scanty, to document the genetic uniqueness in these cattle.

RESULTS

We analyzed 40 Yakutian cattle with strict quality control for ROH detection yielding 683 homozygous segments, averaging 17 per individual with an average length of 9 Mb. ROH regions were found to be involved in important pathways pertaining to cold adaptation. Autozygosity ranged from 1% to 12% of the genome, with a relatively low average inbreeding coefficient (FROH) of 0.057, as compared to other breeds. Also, the different diversity indicators, namely, principal component analysis, heterozygosity, and effective population size analysis, revealed the prevalence of genetic diversity within the breed.

CONCLUSION

Our findings on ROH are the first of its kind in Yakutian cattle that support their adaptability to colder environments, as evidenced by low inbreeding and high genetic diversity.

Integration of ssGWAS and ROH analyses for uncovering genetic variants associated with reproduction traits in Large White pigs

The low heritability of reproduction traits such as total number born (TNB), number born alive (NBA) and adjusted litter weight until 21 days at weaning (ALW) poses a challenge for genetic improvement. In this study, we aimed to identify genetic variants that influence these traits and evaluate the accuracy of genomic selection (GS) using these variants as genomic features. We performed single-step genome-wide association studies (ssGWAS) on 17 823 Large White (LW) pigs, of which 2770 were genotyped by 50K single nucleotide polymorphism (SNP) chips. Additionally, we analyzed runs of homozygosity (ROH) in the population and tested their effects on the traits. The genomic feature best linear unbiased prediction (GFBLUP) was then carried out in an independent population of 350 LW pigs using identified trait-related SNP subsets as genomic features. As a result, our findings identified five, one and four SNP windows that explaining more than 1% of genetic variance for ALW, TNB, and NBA, respectively and discovered 358 hotspots and nine ROH islands. The ROH SSC1:21814570-27186456 and SSC11:7220366-14276394 were found to be significantly associated with ALW and NBA, respectively. We assessed the genomic estimated breeding value accuracy through 20 replicates of five-fold cross-validation. Our findings demonstrate that GFBLUP, incorporating SNPs located in effective ROH (p-value < 0.05) as genomic features, might enhance GS accuracy for ALW compared with GBLUP. Additionally, using SNPs explaining more than 0.1% of the genetic variance in ssGWAS for NBA as genomic features might improve the GS accuracy, too. However, it is important to note that the incorporation of inappropriate genomic features can significantly reduce GS accuracy. In conclusion, our findings provide valuable insights into the genetic mechanisms of reproductive traits in pigs and suggest that the ssGWAS and ROH have the potential to enhance the accuracy of GS for reproductive traits in LW pigs.

Genomic characterisation and diversity assessment of eight endangered Belgian sheep breeds

Assessing the genetic diversity of local breeds is essential for conserving these unique breeds, which may possess unique traits. This study provides the genomic characterisation of eight indigenous sheep breeds in Belgium based on pedigree and single nucleotide polymorphism (SNP) analysis. A total of 687 sheep were genotyped and were subjected to a rigorous quality control, resulting in a set of 45 978 autosomal SNPs. Pedigree analysis showed breed-average inbreeding estimates between 3.3% and 11.3%. The genomic analysis included an assessment of runs of homozygosity (ROH) to examine the genomic inbreeding coefficient, with breed-average inbreeding coefficients estimated between 4.1% and 8.5%. Runs of homozygosity islands were identified in six of the eight breeds studied, with some exhibiting an incidence of up to 58%. Interestingly, several ROH islands overlapped with other breeds included in this study, as well as with international sheep breeds. Pedigree-based effective population sizes were estimated below 100 for all breeds, whereas genomic-based effective population sizes were below 24, indicating that these eight local sheep breeds are endangered. Principal component analysis, admixture analyses, and Fst computations were used to study the population structure and genetic differences. A neighbour-joining tree using 95 international sheep breeds positioned the eight local breeds in the group of milksheep, Texel sheep and the Scandinavian breeds. Additionally, the investigation of paternal oY1 genotypes revealed diverse lineage origins within the Belgian sheep population. This study refines and deepens our knowledge about the local sheep breeds in Belgium, thereby improving their management and conservation. Moreover, as these breeds are linked to other international breeds, these insights are significant for the global scientific community.

Runs of homozygosity analysis for selection signatures in the Yellow Korean native chicken

OBJECTIVE

Yellow Korean native chicken (KNC-Y) is one of the five pure Korean indigenous chicken breeds that were restored through a government project in 1992. KNC-Y is recognized for its superior egg production performance compared to other KNC lines. In this study, we performed runs of homozygosity (ROH) analysis to discover selection signatures associated with egg production traits in the KNC-Y population.

METHODS

A total of 675 DNA samples from KNC-Y were genotyped to generate single nucleotide polymorphism (SNP) data using custom 60K Affymetrix SNP chips. ROH analysis was performed using PLINK software, with predefined parameters set for the analysis. The threshold of ROH island was defined as the top 1% frequency of SNPs withing the ROH among the population.

RESULTS

In the KNC-Y population, a total of 29,958 runs of homozygosity (ROH) fragments were identified. The average total length of ROH was 120.84 Mb, with each ROH fragment having an average length of 2.71 Mb. The calculated ROH-based inbreeding coefficient (FROH) was 0.13. Furthermore, we revealed the presence of ROH islands on chromosomes 1, 2, 4, 5, 7, 8, and 11. Within the identified regions, a total of 111 genes were annotated, and among them were genes related to economic traits, including PRMT3, ANO5, HDAC4, LSS, PLA2G4A, and PTGS2. Most of the overlapping quantitative trait locus regions with ROH islands were found to be associated with production traits.

CONCLUSION

This study conducted a comprehensive analysis of ROH in the KNC-Y population. Notably, among the findings, the PTGS2 gene is believed to play a crucial role in influencing the laying performance of KNC-Y.

Assessment of genetic diversity, inbreeding and collection completeness of Jersey bulls in the US National Animal Germplasm Program

Genomic selection and extensive use of a few elite bulls through artificial insemination are leading to reduced genetic diversity in Jersey cattle. Conservation of genetic diversity through gene banks can protect a breed's genetic diversity and genetic gain, ensuring continued genetic advancement in the future. The availability of genomic information in the US National Animal Germplasm Program (NAGP) facilitates characterization of Jersey bulls in the germplasm collection. Therefore, in this study, we compared the genetic diversity and inbreeding between Jersey bulls in the NAGP and the national cooperator database (NCD). The NCD is maintained and curated by the Council on Dairy Cattle Breeding (CDCB). We found the genetic diversity to be marginally higher in NAGP (Ho = 0.34 ± 0.17) relative to the NCD population (Ho = 0.33 ± 0.16). The average pedigree and genomic inbreeding (FPED, FGRM, FROH > 2Mb) were similar between the groups, with estimates of 7.6% with FPED, 11.07% with FGRM and 20.13% with FROH > 2Mb. An increasing trend in inbreeding was detected, and a significantly higher level of inbreeding was estimated among the older bulls in the NAGP collection, suggesting an overrepresentation of the genetics from elite bulls. Results from principal component analyses (PCA) provided evidence that the NAGP collection is representative of the genetic variation found in the NCD population and a broad majority of the loci segregating (98.2%) in the NCD population were also segregating in the NAGP. Ward's clustering was used to assess collection completeness of Jerseys in the NAGP by comparison with top 1000 sires of bulls, top 1000 sires of cow, and bulls with high Lifetime Net Merit (NM$). All the clusters were represented in the NAGP suggesting that most of the genetic diversity in the US Jersey population is represented in the NAGP and confirmed the PCA results. The decade of birth was the major driver grouping bulls into clusters, suggesting the importance of selection over time. Selection signature analysis between the historic bulls in the NAGP with the newer bulls, born in the decade after implementation of genomic selection, identified selection for milk production, fat and protein yield, fertility, health, and reproductive traits. Cluster analysis revealed that the NAGP has captured allele frequency changes over time associated with selection, validating the strategy of repeated sampling and suggests that the continuation of a repeated sampling policy is essential for the germplasm collection to maintain its future utility. While NAGP should continue to collect bulls that have large influence on the population due to selection, care should be taken to include the entire breadth of bulls, including low merit bulls.

Genomic and Conventional Inbreeding Coefficient Estimation Using Different Estimator Models in Korean Duroc, Landrace, and Yorkshire Breeds Using 70K Porcine SNP BeadChip

The purpose of this study was to estimate the homozygosity distribution and compute genomic and conventional inbreeding coefficients in three genetically diverse pig breed populations. The genomic and pedigree data of Duroc (1586), Landrace (2256), and Yorkshire (3646) were analyzed. We estimated and compared various genomic and pedigree inbreeding coefficients using different models and approaches. A total of 709,384 ROH segments in Duroc, 816,898 in Landrace, and 1,401,781 in Yorkshire, with average lengths of 53.59 Mb, 56.21 Mb, and 53.46 Mb, respectively, were identified. Relatively, the Yorkshire breed had the shortest ROH segments, whereas the Landrace breed had the longest mean ROH segments. Sus scrofa chromosome 1 (SSC1) had the highest chromosomal coverage by ROH across all breeds. Across breeds, an absolute correlation (1.0) was seen between FROH total and FROH1-2Mb, showing that short ROH were the primary contributors to overall FROH values. The overall association between genomic and conventional inbreeding was weak, with values ranging from 0.058 to 0.140. In contrast, total genomic inbreeding (FROH) and ROH classes showed a strong association, ranging from 0.663 to 1.00, across the genotypes. The results of genomic and conventional inbreeding estimates improve our understanding of the genetic diversity among genotypes.

Analyzing Runs of Homozygosity Reveals Patterns of Selection in German Brown Cattle

An increasing trend in ancestral and classical inbreeding coefficients as well as inbreeding depression for longevity were found in the German Brown population. In addition, the proportion of US Brown Swiss genes is steadily increasing in German Browns. Therefore, the aim of the present study was to analyze the presence and genomic localization of runs of homozygosity (ROH) in order to evaluate their associations with the proportion of US Brown Swiss genes and survival rates of cows to higher lactations. Genotype data were sampled in 2364 German Browns from 258 herds. The final data set included 49,693 autosomal SNPs. We identified on average 35.996 ± 7.498 ROH per individual with a mean length of 8.323 ± 1.181 Mb. The genomic inbreeding coefficient FROH was 0.122 ± 0.032 and it decreased to 0.074, 0.031 and 0.006, when genomic homozygous segments > 8 Mb (FROH>8), >16 Mb (FROH>16) and >32 Mb (FROH>32) were considered. New inbreeding showed the highest correlation with FROH>32, whereas ancestral inbreeding coefficients had the lowest correlations with FROH>32. The correlation between the classical inbreeding coefficient and FROH was 0.572. We found significantly lower FROH, FROH>4, FROH>8 and FIS for US Brown Swiss proportions <60% compared to >80%. Cows surviving to the 2nd, 4th, 6th, 8th, and 10th lactation had lower genomic inbreeding for FROH and up to FROH>32, which was due to a lower number of ROH and a shorter average length of ROH. The strongest ROH island and consensus ROH shared by 50% of the animals was found on BTA 6 at 85-88 Mb. The genes located in this genomic region were associated with longevity (NPFFR2 and ADAMTS3), udder health and morphology (SLC4A4, NPFFR2, GC and RASSF6), milk production, milk protein percentage, coagulation properties of milk and milking speed (CSN3). On BTA 2, a ROH island was detected only in animals with <60% US Brown Swiss genes. Genes within this region are predominantly important for dual-purpose cattle breeds including Original Browns. For cows reaching more than 9 lactations, an exclusive ROH island was identified on BTA 7 with genes assumed to be associated with longevity. The analysis indicated that genomic homozygous regions important for Original Browns are still present and also ROH containing genes affecting longevity may have been identified. The breeding of German Browns should prevent any further increase in genomic inbreeding and run a breeding program with balanced weights on production, robustness and longevity.

Genomic analysis of inbreeding and coancestry in Nordic Jersey and Holstein dairy cattle populations

In recent years, genomic selection (GS) has accelerated genetic gain in dairy cattle breeds worldwide. Despite the evident genetic progress, several dairy populations have also encountered challenges such as heightened inbreeding rates and reduced effective population sizes. The challenge has been to find a balance between achieving substantial genetic gain while managing genetic diversity within the population, thereby mitigating the negative effects of inbreeding depression. This study aims to elucidate the impact of GS on pedigree and genomic rates of inbreeding (ΔF) and coancestry (ΔC) in Nordic Jersey (NJ) and Holstein (NH) cattle populations. Furthermore, key genetic metrics, including the generation interval (L), effective population size (Ne), and future effective population size (FNe) were assessed between 2 time periods, before and after GS, and across distinct animal cohorts in both breeds: females, bulls, and approved semen-producing bulls (AI-sires). Analysis of ΔF and ΔC revealed distinct trends across the studied periods and animal groups. Notably, there was a consistent increase in yearly ΔF for most animal groups in both breeds. An exception was observed in NH AI-sires, which demonstrated a slight decrease in yearly ΔF. Moreover, NJ displayed minimal changes in yearly ΔC between the periods, whereas NH exhibited elevated ΔC values across all animal groups. Particularly striking was the substantial increase in yearly ΔC within the NH female population, surging from 0.02% to 0.39% between the periods. Implementation of GS resulted in a reduction of the generation interval across all animal cohorts in both NJ and NH breeds. However, the extent of reduction was more pronounced in males compared with females. This reduction in generation interval influenced generational changes in ΔF and ΔC. Bulls and AI-sires of both breeds exhibited reduced generational ΔF between periods, in contrast to females that demonstrated an opposing pattern. Between the periods, NJ maintained a relatively stable Ne (29.4 before and 30.3 after GS), whereas NH experienced a notable decline from 54.3 to 42.8. Female groups in both breeds displayed a negative Ne trend, whereas males demonstrated either neutral or positive Ne developments. Regarding FNe, NJ exhibited positive FNe development with an increase from 40.7 to 57.2. The opposite was observed in NH, where FNe decreased from 198.8 to 42.7. In summary, it was evident that the genomic methods could detect differences between the populations and changes in ΔF and ΔC more efficiently than pedigree methods. Implementation of GS yielded positive outcomes within the NJ population regarding the rate of coancestry but the opposite was observed with NH. Moreover, analysis of ΔC data hints at the potential to decrease future ΔF through informed mating strategies. Conversely, NH faces more pressing concerns, even though ΔF remains comparatively modest in contrast to what has been observed in other Holstein populations. These findings underscore the necessity of genomic control of inbreeding and coancestry with strategic changes in the Nordic breeding schemes for dairy to ensure long-term sustainability in the forthcoming years.

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.

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

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.

Detection and evaluation of parameters influencing the identification of heterozygous-enriched regions in Holstein cattle based on SNP chip or whole-genome sequence data

BACKGROUND

A heterozygous-enriched region (HER) is a genomic region with high variability generated by factors such as balancing selection, introgression, and admixture processes. In this study, we evaluated the genomic background of HERs and the impact of different parameters (i.e., minimum number of SNPs in a HER, maximum distance between two consecutive SNPs, minimum length of a HER, maximum number of homozygous allowed in a HER) and scenarios [i.e., different SNP panel densities and whole-genome sequence (WGS)] on the detection of HERs. We also compared HERs characterized in Holstein cattle with those identified in Angus, Jersey, and Norwegian Red cattle using WGS data.

RESULTS

The parameters used for the identification of HERs significantly impact their detection. The maximum distance between two consecutive SNPs did not impact HERs detection as the same average of HERs (269.31 ± 787.00) was observed across scenarios. However, the minimum number of markers, maximum homozygous markers allowed inside a HER, and the minimum length size impacted HERs detection. For the minimum length size, the 10 Kb scenario showed the highest average number of HERs (1,364.69 ± 1,483.64). The number of HERs decreased as the minimum number of markers increased (621.31 ± 1,271.83 to 6.08 ± 21.94), and an opposite pattern was observed for the maximum homozygous markers allowed inside a HER (54.47 ± 195.51 to 494.89 ± 1,169.35). Forty-five HER islands located in 23 chromosomes with high Tajima's D values and differential among the observed and estimated heterozygosity were detected in all evaluated scenarios, indicating their ability to potentially detect regions under balancing selection. In total, 3,440 markers and 28 genes previously related to fertility (e.g., TP63, ZSCAN23, NEK5, ARHGAP44), immunity (e.g., TP63, IGC, ARHGAP44), residual feed intake (e.g., MAYO9A), stress sensitivity (e.g., SERPINA6), and milk fat percentage (e.g., NOL4) were identified. When comparing HER islands among breeds, there were substantial overlaps between Holstein with Angus (95.3%), Jersey (94.3%), and Norwegian Red cattle (97.1%), indicating conserved HER across taurine breeds.

CONCLUSIONS

The detection of HERs varied according to the parameters used, but some HERs were consistently identified across all scenarios. Heterozygous genotypes observed across generations and breeds appear to be conserved in HERs. The results presented could serve as a guide for defining HERs detection parameters and further investigating their biological roles in future studies.

Identification of Heilongjiang crossbred beef cattle pedigrees and reveals functional genes related to economic traits based on whole-genome SNP data

Introduction: To enhance the beef cattle industry, Heilongjiang Province has developed a new Crossbred beef cattle variety through crossbreeding with exotic commercial breeds. This new variety exhibits relatively excellent meat quality, and efficient reproductive performance, catering to market demands. Method: This study employed whole genome resequencing technology to analyze the genetic pedigree and diversity of 19 Heilongjiang Crossbred beef cattle, alongside 59 published genomes from East Asian, Eurasian, and European taurine cattle as controls. In addition, genes related to production traits were also searched by identifying Runs of Homozygosity (ROH) islands and important fragments from ancestors. Results: A total of 14,427,729 biallelic SNPs were discovered, with the majority located in intergenic and intron regions and a small percentage in exon regions, impacting protein function. Population genetic analyses including Principal Component Analysis (PCA), Neighbor-Joining (NJ) tree, and ADMIXTURE identified Angus, Holstein, and Mishima as the main ancestors of Crossbred beef cattle. In genetic diversity analysis, nucleotide diversity, linkage disequilibrium, and inbreeding coefficient analysis reveal that the genetic diversity of Crossbred beef cattle is at a moderate level, and a higher inbreeding coefficient indicates the need for careful breeding management. In addition, some genes related to economic traits are identified through the identification of Runs of Homozygosity (ROH) islands and important fragments from ancestors. Conclusion: This comprehensive genomic characterization supports the targeted improvement of economically important traits in Crossbred beef cattle, facilitating advanced breeding strategies.

Genomic insights into adaptation and inbreeding among Sub-Saharan African cattle from pastoral and agropastoral systems

Background

In Sub-Saharan Africa (SSA), cattle are crucial for socioeconomic stability yet face numerous environmental stressors such as diseases, parasites, and extreme heat within pastoral and agropastoral systems. Despite their significance, gaps remain in understanding how genetic diversity and inbreeding influence traits essential for disease resistance and environmental adaptability. This study examines the genomic adaptations that enable SSA cattle to thrive under these conditions and assesses the impact of inbreeding on such adaptive traits.

Methods

We analyzed genomic data from 113 cattle across four breeds-Kuri, N'dama, Zebu-Fulani, and Zebu-Bororo-employing Runs of Homozygosity (ROH) and Integrated Haplotype Score (iHS) analyses to identify historical and recent genetic selections. Strict quality controls using PLINK software ensured accurate genomic pattern identification related to adaptation and inbreeding.

Results

ROH analysis revealed islands with genes such as RSAD2, CMPK2, and NOTCH1, which are involved in immune response and cellular stress management, highlighting regions of historical selection that have likely provided adaptive advantages in overcoming environmental and pathogenic stresses. In contrast, iHS analysis identified genes under recent selection like HIPK1, involved in stress response regulation, and EPHA5, which plays a crucial role in neural development and synaptic functions, potentially equipping these breeds with novel adaptations to ongoing and emergent environmental challenges.

Conclusion

This research confirms that selective pressures inherent in pastoral and agropastoral systems profoundly influence the genetic structure of SSA cattle. By delineating the genetic bases of key adaptive traits, our study offers crucial insights for targeted breeding programs to enhance cattle resilience and productivity. These findings provide a valuable framework for future genetic improvements and conservation strategies, crucial for sustainable livestock management and economic stability in SSA.

Analysis of Runs of Homozygosity in Aberdeen Angus Cattle

A large number of cattle breeds have marked phenotypic differences. They are valuable models for studying genome evolution. ROH analysis can facilitate the discovery of genomic regions that may explain phenotypic differences between breeds affecting traits of economic importance. This paper investigates genome-wide ROH of 189 Aberdeen Angus bulls using the Illumina Bovine GGP HD Beadchip150K to structurally and functionally annotate genes located within or in close ROH of the Aberdeen Angus cattle genome. The method of sequential SNP detection was used to determine the ROH. Based on this parameter, two ROH classes were allocated. The total length of all ROH islands was 11,493 Mb. As a result of studying the genomic architecture of the experimental population of Aberdeen Angus bulls, nine ROH islands and 255 SNPs were identified. Thirteen of these overlapped with regions bearing 'selection imprints' previously identified in other breeds of cattle, and five of these regions were identified in other Aberdeen Angus populations. The total length of the ROH islands was 11,493 Mb. The size of individual islands ranged from 0.038 to 1.812 Mb. Structural annotation showed the presence of 87 genes within the identified ROH islets.

Using high-density SNP data to unravel the origin of the Franches-Montagnes horse breed

BACKGROUND

The Franches-Montagnes (FM) is the last native horse breed of Switzerland, established at the end of the 19th century by cross-breeding local mares with Anglo-Norman stallions. We collected high-density SNP genotype data (Axiom™ 670 K Equine genotyping array) from 522 FM horses, including 44 old-type horses (OF), 514 European Warmblood horses (WB) from Sweden and Switzerland (including a stallion used for cross-breeding in 1990), 136 purebred Arabians (AR), 32 Shagya Arabians (SA), and 64 Thoroughbred (TB) horses, as introgressed WB stallions showed TB origin in their pedigrees. The aim of the study was to ascertain fine-scale population structures of the FM breed, including estimation of individual admixture levels and genomic inbreeding (FROH) by means of Runs of Homozygosity.

RESULTS

To assess fine-scale population structures within the FM breed, we applied a three-step approach, which combined admixture, genetic contribution, and FROH of individuals into a high-resolution network visualization. Based on this approach, we were able to demonstrate that population substructures, as detected by model-based clustering, can be either associated with a different genetic origin or with the progeny of most influential sires. Within the FM breed, admixed horses explained most of the genetic variance of the current breeding population, while OF horses only accounted for a small proportion of the variance. Furthermore, we illustrated that FM horses showed high TB admixture levels and we identified inconsistencies in the origin of FM horses descending from the Arabian stallion Doktryner. With the exception of WB, FM horses were less inbred compared to the other breeds. However, the relatively few but long ROH segments suggested diversity loss in both FM subpopulations. Genes located in FM- and OF-specific ROH islands had known functions involved in conformation and behaviour, two traits that are highly valued by breeders.

CONCLUSIONS

The FM remains the last native Swiss breed, clearly distinguishable from other historically introgressed breeds, but it suffered bottlenecks due to intensive selection of stallions, restrictive mating choices based on arbitrary definitions of pure breeding, and selection of rare coat colours. To preserve the genetic diversity of FM horses, future conservation managements strategies should involve a well-balanced selection of stallions (e.g., by integrating OF stallions in the FM breeding population) and avoid selection for rare coat colours.

Genetic characteristics and selection signatures between Southern Chinese local and commercial chickens

The introduction of exotic breeds and the cultivation of new lines by breeding companies have posed challenges to native chickens in South China, including loss of breed characteristics, decreased genetic diversity, and declining purity. Understanding the population genetic structure and genetic diversity of native chickens in South China is crucial for further advancements in breeding efforts. In this study, we analyzed the population genetic structure and genetic diversity of 321 individuals from 10 different breeds in South China. By comparing commercial chickens with native ones, we identified selection signatures occurring between local chickens and commercial breeds. The analysis of population genetic structure revealed that the native chicken populations in South China exhibited a considerable level of genetic diversity. Moreover, the commercial lines of Xiaobai chicken and Huangma chicken displayed even higher levels of genetic diversity, which distinguished them from other native varieties at the clustering level. However, certain individuals within these commercial varieties showed a discernible genetic relationship with the native populations. Notably, both commercial varieties also retained a significant degree of genetic similarity to their respective native counterparts. In order to investigate the genomic changes occurring during the commercialization of native chickens, we employed 4 methods (Fst, ROD, XPCLR, and XPEHH) to identify potential candidate regions displaying selective signatures in Southern Chinese native chicken population. A total of 168 (identified by Fst and ROD) and 86 (identified by XPCLR and XPEHH) overlapping genes were discovered. Functional annotation analysis revealed that these genes may be associated with reproduction and growth (SAMSN1, HYLS1, ROBO3, FGF14, PRSS23), musculoskeletal development (DNER, MYBPC1, DGKB, ORC1, KLF10), disease resistance and environmental adaptability (PUS3, CRB2, CALD1, USP15, SGCD, LTBP1), as well as egg production (ADGRB3, ACSF3). Overall, native chickens in South China harbor numerous selective sweep regions compared to commercial chickens, enriching valuable genomic resources for future genetic research and breeding conservation.

Genomic Characterization of Local Croatian Sheep Breeds-Effective Population Size, Inbreeding & Signatures of Selection

The Istrian (IS) and the Pag sheep (PS) are local Croatian breeds which provide significant income for the regional economy and have a cultural and traditional importance for the inhabitants. The aim of this study was to estimate some important population specific genetic parameters in IS (N = 1293) and PS (N = 2637) based on genome wide SNPs. Estimates of linkage disequilibrium effective population size (Ne) evidenced more genetic variability in PS (Ne = 838) compared to IS (Ne = 197), regardless of historical time (both recent and ancient genetic variability). The discrepancy in the recent genetic variability between these breeds was additionally confirmed by the estimates of genomic inbreeding (FROH), which was estimated to be notably higher in IS (FROH>2 = 0.062) than in PS (FROH>2 = 0.029). The average FROH2-4, FROH4-8, FROH8-16, and FROH>16 were 0.26, 1.65, 2.14, and 3.72 for IS and 0.22, 0.61, 0.75, and 1.58 for PS, thus evidencing a high contribution of recent inbreeding in the overall inbreeding. One ROH island with > 30% of SNP incidence in ROHs was detected in IS (OAR6; 34,253,440-38,238,124 bp) while there was no ROH islands detected in PS. Seven genes (CCSER1, HERC3, LCORL, NAP1L5, PKD2, PYURF, and SPP1) involved in growth, feed intake, milk production, immune responses, and resistance were associated with the found autozygosity. The results of this study represent the first comprehensive insight into genomic variability of these two Croatian local sheep breeds and will serve as a baseline for setting up the most promising strategy of genomic Optimum Contribution Selection.

Unraveling genomic diversity and positive selection signatures of Qaidam cattle through whole-genome re-sequencing

Qaidam cattle are a typical Chinese native breed inhabiting northwest China. They bear the characteristics of high cold and roughage tolerance, low-oxygen adaptability and good meat quality. To analyze the genetic diversity of Qaidam cattle, 60 samples were sequenced using whole-genome resequencing technology, along with 192 published sets of whole-genome sequencing data of Indian indicine cattle, Chinese indicine cattle, North Chinese cattle breeds, East Asian taurine cattle, Eurasian taurine cattle and European taurine cattle as controls. It was found that Qaidam cattle have rich genetic diversity in Bos taurus, but the degree of inbreeding is also high, which needs further protection. The phylogenetic analysis, principal component analysis and ancestral component analysis showed that Qaidam cattle mainly originated from East Asian taurine cattle. Qaidam cattle had a closer genetic relationship with the North Chinese cattle breeds and the least differentiation from Mongolian cattle. Annotating the selection signals obtained by composite likelihood ratio, nucleotide diversity analysis, integrated haplotype score, genetic differentiation index, genetic diversity ratio and cross-population extended haplotype homozygosity methods, several genes associated with immunity, reproduction, meat, milk, growth and adaptation showed strong selection signals. In general, this study provides genetic evidence for understanding the germplasm characteristics of Qaidam cattle. At the same time, it lays a foundation for the scientific and reasonable protection and utilization of genetic resources of Chinese local cattle breeds, which has great theoretical and practical significance.

Analysis of 206 whole-genome resequencing reveals selection signatures associated with breed-specific traits in Hu sheep

As an invaluable Chinese sheep germplasm resource, Hu sheep are renowned for their high fertility and beautiful wavy lambskins. Their distinctive characteristics have evolved over time through a combination of artificial and natural selection. Identifying selection signatures in Hu sheep can provide a straightforward insight into the mechanism of selection and further uncover the candidate genes associated with breed-specific traits subject to selection. Here, we conducted whole-genome resequencing on 206 Hu sheep individuals, each with an approximate 6-fold depth of coverage. And then we employed three complementary approaches, including composite likelihood ratio, integrated haplotype homozygosity score and the detection of runs of homozygosity, to detect selection signatures. In total, 10 candidate genomic regions displaying selection signatures were simultaneously identified by multiple methods, spanning 88.54 Mb. After annotating, these genomic regions harbored collectively 92 unique genes. Interestingly, 32 candidate genes associated with reproduction were distributed in nine genomic regions detected. Out of them, two stood out as star candidates: BMPR1B and GNRH2, both of which have documented associations with fertility, and a HOXA gene cluster (HOXA1-5, HOXA9, HOXA10, HOXA11 and HOXA13) had also been linked to fertility. Additionally, we identified other genes that are related to hair follicle development (LAMTOR3, EEF1A2), ear size (HOXA1, KCNQ2), fat tail formation (HOXA10, HOXA11), growth and development (FAF1, CCNDBP1, GJB2, GJA3), fat deposition (ACOXL, JAZF1, HOXA3, HOXA4, HOXA5, EBF4), immune (UBR1, FASTKD5) and feed intake (DAPP1, RNF17, NPBWR2). Our results offer novel insights into the genetic mechanisms underlying the selection of breed-specific traits in Hu sheep and provide a reference for sheep genetic improvement programs.

Tracing the Dynamical Genetic Diversity Changes of Russian Livni Pigs during the Last 50 Years with the Museum, Old, and Modern Samples

The pig industry is usually considered an intensive livestock industry, mainly supported by hybrid breeding between commercial pig breeds. However, people's pursuit of a more natural environment and higher meat quality has led to an increasing demand for eco-friendly and diverse pig feeding systems. Therefore, the importance of rearing and conserving local pig breeds is increasing. The Livni pig is a local breed with good adaptability to the environmental and fodder conditions in central Russia. In this study, we aimed to analyze the genetic diversity and population structure of Livni pigs using whole-genome single nucleotide polymorphism (SNP) data. We utilized the Porcine GGP HD BeadChip on genotype samples from old (n = 32, 2004) and modern (n = 32, 2019) populations of Livni pigs. For the museum samples of Livni pigs (n = 3), we extracted DNA from their teeth, performed genomic sequencing, and obtained SNP genotypes from the whole-genome sequences. SNP genotypes of Landrace (n = 32) and Large White (n = 32) pigs were included for comparative analysis. We observed that the allelic richness of Livni pigs was higher than those of Landrace and Large White pigs (AR = 1.775-1.798 vs. 1.703 and 1.668, respectively). The effective population size estimates (NE5 = 108 for Livni pigs, NE5 = 59 for Landrace and Large White pigs) confirmed their genetic diversity tendency. This was further supported by the length and number of runs of homozygosity, as well as the genomic inbreeding coefficient (almost twofold lower in Livni pigs compared to Landrace and Large White pigs). These findings suggest that the Livni pig population exhibits higher genetic diversity and experiences lower selection pressure compared to commercial pig populations. Furthermore, both principal component and network tree analyses demonstrated a clear differentiation between Livni pigs and transboundary commercial pigs. The TreeMix results indicated gene flow from Landrace ancestors to Livni pigs (2019) and from Large White ancestors to Livni pigs (2004), which was consistent with their respective historical breeding backgrounds. The comparative analysis of museum, old, and modern Livni pigs indicated that the modern Livni pig populations have preserved their historical genomic components, suggesting their potential suitability for future design selection programs.

Whole-genome resequencing reveals genomic variation and dynamics in Ethiopian indigenous goats

Ethiopia has about 52 million indigenous goats with marked phenotypic variability, which is the outcome of natural and artificial selection. Here, we obtained whole-genome sequence data of three Ethiopian indigenous goat populations (Arab, Fellata, and Oromo) from northwestern Ethiopia and analyzed their genome-wide genetic diversity, population structure, and signatures of selection. We included genotype data from four other Ethiopian goat populations (Abergelle, Keffa, Gumuz, and Woyto-Guji) and goats from Asia; Europe; and eastern, southern, western, and northern Africa to investigate the genetic predisposition of the three Ethiopian populations and performed comparative genomic analysis. Genetic diversity analysis showed that Fellata goats exhibited the lowest heterozygosity values (Ho = 0.288 ± 0.005 and He = 0.334 ± 0.0001). The highest values were observed in Arab goats (Ho = 0.310 ± 0.010 and He = 0.347 ± 4.35e-05). A higher inbreeding coefficient (FROH = 0.137 ± 0.016) was recorded for Fellata goats than the 0.105 ± 0.030 recorded for Arab and the 0.112 ± 0.034 recorded for Oromo goats. This indicates that the Fellata goat population should be prioritized in future conservation activities. The three goat populations showed the majority (∼63%) of runs of homozygosity in the shorter (100-150 Kb) length category, illustrating ancient inbreeding and/or small founder effects. Population relationship and structure analysis separated the Ethiopian indigenous goats into two distinct genetic clusters lacking phylogeographic structure. Arab, Fellata, Oromo, Abergelle, and Keffa represented one genetic cluster. Gumuz and Woyto-Guji formed a separate cluster and shared a common genetic background with the Kenyan Boran goat. Genome-wide selection signature analysis identified nine strongest regions spanning 163 genes influencing adaptation to arid and semi-arid environments (HOXC12, HOXC13, HOXC4, HOXC6, and HOXC9, MAPK8IP2), immune response (IL18, TYK2, ICAM3, ADGRG1, and ADGRG3), and production and reproduction (RARG and DNMT1). Our results provide insights into a thorough understanding of genetic architecture underlying selection signatures in Ethiopian indigenous goats in a semi-arid tropical environment and deliver valuable information for goat genetic improvement, conservation strategy, genome-wide association study, and marker-assisted breeding.

A High-Coverage Mesolithic Aurochs Genome and Effective Leveraging of Ancient Cattle Genomes Using Whole Genome Imputation

Ancient genomic analyses are often restricted to utilizing pseudohaploid data due to low genome coverage. Leveraging low-coverage data by imputation to calculate phased diploid genotypes that enables haplotype-based interrogation and single nucleotide polymorphism (SNP) calling at unsequenced positions is highly desirable. This has not been investigated for ancient cattle genomes despite these being compelling subjects for archeological, evolutionary, and economic reasons. Here, we test this approach by sequencing a Mesolithic European aurochs (18.49×; 9,852 to 9,376 calBCE) and an Early Medieval European cow (18.69×; 427 to 580 calCE) and combine these with published individuals: two ancient and three modern. We downsample these genomes (0.25×, 0.5×, 1.0×, and 2.0×) and impute diploid genotypes, utilizing a reference panel of 171 published modern cattle genomes that we curated for 21.7 million (Mn) phased SNPs. We recover high densities of correct calls with an accuracy of >99.1% at variant sites for the lowest downsample depth of 0.25×, increasing to >99.5% for 2.0× (transversions only, minor allele frequency [MAF] ≥ 2.5%). The recovery of SNPs correlates with coverage; on average, 58% of sites are recovered for 0.25× increasing to 87% for 2.0×, utilizing an average of 3.5 million (Mn) transversions (MAF ≥2.5%), even in the aurochs, despite the highest temporal distance from the modern reference panel. Our imputed genomes behave similarly to directly called data in allele frequency-based analyses, for example consistently identifying runs of homozygosity >2 Mb, including a long homozygous region in the Mesolithic European aurochs.

Genetic legacy and adaptive signatures: investigating the history, diversity, and selection signatures in Rendena cattle resilient to eighteenth century rinderpest epidemics

BACKGROUND

Rendena is a dual-purpose cattle breed, which is primarily found in the Italian Alps and the eastern areas of the Po valley, and recognized for its longevity, fertility, disease resistance and adaptability to steep Alpine pastures. It is categorized as 'vulnerable to extinction' with only 6057 registered animals in 2022, yet no comprehensive analyses of its molecular diversity have been performed to date. The aim of this study was to analyse the origin, genetic diversity, and genomic signatures of selection in Rendena cattle using data from samples collected in 2000 and 2018, and shed light on the breed's evolution and conservation needs.

RESULTS

Genetic analysis revealed that the Rendena breed shares genetic components with various Alpine and Po valley breeds, with a marked genetic proximity to the Original Braunvieh breed, reflecting historical restocking efforts across the region. The breed shows signatures of selection related to both milk and meat production, environmental adaptation and immune response, the latter being possibly the result of multiple rinderpest epidemics that swept across the Alps in the eighteenth century. An analysis of the Rendena cattle population spanning 18 years showed an increase in the mean level of inbreeding over time, which is confirmed by the mean number of runs of homozygosity per individual, which was larger in the 2018 sample.

CONCLUSIONS

The Rendena breed, while sharing a common origin with Brown Swiss, has developed distinct traits that enable it to thrive in the Alpine environment and make it highly valued by local farmers. Preserving these adaptive features is essential, not only for maintaining genetic diversity and enhancing the ability of this traditional animal husbandry to adapt to changing environments, but also for guaranteeing the resilience and sustainability of both this livestock system and the livelihoods within the Rendena valley.

Exploring genomic inbreeding and selection signatures in a commercial Brangus herd through functional annotation

Composite breeds, including Brangus, are widely utilized in subtropical and tropical regions to harness the advantages of both Bos t. taurus and Bos t. indicus breeds. The formation and subsequent selection of composite breeds may result in discernible signatures of selection and shifts in genomic population structure. The objectives of this study were to 1) assess genomic inbreeding, 2) identify signatures of selection, 3) assign functional roles to these signatures in a commercial Brangus herd, and 4) contrast signatures of selection between selected and non-selected cattle from the same year. A total of 4035 commercial Brangus cattle were genotyped using the GGP-F250K array. Runs of Homozygosity (ROH) were used to identify signatures of selection and calculate genomic inbreeding. Quantitative trait loci (QTL) enrichment analysis and literature search identified phenotypic traits linked to ROH islands. Genomic inbreeding averaged 5%, primarily stemming from ancestors five or more generations back. A total of nine ROH islands were identified, QTL enrichment analysis revealed traits related to growth, milk composition, carcass, reproductive, and meat quality traits. Notably, the ROH island on BTA14 encompasses the pleiomorphic adenoma (PLAG1) gene, which has been linked to growth, carcass, and reproductive traits. Moreover, ROH islands associated with milk yield and composition were more pronounced in selected replacement heifers of the population, underscoring the importance of milk traits in cow-calf production. In summary, our research sheds light on the changing genetic landscape of the Brangus breed due to selection pressures and reveals key genomic regions impacting production traits.

Selective signatures in composite MONTANA TROPICAL beef cattle reveal potential genomic regions for tropical adaptation

Genomic regions related to tropical adaptability are of paramount importance for animal breeding nowadays, especially in the context of global climate change. Moreover, understanding the genomic architecture of these regions may be very relevant for aiding breeding programs in choosing the best selection scheme for tropical adaptation and/or implementing a crossbreeding scheme. The composite MONTANA TROPICAL® population was developed by crossing cattle of four different biological types to improve production in harsh environments. Pedigree and genotype data (51962 SNPs) from 3215 MONTANA TROPICAL® cattle were used to i) characterize the population structure; ii) identify signatures of selection with complementary approaches, i.e. Integrated Haplotype Score (iHS) and Runs of Homozygosity (ROH); and iii) understand genes and traits related to each selected region. The population structure based on principal components had a weak relationship with the genetic contribution of the different biological types. Clustering analyses (ADMIXTURE) showed different clusters according to the number of generations within the composite population. Considering results of both selection signatures approaches, we identified only one consensus region on chromosome 20 (35399405-40329703 bp). Genes in this region are related to immune function, regulation of epithelial cell differentiation, and cell response to ionizing radiation. This region harbors the slick locus which is related to slick hair and epidermis anatomy, both of which are related to heat stress adaptation. Also, QTLs in this region were related to feed intake, milk yield, mastitis, reproduction, and slick hair coat. The signatures of selection detected here arose in a few generations after crossbreeding between contrasting breeds. Therefore, it shows how important this genomic region may be for these animals to thrive in tropical conditions. Further investigations on sequencing this region can identify candidate genes for animal breeding and/or gene editing to tackle the challenges of climate change.

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.

Genetic Variation and Composition of Two Commercial Estonian Dairy Cattle Breeds Assessed by SNP Data

The aims of this study were to assess the genomic relatedness of Estonian and selected European dairy cattle breeds and to examine the within-breed diversity of two Estonian dairy breeds using genome-wide SNP data. This study was based on a genotyped heifer population of the Estonian Red (ER) and Estonian Holstein (EH) breeds, including about 10% of all female cattle born in 2017-2020 (sample sizes n = 215 and n = 2265, respectively). The within-breed variation study focused on the level of inbreeding using the ROH-based inbreeding coefficient. The genomic relatedness analyses were carried out among two Estonian and nine European breeds from the WIDDE database. Admixture analysis revealed the heterogeneity of ER cattle with a mixed pattern showing several ancestral populations containing a relatively low proportion (1.5-37.0%) of each of the reference populations used. There was a higher FROH in EH (FROH = 0.115) than in ER (FROH = 0.044). Compared to ER, the long ROHs of EH indicated more closely related parents. The paternal origin of the genetic material used in breeding had a low effect on the inbreeding level. However, among EH, the highest genomic inbreeding was estimated in daughters of USA-born sires.

Resequencing Analyses Revealed Genetic Diversity and Selection Signatures during Rabbit Breeding and Improvement

Domestication has shaped the diverse characteristics of rabbits, including coat color, fur structure, body size, and various physiological traits. Utilizing whole-genome resequencing (DNBSEQ-T7), we analyzed the genetic diversity, population structure, and genomic selection across 180 rabbits from 17 distinct breeds to uncover the genetic basis of these traits. We conducted whole-genome sequencing on 17 rabbit breeds, identifying 17,430,184 high-quality SNPs and analyzing genomic diversity, patterns of genomic variation, population structure, and selection signatures related to coat color, coat structure, long hair, body size, reproductive capacity, and disease resistance. Through PCA and NJ tree analyses, distinct clusters emerged among Chinese indigenous rabbits, suggesting varied origins and domestication histories. Selective sweep testing pinpointed regions and genes linked to domestication and key morphological and economic traits, including those affecting coat color (TYR, ASIP), structure (LIPH), body size (INSIG2, GLI3), fertility (EDNRA, SRD5A2), heat stress adaptation (PLCB1), and immune response (SEC31A, CD86, LAP3). Our study identified key genomic signatures of selection related to traits such as coat color, fur structure, body size, and fertility; these findings highlight the genetic basis underlying phenotypic diversification in rabbits and have implications for breeding programs aiming to improve productive, reproductive, and adaptive traits. The detected genomic signatures of selection also provide insights into rabbit domestication and can aid conservation efforts for indigenous breeds.

Genome-wide analysis of population structure, effective population size and inbreeding in Iranian and exotic horses

Population structure and genetic diversity are the key parameters to study the breeding history of animals. This research aimed to provide a characterization of the population structure and to compare the effective population size (Ne), LD decay, genetic diversity, and genomic inbreeding in Iranian native Caspian (n = 38), Turkmen (n = 24) and Kurdish (n = 29) breeds and some other exotic horses consisting of Arabian (n = 24), Fell pony (n = 21) and Akhal-Teke (n = 20). A variety of statistical population analysis techniques, such as principal component analysis (PCA), discriminant analysis of principal component (DAPC) and model-based method (STRUCTURE) were employed. The results of the population analysis clearly demonstrated a distinct separation of native and exotic horse breeds and clarified the relationships between studied breeds. The effective population size (Ne) for the last six generations was estimated 54, 49, 37, 35, 27 and 26 for the Caspian, Kurdish, Arabian, Turkmen, Akhal-Teke and Fell pony breeds, respectively. The Caspian breed showed the lowest LD with an average r2 value of 0.079, while the highest was observed in Fell pony (0.148). The highest and lowest average observed heterozygosity were found in the Kurdish breeds (0.346) and Fell pony (0.290) breeds, respectively. The lowest genomic inbreeding coefficient based on run of homozygosity (FROH) and excess of homozygosity (FHOM) was in the Caspian and Kurdish breeds, respectively, while based on genomic relationship matrix) FGRM) and correlation between uniting gametes) FUNI) the lowest genomic inbreeding coefficient was found in the Kurdish breed. The estimation of genomic inbreeding rates in the six breeds revealed that FROH yielded lower estimates compared to the other three methods. Additionally, the Iranian breeds displayed lower levels of inbreeding compared to the exotic breeds. Overall, the findings of this study provide valuable insights for the development of effective breeding management strategies aimed at preserving these horse breeds.

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

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).

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight

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.

Assessment of runs of homozygosity, heterozygosity-rich regions and genomic inbreeding estimates in a subpopulation of Guzerá (Bos indicus) dual-purpose cattle

For decades, inbreeding in cattle has been evaluated using pedigree information. Nowadays, inbreeding coefficients can be obtained using genomic information such as runs of homozygosity (ROH). The aims of this study were to quantify ROH and heterozygosity-rich regions (HRR) in a subpopulation of Guzerá dual-purpose cattle, to examine ROH and HRR islands, and to compare inbreeding coefficients obtained by ROH with alternative genomic inbreeding coefficients. A subpopulation of 1733 Guzerá animals genotyped for 50k SNPs was used to obtain the ROH and HRR segments. Inbreeding coefficients by ROH (FROH ), by genomic relationship matrix based on VanRaden's method 1 using reference allele frequency in the population (FGRM ), by genomic relationship matrix based on VanRaden's method 1 using allele frequency fixed in 0.5 (FGRM_0.5 ), and by the proportion of homozygous loci (FHOM ) were calculated. A total of 15,660 ROH were identified, and the chromosome with the highest number of ROH was BTA6. A total of 4843 HRRs were identified, and the chromosome with the highest number of HRRs was BTA23. No ROH and HRR islands were identified according to established criteria, but the regions closest to the definition of an island were examined from 64 to 67 Mb of BTA6, from 36 to 37 Mb of BTA2 and from 0.50 to 1.25 Mb of BTA23. The genes identified in ROH islands have previously been associated with dairy and beef traits, while genes identified on HRR islands have previously been associated with reproductive traits and disease resistance. FROH was equal to 0.095 ± 0.084, and its Spearman correlation with FGRM was low (0.44) and moderate-high with FHOM (0.79) and with FGRM_0.5 (0.80). The inbreeding coefficients determined by ROH were higher than other cattle breeds' and higher than pedigree-based inbreeding in the Guzerá breed obtained in previous studies. It is recommended that future studies investigate the effects of inbreeding determined by ROH on the traits under selection in the subpopulation studied.

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.

Genome-wide detection of positive and balancing signatures of selection shared by four domesticated rainbow trout populations (Oncorhynchus mykiss)

BACKGROUND

Evolutionary processes leave footprints along the genome over time. Highly homozygous regions may correspond to positive selection of favorable alleles, while maintenance of heterozygous regions may be due to balancing selection phenomena. We analyzed data from 176 fish from four disconnected domestic rainbow trout populations that were genotyped using a high-density Axiom Trout genotyping 665K single nucleotide polymorphism array, including 20 from the US and 156 from three French lines. Using methods based on runs of homozygosity and extended haplotype homozygosity, we detected signatures of selection in these four populations.

RESULTS

Nine genomic regions that included 253 genes were identified as being under positive selection in all four populations Most were located on chromosome 2 but also on chromosomes 12, 15, 16, and 20. In addition, four heterozygous regions that contain 29 genes that are putatively under balancing selection were also shared by the four populations. These were located on chromosomes 10, 13, and 19. Regardless of the homozygous or heterozygous nature of the regions, in each region, we detected several genes that are highly conserved among vertebrates due to their critical roles in cellular and nuclear organization, embryonic development, or immunity. We identified new candidate genes involved in rainbow trout fitness, as well as 17 genes that were previously identified to be under positive selection, 10 of which in other fishes (auts2, atp1b3, zp4, znf135, igf-1α, brd2, col9a2, mrap2, pbx1, and emilin-3).

CONCLUSIONS

Using material from disconnected populations of different origins allowed us to draw a genome-wide map of signatures of positive selection that are shared between these rainbow trout populations, and to identify several regions that are putatively under balancing selection. These results provide a valuable resource for future investigations of the dynamics of genetic diversity and genome evolution during domestication.

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.

The Addis Ababa Lions: Whole-Genome Sequencing of a Rare and Precious Population

Lions are widely known as charismatic predators that once roamed across the globe, but their populations have been greatly affected by environmental factors and human activities over the last 150 yr. Of particular interest is the Addis Ababa lion population, which has been maintained in captivity at around 20 individuals for over 75 yr, while many wild African lion populations have become extinct. In order to understand the molecular features of this unique population, we conducted a whole-genome sequencing study on 15 Addis Ababa lions and detected 4.5 million distinct genomic variants compared with the reference African lion genome. Using functional annotation, we identified several genes with mutations that potentially impact various traits such as mane color, body size, reproduction, gastrointestinal functions, cardiovascular processes, and sensory perception. These findings offer valuable insights into the genetics of this threatened lion population.

Relatedness and genomic inbreeding in a sample of Timor ponies

Timor ponies (TP) were first shipped to Australia in the early 1800s and were highly valued as transport and pack animals, which resulted in TPs contributing to the development of Australian horse breeds. Today, while the exact number of TPs in Australia is currently unknown, there has been recent interest in establishing a domestic breeding program for Australian TPs. The aim of this study was to evaluate the relatedness of a sample of TPs, as well as provide estimates of genomic inbreeding levels to better inform the feasibility of using these animals as founders for a domestic breeding program. Hair samples from each horse were genotyped using the Illumina 80K Infinium Equine genotyping array and data were analysed using PLINK v1.90b7, KING 2.3.2 and R v4.3.1. The results illustrate that there are distantly related and minimally inbred horses within the sampled TPs. Lengths of the ROH segments also indicated that recent inbreeding events are likely to only have occurred in a third of the horses. Overall, these results are promising for the success of a domestic TP breeding program; however, considering the low number of domestic TPs known to reside in Australia, there would certainly still be substantial benefits to incorporating additional TPs either directly from Timor or from areas in Australia that are believed to contain wild descendants of TPs.

Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle

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.

Comparative peculiarities of genomic diversity in Gallus gallus domesticus chickens with decorative plumage: the muffs and beard phenotype

Throughout history, humans have been attempting to develop the ornamental features of domestic animals in addition to their productive qualities. Many chicken breeds have developed tufts of elongated feathers that jut out from the sides and bottom of the beak, leading to the phenotype known as muffs and beard. It is an incomplete autosomal dominant phenotype determined by the Mb locus localised on chromosome GGA27. This project aimed to analyse the genetic diversity of chicken breeds using full genomic genotyping with the Chicken 60K BeadChip. A total of 53,313 Single Nucleotide Polymorphisms were analysed. DNA was obtained from breeds with the muffs and beard as a marker phenotype: Faverolles (n = 20), Ukrainian Muffed (n = 18), Orloff (n = 20), Novopavlov White (n = 20), and Novopavlov Coloured (n = 15). The Russian White (n = 20) was selected as an alternative breed without the muffs and beard phenotype. The chickens are owned by the Centre of Collective Use "Genetic Collection of Rare and Endangered Breeds of Chickens" (St. Petersburg region, Pushkin), and are also included in the Core Shared Research Facility (CSRF) and/or Large-Scale Research Facility (LSRF). Multidimensional scaling revealed that the Novopavlov White and the Novopavlov Coloured populations formed a separate group. The Ukrainian Muffed and the Orloff have also been combined into a separate group. Based on cluster analysis, with the cross-validation error and the most probable number of clusters K = 4 taken into account, the Orloff was singled out as a separate group. The Ukrainian Muffed exhibited a notable similarity with the Orloff under the same conditions. At K = 5, the populations of the Novopavlov White and the Novopavlov Coloured diverged. Only at K = 6, a distinct and separate cluster was formed by the Ukrainian Muffed. The Russian White had the greatest number of short (1-2 Mb) homozygous regions. If the HOXB8 gene is located between 3.402 and 3.404 Mb on chromosome GGA27, homozygous regions are rarely found in the chickens with the muffs and beard phenotype. Scanning the chicken genome with the Chicken 60K BeadChip provided enough information about the genetic diversity of the chicken breeds for the peculiarities of the development of the ornamental muffs and beard phenotypes in them to be understood. For example, Phoenix bantams, whose tail feathers grow throughout their lives, require greater consideration of husbandry conditions.

Whole genome sequencing identified genomic diversity and candidated genes associated with economic traits in Northeasern Merino in China

Introduction: Northeast Merino (NMS) is a breed developed in Northeast China during the 1960s for wool and meat production. It exhibits excellent traits such as high wool yield, superior meat quality, rapid growth rate, robust disease resistance, and adaptability to cold climates. However, no studies have used whole-genome sequencing data to investigate the superior traits of NMS. Methods: In this study, we investigated the population structure, genetic diversity, and selection signals of NMS using whole-genome sequencing data from 20 individuals. Two methods (integrated haplotype score and composite likelihood ratio) were used for selection signal analysis, and the Fixation Index was used to explore the selection signals of NMS and the other two breeds, Mongolian sheep and South African meat Merino. Results: The results showed that NMS had low inbreeding levels, high genomic diversity, and a pedigree of both Merino breeds and Chinese local breeds. A total length of 14.09 Mb genomic region containing 287 genes was detected using the two methods. Further exploration of the functions of these genes revealed that they are mainly concentrated in wool production performance (IRF2BP2, MAP3K7, and WNT3), meat production performance (NDUFA9, SETBP1, ZBTB38, and FTO), cold resistance (DNAJC13, LPGAT1, and PRDM16), and immune response (PRDM2, GALNT8, and HCAR2). The selection signals of NMS and the other two breeds annotated 87 and 23 genes, respectively. These genes were also mainly focused on wool and meat production performance. Conclusion: These results provide a basis for further breeding improvement, comprehensive use of this breed, and a reference for research on other breeds.

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.

Effect of genotyping density on the detection of runs of homozygosity and heterozygosity in cattle

Runs of homozygosity (ROHom) are contiguous stretches of homozygous regions of the genome. In contrast, runs of heterozygosity (ROHet) are heterozygosity-rich regions. The detection of these two types of genomic regions (ROHom and ROHet) is influenced by the parameters involved in their identification and the number of available single-nucleotide polymorphisms (SNPs). The present study aimed to test the effect of chip density in detecting ROHom and ROHet in the Italian Simmental cattle breed. A sample of 897 animals were genotyped at low density (50k SNP; 397 individuals), medium density (140k SNP; 348 individuals), or high density (800k SNP; 152 individuals). The number of ROHom and ROHet per animal (nROHom and nROHet, respectively) and their average length were calculated. ROHom or ROHet shared by more than one animal and the number of times a particular SNP was inside a run were also computed (SNPROHom and SNPROHet). As the chip density increased, the nROHom increased, whereas their average length decreased. In contrast, the nROHet decreased and the average length increased as the chip density increased. The most repeated ROHom harbored no genes, whereas in the most repeated ROHet four genes (SNRPN, SNURF, UBE3A, and ATP10A) previously associated with reproductive traits were found. Across the 3 datasets, 31 SNP, located on Bos taurus autosome (BTA) 6, and 37 SNP (located on BTA21) exceeded the 99th percentile in the distribution of the SNPROHom and SNPROHet, respectively. The genomic region on BTA6 mapped the SLIT2, PACRGL, and KCNIP4 genes, whereas 19 and 18 genes were mapped on BTA16 and BTA21, respectively. Interestingly, most of genes found through the ROHet analysis were previously reported to be related to health, reproduction, and fitness traits. The results of the present study confirm that the detection of ROHom is more reliable when the chip density increases, whereas the ROHet trend seems to be the opposite. Genes and quantitative trait loci (QTL) mapped in the highlighted regions confirm that ROHet can be due to balancing selection, thus related to fitness traits, health, and reproduction, whereas ROHom are mainly involved in production traits. The results of the present study strengthened the usefulness of these parameters in analyzing the genomes of livestock and their biological meaning.

Genome-wide detection of selection signatures in Jianli pigs reveals novel cis-regulatory haplotype in EDNRB associated with two-end black coat color

BACKGROUND

Jianli pig, a renowned indigenous breed in China, has the characteristics of a two-end black (TEB) coat color, excellent meat quality, strong adaptability and increased prolificacy. However, there is limited information available regarding the genetic diversity, population structure and genomic regions under selection of Jianli pig. On the other hand, the genetic mechanism of TEB coat color has remained largely unknown.

RESULTS

In this study, the whole genome resequencing of 30 Jianli pigs within a context of 153 individuals representing 13 diverse breeds was performed. The population structure analysis revealed that Jianli pigs have close genetic relationships with the Tongcheng pig breed, their geographical neighbors. Three methods (observed heterozygosity, expected heterozygosity, and runs of homozygosity) implied a relatively high level of genetic diversity and, a low inbreeding coefficient in Jianli compared with other pigs. We used Fst and XP-EHH to detect the selection signatures in Jianli pigs compared with Asian wild boar. A total of 451 candidate genes influencing meat quality (CREBBP, ADCY9, EEPD1 and HDAC9), reproduction (ESR1 and FANCA), and coat color (EDNRB, MITF and MC1R), were detected by gene annotation analysis. Finally, to fine-map the genomic region for the two-end black (TEB) coat color phenotype in Jianli pigs, we performed three signature selection methods between the TEB coat color and no-TEB coat color pig breeds. The current study, further confirmed that the EDNRB gene is a candidate gene for TEB color phenotype found in Chinese pigs, including Jinhua pigs, and the haplotype harboring 25 SNPs in the EDNRB gene may promote the formation of TEB coat color. Further ATAC-seq and luciferase reporter assays of these regions suggest that the 25-SNPs region was a strong candidate causative mutation that regulates the TEB coat color phenotype by altering enhancer function.

CONCLUSION

Our results advanced the understanding of the genetic mechanism behind artificial selection, and provided further resources for the protection and breeding improvement of Jianli pigs.

The novel IMAGE001 genotyping array as a valuable alternative for genetic diversity screening in chicken: a demonstration in a local chicken breed in Belgium

Screening for genetic diversity in livestock species breeds is of utmost importance, especially for local, small populations that are at the risk of extinction. Luckily, recent developments in technology increase access to genotyping, also for numerically small breeds. One of these new technologies is the IMAGE001 single nucleotide polymorphism genotyping array that includes markers for 6 different species (cow, pig, sheep, chicken, horse and goat). For our current study, we studied the Turkey-headed Malines chicken, a local chicken breed in Belgium, for the first time. A total of 110 animals were genotyped, together with 29 samples from 4 supposedly related breeds. The genotypes were used to assess the genetic diversity of this local breed. Our analysis revealed an average inbreeding coefficient of 0.20 through runs of homozygosity analysis, and effective population size estimation based on linkage disequilibrium indicated a low genetic diversity (Ne = 34). Moreover, a principal component analysis and a genetic differentiation study (FST) were performed using these marker data to position the Turkey-headed Malines relative to the 4 other indigenous Belgian chicken breeds. Finally, we discussed the practical implications of the overlap between the IMAGE001 array and other existing chicken genotyping arrays. This study is the first use of the novel IMAGE001 array to evaluate a local chicken breed, and demonstrates it as a viable option for genomic characterization a breed. Moreover, with this research, we are able to provide a good basis for further evaluation of the Belgian chicken heritage.

Genome-Wide Association Studies and Runs of Homozygosity to Identify Reproduction-Related Genes in Yorkshire Pig Population

Reproductive traits hold considerable economic importance in pig breeding and production. However, candidate genes underpinning the reproductive traits are still poorly identified. In the present study, we executed a genome-wide association study (GWAS) and runs of homozygosity (ROH) analysis using the PorcineSNP50 BeadChip array for 585 Yorkshire pigs. Results from the GWAS identified two genome-wide significant and eighteen suggestive significant single nucleotide polymorphisms (SNPs) associated with seven reproductive traits. Furthermore, we identified candidate genes, including ELMO1, AOAH, INSIG2, NUP205, LYPLAL1, RPL34, LIPH, RNF7, GRK7, ETV5, FYN, and SLC30A5, which were chosen due to adjoining significant SNPs and their functions in immunity, fertilization, embryonic development, and sperm quality. Several genes were found in ROH islands associated with spermatozoa, development of the fetus, mature eggs, and litter size, including INSL6, TAF4B, E2F7, RTL1, CDKN1C, and GDF9. This study will provide insight into the genetic basis for pig reproductive traits, facilitating reproduction improvement using the marker-based selection methods.

Assessing genomic diversity and selective pressures in Bashan cattle by whole-genome sequencing data

Specific ecological environments and domestication have continuously influenced the physiological characteristics of Chinese indigenous cattle. Among them, Bashan cattle belongs to one of the indigenous breeds. However, the genomic diversity of Bashan cattle is still unknown. Published whole-genome sequencing (WGS) data of 13 Bashan cattle and 48 worldwide cattle were used to investigate the genetic composition and selection characteristics of Bashan cattle. The population structure analysis revealed that Bashan cattle harbored ancestries with East Asian taurine and Chinese indicine. Genetic diversity analysis implied the relatively high genomic diversity in Bashan cattle. Through the identification of containing >5 nsSNPs or frameshift mutations genes in Bashan cattle, a large number of pathways related to sensory perception were discovered. CLR, θπ ratio, FST, and XP-EHH methods were used to detect the candidate signatures of positive selection in Bashan cattle. Among the identified genes, most of the enriched signal pathways were related to environmental information processing, biological systems, and metabolism. We mainly reported genes related to the nervous system (HCN1, KATNA1, FSTL1, GRIK2, and CPLX2), immune (CD244, SLAMF1, LY9, and CD48), and reproduction (AKR1C1, AKR1C3, AKR1C4, and TUSC3). Our findings will be significant in understanding the molecular basis underlying phenotypic variation of breed-related traits and improving productivity in Bashan cattle.