Dataset Information

Runs of homozygosity and analysis of inbreeding depression.

ABSTRACT: Pedigree information was traditionally used to assess inbreeding. The availability of high-density marker panels provides an alternative to assess inbreeding, particularly in the presence of incomplete and error-prone pedigrees. Assessment of autozygosity across chromosomal segments using runs of homozygosity (ROH) has emerged as a valuable tool to estimate inbreeding due to its general flexibility and ability to quantify the chromosomal contribution to genome-wide inbreeding. Unfortunately, the identification of ROH segments is sensitive to the parameters used during the search process. These parameters are heuristically set, leading to significant variation in the results. The minimum length required to identify an ROH segment has major effects on the estimation of inbreeding and inbreeding depression, yet it is arbitrarily set. To overcome this limitation, a search algorithm to approximate mutation enrichment was developed to determine the minimum length of ROH segments. It consists of finding genome segments with significant effect differences in trait means between animals with high and low burdens of autozygous intervals with a specific length. The minimum length could be determined heuristically as the smallest interval at which a significant signal is detected. The proposed method was tested in an inbred Hereford cattle population genotyped for 30,220 SNPs. Phenotypes recorded for six traits were used for the approximation of mutation loads. The estimated minimum length was around 1 Mb for yearling weight (YW) and average daily gain (ADG) and 4 Mb for birth weight and weaning weight. These trait-specific thresholds estimated using the proposed method could be attributed to a trait-dependent effect of homozygosity. The detection of significant inbreeding effects was well aligned with the estimated thresholds, especially for YW and ADG. Although highly deleterious alleles are expected to be more frequent in recent inbreeding (long ROH), short ROH segments (<5 Mb) could contain a large number of less deleterious mutations with substantial joint effects on some traits (YW and ADG). Our results highlight the importance of accurate estimation of the ROH-based inbreeding and the necessity to consider a trait-specific minimum length threshold for the identification of ROH segments in inbreeding depression analyses. These thresholds could be determined using the proposed method provided the availability of phenotypic information.

SUBMITTER: Sumreddee P

PROVIDER: S-EPMC7846130 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BACKGROUND:Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (FROH), genomic relationship matrix approach (FGRM) and based on the observed versus expected number of homozygous genotypes (FHOM), and from pedigree-based coefficient (FPED). RESULTS:ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH1-2 Mb) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of FROH, FGRM, FHOM, and FPED approaches. FPED estimates ranged from 0.00 to 0.327 and FROH from 0.001 to 0.201. Low to moderate correlations were observed between FPED-FROH and FGRM-FROH, with values ranging from -0.11 to 0.51. Low to high correlations were observed between FROH-FHOM and moderate between FPED-FHOM and FGRM-FHOM. Correlations between FROH from different lengths and FPED gradually increased with ROH length. CONCLUSIONS:Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low FPED-FROH correlations for small segments indicate that FPED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that FROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.

Project description:BackgroundConsecutive homozygous fragments of a genome inherited by offspring from a common ancestor are known as runs of homozygosity (ROH). ROH can be used to calculate genomic inbreeding and to identify genomic regions that are potentially under historical selection pressure. The dataset of our study consisted of 254 Azeri (AZ) and 115 Khuzestani (KHZ) river buffalo genotyped for ~?65,000 SNPs for the following two purposes: 1) to estimate and compare inbreeding calculated using ROH (FROH), excess of homozygosity (FHOM), correlation between uniting gametes (FUNI), and diagonal elements of the genomic relationship matrix (FGRM); 2) to identify frequently occurring ROH (i.e. ROH islands) for our selection signature and gene enrichment studies.ResultsIn this study, 9102 ROH were identified, with an average number of 21.2?±?13.1 and 33.2?±?15.9 segments per animal in AZ and KHZ breeds, respectively. On average in AZ, 4.35% (108.8?±?120.3?Mb), and in KHZ, 5.96% (149.1?±?107.7?Mb) of the genome was autozygous. The estimated inbreeding values based on FHOM, FUNI and FGRM were higher in AZ than they were in KHZ, which was in contrast to the FROH estimates. We identified 11 ROH islands (four in AZ and seven in KHZ). In the KHZ breed, the genes located in ROH islands were enriched for multiple Gene Ontology (GO) terms (P ? 0.05). The genes located in ROH islands were associated with diverse biological functions and traits such as body size and muscle development (BMP2), immune response (CYP27B1), milk production and components (MARS, ADRA1A, and KCTD16), coat colour and pigmentation (PMEL and MYO1A), reproductive traits (INHBC, INHBE, STAT6 and PCNA), and bone development (SUOX).ConclusionThe calculated FROH was in line with expected higher inbreeding in KHZ than in AZ because of the smaller effective population size of KHZ. Thus, we find that FROH can be used as a robust estimate of genomic inbreeding. Further, the majority of ROH peaks were overlapped with or in close proximity to the previously reported genomic regions with signatures of selection. This tells us that it is likely that the genes in the ROH islands have been subject to artificial or natural selection.

Dataset Information

Runs of homozygosity and analysis of inbreeding depression.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets