Genomic Inbreeding Coefficients Research Articles

Development of a haplotype reference panel and analysis of genetic structure in Golden Pompano (Trachinotus ovatus)

Trachinotus ovatus, a valuable marine fish, is primarily farmed in the South China Sea. Intensive aquaculture can lead to genetic uniformity and disease risks. Additionally, the high cost of genome sequencing restricts knowledge of genetic variation, hindering genetic analysis and breeding advancements. This study investigated the genetic diversity and structure of Trachinotus ovatus populations within the South China Sea aquaculture context, with the goal of strengthening sustainability. Additionally, the study also enhanced breeding strategies by developing a high-quality haplotype reference panel and an efficient imputation process. Through high-throughput sequencing (depth＞20X), 173 samples from three key breeding sites (SanYa, Hainan; Zhanjiang, Guangdong; and Beihai, Guangxi) were analyzed, revealing substantial genetic diversity and distinct population structures. Notably, the Beihai (BH) population exhibited significant stratification, suggesting regional genetic uniqueness. Low levels of genomic inbreeding coefficient and the identification of unique SNP loci across populations underscored the genetic richness within these Trachinotus ovatus communities. A major constructed of this work is the creation of the first high-quality haplotype reference panel for the species, demonstrating genotype imputation's viability as a cost-effective alternative for large-scale genetic studies in aquaculture. In particular, a Shapeit reference panel with 5476,550 SNPs was created, excelling in accuracy across various depths, notably achieving up to 95.8 % imputation accuracy at 1X sequencing depth. This high accuracy across populations shows the importance of genetic diversity in the panel. This research not only advances the understanding of Trachinotus ovatus genetic foundations but also provides crucial insights for conservation efforts and genetic analysis. By establishing a robust haplotype reference panel for Trachinotus ovatus helps lower the costs of future genetic research and variation studies, this study lays the foundation for sustainable aquaculture practices and breeding programs.

Genetic Diversity, Selection Signatures, and Genome-Wide Association Study Identify Candidate Genes Related to Litter Size in Hu Sheep.

Hu sheep is a renowned prolific local sheep breed in China, widely distributed across the country due to its excellent reproductive performance. Deciphering the molecular mechanisms underlying the high fecundity of Hu sheep is crucial for improving the litter size of ewes. In this study, we genotyped 830 female Hu sheep using the Illumina OvineSNP50 BeadChip and performed genetic diversity analysis, selection signature detection, and a genome-wide association study (GWAS) for litter size. Our results revealed that the Hu sheep population exhibits relatively high genetic diversity. A total of 4927 runs of homozygosity (ROH) segments were detected, with the majority (74.73%) being short in length. Different genomic inbreeding coefficients (FROH, FHOM, FGRM, and FUNI) ranged from -0.0060 to 0.0126, showing low levels of inbreeding in this population. Additionally, we identified 91 candidate genomic regions through three complementary selection signature methods, including ROH, composite likelihood ratio (CLR), and integrated haplotype score (iHS), and annotated 189 protein-coding genes. Moreover, we observed two significant SNPs related to the litter size of Hu sheep using GWAS analysis based on a repeatability model. Integrating the selection signatures and the GWAS results, we identified 15 candidate genes associated with litter size, among which BMPR1B and UNC5C were particularly noteworthy. These findings provide valuable insights for improving the reproductive performance and breeding of high-fecundity lines of Hu sheep.

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

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.

Genetic diversity and signatures of selection in Icelandic horses and Exmoor ponies

BackgroundThe Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds.ResultsHO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20).ConclusionsEvaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

Whole-genome resequencing identifies candidate genes associated with heat adaptation in chickens

The wide distribution and diverse varieties of chickens make them important models for studying genetic adaptation. The aim of this study was to identify genes that alter heat adaptation in commercial chicken breeds by comparing genetic differences between tropical and cold-resistant chickens. We analyzed whole-genome resequencing data of 186 chickens across various regions in Asia, including the following breeds: Bian chickens (B), Dagu chickens (DG), Beijing-You chickens (BY), and Gallus gallus jabouillei from China; Gallus gallus murghi from India; Vietnam native chickens (VN); Thailand native chickens (TN) and Gallus gallus spadiceus from Thailand; and Indonesia native chickens (IN), Gallus gallus gallus, and Gallus gallus bankiva from Indonesia. In total, 5,454,765 SNPs were identified for further analyses. Population genetic structure analysis revealed that each local chicken breed had undergone independent evolution. Additionally, when K = 5, B, BY, and DG chickens shared a common ancestor and exhibited high levels of inbreeding, suggesting that northern cold-resistant chickens are likely the result of artificial selection. In contrast, the runs of homozygosity (ROH) and the ROH-based genomic inbreeding coefficient (FROH) results for IN, TN, and VN chickens showed low levels of inbreeding. Low population differentiation index values indicated low differentiation levels, suggesting low genetic diversity in tropical chickens, implying increased vulnerability to environmental changes, decreased adaptability, and disease resistance. Whole-genome selection sweep analysis revealed 69 candidate genes, including LGR4, G6PC, and NBR1, between tropical and cold-resistant chickens. The genes were further subjected to GO and KEGG enrichment analyses, revealing that most of the genes were primarily enriched in biological synthesis processes, metabolic processes, central nervous system development, ion transmembrane transport, and the Wnt signaling pathway. Our study identified heat adaptation genes and their functions in chickens that primarily affect chickens in high-temperature environments through metabolic pathways. These heat-resistance genes provide a theoretical basis for improving the heat-adaptation capacity of commercial chicken breeds.

Development of Single-Nucleotide Polymorphism Markers and Population Genetic Analysis of the Hadal Amphipod Alicella gigantea across the Mariana and New Britain Trenches

Alicella gigantea, the largest amphipod scavengers found to date, play key roles in the food web of the hadal ecosystem. However, the genetic structure of A. gigantea populations among different trenches has not been reported yet. In this study, SNP (single-nucleotide polymorphism) markers were developed for three A. gigantea geographic populations collected from the southern Mariana Trench (SMT), the central New Britain Trench (CNBT), and the eastern New Britain Trench (ENBT), based on the SLAF-seq (specific locus amplified fragment sequencing) technology. A total of 570,168 filtered SNPs were screened out for subsequent population genetic analysis. Results showed that the inbreeding levels across the three geographic populations were relatively low, and the genomic inbreeding coefficients of the three populations were similar in magnitude. Based on the results of phylogenetic analysis, population structure analysis, and principal component analysis, it is believed that the three A. gigantea geographic populations belong to the same population, and the kinship relationship between the ENBT and CNBT populations is close. Moreover, the differential candidate adaptive sites on the SNPs suggest that there may be variations in metabolic rates among the three geographic populations, possibly linked to differences in food availability and sources in different trenches, ultimately resulting in different survival strategies in A. gigantea populations within distinct trenches. Compared with the Mariana Trench, the New Britain Trench has a richer organic matter input, and it is speculated that the A. gigantea Mariana Trench population may adopt a lower metabolic rate to cope with the harsher environment of nutrient deficiency.

Comprehensive analysis of runs of homozygosity and heterozygosity in Holstein cattle on the basis of medium and high density SNP panels and large population sample

Abstract This study reports runs of homozygosity (ROH) and heterozygosity (ROHet) distributed in a large population of Holstein cattle on the basis of two microarrays of medium (50k; 2163 animals; 54 609 SNPs) and high single nucleotide polymorphism (SNP) density (HD; 600 animals; 777 692 SNPs). To assess the inbreeding values of Holstein cattle, the ROH-based genomic inbreeding coefficient (FROH) was calculated. The comparison of SNP panels suggested that FROH values above 4 Mb should be considered for panels of medium densities as a relatively reliable measure of inbreeding. Moreover, ROH hotspots and coldspots were identified and compared between the HD and 50k SNP panels and were carefully examined for association with production and functional traits. The obtained results pinpointed genomic regions presumably under selection pressure in Holstein cattle. The regions overlapped with a large number of genes, including GHR, GBF1, SUMF1, CCL28, NIM1K, U6, BTRC and FABP1, many of which are involved in important Holstein cattle characteristics. We also found that some ROH hotspots and coldspots identified with the HD panel were not detected with the 50k panel, mainly because of insufficient SNP density in certain genomic regions. This suggests that using medium-density panels might not be the best choice when precise identification of ROH patterns is the main goal. In summary, in this work, we confirmed that a high-density SNP panel compared to a medium-density SNP panel allows for more precise identification of ROH patterns, especially in the case of short ROH that could be associated with ancestral inbreeding.

Genomic inbreeding coefficients using imputation genotypes: Assessing the effect of ancestral genotyping in Holstein-Friesian dairy cows

The objective of this study was to assess the effect of using or not the genotypes of the parents of a cow for imputing single nucleotide polymorphisms (SNP), on the estimation of genomic inbreeding coefficients of cows. Imputation (i.e., genotyped plus imputed) genotypes from 68,127 Italian Holstein dairy cows registered in the Italian National Association of Holstein, Brown and Jersey Breeders (ANAFIBJ) were analyzed. Cows were genotyped with the HD Illumina Infinium BovineHD BeadChip and GeneSeek Genomic Profiler HD-150K, and the MD GeneSeek Genomic Profiler 3, GeneSeek Genomic Profiler 4, GeneSeek MD and the Labogena MD. To assess differences among estimators genomic inbreeding coefficients were estimated with 4 PLINK v1.9 estimators (F, Fhat1, 2, 3), 2 genomic relationship matrix (grm) based estimators (Fgrm and Fgrm2; with the latter including also pedigree information) and one estimator of runs of homozygosity (ROH; FROH). Assuming that the correct genomic inbreeding coefficients should be those estimated from genotyped SNP, a comparison of the genomic inbreeding coefficients estimated either with the genotyped SNP or the SNP after imputation was made. Information on the presence or absence of genotypic information from sire, dam and maternal grandsire during the imputation was investigated. Genomic inbreeding coefficients estimated with genotyped SNP or SNP after imputation were consistent for F, Fhat3, Fgrm2 and FROH, when at least one of the parents was genotyped. Biased (mainly higher) genomic inbreeding coefficients of imputation SNP were observed in cows that were genotyped with MD SNP panels whose SNP were poorly represented in the selected imputation SNP data set and also did not have their parents genotyped compared with what expected based on actual genotype data. For cows genotyped with MD the estimators Fhat1, Fhat2 and Fgrm provided higher genomic inbreeding coefficients of imputation SNP even with both parents and the maternal grandsire genotyped. Overall, FROH was the most robust estimator, followed by F and Fhat3. Our findings suggest that SNP selection, parental genotyping and estimator should be considered for designing imputation strategies in dairy cattle for estimating genomic inbreeding with imputation SNP. For computing genomic inbreeding coefficients, it is recommendable to have at least one parent genotyped and use an ROH based estimator.

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.

Reproductive traits in Pura Raza Española mares manifest inbreeding depression from low levels of homozygosity.

Inbreeding depression is a genetic phenomenon associated with the loss of fitness and mean phenotypic performance due to mating between relatives. Historically, inbreeding coefficients have been estimated from pedigree information. However, the onset of genomic selection programs provides large datasets of individuals genotyped using SNP arrays, enabling more precise assessment of an individual's genomic-level inbreeding using genomic data. One of the traits most sensitive to issues stemming from increased inbreeding is reproduction. This is particularly important in equine, in which fertility is only moderate compared to other livestock species. To explore this further, we evaluated the effect of inbreeding on five reproductive traits (age at first foaling (AFF), average interval between foalings (AIF), total number of foalings (NF), productive life (PL) and reproductive efficiency (RE)) in Pura Raza Español mares using genomic data. Residual predicted phenotypes were obtained by purging these traits through the REML (wgResidual ) and ssGREML (gResidual ) approaches in reproductive data of 29,847 PRE mares using the BLUPF90+ program. Next, we used pedigree-based (Fped ) and ROH-based genomic (FROH ) inbreeding coefficients derived from 1018 animals genotyped with 61,271 SNPs to estimate the inbreeding depression (linear regression). Our results indicated significant levels of inbreeding depression for all reproductive traits, with the exception of the AIF trait when Fped was used. However, all traits were negatively affected by the increase in genomic inbreeding, and FROH was found to capture more inbreeding depression than Fped . Likewise, REML models (ssGREML) using genomic data for estimated predicted residual phenotypes resulted in higher variance explained by the model compared with the models not using genomics (REML). Finally, a segmented regression analysis was conducted to evaluate the effect of inbreeding depression, revealing that the levels of genealogical and genomic homozygosity do not manifest uniformly in reproductive traits. In contrast, the levels of inbreeding depression ranged from low to high as homozygosity increased. This analysis also showed that reproductive traits are very sensitive to inbreeding depression, even with relatively low levels of homozygosity.

Whole-Genome Resequencing Reveals Genetic Diversity and Wool Trait-Related Genes in Liangshan Semi-Fine-Wool Sheep

Understanding the genetic makeup of local sheep breeds is essential for their scientific conservation and sustainable utilization. The Liangshan semi-fine-wool sheep (LSS), a Chinese semi-fine-wool breed renowned for its soft wool, was analyzed using whole-genome sequencing data including 35 LSS, 84 sheep from other domestic breeds, and 20 Asiatic mouflons. We investigated the genetic composition of LSS by conducting analyses of the population structure, runs of homozygosity, genomic inbreeding coefficients, and selection signature. Our findings indicated that LSS shares greater genetic similarity with Border Leicester and Romney sheep than with Tibetan (TIB), Yunnan (YNS), and Chinese Merino sheep. Genomic analysis indicated low to moderate inbreeding coefficients, ranging from 0.014 to 0.154. In identifying selection signals across the LSS genome, we pinpointed 195 candidate regions housing 74 annotated genes (e.g., IRF2BP2, BVES, and ALOX5). We also found the overlaps between the candidate regions and several known quantitative trait loci related to wool traits, such as the wool staple length and wool fiber diameter. A selective sweep region, marked by the highest value of cross-population extended haplotype homozygosity, encompassed IRF2BP2—an influential candidate gene affecting fleece fiber traits. Furthermore, notable differences in genotype frequency at a mutation site (c.1051 + 46T &gt; C, Chr25: 6,784,190 bp) within IRF2BP2 were observed between LSS and TIB and YNS sheep (Fisher’s exact test, p &lt; 2.2 × 10−16). Taken together, these findings offer insights crucial for the conservation and breeding enhancement of LSS.

Genome-wide mining of diversity and evolutionary signatures revealed selective hotspots in Indian Sahiwal cattle

The Sahiwal cattle breed is the best indigenous dairy cattle breed, and it plays a pivotal role in the Indian dairy industry. This is due to its exceptional milk-producing potential, adaptability to local tropical conditions, and its resilience to ticks and diseases. The study aimed to identify selective sweeps and estimate intrapopulation genetic diversity parameters in Sahiwal cattle using ddRAD sequencing-based genotyping data from 82 individuals. After applying filtering criteria, 78,193 high-quality SNPs remained for further analysis. The population exhibited an average minor allele frequency of 0.221 ± 0.119. Genetic diversity metrics, including observed (0.597 ± 0.196) and expected heterozygosity (0.433 ± 0.096), nucleotide diversity (0.327 ± 0.114), the proportion of polymorphic SNPs (0.726), and allelic richness (1.323 ± 0.134), indicated ample genomic diversity within the breed. Furthermore, an effective population size of 74 was observed in the most recent generation. The overall mean linkage disequilibrium (r2) for pairwise SNPs was 0.269 ± 0.057. Moreover, a greater proportion of short Runs of Homozygosity (ROH) segments were observed suggesting that there may be low levels of recent inbreeding in this population. The genomic inbreeding coefficients, computed using different inbreeding estimates (FHOM, FUNI, FROH, and FGROM), ranged from −0.0289 to 0.0725. Subsequently, we found 146 regions undergoing selective sweeps using five distinct statistical tests: Tajima's D, CLR, |iHS|, |iHH12|, and ROH. These regions, located in non-overlapping 500 kb windows, were mapped and revealed various protein-coding genes associated with enhanced immune systems and disease resistance (IFNL3, IRF8, BLK), as well as production traits (NRXN1, PLCE1, GHR). Notably, we identified interleukin 2 (IL2) on Chr17: 35217075–35223276 as a gene linked to tick resistance and uncovered a cluster of genes (HSPA8, UBASH3B, ADAMTS18, CRTAM) associated with heat stress. These findings indicate the evolutionary impact of natural and artificial selection on the environmental adaptation of the Sahiwal cattle population.

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 250Kb, 120Kb/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.

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.

Distribution of Runs of Homozygosity and Their Relationship with Candidate Genes for Productivity in Kazakh Meat-Wool Sheep Breed.

Increasing the fertility of sheep remains one of the crucial issues of modern sheep breeding. The Kazakh meat-wool sheep is an excellent breed with high meat and wool productivity and well adapted to harsh conditions. Nowadays, runs of homozygosity (ROHs) are considered a suitable approach for studying the genetic characteristics of farm animals. The aims of the study were to analyze the distribution of ROHs, describe autozygosity, and detect genomic regions with high ROH islands. In this study, we genotyped a total of 281 Kazakh meat-wool sheep using the Illumina iScan® system (EquipNet, Canton, MA, USA) via Ovine SNP50 BeadChip array. As a results, a total of 15,069 ROHs were found in the three Kazakh meat-wool sheep populations. The mean number of ROH per animal across populations varied from 40.3 (POP1) to 42.2 (POP2) in the category 1+ Mb. Furthermore, the number of ROH per animal in ROH1-2 Mb were much higher than ROH2-4 Mb and ROH8-16 Mb in the three sheep populations. Most of individuals had small number of ROH>16 Mb. The highest and lowest genomic inbreeding coefficient values were observed in POP2 and POP3, respectively. The estimated FROH presented the impact that recent inbreeding has had in all sheep populations. Furthermore, a set of interesting candidate genes (BMP2, BMPR2, BMPRIB, CLOCK, KDM2B, TIAM1, TASP1, MYBPC1, MYOM1, and CACNA2D1), which are related to the productive traits, were found. Collectively, these findings will contribute to the breeding and conservation strategies of the Kazakh meat-wool sheep breed.

Genome-wide runs of homozygosity in Indian goat breeds reared under small holder production system

Small holder goat production system is often constrained by high inbreeding due to small herd size and uncontrolled breeding. Runs of homozygosity (ROH), the long stretches of homozygous genomic DNA is the most reliable measure of autozygosity. The objective of present study was to investigate the number and coverage of ROH in the genome and to estimate the genomic inbreeding coefficient (FROH) in two divergent Indian goat breeds viz; Attappady Black and Malabari goats farmed under small holder production system. The genome-wide single nucleotide polymorphism (SNP) marker data obtained by Illumina SNP50 goat BeadChip genotyping of 48 goats (24 each of Attappady Black and Malabari goats) sampled from farmers flock in their respective breeding tract was utilised for ROH detection. The ROH of > 1 Mb was defined as tracts of homozygous genotypes with at least 15 SNPs and less than five missing SNPs. The ROH detection was carried out with PLINK 1.9 version. The proportion of goats lacking ROH (per cent) was 66.67 in Attappady Black and 39.13 in Malabari goats. The mean number and length of ROH of 1 Mb or more was higher in Malabari goats (2.1304 ± 0.6369 and 6.0016 ± 0.5271 Mb) compared to Attappady Black goats (0.4583 ± 0.1343 and 4.0197 ± 0.3755 Mb). Mean ± S.E of sum of ROH coverage of genome per individual (SROH) was 1.8423 ± 0.05976 and 8.8236 ± 1.2785 Mb in Attappady Black and Malabari goats, respectively and differed among the two breeds significantly (p ≤ 0.05). Analysis of the distribution of ROH according to their length showed that majority of the detected ROH were < 8 Mb in length in both breeds with no ROH of > 8 Mb in Attappady Black goats. Mean ± S.E of genomic inbreeding coefficient (FROH≥1 Mb) in per cent was 0.519 ± 0.19 and 0.07480 ± 0.0243 in Malabari goats and Attappady Black goats, respectively. The pattern of ROH estimated in the present analysis indicates the inbreeding history of the goats farmed under small holder production system. Results of present study indicate that both breeds are less consanguineous and showed no evidence of recent inbreeding.

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

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

Genome-wide diversity and admixture of five indigenous cattle populations from the Tigray region of northern Ethiopia.

The Tigray region, where we found around eight per cent of the indigenous cattle population of Ethiopia, is considered as the historic centre of the country, with the ancient pre-Aksumite and Aksumite civilisations in contact with the civilisations of the Fertile Crescent and the Indian subcontinent. Here, we used whole genome sequencing data to characterise the genomic diversity, relatedness, and admixture of five cattle populations (Abergelle, Arado, Begait, Erob, and Raya) indigenous to the Tigray region of Ethiopia. We detected 28 to 29 million SNPs and 2.7 to 2.9 million indels in each population, of which 7% of SNPs and 34% of indels were novel. Functional annotation of the variants showed around 0.01% SNPs and 0.22%-0.27% indels in coding regions. Enrichment analysis of genes overlapping missense private SNPs revealed 20 significant GO terms and KEGG pathways that were shared by or specific to breeds. They included important genes associated with morphology (SCN4A, TAS1R2 and KCNG4), milk yield (GABRG1), meat quality (MMRN2, VWC2), feed efficiency (PCDH8 and SLC26A3), immune response (LAMC1, PCDH18, CELSR1, TLR6 and ITGA5), heat resistance (NPFFR1 and HTR7) and genes belonging to the olfactory gene family, which may be related to adaptation to harsh environments. Tigray indigenous cattle are very diverse. Their genome-wide average nucleotide diversity ranged from 0.0035 to 0.0036. The number of heterozygous SNPs was about 0.6-0.7 times higher than homozygous ones. The within-breed average number of ROHs ranged from 777.82 to 1000.45, with the average sum of the length of ROHs ranging from 122.01Mbp to 163.88Mbp. The genomic inbreeding coefficients differed among animals and breeds, reaching up to 10% in some Begait and Raya animals. Tigray indigenous cattle shared a common ancestry with Asian indicine (85.6%-88.7%) and African taurine (11.3%-14.1%) cattle, with very small, if any, European taurine introgression. This study identified high within-breed genetic diversity representing an opportunity for breeding improvement programs and, also, significant novel variants that could increase the number of known cattle variants, an important contribution to the knowledge of domestic cattle genetic diversity.