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.

Genomic inbreeding estimation, runs of homozygosity, and heterozygosity-enriched regions uncover signals of selection in the Quarter Horse racing line

With the advent of genomics, significant progress has been made in the genetic improvement of livestock species, particularly through increased accuracy in predicting breeding values for selecting superior animals and the possibility of performing a high-resolution genetic scan throughout the genome of an individual. The main objectives of this study were to estimate the individual genomic inbreeding coefficient based on runs of homozygosity (FROH ), to identify and characterize runs of homozygosity and heterozygosity (ROH and ROHet, respectively; length and distribution) throughout the genome, and to map selection signatures in relevant chromosomal regions in the Quarter Horse racing line. A total of 336 animals registered with the Brazilian Association of Quarter Horse Breeders (ABQM) were genotyped. One hundred and twelve animals were genotyped using the Equine SNP50 BeadChip (Illumina, USA), with 54,602 single nucleotide polymorphisms (SNPs; 54K). The remaining 224 samples were genotyped using the Equine SNP70 BeadChip (Illumina, USA) with 65,157 SNPs (65K). To ensure data quality, we excluded animals with a call rate below 0.9. We also excluded SNPs located on non-autosomal chromosomes, as well as those with a call rate below 0.9 or a p-value below 1 × 10-5 for Hardy-Weinberg equilibrium. The results indicate moderate to high genomic inbreeding, with 46,594 ROH and 16,101 ROHet detected. In total, 30 and 14 candidate genes overlap with ROH and ROHet regions, respectively. The ROH islands showed genes linked to crucial biological processes, such as cell differentiation (CTBP1, WNT5B, and TMEM120B), regulation of glucose metabolic process (MAEA and NKX1-1), heme transport (PGRMC2), and negative regulation of calcium ion import (VDAC1). In ROHet, the islands showed genes related to respiratory capacity (OR7D19, OR7D4G, OR7D4E, and OR7D4J) and muscle repair (EGFR and BCL9). These findings could aid in selecting animals with greater regenerative capacity and developing treatments for muscle disorders in the QH breed. This study serves as a foundation for future research on equine breeds. It can contribute to developing reproductive strategies in animal breeding programs to improve and preserve the Quarter Horse breed.

Influence of autozygosity on common disease risk across the phenotypic spectrum

Autozygosity is associated with rare Mendelian disorders and clinically relevant quantitative traits. We investigated associations between the fraction of the genome in runs of homozygosity (FROH) and common diseases in Genes & Health (n = 23,978 British South Asians), UK Biobank (n = 397,184), and 23andMe. We show that restricting analysis to offspring of first cousins is an effective way of reducing confounding due to social/environmental correlates of FROH. Within this group in G&H+UK Biobank, we found experiment-wide significant associations between FROH and twelve common diseases. We replicated associations with type 2 diabetes (T2D) and post-traumatic stress disorder via within-sibling analysis in 23andMe (median n = 480,282). We estimated that autozygosity due to consanguinity accounts for 5%-18% of T2D cases among British Pakistanis. Our work highlights the possibility of widespread non-additive genetic effects on common diseases and has important implications for global populations with high rates of consanguinity.

Patterns of reproduction and autozygosity distinguish the breeding from non-breeding gray wolves of Yellowstone National Park

For species of management concern, accurate estimates of inbreeding and associated consequences on reproduction are crucial for predicting their future viability. However, few studies have partitioned this aspect of genetic viability with respect to reproduction in a group-living social mammal. We investigated the contributions of foundation stock lineages, putative fitness consequences of inbreeding, and genetic diversity of the breeding versus non-reproductive segment of the Yellowstone National Park gray wolf population. Our dataset spans 25 years and seven generations since reintroduction, encompassing 152 nuclear families and 329 litters. We found over 87% of the pedigree foundation genomes persisted and report influxes of allelic diversity from two translocated wolves from a divergent source in Montana. As expected for group-living species, mean kinship significantly increased over time but with minimal loss of observed heterozygosity. Strikingly, the reproductive portion of the population carried a significantly lower genome-wide inbreeding coefficients, autozygosity, and more rapid decay for linkage disequilibrium relative to the non-breeding population. Breeding wolves had significantly longer lifespans and lower inbreeding coefficients than non-breeding wolves. Our model revealed that the number of litters was negatively significantly associated with heterozygosity (R=-0.11). Our findings highlight genetic contributions to fitness, and the importance of the reproductively active individuals in a population to counteract loss of genetic variation in a wild, free-ranging social carnivore. It is crucial for managers to mitigate factors that significantly reduce effective population size and genetic connectivity, which supports the dispersion of genetic variation that aids in rapid evolutionary responses to environmental challenges.

Declining autozygosity over time: An exploration in over 1 million individuals from three diverse cohorts

Previous studies have hypothesized that autozygosity is decreasing over generational time. However, these studies were limited to relatively small samples (n < 11,000) lacking in diversity, which may limit the generalizability of their findings. We present data that partially support this hypothesis from three large cohorts of diverse ancestries, two from the US (All of Us, n = 82,474; the Million Veteran Program, n = 622,497) and one from the UK (UK Biobank, n = 380,899). Our results from a mixed-effect meta-analysis demonstrate an overall trend of decreasing autozygosity over generational time (meta-analyzed slope = -0.029, SE = 0.009, p = 6.03e-4). On the basis of our estimates, we would predict FROH to decline 0.29% for every 20-year increase in birth year. We determined that a model including an ancestry-by-country interaction term fit the data best, indicating that ancestry differences in this trend differ by country. We found further evidence to suggest a difference between the US and UK cohorts by meta-analyzing within country, observing a significant negative estimate in the US cohorts (meta-analyzed slope = -0.058, SE = 0.015, p = 1.50e-4) but a non-significant estimate in the UK (meta-analyzed slope = -0.001, SE = 0.008, p = 0.945). The association between autozygosity and birth year was substantially attenuated when accounting for educational attainment and income (meta-analyzed slope = -0.011, SE = 0.008, p = 0.167), suggesting they may partially account for decreasing autozygosity over time. Overall, we demonstrate decreasing autozygosity over time in a large, modern sample and speculate that this trend can be attributed to increases in urbanization and panmixia and differences in sociodemographic processes lead to country-specific differences in the rate of decline.

Assessment and Distribution of Runs of Homozygosity in Horse Breeds Representing Different Utility Types

The present study reports runs of homozygosity (ROH) distribution in the genomes of six horse breeds (571 horses in total) representing three horse types (primitive, light, and draft horses) based on the 65k Equine BeadChip assay. Of major interest was the length, quantity, and frequency of ROH characteristics, as well as differences between horse breeds and types. Noticeable differences in the number, length and distribution of ROH between breeds were observed, as well as in genomic inbreeding coefficients. We also identified regions of the genome characterized by high ROH coverage, known as ROH islands, which may be signals of recent selection events. Eight to fourteen ROH islands were identified per breed, which spanned multiple genes. Many were involved in important horse breed characteristics, including WFIKNN2, CACNA1G, STXBP4, NOG, FAM184B, QDPR, LCORL, and the zinc finger protein family. Regions of the genome with zero ROH occurrences were also of major interest in specific populations. Depending on the breed, we detected between 2 to 57 no-ROH regions and identified 27 genes in these regions that were common for five breeds. These genes were involved in, e.g., muscle contractility (CACNA1A) and muscle development (miR-23, miR-24, miR-27). To sum up, the obtained results can be furthered analyzed in the topic of identification of markers unique for specific horse breed characteristics.

Unravelling the effect of environment on the genome of Sarda breed ewes using Runs of Homozygosity

Natural adaptation and artificial selection have shaped the genome of modern livestock breeds. Among SNP-based metrics that are used to detect signatures of selection at genome-wide level, runs of homozygosity (ROH) are getting increasing popularity. In this paper, ROH distribution and features of a sample of 823 Sarda breed ewes farmed at different levels of altitude are analysed to investigate the effect of the environment on the patterns of homozygosity. A total of 46,829 (33,087 unique) ROH were detected. OAR2 exhibited the largest average number of ROH per animal. The most frequent ROH (OAR27, 38.9-44.2 Mb) was shared by 327. ROH length was statistically affected (p < 0.001) by both the altitude and temperature of the place where the flock was located. The highest probability of a SNP falling in a ROH was observed for hill ewes, whereas the smallest one for mountain. A total of 457 SNP exceeded the 99th percentile of the ROH count per SNP distribution and were considered significant. These markers mapped in eight chromosomes and they clustered into 17 ROH islands, where 80 candidate genes were mapped. Results of this study highlighted differences in the ROH distribution and features among sheep farmed in flocks located at different levels of altitude, confirming the role of environmental adaptability in shaping the genome of this breed.

More Than a Moggy; A Population Genetics Analysis of the United Kingdom's Non-Pedigree Cats

The domestic cat is one of the most popular pets in the world. It is estimated that 89–92% of domestic cats in the UK are non-pedigree Domestic shorthair (DSH), Domestic longhair (DLH), or Domestic semi-longhair cats (DSLH). Despite their popularity, little is known of the UK non-pedigree cats’ population structure and breeding dynamics. Using a custom designed single nucleotide variant (SNV) array, this study investigated the population genetics of 1344 UK cats. Principal components analysis (PCA) and fastSTRUCTURE analysis verified that the UK’s DSH, DLH, and DSLH cats are random-bred, rather than admixed, mix breed, or crossbred. In contrast to pedigree cats, the linkage disequilibrium of these random-bred cats was least extensive and decayed rapidly. Homozygosity by descent (HBD) analysis showed the majority of non-pedigree cats had proportionally less of their genome in HBD segments compared to pedigree cats, and that these segments were older. Together, these findings suggest that the DSH, DLH, and DSLH cats should be considered as a population of random-bred cats rather than a crossbred or pedigree-admixed cat. Unexpectedly, 19% of random-bred cat genomes displayed a higher proportion of HBD segments associated with more recent inbreeding events. Therefore, while non-pedigree cats as a whole are genetically diverse, they are not impervious to inbreeding and its health risks.

Human inbreeding has decreased in time through the Holocene

The history of human inbreeding is controversial.¹ In particular, how the development of sedentary and/or agricultural societies may have influenced overall inbreeding levels, relative to those of hunter-gatherer communities, is unclear.^2-5 Here, we present an approach for reliable estimation of runs of homozygosity (ROHs) in genomes with ≥3× mean sequence coverage across >1 million SNPs and apply this to 411 ancient Eurasian genomes from the last 15,000 years.^5-34 We show that the frequency of inbreeding, as measured by ROHs, has decreased over time. The strongest effect is associated with the Neolithic transition, but the trend has since continued, indicating a population size effect on inbreeding prevalence. We further show that most inbreeding in our historical sample can be attributed to small population size instead of consanguinity. Cases of high consanguinity were rare and only observed among members of farming societies in our sample. Despite the lack of evidence for common consanguinity in our ancient sample, consanguineous traditions are today prevalent in various modern-day Eurasian societies,^1,35-37 suggesting that such practices may have become widespread within the last few millennia.

Analysis of Homozygous-by-Descent (HBD) Segments for Purebred and Crossbred Pigs in Russia

Intensive selection raises the efficiency of pig farming considerably, but it also promotes the accumulation of homozygosity, which can lead to an increase in inbreeding and the accumulation of deleterious variation. The analysis of segments homozygous-by-descent (HBD) and non-HBD segments in purebred and crossbred pigs is of great interest. Research was carried out on 657 pigs, of which there were Large White (LW, n = 280), Landrace (LR, n = 218) and F1 female (♂LR × ♀LW) (F1, n = 159). Genotyping was performed using the GeneSeek^® GGP Porcine HD Genomic Profiler v1 (Illumina Inc., USA). To identify HBD segments and estimate autozygosity (inbreeding coefficient), we used the multiple HBD classes model. LW pigs exhibited 50,420 HBD segments, an average of 180 per animal; LR pigs exhibited 33,586 HBD segments, an average of 154 per animal; F1 pigs exhibited 21,068 HBD segments, an average of 132 per animal. The longest HBD segments in LW were presented in SSC1, SSC13 and SSC15; in LR, in SSC1; and in F1, in SSC15. In these segments, 3898 SNPs localized in 1252 genes were identified. These areas overlap with 441 QTLs (SSC1—238 QTLs; SSC13—101 QTLs; and SSC15—102 QTLs), including 174 QTLs for meat and carcass traits (84 QTLs—fatness), 127 QTLs for reproduction traits (100 QTLs—litter traits), 101 for production traits (69 QTLs—growth and 30 QTLs—feed intake), 21 QTLs for exterior traits (9 QTLs—conformation) and 18 QTLs for health traits (13 QTLs—blood parameters). Thirty SNPs were missense variants. Whilst estimating the potential for deleterious variation, six SNPs localized in the NEDD4, SEC11C, DCP1A, CCT8, PKP4 and TENM3 genes were identified, which may show deleterious variation. A high frequency of potential deleterious variation was noted for LR in DCP1A, and for LW in TENM3 and PKP4. In all cases, the genotype frequencies in F1 were intermediate between LR and LW. The findings presented in our work show the promise of genome scanning for HBD as a strategy for studying population history, identifying genomic regions and genes associated with important economic traits, as well as deleterious variation.

Ethiopian indigenous goats offer insights into past and recent demographic dynamics and local adaptation in sub-Saharan African goats

Knowledge on how adaptive evolution and human socio-cultural and economic interests shaped livestock genomes particularly in sub-Saharan Africa remains limited. Ethiopia is in a geographic region that has been critical in the history of African agriculture with ancient and diverse human ethnicity and bio-climatic conditions. Using 52K genome-wide data analysed in 646 individuals from 13 Ethiopian indigenous goat populations, we observed high levels of genetic variation. Although runs of homozygosity (ROH) were ubiquitous genome-wide, there were clear differences in patterns of ROH length and abundance and in effective population sizes illustrating differences in genome homozygosity, evolutionary history, and management. Phylogenetic analysis incorporating patterns of genetic differentiation and gene flow with ancestry modelling highlighted past and recent intermixing and possible two deep ancient genetic ancestries that could have been brought by humans with the first introduction of goats in Africa. We observed four strong selection signatures that were specific to Arsi-Bale and Nubian goats. These signatures overlapped genomic regions with genes associated with morphological, adaptation, reproduction and production traits due possibly to selection under environmental constraints and/or human preferences. The regions also overlapped uncharacterized genes, calling for a comprehensive annotation of the goat genome. Our results provide insights into mechanisms leading to genome variation and differentiation in sub-Saharan Africa indigenous goats.

Diversity analysis, runs of homozygosity and genomic inbreeding reveal recent selection in Blanco Orejinegro cattle

Blanco Orejinegro (BON) cattle have 500 years of adaptation to the Colombian tropic, but little is known about their genetic history. Our aim was to estimate levels of linkage disequilibrium (LD), effective population size (Ne), genomic inbreeding for runs of homozygosity (F_ROH ), genomic relation matrix (F_GRM ), excess of homozygotes (F_HOM ) and pedigree information (F_PEDCOMP ) and to characterize the runs of homozygosity (ROH), searching for selection signatures. A total of 419 BON animals were genotyped, 70 with a 150K chip and 349 with a 50K chip. Next, an imputation to 50K was performed, and, after editing, databases of 40K were obtained. The PLINK v1.90 and R programs were used to estimate LD, ROH, F_ROH and F_HOM . The SNeP v1.1 program was used to obtain Ne, and PreGSf90 was used to elaborate the scaled G matrix. The MTDFNRM program was used to estimate F_PEDCOMP . The LD mean as r² at 1 Mb was 0.21 (r²  > 0.30 at a distance of 96.72kb), and Ne was 123 ± 1. A total of 7,652 homozygous segments were obtained, with a mean of 18.35 ± 0.55 ROH/animal. Most of the genome was covered by long ROHs (ROH_{>8 Mb}  = 4.86%), indicating significant recent inbreeding. The average inbreeding coefficient for F_PEDCOM , F_GRM , F_HOM and F_ROH was 4.41%, 4.18%, 5.58% and 6.78%, respectively. The highest correlation was observed between F_HOM and F_ROH (0.95). ROH hotspots/islands were defined using the extreme values of a box plot that was generated, and correspond to QTLs related to milk yield (55.11%), external appearance (13.47%), production (13.30%), reproduction (8.15%), health (5.24%) and meat carcass (4.74%).

Runs of homozygosity provide a genome landscape picture of inbreeding and genetic history of European autochthonous and commercial pig breeds

ROHs are long stretches of DNA homozygous at each polymorphic position. The proportion of genome covered by ROHs and their length are indicators of the level and origin of inbreeding. Frequent common ROHs within the same population define ROH islands and indicate hotspots of selection. In this work, we investigated ROHs in a total of 1131 pigs from 20 European local pig breeds and in three cosmopolitan breeds, genotyped with the GGP Porcine HD Genomic Profiler. plink software was used to identify ROHs. Size classes and genomic inbreeding parameters were evaluated. ROH islands were defined by evaluating different thresholds of homozygous SNP frequency. A functional overview of breed-specific ROH islands was obtained via over-representation analyses of GO biological processes. Mora Romagnola and Turopolje breeds had the largest proportions of genome covered with ROH (~1003 and ~955 Mb respectively), whereas Nero Siciliano and Sarda breeds had the lowest proportions (~207 and 247 Mb respectively). The highest proportion of long ROH (>16 Mb) was in Apulo-Calabrese, Mora Romagnola and Casertana. The largest number of ROH islands was identified in the Italian Landrace (n = 32), Cinta Senese (n = 26) and Lithuanian White Old Type (n = 22) breeds. Several ROH islands were in regions encompassing genes known to affect morphological traits. Comparative ROH structure analysis among breeds indicated the similar genetic structure of local breeds across Europe. This study contributed to understanding of the genetic history of the investigated pig breeds and provided information to manage these pig genetic resources.

Recent Consanguinity and Outbred Autozygosity Are Associated With Increased Risk of Late-Onset Alzheimer's Disease

Prior work in late-onset Alzheimer's disease (LOAD) has resulted in discrepant findings as to whether recent consanguinity and outbred autozygosity are associated with LOAD risk. In the current study, we tested the association between consanguinity and outbred autozygosity with LOAD in the largest such analysis to date, in which 20 LOAD GWAS datasets were retrieved through public databases. Our analyses were restricted to eight distinct ethnic groups: African-Caribbean, Ashkenazi-Jewish European, European-Caribbean, French-Canadian, Finnish European, North-Western European, South-Eastern European, and Yoruba African for a total of 21,492 unrelated subjects (11,196 LOAD and 10,296 controls). Recent consanguinity determination was performed using FSuite v1.0.3, according to subjects' ancestral background. The level of autozygosity in the outbred population was assessed by calculating inbreeding estimates based on the proportion (FROH) and the number (NROH) of runs of homozygosity (ROHs). We analyzed all eight ethnic groups using a fixed-effect meta-analysis, which showed a significant association of recent consanguinity with LOAD (N = 21,481; OR = 1.262, P = 3.6 × 10-4), independently of APOE ∗4 (N = 21,468, OR = 1.237, P = 0.002), and years of education (N = 9,257; OR = 1.274, P = 0.020). Autozygosity in the outbred population was also associated with an increased risk of LOAD, both for F ROH (N = 20,237; OR = 1.204, P = 0.030) and N ROH metrics (N = 20,237; OR = 1.019, P = 0.006), independently of APOE ∗4 [(F ROH, N = 20,225; OR = 1.222, P = 0.029) (N ROH, N = 20,225; OR = 1.019, P = 0.007)]. By leveraging the Alzheimer's Disease Sequencing Project (ADSP) whole-exome sequencing (WES) data, we determined that LOAD subjects do not show an enrichment of rare, risk-enhancing minor homozygote variants compared to the control population. A two-stage recessive GWAS using ADSP data from 201 consanguineous subjects in the discovery phase followed by validation in 10,469 subjects led to the identification of RPH3AL p.A303V (rs117190076) as a rare minor homozygote variant increasing the risk of LOAD [discovery: Genotype Relative Risk (GRR) = 46, P = 2.16 × 10-6; validation: GRR = 1.9, P = 8.0 × 10-4]. These results confirm that recent consanguinity and autozygosity in the outbred population increase risk for LOAD. Subsequent work, with increased samples sizes of consanguineous subjects, should accelerate the discovery of non-additive genetic effects in LOAD.

Genetic Background and Inbreeding Depression in Romosinuano Cattle Breed in Mexico

Simple Summary

The objective of this study was to evaluate the genetic background and inbreeding depression in the Mexican Romosinuano cattle using pedigree and genomic information. Inbreeding was estimated using pedigree (FPED) and genomic information based on the genomic relationship matrix (FGRM) and runs of homozygosity (FROH). Linkage disequilibrium (LD) was evaluated using the correlation between pairs of loci, and the effective population size (Ne) was calculated based on LD and pedigree information. The pedigree file consisted of 4875 animals; 71 had genotypes. LD decreased with the increase in distance between markers, and Ne estimated using genomic information decreased from 610 to 72 animals (from 109 to 1 generation ago), the Ne estimated using pedigree information was 86.44. The number of runs of homozygosity per animal ranged between 18 and 102 segments with an average of 55. The average inbreeding was 2.98 ± 2.81, 2.98 ± 4.01, and 7.28 ± 3.68% for FPED, FGRM, and FROH, respectively. A 1% increase in inbreeding decreased birth weight by 0.103 kg and weaning weight by 0.685 kg. A strategy such as optimum genetic contributions to maximize selection response and manage the long-term genetic variability and inbreeding could lead to sustainable breeding programs for the Mexican Romosinuano cattle breed.

Abstract

The ultimate goal of genetic selection is to improve genetic progress by increasing favorable alleles in the population. However, with selection, homozygosity, and potentially harmful recessive alleles can accumulate, deteriorating genetic variability and hampering continued genetic progress. Such potential adverse side effects of selection are of particular interest in populations with a small effective population size like the Romosinuano beef cattle in Mexico. The objective of this study was to evaluate the genetic background and inbreeding depression in Mexican Romosinuano cattle using pedigree and genomic information. Inbreeding was estimated using pedigree (FPED) and genomic information based on the genomic relationship matrix (FGRM) and runs of homozygosity (FROH) of different length classes. Linkage disequilibrium (LD) was evaluated using the correlation between pairs of loci, and the effective population size (Ne) was calculated based on LD and pedigree information. The pedigree file consisted of 4875 animals born between 1950 and 2019, of which 71 had genotypes. LD decreased with the increase in distance between markers, and Ne estimated using genomic information decreased from 610 to 72 animals (from 109 to 1 generation ago), the Ne estimated using pedigree information was 86.44. The reduction in effective population size implies the existence of genetic bottlenecks and the decline of genetic diversity due to the intensive use of few individuals as parents of the next generations. The number of runs of homozygosity per animal ranged between 18 and 102 segments with an average of 55. The shortest and longest segments were 1.0 and 36.0 Mb long, respectively, reflecting ancient and recent inbreeding. The average inbreeding was 2.98 ± 2.81, 2.98 ± 4.01, and 7.28 ± 3.68% for FPED, FGRM, and FROH, respectively. The correlation between FPED and FGRM was −0.25, and the correlations among FPED and FROH of different length classes were low (from 0.16 to 0.31). The correlations between FGRM and FROH of different length classes were moderate (from 0.44 to 0.58), indicating better agreement. A 1% increase in population inbreeding decreased birth weight by 0.103 kg and weaning weight by 0.685 kg. A strategy such as optimum genetic contributions to maximize selection response and manage the long-term genetic variability and inbreeding could lead to more sustainable breeding programs for the Mexican Romosinuano beef cattle breed.

Exploiting the Autozygome to Support Previously Published Mendelian Gene-Disease Associations: An Update

There is a growing interest in standardizing gene-disease associations for the purpose of facilitating the proper classification of variants in the context of Mendelian diseases. One key line of evidence is the independent observation of pathogenic variants in unrelated individuals with similar phenotypes. Here, we expand on our previous effort to exploit the power of autozygosity to produce homozygous pathogenic variants that are otherwise very difficult to encounter in the homozygous state due to their rarity. The identification of such variants in genes with only tentative associations to Mendelian diseases can add to the existing evidence when observed in the context of compatible phenotypes. In this study, we report 20 homozygous variants in 18 genes (ADAMTS18, ARNT2, ASTN1, C3, DMBX1, DUT, GABRB3, GM2A, KIF12, LOXL3, NUP160, PTRHD1, RAP1GDS1, RHOBTB2, SIGMAR1, SPAST, TENM3, and WASHC5) that satisfy the ACMG classification for pathogenic/likely pathogenic if the involved genes had confirmed rather than tentative links to diseases. These variants were selected because they were truncating, founder with compelling segregation or supported by robust functional assays as with the DUT variant that we present its validation using yeast model. Our findings support the previously reported disease associations for these genes and represent a step toward their confirmation.

Genome-wide analysis of homozygosity regions in european simmental bulls

The study of Runs of Homozygosity (ROH) is a useful approach for the characterization of the genome of livestock populations. Due to their high relationship with autozygosity, ROH allow to make inference about population genetic history, to estimate the level of inbreeding, to assess within breed heterogeneity and to detect the footprints of selection on livestock genomes. Aim of this study was to investigate the distribution of runs of homozygosity in bulls belonging to five European Simmental populations and to assess the relationship between three production traits (milk yield, fat and protein contents) and autozygosity. ROH count, distribution and ROH-based coefficient of inbreeding (F_ROH ) were calculated for 3,845 Simmental bulls of five different European countries: Austria (AT), Switzerland (CH), Czech Republic (CZ), Germany (DE) and Italy (IT). Average values of ROH number per animal, and total genome length covered by ROH were 77.8 ± 20.7 and 205 ± 74.4 Mb, respectively. Bulls from AT, DE and IT exhibited similar ROH characteristics. Swiss animals showed the highest (12.6%), while CZ the lowest (4.6%) F_ROH coefficient. The relationship between ROH occurrence and milk production traits was investigated through a genome-wide ROH-traits association analysis (GWRA). A total of 34 regions previously associated with milk traits (yield and/or composition) were identified by GWRA. Results of the present research highlight a mixed genetic background in the 5 European Simmental populations, with the possible presence of three subgroups. Moreover, a strong relationship between autozygosity and production traits has been detected.

Estimates of Autozygosity Through Runs of Homozygosity in Farmed Coho Salmon

The characterization of runs of homozygosity (ROH), using high-density single nucleotide polymorphisms (SNPs) allows inferences to be made about the past demographic history of animal populations and the genomic ROH has become a common approach to characterize the inbreeding. We aimed to analyze and characterize ROH patterns and compare different genomic and pedigree-based methods to estimate the inbreeding coefficient in two pure lines (POP A and B) and one recently admixed line (POP C) of coho salmon (Oncorhynchus kisutch) breeding nuclei, genotyped using a 200 K Affymetrix Axiom^® myDesign Custom SNP Array. A large number and greater mean length of ROH were found for the two “pure” lines and the recently admixed line (POP C) showed the lowest number and smaller mean length of ROH. The ROH analysis for different length classes suggests that all three coho salmon lines the genome is largely composed of a high number of short segments (<4 Mb), and for POP C no segment >16 Mb was found. A high variable number of ROH, mean length and inbreeding values across chromosomes; positively the consequence of artificial selection. Pedigree-based inbreeding values tended to underestimate genomic-based inbreeding levels, which in turn varied depending on the method used for estimation. The high positive correlations between different genomic-based inbreeding coefficients suggest that they are consistent and may be more accurate than pedigree-based methods, given that they capture information from past and more recent demographic events, even when there are no pedigree records available.

A detailed characteristics of bias associated with long runs of homozygosity identification based on medium density SNP microarrays

In the present study, runs of homozygosity (ROH) detected with the use of a standard bovine 54k single nucleotide polymorphism (SNP) genotyping assay and two different ROH detection approaches, based on 50 (M1) or 15 (M2) consecutive SNPs, were compared with results of whole genome sequencing. Both microarray-based methods accurately recognised medium-sized ROH, however, it was found that M2 method seemed to better than M1 identify short ROH, but highly overestimated their number, leading to numerous false positive calls. Moreover, long ROH identified with microarray data tended to break into shorter segments in sequencing data because of the presence of regions with high heterozygosity within the ROH sequences. This may indicate, that these long ROH are formed by closely positioned shorter homozygous segments that may be of older origin or may be created by two similar but not identical haplotypes, showing minor internal recombination signs. Such finding also suggests that at least some of the results of previous studies in regard to long ROH may be biased leading to inaccurate estimations of genomes autozygosity via ROH classification into length categories.

Genome-Wide Runs of Homozygosity Revealed Selection Signatures in Bos indicus

Genome-wide runs of homozygosity (ROH) are suitable for understanding population history, calculating genomic inbreeding, deciphering genetic architecture of complex traits and diseases as well as identifying genes linked with agro-economic traits. Autozygosity and ROH islands, genomic regions with elevated ROH frequencies, were characterized in 112 animals of seven Indian native cattle breeds (B. indicus) using BovineHD BeadChip. In total, 4138 ROH were detected. The average number of ROH per animal was maximum in draft breed, Kangayam (63.62 ± 22.71) and minimum in dairy breed, Sahiwal (24.62 ± 11.03). The mean ROH length was maximum in Vechur (6.97 Mb) and minimum in Hariana (4.04 Mb). Kangayam revealed the highest ROH based inbreeding (F_{ROH > 1Mb} = 0.113 ± 0.059), whereas Hariana (F_{ROH > 1Mb} = 0.042 ± 0.031) and Sahiwal (F_{ROH > 1Mb} = 0.043 ± 0.048) showed the lowest. The high standard deviation observed in each breed highlights a considerable variability in autozygosity. Out of the total autozygous segments observed in each breed except Vechur, > 80% were of short length (< 8 Mb) and contributed almost 50% of the genome proportion under ROH. However, in Vechur cattle, long ROH contributed 75% of the genome proportion under ROH. ROH patterns revealed Hariana and Sahiwal breeds as less consanguineous, while recent inbreeding was apparent in Vechur. Maximum autozygosity observed in Kangayam is attributable to both recent and ancient inbreeding. The ROH islands were harbouring higher proportion of QTLs for production traits (20.68% vs. 14.64%; P≤ 0.05) but lower for reproductive traits (11.49% vs. 15.76%; P≤ 0.05) in dairy breeds compared to draft breed. In draft cattle, genes associated with resistant to diseases/higher immunity (LYZL1, SVIL, and GPX4) and stress tolerant (CCT4) were identified in ROH islands; while in dairy breeds, for milk production (PTGFR, CSN1S1, CSN2, CSN1S2, and CSN3). Significant difference in ROH islands among large and short statured breeds was observed at chromosome 3 and 5 involving genes like PTGFR and HMGA2 responsible for milk production and stature, respectively. PCA analysis on consensus ROH regions revealed distinct clustering of dairy, draft and short stature cattle breeds.

Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (FROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (FHOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (FPED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. FHOM was highly correlated (0.671 to 0.985) with FROH measures in all breeds. Apulo-Calabrese and Casertana had the highest FROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (FROH1) to 16 Mbp (FROH16)), whereas the lowest FROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). FROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all FROH measures and FPED (from 0.514 to 0.523) and between FHOM and FPED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between FROH and FPED. This might be due to the imperfect measure of FPED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.

A Comprehensive Analysis of Runs of Homozygosity of Eleven Cattle Breeds Representing Different Production Types

Simple Summary

Runs of homozygosity (ROH) regions are known to be common in the genomes of cattle and have become a subject of interest of various research in recent years. ROH can be used as a valuable tool to estimate inbreeding, which needs to be controlled in livestock populations. Moreover, analysis of ROH is considered to be an effective method of identifying genome regions that are a subject of selection pressure, which may help in understanding the genetic aspects of production traits under selection. In this study, we analyzed ROH characteristics of 11 cattle breeds, both commercial and native, maintained in Poland. We presented distinct differences in the length, quantity and frequency of ROH between the analyzed breeds as well as in the levels of genomic inbreeding. Higher levels of inbreeding were characteristic for commercial breeds, especially beef breeds. In addition, within ROH islands, we observed a number of genes with a confirmed influence on the level of production traits. The presented results and identified genes can be a basis for further research focused on the identification of genes and markers essential in the determination of the most important production traits in cattle.

Abstract

In the presented research, BovineSNP50 microarrays (Illumina) were applied to determine runs of homozygosity in the genomes of 11 cattle breeds maintained in Poland. These cattle breeds represent three basic utility types: milk, meat and dual purpose. Analysis of runs of homozygosity allowed the evaluation of the level of autozygosity within each breed in order to calculate the genomic inbreeding coefficient (F_ROH), as well as to identify regions of the genome with a high frequency of ROH occurrence, which may reflect traces of directional selectin left in their genomes. Visible differences in the length and distribution of runs of homozygosity in the genomes of the analyzed cattle breeds have been observed. The highest mean number and mean sums of lengths of runs of homozygosity were characteristic for Hereford cattle and intermediate for the Holstein-Friesian Black-and-White variety, Holstein-Friesian Red-and-White variety, Simmental, Limousin, Montbeliarde and Charolais breeds. However, lower values were observed for cattle of conserved breeds. Moreover, the selected livestock differed in the level of inbreeding estimated using the F_ROH coefficient. In regions of the genome with a high frequency of ROH occurrence, which may reflect the impact of directional selection, a number of genes were observed that can be potentially related to the production traits which are under selection pressure for specific production types. The most important detected genes were GHR, MSTN, DGAT1, FABP4, and TRH, with a known influence on the milk and meat traits of the studied cattle breeds.

Autozygosity and Genetic Differentiation of Landrace and Large White Pigs as Revealed by the Genetic Analyses of Crossbreds

Genomic information from crossbreds is routinely generated for genomic evaluations. The objective of this study is to investigate autozygosity and genetic differentiation in Landrace by Large-White breeds by using the genotypic information of SNP arrays in 1,173 crossbreds. A maximum likelihood approach was developed to estimate the probability of autozygosity (F_L). Regions of differentiation between breeds were investigated using F_ST and the difference in allele frequencies between the two parental breeds (릌Δ) at each single-nucleotide polymorphism (SNP) position. A maximum likelihood approach was proposed to estimate allele frequencies in the parental populations. The average length of runs of homozygosity (ROH) across the genome was 3.91, 2.3, and 0.7 Mb for segments with at least 25, 15, and 5 SNPs, respectively. Average age to coalesce was 46, 414, and 388 years for segments with at least 25, 15, and 5 SNPs, respectively. The probability of autozygosity was not uniform along the crossbred genome, being higher at the center for most chromosomes. The correlation between autozygosity and distance to the closest telomere was positive and significant in most chromosomes, which could be attributed to the higher recombination rate near telomeres. We also report a relatively high negative correlation between probability of recombination (from a published map) and probability of autozygosity. It supports that structural characteristics of the chromosomes related to recombination rate determine autozygosity at each chromosomal position of the pig genome. The average is Δ across the genome was 0.17 (SD = 0.16). After testing for differences in allele frequencies between the parental breeds, there were 4,184 SNPs with a likelihood ratio test, LRT ≥ 32.02. The average F_ST across the genome was 0.038 (SD = 0.059). There were 2,949 SNPs with F_ST > 0.125. The correlation between estimates of F_L and estimates of F_ST across the genome was -0.10 (SE = 0.006). Analysis of the gene content of the genomic regions with the 2000 SNPs with highest LRT for F_L and high F_ST showed overrepresentation of genes with a regulatory function. Genes with biological functions associated with production, such as tissue development, anatomical structure, and animal organ development, were also overrepresented in regions with a high F_ST.

Relationship between Deleterious Variation, Genomic Autozygosity, and Disease Risk: Insights from The 1000 Genomes Project

Genomic regions of autozygosity (ROAs) represent segments of individual genomes that are homozygous for haplotypes inherited identical-by-descent (IBD) from a common ancestor. ROAs are nonuniformly distributed across the genome, and increased ROA levels are a reported risk factor for numerous complex diseases. Previously, we hypothesized that long ROAs are enriched for deleterious homozygotes as a result of young haplotypes with recent deleterious mutations-relatively untouched by purifying selection-being paired IBD as a consequence of recent parental relatedness, a pattern supported by ROA and whole-exome sequence data on 27 individuals. Here, we significantly bolster support for our hypothesis and expand upon our original analyses using ROA and whole-genome sequence data on 2,436 individuals from The 1000 Genomes Project. Considering CADD deleteriousness scores, we reaffirm our previous observation that long ROAs are enriched for damaging homozygotes worldwide. We show that strongly damaging homozygotes experience greater enrichment than weaker damaging homozygotes, while overall enrichment varies appreciably among populations. Mendelian disease genes and those encoding FDA-approved drug targets have significantly increased rates of gain in damaging homozygotes with increasing ROA coverage relative to all other genes. In genes implicated in eight complex phenotypes for which ROA levels have been identified as a risk factor, rates of gain in damaging homozygotes vary across phenotypes and populations but frequently differ significantly from non-disease genes. These findings highlight the potential confounding effects of population background in the assessment of associations between ROA levels and complex disease risk, which might underlie reported inconsistencies in ROA-phenotype associations.

Runs of homozygosity: current knowledge and applications in livestock

This review presents a broader approach to the implementation and study of runs of homozygosity (ROH) in animal populations, focusing on identifying and characterizing ROH and their practical implications. ROH are continuous homozygous segments that are common in individuals and populations. The ability of these homozygous segments to give insight into a population's genetic events makes them a useful tool that can provide information about the demographic evolution of a population over time. Furthermore, ROH provide useful information about the genetic relatedness among individuals, helping to minimize the inbreeding rate and also helping to expose deleterious variants in the genome. The frequency, size and distribution of ROH in the genome are influenced by factors such as natural and artificial selection, recombination, linkage disequilibrium, population structure, mutation rate and inbreeding level. Calculating the inbreeding coefficient from molecular information from ROH (F_ROH ) is more accurate for estimating autozygosity and for detecting both past and more recent inbreeding effects than are estimates from pedigree data (F_PED ). The better results of F_ROH suggest that F_ROH can be used to infer information about the history and inbreeding levels of a population in the absence of genealogical information. The selection of superior animals has produced large phenotypic changes and has reshaped the ROH patterns in various regions of the genome. Additionally, selection increases homozygosity around the target locus, and deleterious variants are seen to occur more frequently in ROH regions. Studies involving ROH are increasingly common and provide valuable information about how the genome's architecture can disclose a population's genetic background. By revealing the molecular changes in populations over time, genome-wide information is crucial to understanding antecedent genome architecture and, therefore, to maintaining diversity and fitness in endangered livestock breeds.

Importance of Genetic Studies in Consanguineous Populations for the Characterization of Novel Human Gene Functions

Consanguineous offspring have elevated levels of homozygosity. Autozygous stretches within their genome are likely to harbour loss of function (LoF) mutations which will lead to complete inactivation or dysfunction of genes. Studying consanguineous offspring with clinical phenotypes has been very useful for identifying disease causal mutations. However, at present, most of the genes in the human genome have no disorder associated with them or have unknown function. This is presumably mostly due to the fact that homozygous LoF variants are not observed in outbred populations which are the main focus of large sequencing projects. However, another reason may be that many genes in the genome—even when completely “knocked out,” do not cause a distinct or defined phenotype. Here, we discuss the benefits and implications of studying consanguineous populations, as opposed to the traditional approach of analysing a subset of consanguineous families or individuals with disease. We suggest that studying consanguineous populations “as a whole” can speed up the characterisation of novel gene functions as well as indicating nonessential genes and/or regions in the human genome. We also suggest designing a single nucleotide variant (SNV) array to make the process more efficient.

Identification of three novel ECEL1 mutations in three families with distal arthrogryposis type 5D

Arthrogryposis refers to congenital contracture in at least two different body parts. When distal joints are primarily involved, the term distal arthrogryposis (DA) is used. The recognition of clinically distinct subtypes of DA has proven very useful in mapping the disease genes for this genetically heterogeneous condition. DA5D is characterized by ocular involvement usually in the form of ptosis and incomitant strabismus, but extraocular manifestations have also been reported. In a multiplex consanguineous family with DA5D, we combined autozygosity mapping and exome sequencing to identify a novel mutation in ECEL1. This was followed by targeted sequencing of this gene in another two extended consanguineous family with the same phenotype, which revealed two additional novel homozygous mutations. Our results support the recent identification of mutations in ECEL1 as a disease gene in DA5D and expand the clinical and allelic spectrum of this condition.

Estimating the human mutation rate using autozygosity in a founder population

Knowledge of the rate and pattern of new mutation is critical to the understanding of human disease and evolution. We used extensive autozygosity in a genealogically well-defined population of Hutterites to estimate the human sequence mutation rate over multiple generations. We sequenced whole genomes from 5 parent-offspring trios and identified 44 segments of autozygosity. Using the number of meioses separating each pair of autozygous alleles and the 72 validated heterozygous single-nucleotide variants (SNVs) from 512 Mb of autozygous DNA, we obtained an SNV mutation rate of 1.20 × 10(-8) (95% confidence interval 0.89-1.43 × 10(-8)) mutations per base pair per generation. The mutation rate for bases within CpG dinucleotides (9.72 × 10(-8)) was 9.5-fold that of non-CpG bases, and there was strong evidence (P = 2.67 × 10(-4)) for a paternal bias in the origin of new mutations (85% paternal). We observed a non-uniform distribution of heterozygous SNVs (both newly identified and known) in the autozygous segments (P = 0.001), which is suggestive of mutational hotspots or sites of long-range gene conversion.

Pedigree and marker information requirements to monitor genetic variability

There are several measures available to describe the genetic variability of populations. The average inbreeding coefficient of a population based on pedigree information is a frequently chosen option. Due to the developments in molecular genetics it is also possible to calculate inbreeding coefficients based on genetic marker information. A simulation study was carried out involving ten sires and 50 dams. The animals were mated over a period of 20 discrete generations. The population size was kept constant. Different situations with regard to the level of polymorphism and initial allele frequencies and mating scheme (random mating, avoidance of full sib mating, avoidance of full sib and half sib mating) were considered. Pedigree inbreeding coefficients of the last generation using full pedigree or 10, 5 and 2 generations of the pedigree were calculated. Marker inbreeding coefficients based on different sets of microsatellite loci were also investigated. Under random mating, pedigree-inbreeding coefficients are clearly more closely related to true autozygosity (i.e., the actual proportion of loci with alleles identical by descent) than marker-inbreeding coefficients. If mating is not random, the demands on the quality and quantity of pedigree records increase. Greater attention must be paid to the correct parentage of the animals.

Lack of involvement of known DNA methyltransferases in familial hydatidiform mole implies the involvement of other factors in establishment of imprinting in the human female germline

BACKGROUND

Differential methylation of the two alleles is a hallmark of imprinted genes. Correspondingly, loss of DNA methyltransferase function results in aberrant imprinting and abnormal post-fertilization development. In the mouse, mutations of the oocyte-specific isoform of the DNA methyltransferase Dnmt1 (Dnmt1o) and of the methyltransferase-like Dnmt3L gene result in specific failures of imprint establishment or maintenance, at multiple loci. We have previously shown in humans that an analogous inherited failure to establish imprinting at multiple loci in the female germline underlies a rare phenotype of recurrent hydatidiform mole.

RESULTS

We have identified a human homologue of the murine Dnmt1o and assessed its pattern of expression. Human DNMT1o mRNA is detectable in mature oocytes and early fertilized embryos but not in any somatic tissues analysed. The somatic isoform of DNMT1 mRNA, in contrast, is not detectable in human oocytes. In the previously-described family with multi-locus imprinting failure, mutation of DNMT1o and of the other known members of this gene family has been excluded.

CONCLUSIONS

Mutation of the known DNMT genes does not underlie familial hydatidiform mole, at least in the family under study. This suggests that trans-acting factors other than the known methyltransferases are required for imprint establishment in humans, a concept that has indirect support from recent biochemical studies of DNMT3L.