利用高密度SNP标记分析中国荷斯坦牛基因组近交

畜禽育种中传统上利用系谱信息评估群体近交程度?近年来随着高通量单核苷酸多态(single nucleotide polymorphism, SNP)检测成本降低,使利用基因组信息分析真实的基因组近交程度成为可能?本研究利用牛54 K SNP 芯片数据统计了北京地区2107头荷斯坦牛基因组上的长纯合片段(runs of homozygosity, ROH)的频率和分布,计算了2种基因组近交系数,即染色体上ROH的长度占基因组总长度的比例(Froh)及个体所有标记基因型中纯合子所占比例,即基因组纯合度(Fhom),进而分析了两种基因组近交系数之间的相关性以及基因组近交与系谱近交系数之间的相关性?结果表明,共检测到44 676个ROH片段,其长度主要分布在1~10 Mb之间?不同长度的ROH散布于个体基因组内,短ROH较长ROH更为常见?ROH在染色体上并非均匀分布,ROH频率最高的区域为10号染色体中部?两种基因组近交系数之间相关性很高(91%以上),但基因组近交与系谱近交之间的相关性较低(低于50%)?系谱完整性是影响基因组近交与系谱近交结果一致的重要因素,基因组近交系数能够反映个体真实的近交,本研究为评估群体近交水平提供了有力工具?     

Regions of Homozygosity in the Porcine Genome: Consequence of Demography and the Recombination Landscape

Inbreeding has long been recognized as a primary cause of fitness reduction in both wild and domesticated populations. Consanguineous matings cause inheritance of haplotypes that are identical by descent (IBD) and result in homozygous stretches along the genome of the offspring. Size and position of regions of homozygosity (ROHs) are expected to correlate with genomic features such as GC content and recombination rate, but also direction of selection. Thus, ROHs should be non-randomly distributed across the genome. Therefore, demographic history may not fully predict the effects of inbreeding. The porcine genome has a relatively heterogeneous distribution of recombination rate, making Sus scrofa an excellent model to study the influence of both recombination landscape and demography on genomic variation. This study utilizes next-generation sequencing data for the analysis of genomic ROH patterns, using a comparative sliding window approach. We present an in-depth study of genomic variation based on three different parameters: nucleotide diversity outside ROHs, the number of ROHs in the genome, and the average ROH size. We identified an abundance of ROHs in all genomes of multiple pigs from commercial breeds and wild populations from Eurasia. Size and number of ROHs are in agreement with known demography of the populations, with population bottlenecks highly increasing ROH occurrence. Nucleotide diversity outside ROHs is high in populations derived from a large ancient population, regardless of current population size. In addition, we show an unequal genomic ROH distribution, with strong correlations of ROH size and abundance with recombination rate and GC content. Global gene content does not correlate with ROH frequency, but some ROH hotspots do contain positive selected genes in commercial lines and wild populations. This study highlights the importance of the influence of demography and recombination on homozygosity in the genome to understand the effects of inbreeding.     

Runs of Homozygosity: Association with Coronary Artery Disease and Gene Expression in Monocytes and Macrophages

Runs of homozygosity (ROHs) are recognized signature of recessive inheritance. Contributions of ROHs to the genetic architecture of coronary artery disease and regulation of gene expression in cells relevant to atherosclerosis are not known. Our combined analysis of 24,320 individuals from 11 populations of white European ethnicity showed an association between coronary artery disease and both the count and the size of ROHs. Individuals with coronary artery disease had approximately 0.63 (95% CI: 0.4–0.8) excess of ROHs when compared to coronary-artery-disease-free control subjects (p = 1.49 × 10⁻⁹). The average total length of ROHs was approximately 1,046.92 (95% CI: 634.4–1,459.5) kb greater in individuals with coronary artery disease than control subjects (p = 6.61 × 10⁻⁷). None of the identified individual ROHs was associated with coronary artery disease after correction for multiple testing. However, in aggregate burden analysis, ROHs favoring increased risk of coronary artery disease were much more common than those showing the opposite direction of association with coronary artery disease (p = 2.69 × 10⁻³³). Individual ROHs showed significant associations with monocyte and macrophage expression of genes in their close proximity—subjects with several individual ROHs showed significant differences in the expression of 44 mRNAs in monocytes and 17 mRNAs in macrophages when compared to subjects without those ROHs. This study provides evidence for an excess of homozygosity in coronary artery disease in outbred populations and suggest the potential biological relevance of ROHs in cells of importance to the pathogenesis of atherosclerosis.     

Runs of homozygosity have utility in mammalian conservation and evolutionary studies

Runs of homozygosity (ROHs) arise due the transmission from parents to offspring of segments that are either identical by decent (IBD) or identical by state (IBS). The former is due to consanguineous matings whereas the latter is due to demographic processes. ROHs reduce individual nucleotide diversity (θ) as a function of homozygosity, and thus ROH distributions and θ are expected to vary among species because inbreeding levels, recombination rates, and demographic histories vary widely. To help interpret genetic diversity within and among species, we utilized genome sequence data from 78 mammalian species to compare θ and ROH burden (i.e., number and length of ROHs in the genome) among groups of mammals to assess genomic signatures of inbreeding. We compared θ and ROHs: (i) among threatened and non-threatened mammals to determine the significance of contemporary conservation status; (ii) among carnivorous and non-carnivorous mammals to determine the relevance of trophic effects; (iii) relative to body size because mutation rates generally vary with body mass; and (iv) across mammals from different latitudes to test for gradients in genomic diversity (e.g., due to effects of historic climatic regimes). Our results illustrate the considerable variance in genomic diversity across mammals, and that trophic level, body mass, and latitude have significant effects on θ and ROH burden. However, conservation status was not a reliable indicator of genomic diversity. We argue that genetic or genomic diversity should be an explicit component of conservation status, as such diversity is critical to the long-term sustainability of populations, and anticipate that ROHs will become more commonly used to estimate inbreeding in wild animals.     

Effect of reduced genomic representation on using runs of homozygosity for inbreeding characterization

Genomic measures of inbreeding based on identical-by-descent (IBD) segments are increasingly used to measure inbreeding and mostly estimated on SNP arrays and whole-genome sequencing (WGS) data. However, some softwares recurrently used for their estimation assume that genomic positions which have not been genotyped are nonvariant. This might be true for WGS data, but not for reduced genomic representations and can lead to spurious IBD segments estimation. In this project, we simulated the outputs of WGS, two SNP arrays of different sizes and RAD-sequencing for three populations with different sizes and histories. We compare the results of IBD segments estimation with two softwares: runs of homozygosity (ROHs) estimated with PLINK and homozygous-by-descent (HBD) segments estimated with RZooRoH. We demonstrate that to obtain meaningful estimates of inbreeding, RZooRoH requires a SNPs density 11 times smaller compared to PLINK: ranks of inbreeding coefficients were conserved among individuals above 22 SNPs/Mb for PLINK and 2 SNPs/Mb for RZooRoH. We also show that in populations with simple demographic histories, distribution of ROHs and HBD segments are correctly estimated with both SNP arrays and WGS. PLINK correctly estimated distribution of ROHs with SNP densities above 22 SNPs/Mb, while RZooRoH correctly estimated distribution of HBD segments with SNPs densities above 11 SNPs/Mb. However, in a population with a more complex demographic history, RZooRoH resulted in better distribution of IBD segments estimation compared to PLINK even with WGS data. Consequently, we advise researchers to use either methods relying on excess homozygosity averaged across SNPs or model-based HBD segments calling methods for inbreeding estimations.     

A genome-wide homozygosity association study identifies runs of homozygosity associated with rheumatoid arthritis in the human major histocompatibility complex

Rheumatoid arthritis (RA) is a chronic inflammatory disorder with a polygenic mode of inheritance. This study examined the hypothesis that runs of homozygosity (ROHs) play a recessive-acting role in the underlying RA genetic mechanism and identified RA-associated ROHs. Ours is the first genome-wide homozygosity association study for RA and characterized the ROH patterns associated with RA in the genomes of 2,000 RA patients and 3,000 normal controls of the Wellcome Trust Case Control Consortium. Genome scans consistently pinpointed two regions within the human major histocompatibility complex region containing RA-associated ROHs. The first region is from 32,451,664 bp to 32,846,093 bp (-log10(p)>22.6591). RA-susceptibility genes, such as HLA-DRB1, are contained in this region. The second region ranges from 32,933,485 bp to 33,585,118 bp (-log10(p)>8.3644) and contains other HLA-DPA1 and HLA-DPB1 genes. These two regions are physically close but are located in different blocks of linkage disequilibrium, and ～40% of the RA patients' genomes carry these ROHs in the two regions. By analyzing homozygote intensities, an ROH that is anchored by the single nucleotide polymorphism rs2027852 and flanked by HLA-DRB6 and HLA-DRB1 was found associated with increased risk for RA. The presence of this risky ROH provides a 62% accuracy to predict RA disease status. An independent genomic dataset from 868 RA patients and 1,194 control subjects of the North American Rheumatoid Arthritis Consortium successfully validated the results obtained using the Wellcome Trust Case Control Consortium data. In conclusion, this genome-wide homozygosity association study provides an alternative to allelic association mapping for the identification of recessive variants responsible for RA. The identified RA-associated ROHs uncover recessive components and missing heritability associated with RA and other autoimmune diseases.     

Identification of novel schizophrenia loci by homozygosity mapping using DNA microarray analysis

The recent development of high-resolution DNA microarrays, in which hundreds of thousands of single nucleotide polymorphisms (SNPs) are genotyped, enables the rapid identification of susceptibility genes for complex diseases. Clusters of these SNPs may show runs of homozygosity (ROHs) that can be analyzed for association with disease. An analysis of patients whose parents were first cousins enables the search for autozygous segments in their offspring. Here, using the Affymetrix® Genome-Wide Human SNP Array 5.0 to determine ROHs, we genotyped 9 individuals with schizophrenia (SCZ) whose parents were first cousins. We identified overlapping ROHs on chromosomes 1, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 16, 17, 19, 20, and 21 in at least 3 individuals. Only the locus on chromosome 5 has been reported previously. The ROHs on chromosome 5q23.3–q31.1 include the candidate genes histidine triad nucleotide binding protein 1 (HINT1) and acyl-CoA synthetase long-chain family member 6 (ACSL6). Other overlapping ROHs may contain novel rare recessive variants that affect SCZ specifically in our samples, given the highly heterozygous nature of SCZ. Analysis of patients whose parents are first cousins may provide new insights for the genetic analysis of psychiatric diseases.     

The existence of species rests on a metastable equilibrium between inbreeding and outbreeding. An essay on the close relationship between speciation, inbreeding and recessive mutations

Background

Speciation corresponds to the progressive establishment of reproductive barriers between groups of individuals derived from an ancestral stock. Since Darwin did not believe that reproductive barriers could be selected for, he proposed that most events of speciation would occur through a process of separation and divergence, and this point of view is still shared by most evolutionary biologists today.

Results

I do, however, contend that, if so much speciation occurs, the most likely explanation is that there must be conditions where reproductive barriers can be directly selected for. In other words, situations where it is advantageous for individuals to reproduce preferentially within a small group and reduce their breeding with the rest of the ancestral population. This leads me to propose a model whereby new species arise not by populations splitting into separate branches, but by small inbreeding groups "budding" from an ancestral stock. This would be driven by several advantages of inbreeding, and mainly by advantageous recessive phenotypes, which could only be retained in the context of inbreeding. Reproductive barriers would thus not arise as secondary consequences of divergent evolution in populations isolated from one another, but under the direct selective pressure of ancestral stocks. Many documented cases of speciation in natural populations appear to fit the model proposed, with more speciation occurring in populations with high inbreeding coefficients, and many recessive characters identified as central to the phenomenon of speciation, with these recessive mutations expected to be surrounded by patterns of limited genomic diversity.

Conclusions

Whilst adaptive evolution would correspond to gains of function that would, most of the time, be dominant, this type of speciation by budding would thus be driven by mutations resulting in the advantageous loss of certain functions since recessive mutations very often correspond to the inactivation of a gene. A very important further advantage of inbreeding is that it reduces the accumulation of recessive mutations in genomes. A consequence of the model proposed is that the existence of species would correspond to a metastable equilibrium between inbreeding and outbreeding, with excessive inbreeding promoting speciation, and excessive outbreeding resulting in irreversible accumulation of recessive mutations that could ultimately only lead to extinction.

Reviewer names

Eugene V. Koonin, Patrick Nosil (nominated by Dr Jerzy Jurka), Pierre Pontarotti     

Individual Variation in Inbreeding Depression: The Roles of Inbreeding History and Mutation

We use mutation-selection recursion models to evaluate the relative contributions of mutation and inbreeding history to variation among individuals in inbreeding depression and the ability of experiments to detect associations between individual inbreeding depression and mating system genotypes within populations. Poisson mutation to deleterious additive or recessive alleles generally produces far more variation among individuals in inbreeding depression than variation in history of inbreeding, regardless of selfing rate. Moreover, variation in inbreeding depression can be higher in a completely outcrossing or selfing population than in a mixed-mating population. In an initially random mating population, the spread of a dominant selfing modifier with no pleiotropic effects on male outcross success causes a measurable increase in inbreeding depression variation if its selfing rate is large and inbreeding depression is caused by recessive lethals. This increase is observable during a short period as the modifier spreads rapidly to fixation. If the modifier alters selfing rate only slightly, it fails to spread or causes no measurable increase in inbreeding depression variance. These results suggest that genetic associations between mating loci and inbreeding depression loci could be difficult to demonstrate within populations and observable only transiently during rapid evolution to a substantially new selfing rate.     

Evidence of inbreeding depression on human height

Stature is a classical and highly heritable complex trait, with 80%-90% of variation explained by genetic factors. In recent years, genome-wide association studies (GWAS) have successfully identified many common additive variants influencing human height; however, little attention has been given to the potential role of recessive genetic effects. Here, we investigated genome-wide recessive effects by an analysis of inbreeding depression on adult height in over 35,000 people from 21 different population samples. We found a highly significant inverse association between height and genome-wide homozygosity, equivalent to a height reduction of up to 3 cm in the offspring of first cousins compared with the offspring of unrelated individuals, an effect which remained after controlling for the effects of socio-economic status, an important confounder (χ(2) = 83.89, df = 1; p = 5.2 × 10(-20)). There was, however, a high degree of heterogeneity among populations: whereas the direction of the effect was consistent across most population samples, the effect size differed significantly among populations. It is likely that this reflects true biological heterogeneity: whether or not an effect can be observed will depend on both the variance in homozygosity in the population and the chance inheritance of individual recessive genotypes. These results predict that multiple, rare, recessive variants influence human height. Although this exploratory work focuses on height alone, the methodology developed is generally applicable to heritable quantitative traits (QT), paving the way for an investigation into inbreeding effects, and therefore genetic architecture, on a range of QT of biomedical importance.     

Educational Attainment Influences Levels of Homozygosity through Migration and Assortative Mating

Individuals with a higher education are more likely to migrate, increasing the chance of meeting a spouse with a different ancestral background. In this context, the presence of strong educational assortment can result in greater ancestry differences within more educated spouse pairs, while less educated individuals are more likely to mate with someone with whom they share more ancestry. We examined the association between educational attainment and F _roh (= the proportion of the genome consisting of runs of homozygosity [ROHs]) in ~2,000 subjects of Dutch ancestry. The subjects’ own educational attainment showed a nominally significant negative association with F _roh (p = .045), while the contribution of parental education to offspring F _roh was highly significant (father: p < 10^-5; mother: p = 9×10^-5), with more educated parents having offspring with fewer ROHs. This association was significantly and fully mediated by the physical distance between parental birthplaces (paternal education: p _mediation = 2.4 × 10^-4; maternal education: p _mediation = 2.3 × 10^-4), which itself was also significantly associated with F _roh (p = 9 × 10^-5). Ancestry-informative principal components from the offspring showed a significantly decreasing association with geography as parental education increased, consistent with the significantly higher migration rates among more educated parents. Parental education also showed a high spouse correlation (Spearman’s ρ = .66, p = 3 × 10^-262). We show that less educated parents are less likely to mate with the more mobile parents with a higher education, creating systematic differences in homozygosity due to ancestry differences not directly captured by ancestry-informative principal components (PCs). Understanding how behaviors influence the genomic structure of a population is highly valuable for studies on the genetic etiology of behavioral, cognitive, and social traits.     

Age-specific differences in blood pressure among inbred and non-inbred north Indian children

Genetic and inbreeding influences on systolic blood pressure (SBP), diastolic blood pressure (DBP) and mean arterial blood pressure (MBP) were examined among 3015 children (1527 males and 1488 females) from the Aligarh district, Uttar Pradesh in north India. The subjects included offspring of first cousins, first cousins once removed, second cousins and unrelated spouses from the same population. The measurements of the inbred children were compared with those of their non-inbred relatives in at least 80% of the cases (matched controls). Two unique findings emerge from this study. First a consistent increase in mean values of SBP, DBP and MBP with increasing inbreeding coefficients have been observed among all age groups, including both the sexes. The results suggest that the hypothesis for a recessive gene or genes could be held responsible for higher BP. Secondly, the effects of inbreeding on mean blood pressure among children and adults may not necessarily be in the same direction. It can be said, therefore, that studies on inbreeding effects using matched controls may provide more direct information regarding the genetics of blood pressure, which has been considerably underestimated in earlier studies.     

Copy-Number Disorders Are a Common Cause of Congenital Kidney Malformations

We examined the burden of large, rare, copy-number variants (CNVs) in 192 individuals with renal hypodysplasia (RHD) and replicated findings in 330 RHD cases from two independent cohorts. CNV distribution was significantly skewed toward larger gene-disrupting events in RHD cases compared to 4,733 ethnicity-matched controls (p = 4.8 × 10⁻¹¹). This excess was attributable to known and novel (i.e., not present in any database or in the literature) genomic disorders. All together, 55/522 (10.5%) RHD cases harbored 34 distinct known genomic disorders, which were detected in only 0.2% of 13,839 population controls (p = 1.2 × 10⁻⁵⁸). Another 32 (6.1%) RHD cases harbored large gene-disrupting CNVs that were absent from or extremely rare in the 13,839 population controls, identifying 38 potential novel or rare genomic disorders for this trait. Deletions at the HNF1B locus and the DiGeorge/velocardiofacial locus were most frequent. However, the majority of disorders were detected in a single individual. Genomic disorders were detected in 22.5% of individuals with multiple malformations and 14.5% of individuals with isolated urinary-tract defects; 14 individuals harbored two or more diagnostic or rare CNVs. Strikingly, the majority of the known CNV disorders detected in the RHD cohort have previous associations with developmental delay or neuropsychiatric diseases. Up to 16.6% of individuals with kidney malformations had a molecular diagnosis attributable to a copy-number disorder, suggesting kidney malformations as a sentinel manifestation of pathogenic genomic imbalances. A search for pathogenic CNVs should be considered in this population for the diagnosis of their specific genomic disorders and for the evaluation of the potential for developmental delay.     

Genome-wide identification of runs of homozygosity islands and associated genes in local dairy cattle breeds

Runs of homozygosity (ROH) are widely used as predictors of whole-genome inbreeding levels in cattle. They identify regions that have an unfavorable effect on a phenotype when homozygous, but also identify the genes associated with traits of economic interest present in these regions. Here, the distribution of ROH islands and enriched genes within these regions in four dairy cattle breeds were investigated. Cinisara (71), Modicana (72), Reggiana (168) and Italian Holstein (96) individuals were genotyped using the 50K v2 Illumina BeadChip. The genomic regions most commonly associated with ROHs were identified by selecting the top 1% of the single nucleotide polymorphisms (SNPs) most commonly observed in the ROH of each breed. In total, 11 genomic regions were identified in Cinisara and Italian Holstein, and eight in Modicana and Reggiana, indicating an increased ROH frequency level. Generally, ROH islands differed between breeds. The most homozygous region (>45% of individuals with ROH) was found in Modicana on chromosome 6 within a quantitative trail locus affecting milk fat and protein concentrations. We identified between 126 and 347 genes within ROH islands, which are involved in multiple signaling and signal transduction pathways in a wide variety of biological processes. The gene ontology enrichment provided information on possible molecular functions, biological processes and cellular components under selection related to milk production, reproduction, immune response and resistance/susceptibility to infection and diseases. Thus, scanning the genome for ROH could be an alternative strategy to detect genomic regions and genes related to important economic traits.     

Regions of homozygosity and their impact on complex diseases and traits

Regions of homozygosity (ROHs) are more abundant in the human genome than previously thought. These regions are without heterozygosity, i.e. all the genetic variations within the regions have two identical alleles. At present there are no standardized criteria for defining the ROHs resulting in the different studies using their own criteria in the analysis of homozygosity. Compared to the era of genotyping microsatellite markers, the advent of high-density single nucleotide polymorphism genotyping arrays has provided an unparalleled opportunity to comprehensively detect these regions in the whole genome in different populations. Several studies have identified ROHs which were associated with complex phenotypes such as schizophrenia, late-onset of Alzheimer’s disease and height. Collectively, these studies have conclusively shown the abundance of ROHs larger than 1 Mb in outbred populations. The homozygosity association approach holds great promise in identifying genetic susceptibility loci harboring recessive variants for complex diseases and traits.     

Inbreeding depression for various traits in two cultured populations of the American bay scallop, Argopecten irradians irradians Lamarck (1819) introduced into China

This paper examines the effect of inbreeding level of population on the magnitude of inbreeding depression expressed by comparing them between two cultured populations (A and B) in the hermaphroditic animal of the bay scallop Argopecten irradians irradians. Population A is expected to have less genetic variations and higher inbreeding level due to longer cultured history (20 generations) and less “ancestral” individuals (26 individuals) than population B due to shorter cultured history (4 generations) and more “ancestral” individuals (406 individuals). Two groups within each population were produced, one using self-fertilization and one using mass-mating within the same population. Selfed offspring (AS and BS) from two populations both had lower fitness components than their mass-mated counterparts (AM and BM) and exhibited inbreeding depression for all examined traits, e.g. lower hatching, less viability and slower growth, indicating that inbreeding depression is a common feature in this animal. Fitness components in all traits of offspring from population A significantly differed those from population B and the magnitude of inbreeding depression for all traits in population A with higher inbreeding level was significantly smaller than that in population B with lower inbreeding level, indicating that both fitness components and magnitude of inbreeding depression were significantly affected by inbreeding level of populations and genetic load harbored in population A may be partially purged through inbreeding. Moreover, the magnitude of inbreeding depression in the two populations both varied among traits and life history stages. The present results support the partial-dominance hypothesis of inbreeding depression.     

Blood cis-eQTL Analysis Fails to Identify Novel Association Signals among Sub-Threshold Candidates from Genome-Wide Association Studies in Restless Legs Syndrome

Restless legs syndrome (RLS) is a common neurologic disorder characterized by nightly dysesthesias affecting the legs primarily during periods of rest and relieved by movement. RLS is a complex genetic disease and susceptibility factors in six genomic regions have been identified by means of genome-wide association studies (GWAS). For some complex genetic traits, expression quantitative trait loci (eQTLs) are enriched among trait-associated single nucleotide polymorphisms (SNPs). With the aim of identifying new genetic susceptibility factors for RLS, we assessed the 332 best-associated SNPs from the genome-wide phase of the to date largest RLS GWAS for cis-eQTL effects in peripheral blood from individuals of European descent. In 740 individuals belonging to the KORA general population cohort, 52 cis-eQTLs with p_nominal<10⁻³ were identified, while in 976 individuals belonging to the SHIP-TREND general population study 53 cis-eQTLs with p_nominal<10⁻³ were present. 23 of these cis-eQTLs overlapped between the two cohorts. Subsequently, the twelve of the 23 cis-eQTL SNPs, which were not located at an already published RLS-associated locus, were tested for association in 2449 RLS cases and 1462 controls. The top SNP, located in the DET1 gene, was nominally significant (p<0.05) but did not withstand correction for multiple testing (p = 0.42). Although a similar approach has been used successfully with regard to other complex diseases, we were unable to identify new genetic susceptibility factor for RLS by adding this novel level of functional assessment to RLS GWAS data.     

The pattern of runs of homozygosity and genomic inbreeding in world-wide sheep populations

Genome-wide pattern of runs of homozygosity (ROH) across ovine genome can provide a useful resource for studying diversity and demography history in sheep. We analyzed 50 k SNPs chip data of 2536 animals to identify pattern, distribution and level of ROHs in 68 global sheep populations. A total of 60,301 ROHs were detected in all breeds. The majority of the detected ROHs were <16 Mb and the average total number of ROHs per individual was 23.8 ± 13.8. The ROHs greater than 1 Mb covered on average 8.2% of the sheep autosomes, 1% of which was related to the ROHs with 1–4 Mb of length. The mean sum of ROH length in two-thirds of the populations was less than 250 Mb ranging from 21.7 to near 570 Mb. The level of genomic inbreeding was relatively low. The average of the inbreeding coefficients based on ROH (F_ROH) was 0.09 ± 0.05. It was rising in a stepwise manner with distance from Southwest Asia and maximum values were detected in North European breeds. A total of 465 ROH hotspots were detected in 25 different autosomes which partially surrounding 257 Refseq genes across the genome. Most of the detected genes were related to growth, body weight, meat production and quality, wool production and pigmentation. In conclusion, our analysis showed that the sheep genome, compared with other livestock species such as cattle and pig, displays low levels of homozygosity and appropriate genetic diversity for selection response and genetic merit gain.