Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens.

BACKGROUND: Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding. RESULTS: To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively. CONCLUSION: Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.

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.

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

BACKGROUND: The 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. RESULTS: HO 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). CONCLUSIONS: Evaluations 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.

Genomic diversity and demographic history of the endangered Sichuan hill-partridge (Arborophila rufipectus).

Species conservation can be improved by knowledge of genetic diversity and demographic history. The Sichuan hill-partridge (Arborophila rufipectus, SP) is an endangered species endemic to the mountains in southwestern China. However, little is known about this species' genomic variation and demographic history. Here, we present a comprehensive whole-genome analysis of six SP individuals from the Laojunshan National Nature Reserve in Sichuan Province, China. We observe a relatively high genetic diversity and low level of recent inbreeding in the studied SP individuals. This suggests that the current population carries genetic variability that may benefit the long-term survival of this species, and that the present population may be larger than currently recognized. Analyses of demographic history showed that fluctuations in the effective population size of SP are inconsistent with changes of the historical climate. Strikingly, evidence from demographic modeling suggests SPs population decreased dramatically 15,100 years ago after the Last Glacial Maximum, possibly due to refugial isolation and later human interference. These results provide the first detailed and comprehensive genomic insights into genetic diversity, genomic inbreeding levels, and demographic history of the Sichuan hill-partridge, which are crucial for the conservation and management of this endangered species.

Population genomic analysis provides evidence of the past success and future potential of South China tiger captive conservation.

BACKGROUND: Among six extant tiger subspecies, the South China tiger (Panthera tigris amoyensis) once was widely distributed but is now the rarest one and extinct in the wild. All living South China tigers are descendants of only two male and four female wild-caught tigers and they survive solely in zoos after 60 years of effective conservation efforts. Inbreeding depression and hybridization with other tiger subspecies were believed to have occurred within the small, captive South China tiger population. It is therefore urgently needed to examine the genomic landscape of existing genetic variation among the South China tigers. RESULTS: In this study, we assembled a high-quality chromosome-level genome using long-read sequences and re-sequenced 29 high-depth genomes of the South China tigers. By combining and comparing our data with the other 40 genomes of six tiger subspecies, we identified two significantly differentiated genomic lineages among the South China tigers, which harbored some rare genetic variants introgressed from other tiger subspecies and thus maintained a moderate genetic diversity. We noticed that the South China tiger had higher FROH values for longer runs of homozygosity (ROH > 1 Mb), an indication of recent inbreeding/founder events. We also observed that the South China tiger had the least frequent homozygous genotypes of both high- and moderate-impact deleterious mutations, and lower mutation loads than both Amur and Sumatran tigers. Altogether, our analyses indicated an effective genetic purging of deleterious mutations in homozygous states from the South China tiger, following its population contraction with a controlled increase in inbreeding based on its pedigree records. CONCLUSIONS: The identification of two unique founder/genomic lineages coupled with active genetic purging of deleterious mutations in homozygous states and the genomic resources generated in our study pave the way for a genomics-informed conservation, following the real-time monitoring and rational exchange of reproductive South China tigers among zoos.

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.

Comparing pedigree and genomic inbreeding coefficients, and inbreeding depression of reproductive traits in Japanese Black cattle.

BACKGROUND: Pedigree-based inbreeding coefficients have been generally included in statistical models for genetic evaluation of Japanese Black cattle. The use of genomic data is expected to provide precise assessment of inbreeding level and depression. Recently, many measures have been used for genome-based inbreeding coefficients; however, with no consensus on which is the most appropriate. Therefore, we compared the pedigree- ([Formula: see text]) and multiple genome-based inbreeding coefficients, which were calculated from the genomic relationship matrix with observed allele frequencies ([Formula: see text]), correlation between uniting gametes ([Formula: see text]), the observed vs expected number of homozygous genotypes ([Formula: see text]), runs of homozygosity (ROH) segments ([Formula: see text]) and heterozygosity by descent segments ([Formula: see text]). We quantified inbreeding depression from estimating regression coefficients of inbreeding coefficients on three reproductive traits: age at first calving (AFC), calving difficulty (CD) and gestation length (GL) in Japanese Black cattle. RESULTS: The highest correlations with [Formula: see text] were for [Formula: see text] (0.86) and [Formula: see text] (0.85) whereas [Formula: see text] and [Formula: see text] provided weak correlations with [Formula: see text], with range 0.33-0.55. Except for [Formula: see text] and [Formula: see text], there were strong correlations among genome-based inbreeding coefficients ([Formula: see text] 0.94). The estimates of regression coefficients of inbreeding depression for [Formula: see text] was 2.1 for AFC, 0.63 for CD and -1.21 for GL, respectively, but [Formula: see text] had no significant effects on all traits. Genome-based inbreeding coefficients provided larger effects on all reproductive traits than [Formula: see text]. In particular, for CD, all estimated regression coefficients for genome-based inbreeding coefficients were significant, and for GL, that for [Formula: see text] had a significant.. Although there were no significant effects when using overall genome-level inbreeding coefficients for AFC and GL, [Formula: see text] provided significant effects at chromosomal level in four chromosomes for AFC, three chromosomes for CD, and two chromosomes for GL. In addition, similar results were obtained for [Formula: see text]. CONCLUSIONS: Genome-based inbreeding coefficients can capture more phenotypic variation than [Formula: see text]. In particular, [Formula: see text] and [Formula: see text] can be considered good estimators for quantifying inbreeding level and identifying inbreeding depression at the chromosome level. These findings might improve the quantification of inbreeding and breeding programs using genome-based inbreeding coefficients.

Genomic signatures of inbreeding depression for a threatened Aotearoa New Zealand passerine.

For small and isolated populations, the increased chance of mating between related individuals can result in a substantial reduction in individual and population fitness. Despite the increasing availability of genomic data to measure inbreeding accurately across the genome, inbreeding depression studies for threatened species are still scarce due to the difficulty of measuring fitness in the wild. Here, we investigate inbreeding and inbreeding depression for the extensively monitored Tiritiri Matangi island population of a threatened Aotearoa New Zealand passerine, the hihi (Notiomystis cincta). First, using a custom 45 k single nucleotide polymorphism (SNP) array, we explore genomic inbreeding patterns by inferring homozygous segments across the genome. Although all individuals have similar levels of ancient inbreeding, highly inbred individuals are affected by recent inbreeding, which can probably be explained by bottleneck effects such as habitat loss after European arrival and their translocation to the island in the 1990s. Second, we investigate genomic inbreeding effects on fitness, measured as lifetime reproductive success, and its three components, juvenile survival, adult annual survival and annual reproductive success, in 363 hihi. We find that global inbreeding significantly affects juvenile survival but none of the remaining fitness traits. Finally, we employ a genome-wide association approach to test the locus-specific effects of inbreeding on fitness, and identify 13 SNPs significantly associated with lifetime reproductive success. Our findings suggest that inbreeding depression does impact hihi, but at different genomic scales for different traits, and that purging has therefore failed to remove all variants with deleterious effects from this population of conservation concern.

Estimates of genomic inbreeding and identification of candidate regions in Beijing-You chicken populations.

Runs of homozygosity (ROHs) has become an effective method for analysing inbreeding in livestock populations. Moreover, ROHs is well-suited to detect signatures of selection via ROH islands. This study aimed to investigate the occurrence and distribution of ROHs, compare the genomic inbreeding coefficients and identify the genomic regions with high ROH frequencies in different Beijing-You chicken (BY) populations, including a random conservation population (BY_R), a pedigree conservation population (BY_P), and a commercial population obtained from the market (BY_C). Among them, BY_R in 2010 and 2019 were BY_R1 and BY_R2 respectively. A total of 27 916 ROHs were identified. The average number of ROHs per individual across the three BY populations ranged from 213 (BY_P) to 161 (BY_C), and the average length of ROHs ranged from 0.432 Mb (BY_R2) to 0.451 Mb (BY_P). The highest inbreeding coefficient calculated based on ROHs (FROH ) was 0.1 in BY_P, whereas the lowest FROH was 0.0743 in BY_C. In addition, the inbreeding coefficient of BY_R2 (FROH  = 0.0798) was higher than that of BY_R1 (FROH  = 0.0579). Furthermore, the highest proportion of long ROH fragments (>4 Mb) was observed in BY_P and BY_C. This study showed the top 10 ROH islands of each population, and these ROH islands harboured 53 genes, some of which were related to limb development, body size and immune response. These findings contribute to the understanding of genetic diversity and population demography, and might help improve breeding and conservation strategies for BY populations.

A Comprehensive Genomic Analysis of Chinese Indigenous Ningxiang Pigs: Genomic Breed Compositions, Runs of Homozygosity, and Beyond.

Ningxiang pigs are a renowned indigenous pig breed in China, known for their meat quality, disease resistance, and environmental adaptability. In recent decades, consumer demand for meats from indigenous breeds has grown significantly, fueling the selection and crossbreeding of Ningxiang pigs (NXP). The latter has raised concerns about the conservation and sustainable use of Ningxiang pigs as an important genetic resource. To address these concerns, we conducted a comprehensive genomic study using 2242 geographically identified Ningxiang pigs. The estimated genomic breed composition (GBC) suggested 2077 pigs as purebred Ningxiang pigs based on a ≥94% NXP-GBC cut-off. The remaining 165 pigs were claimed to be crosses, including those between Duroc and Ningxiang pigs and between Ningxiang and Shaziling pigs, and non-Ningxiang pigs. Runs of homozygosity (ROH) were identified in the 2077 purebred Ningxiang pigs. The number and length of ROH varied between individuals, with an average of 32.14 ROH per animal and an average total length of 202.4 Mb per animal. Short ROH (1-5 Mb) was the most abundant, representing 66.5% of all ROH and 32.6% of total ROH coverage. The genomic inbreeding estimate was low (0.089) in purebred Ningxiang pigs compared to imported western pig breeds. Nine ROH islands were identified, pinpointing candidate genes and QTLs associated with economic traits of interest, such as reproduction, carcass and growth traits, lipid metabolism, and fat deposition. Further investigation of these ROH islands and candidate genes is anticipated to better understand the genomics of Ningxiang pigs.

Genome-wide run of homozygosity analysis reveals candidate genomic regions associated with environmental adaptations of Tibetan native chickens.

BACKGROUND: In Tibet, the two most important breeds are Tibetan chicken and Lhasa white chicken, and the duo exhibit specific adaptations to the high altitude thereby supplying proteins for humans living in the plateau. These breeds are partly included in the conservation plans because they represent important chicken genetic resources. However, the genetic diversity of these chickens is rarely investigated. Based on whole-genome sequencing data of 113 chickens from 4 populations of Tibetan chicken including Shigatse (SH), Nyemo (NM), Dagze (DZ) and Nyingchi (LZ), as well as Lhasa white (LW) chicken breed, we investigated the genetic diversity of these chicken breeds by genetic differentiation, run of homozygosity (ROH), genomic inbreeding and selection signature analyses. RESULTS: Our results revealed high genetic diversity across the five chicken populations. The linkage disequilibrium decay was highest in LZ, while subtle genetic differentiation was found between LZ and other populations (Fst ranging from 0.05 to 0.10). Furthermore, the highest ROH-based inbreeding estimate (FROH) of 0.11 was observed in LZ. In other populations, the FROH ranged from 0.04 to 0.06. In total, 74, 111, 62, 42 and 54 ROH islands containing SNPs ranked top 1% for concurrency were identified in SH, NM, DZ, LZ and LW, respectively. Genes common to the ROH islands in the five populations included BDNF, CCDC34, LGR4, LIN7C, GLS, LOC101747789, MYO1B, STAT1 and STAT4. This suggested their essential roles in adaptation of the chickens. We also identified a common candidate genomic region harboring AMY2A, NTNG1 and VAV3 genes in all populations. These genes had been implicated in digestion, neurite growth and high-altitude adaptation. CONCLUSIONS: High genetic diversity is observed in Tibetan native chickens. Inbreeding is more intense in the Nyingchi population which is also genetically distant from other chicken populations. Candidate genes in ROH islands are likely to be the drivers of adaptation to high altitude exhibited by the five Tibetan native chicken populations. Our findings contribute to the understanding of genetic diversity offer valuable insights for the genetic mechanism of adaptation, and provide veritable tools that can help in the design and implementation of breeding and conservation strategies for Tibetan native chickens.

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

Genomic data of different resolutions reveal consistent inbreeding estimates but contrasting homozygosity landscapes for the threatened Aotearoa New Zealand hihi.

Inbreeding can lead to a loss of heterozygosity in a population and when combined with genetic drift may reduce the adaptive potential of a species. However, there is uncertainty about whether resequencing data can provide accurate and consistent inbreeding estimates. Here, we performed an in-depth inbreeding analysis for hihi (Notiomystis cincta), an endemic and nationally vulnerable passerine bird of Aotearoa New Zealand. We first focused on subsampling variants from a reference genome male, and found that low-density data sets tend to miss runs of homozygosity (ROH) in some places and overestimate ROH length in others, resulting in contrasting homozygosity landscapes. Low-coverage resequencing and 50 K SNP array densities can yield comparable inbreeding results to high-coverage resequencing approaches, but the results for all data sets are highly dependent on the software settings employed. Second, we extended our analysis to 10 hihi where low-coverage whole genome resequencing, RAD-seq and SNP array genotypes are available. We inferred ROH and individual inbreeding to evaluate the relative effects of sequencing depth versus SNP density on estimating inbreeding coefficients and found that high rates of missingness downwardly bias both the number and length of ROH. In summary, when using genomic data to evaluate inbreeding, studies must consider that ROH estimates are heavily dependent on analysis parameters, data set density and individual sequencing depth.

Effect of recent and ancient inbreeding on production and fertility traits in Canadian Holsteins.

BACKGROUND: Phenotypic performances of livestock animals decline with increasing levels of inbreeding, however, the noticeable decline known as inbreeding depression, may not be due only to the total level of inbreeding, but rather could be distinctly associated with more recent or more ancient inbreeding. Therefore, splitting inbreeding into different age classes could help in assessing detrimental effects of different ages of inbreeding. Hence, this study sought to investigate the effect of recent and ancient inbreeding on production and fertility traits in Canadian Holstein cattle with both pedigree and genomic records. Furthermore, inbreeding coefficients were estimated using traditional pedigree measure (FPED) and genomic measures using segment based (FROH) and marker-by-marker (FGRM) based approaches. RESULTS: Inbreeding depression was found for all production and most fertility traits, for example, every 1% increase in FPED, FROH and FGRM was observed to cause a - 44.71, - 40.48 and - 48.72 kg reduction in 305-day milk yield (MY), respectively. Similarly, an extension in first service to conception (FSTC) of 0.29, 0.24 and 0.31 day in heifers was found for every 1% increase in FPED, FROH and FGRM, respectively. Fertility traits that did not show significant depression were observed to move in an unfavorable direction over time. Splitting both pedigree and genomic inbreeding into age classes resulted in recent age classes showing more detrimental inbreeding effects, while more distant age classes caused more favorable effects. For example, a - 1.56 kg loss in 305-day protein yield (PY) was observed for every 1% increase in the most recent pedigree age class, whereas a 1.33 kg gain was found per 1% increase in the most distant pedigree age class. CONCLUSIONS: Inbreeding depression was observed for production and fertility traits. In general, recent inbreeding had unfavorable effects, while ancestral inbreeding had favorable effects. Given that more negative effects were estimated from recent inbreeding when compared to ancient inbreeding suggests that recent inbreeding should be the primary focus of selection programs. Also, further work to identify specific recent homozygous regions negatively associated with phenotypic traits could be investigated.

The evaluation of genomic homozygosity for Xinjiang inbred population by SNP panels.

Xinjiang inbred cattle is a population which has been highly inbred for 45 years. However, the breed origin of this population cannot be traced back due to the lack of original records. To demonstrate the genetic background of Xinjiang inbred cattle, we analysed the worldwide genomic information of 16 cattle breeds using principal components analysis, and Admixture method. Furthermore, the shared SNP markers of Xinjiang inbred cattle, local Kazakh cattle, Holstein cattle, and Xinjiang Brown cattle were extracted to calculate population genetic parameters and genomic inbreeding indicators in order to evaluate the magnitude of inbreeding in each population. We also evaluated the relationship between inbreeding indicators and body size in the Xinjiang inbred population. Finally, the high frequency runs of homozygosity (ROH) regions for Xinjiang inbred cattle and local Kazakh population were selected for genes and QTL annotations. These results demonstrate that the ancestry proportions of inbreeding breed are similar to those of Kazakh cattle. The genomic homozygosity of Xinjiang inbred cattle is significantly higher than other populations; the inbreeding depression is observed in body size to a certain extent because body size decreased when corresponding homozygosity increased. Totally, six basic bio-pathways and 32 QTL regions that related to bovine economical traits were annotated. Our results provide the insights into breeding strategies, future protection, and utilization plan design for this special genetic material-Xinjiang inbred cattle.

Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits.

BACKGROUND: The aim of this study was to assess genome-wide autozygosity in a Nellore cattle population and to characterize ROH patterns and autozygosity islands that may have occurred due to selection within its lineages. It attempts also to compare estimates of inbreeding calculated from ROH (FROH), genomic relationship matrix (FGRM), and pedigree-based coefficient (FPED). RESULTS: The average number of ROH per animal was 55.15 ± 13.01 with an average size of 3.24 Mb. The Nellore genome is composed mostly by a high number of shorter segments accounting for 78% of all ROH, although the proportion of the genome covered by them was relatively small. The genome autozygosity proportion indicates moderate to high inbreeding levels for classical standards, with an average value of 7.15% (178.70 Mb). The average of FPED and FROH, and their correlations (- 0.05 to 0.26) were low. Estimates of correlation between FGRM-FPED was zero, while the correlation (- 0.01 to - 0.07) between FGRM-FROH decreased as a function of ROH length, except for FROH > 8Mb (- 0.03). Overall, inbreeding coefficients were not high for the genotyped animals. Autozygosity islands were evident across the genome (n = 62) and their genomic location did not largely differ within lineages. Enriched terms (p < 0.01) associated with defense response to bacteria (GO:0042742), immune complex reaction (GO:0045647), pregnancy-associated glycoproteins genes (GO:0030163), and organism growth (GO:0040014) were described within the autozygotic islands. CONCLUSIONS: Low FPED-FROH correlation estimates indicate that FPED is not the most suitable method for capturing ancient inbreeding when the pedigree does not extend back many generations and FROH should be used instead. Enriched terms (p < 0.01) suggest a strong selection for immune response. Non-overlapping islands within the lineages greatly explain the mechanism underlying selection for functionally important traits in Nellore cattle.

Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds.

The availability of dense single nucleotide polymorphism (SNP) assays allows for the determination of autozygous segments based on runs of consecutive homozygous genotypes (ROH). The aim of the present study was to investigate the occurrence and distribution of ROH in 21 Italian sheep breeds using medium-density SNP genotypes in order to characterize autozygosity and identify genomic regions that frequently appeared in ROH within individuals, namely ROH islands. After filtering, the final number of animals and SNPs retained for analyses were 502 and 46 277 respectively. A total of 12 302 ROH were identified. The mean number of ROH per breed ranged from 10.58 (Comisana) to 44.54 (Valle del Belice). The average length of ROH across breeds was 4.55 Mb and ranged from 3.85 Mb (Biellese) to 5.51 Mb (Leccese). Valle del Belice showed the highest value of inbreeding on the basis of ROH (FROH  = 0.099), whereas Comisana showed the lowest (FROH  = 0.016), and high standard deviation values revealed high variability in autozygosity levels within each breed. Differences also existed in the length of ROH. Analysis of the distribution of ROH according to their size showed that, for all breeds, the majority of the detected ROH were <10 Mb in length, with a few long ROH >25 Mb. The levels of ROH that we estimated here reflect the inbreeding history of the investigated sheep breeds. These results also highlight that ancient and recent inbreeding have had an impact on the genome of the Italian sheep breeds and suggest that several animals have experienced recent autozygosity events. Comisana and Bergamasca appeared as the less consanguineous breeds, whereas Barbaresca, Leccese and Valle del Belice showed ROH patterns typically produced by recent inbreeding. Moreover, within the genomic regions most commonly associated with ROH, several candidate genes were detected.

Application of Genomic Estimation Methods of Inbreeding and Population Structure in an Arabian Horse Herd.

Horse breeders rely heavily on pedigrees for identification of ancestry in breeding stock. Inaccurate pedigrees may erroneously assign individuals to false lineages or breed memberships resulting in wrong estimates of inbreeding and coancestry. Moreover, discrepancies in pedigree records can lead breeders seeking to limit inbreeding into making misguided breeding decisions. Genome-wide SNPs provide a quantitative tool to aid in the resolution of lineage assignments and the calculation of genomic measures of relatedness. The aim of this project was to pilot a comparison between pedigree and genomic relatedness and inbreeding measures in a herd of 36 pedigreed Egyptian Arabian horses genotyped using the Equine SNP70 platform (Geneseek, Inc.). Moreover, we sought to estimate the minimum number of markers sufficient for genomic inbreeding calculations. Pedigree inbreeding values were moderately correlated with genomic inbreeding values (r = 0.406), whereas genomic relationships and pedigree relationships have a high correlation (r = 0.77). Although first degree relationships were successfully reconstructed, more distant relationships were difficult to resolve. Multi-dimensional scaling and clustering analysis agreed with within-herd pedigree information. In comparing the herd to a reference sample of United States, Polish, and Egyptian Arabian horses, the herd's historically recorded Egyptian lineage was successfully recovered. We conclude that genomic estimates of inbreeding and relationships are superior to their pedigree counterparts. They can be thus utilized in conservation of valuable lines of livestock, and in breeds at risk for loss of genomic diversity. We also postulate a minimum of 2000 markers in linkage equilibrium to be used for inbreeding estimation.

Using a family-based structure to detect the effects of genomic inbreeding on embryo viability in Holstein cattle.

Recent evidence has suggested that some of the decline in reproductive ability in dairy cattle has been caused by embryonic death. The current study compared expected genomic inbreeding from sire-dam mating pairs to genomic inbreeding from live progeny in an attempt to determine how embryonic inbreeding may affect fertility. A total of 11,484 Holstein cattle with 43,485 SNP markers and pedigree information were available for analysis. A total of 412 sire-dam-progeny trios in which all animals had reliable genotypes were discovered. After removal of trios because of parentage errors, 374 remained for analysis. Additionally, a total of 3,031 animals comprising 3,906 genotyped full-sibling pairs were available for comparison. Expected genomic inbreeding measures were calculated by predicting homozygosity independently per SNP (FPHE) in sire-dam mating pairs and by simulating progeny using phased haplotype information (FROHE and FPHE). Actual genomic inbreeding measures were calculated using the percent homozygosity of all SNP (FPH) and using runs of homozygosity (FROH). Average FPHE values (62.8±0.78%) were slightly lower than FPH (63.1±1.12%), when considering each SNP independently. After phasing haplotypes, FPHE (62.5±0.83%) was again slightly lower than FPH (62.7±1.16%), and FROHE (3.46±1.54%) was slightly lower than FROH (3.53±2.17%). Results suggest increases in expected genomic inbreeding do not explain a large effect on embryo viability at average levels of expected inbreeding. Higher variation in FROH values was present with sire-dam mating pairs exhibiting high FROHE, which may suggest high levels of genomic inbreeding are required for a noticeable effect on overall embryo viability. Genomic inbreeding between full siblings was also compared with moderate correlations (0.47-0.52) present. Overall, expected genomic inbreeding measures were calculated, but results did not suggest a large effect of expected inbreeding on embryo viability.

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.