Detecting autozygosity through runs of homozygosity: a comparison of three autozygosity detection algorithms

Genomic inbreeding analysis reveals resilience and genetic diversity in Indian yak populations

The yak (Bos grunniens), renowned for its adaptability to extreme cold and hypoxic conditions, stands as a remarkable domestic animal crucial for sustaining livelihoods in harsh climates. We conducted a comprehensive analysis of the whole genome sequence data from three distinct Indian yak populations: Arunachali yak (n = 10), Himachali yak (n = 10), and Ladakhi yak (n = 10). The genomic data for Indian yaks were meticulously generated by our laboratory and compared with their Chinese counterpart, the Jinchuan yak (n = 8), for a more nuanced understanding. Our investigation revealed a total of 37,437 runs of homozygosity (ROH) segments in 34 animals representing four distinct yak populations. The Jinchuan yak population exhibited the highest proportion, constituting 80.8 % of total ROHs, predominantly as small segments (<0.1 Mb), accounting for 63 % of the overall ROHs. Further analysis uncovered a significantly higher degree of inbreeding in Chinese yaks compared to their Indian counterparts. The Indian yak populations, in contrast, demonstrated relatively lower and consistent levels of inbreeding. Moreover, we identified ROH hotspots that covered at least 60 % of individuals in our study, indicating their pivotal role in environmental adaptation. A total of five hotspot regions were detected, housing genes such as ENSBGRG00000015023 (WNT2), YIPF4, SPAST, TLN2, and DSG4. These genes are associated with traits including hair follicle initiation, nutrient stress response, microtubule assembly, development of cardiac muscle, hair follicle, and coat color. This observation strongly suggests that there is substantial selection acting on these genes, emphasizing their important role in environmental adaptation among yak populations.

Polygenic risk score portability for common diseases across genetically diverse populations

BACKGROUND

Polygenic risk scores (PRS) derived from European individuals have reduced portability across global populations, limiting their clinical implementation at worldwide scale. Here, we investigate the performance of a wide range of PRS models across four ancestry groups (Africans, Europeans, East Asians, and South Asians) for 14 conditions of high-medical interest.

METHODS

To select the best-performing model per trait, we first compared PRS performances for publicly available scores, and constructed new models using different methods (LDpred2, PRS-CSx and SNPnet). We used 285 K European individuals from the UK Biobank (UKBB) for training and 18 K, including diverse ancestries, for testing. We then evaluated PRS portability for the best models in Europeans and compared their accuracies with respect to the best PRS per ancestry. Finally, we validated the selected PRS models using an independent set of 8,417 individuals from Biobank of the Americas-Genomelink (BbofA-GL); and performed a PRS-Phewas.

RESULTS

We confirmed a decay in PRS performances relative to Europeans when the evaluation was conducted using the best-PRS model for Europeans (51.3% for South Asians, 46.6% for East Asians and 39.4% for Africans). We observed an improvement in the PRS performances when specifically selecting ancestry specific PRS models (phenotype variance increase: 1.62 for Africans, 1.40 for South Asians and 0.96 for East Asians). Additionally, when we selected the optimal model conditional on ancestry for CAD, HDL-C and LDL-C, hypertension, hypothyroidism and T2D, PRS performance for studied populations was more comparable to what was observed in Europeans. Finally, we were able to independently validate tested models for Europeans, and conducted a PRS-Phewas, identifying cross-trait interplay between cardiometabolic conditions, and between immune-mediated components.

CONCLUSION

Our work comprehensively evaluated PRS accuracy across a wide range of phenotypes, reducing the uncertainty with respect to which PRS model to choose and in which ancestry group. This evaluation has let us identify specific conditions where implementing risk-prioritization strategies could have practical utility across diverse ancestral groups, contributing to democratizing the implementation of PRS.

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.

Identification of consensus homozygous regions and their associations with growth and feed efficiency traits in American mink

The recent chromosome-based genome assembly and the newly developed 70K single nucleotide polymorphism (SNP) array for American mink (Neogale vison) facilitate the identification of genetic variants underlying complex traits in this species. The objective of this study was to evaluate the association between consensus runs of homozygosity (ROH) with growth and feed efficiency traits in American mink. A subsample of two mink populations (n = 2,986) were genotyped using the Affymetrix Mink 70K SNP array. The identified ROH segments were included simultaneously, concatenated into consensus regions, and the ROH-based association studies were carried out with linear mixed models considering a genomic relationship matrix for 11 growth and feed efficiency traits implemented in ASReml-R version 4. In total, 298,313 ROH were identified across all individuals, with an average length and coverage of 4.16 Mb and 414.8 Mb, respectively. After merging ROH segments, 196 consensus ROH regions were detected and used for genome-wide ROH-based association analysis. Thirteen consensus ROH regions were significantly (P < 0.01) associated with growth and feed efficiency traits. Several candidate genes within the significant regions are known for their involvement in growth and body size development, including MEF2A, ADAMTS17, POU3F2, and TYRO3. In addition, we found ten consensus ROH regions, defined as ROH islands, with frequencies over 80% of the population. These islands harbored 12 annotated genes, some of which were related to immune system processes such as DTX3L, PARP9, PARP14, CD86, and HCLS1. This is the first study to explore the associations between homozygous regions with growth and feed efficiency traits in American mink. Our findings shed the light on the effects of homozygosity in the mink genome on growth and feed efficiency traits, that can be utilized in developing a sustainable breeding program for mink.

Signatures of positive selection after the introduction of genomic selection in the Finnish Ayrshire population

The Finnish Ayrshire (FAY) belongs to the Nordic Red breeds and is characterized by high milk yield, high milk components, good fertility, and functional conformation. The FAY breeding program is based on genomic selection. Despite the benefits of selection on breeding values, autozygosity in the genome may increase due to selection, and increased autozygosity may cause inbreeding depression in selected traits. However, there is lack of studies concerning selection signatures in the FAY after genomic selection introduction. The aim of this study was to identify signatures of selection in FAY after the introduction of genomic selection. Genomic data included 45,834 SNPs. The genotyped animals were divided into 2 groups: animals born before genomic selection introduction (6,108 cows) and animals born after genomic selection introduction (47,361 cows). We identified the selection signatures using 3 complementary methods: 2 based on identification of selection signatures from runs of homozygosity (ROH) islands and one based on the decay of site-specific extended haplotype between populations at SNP sites (Rsb). In total, we identified 34 ROH islands on chromosomes 1, 3, 6, 8, 12-15, 17, 19, 22, and 26 in FAY animals born before genomic selection (between 1980 and 2011) and 30 ROH islands on chromosomes 1-3, 13-17, 22, and 25-26 in FAY animals born after genomic selection introduction (between 2015 and 2020). We additionally detected 22 ΔROH islands on chromosomes 2-3, 11, 13, 14, 16, 18, 20, and 25-26. Finally, a total of 31 Rsb regions on chromosomes 2, 3, 14, 18, 20, and 25 were identified. Based on the results, genomic selection has favored certain alleles and haplotypes on genomic regions related to traits relevant in the FAY breeding program: milk production, fertility, growth, beef production traits, and feed efficiency. Several genes related to these traits (e.g., PLA2G4A, MECR, CHUK, COX15, RICTOR, SHISA9, and SEMA4G) overlapped or partially overlapped the observed selection signature regions. The association of genotypes within these regions and their effects on traits relevant in the FAY breeding program should be studied and genetic regions undergoing selection monitored in the FAY population.

Does inbreeding contribute to pregnancy loss in Thoroughbred horses?

BACKGROUND

Excessive inbreeding increases the probability of uncovering homozygous recessive genotypes and has been associated with an increased risk of retained placenta and lower semen quality. No genomic analysis has investigated the association between inbreeding levels and pregnancy loss.

OBJECTIVES

To compare genetic inbreeding coefficients (F) of naturally occurring Thoroughbred Early Pregnancy Loss (EPLs), Mid and Late term Pregnancy Loss (MLPL) and Controls. The F value was hypothesised to be higher in cases of pregnancy loss (EPLs and MLPLs) than Controls.

STUDY DESIGN

Observational case-control study.

METHODS

Allantochorion and fetal DNA from EPL (n = 37, gestation age 14-65 days), MLPL (n = 94, gestational age 70 days-24 h post parturition) and Controls (n = 58) were genotyped on the Axiom Equine 670K SNP Genotyping Array. Inbreeding coefficients using Runs of Homozygosity (FROH) were calculated using PLINK software. ROHs were split into size categories to investigate the recency of inbreeding.

RESULTS

MLPLs had significantly higher median number of ROH (188 interquartile range [IQR], 180.8-197.3), length of ROH (3.10, IQR 2.93-3.33), and total number of ROH (590.8, IQR 537.3-632.3), and FROH (0.26, IQR 0.24-0.28) when compared with the Controls and the EPLs (p < 0.05). There was no significant difference in any of the inbreeding indices between the EPLs and Controls. The MLPLs had a significantly higher proportion of long (>10 Mb) ROH (2.5%, IQR 1.6-3.6) than the Controls (1.7%, IQR 0.6-2.5), p = 0.001. No unique ROHs were found in the EPL or MLPL populations.

MAIN LIMITATIONS

SNP-array data does not allow analysis of every base in the sequence.

CONCLUSIONS

This first study of the effect of genomic inbreeding levels on pregnancy loss showed that inbreeding is a contributor to MLPL, but not EPL in the UK Thoroughbred population. Mating choices remain critical, because inbreeding may predispose to MLPL by increasing the risk of homozygosity for specific lethal allele(s).

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.

Dissecting Selective Signatures and Candidate Genes in Grandparent Lines Subject to High Selection Pressure for Broiler Production and in a Local Russian Chicken Breed of Ushanka

Breeding improvements and quantitative trait genetics are essential to the advancement of broiler production. The impact of artificial selection on genomic architecture and the genetic markers sought remains a key area of research. Here, we used whole-genome resequencing data to analyze the genomic architecture, diversity, and selective sweeps in Cornish White (CRW) and Plymouth Rock White (PRW) transboundary breeds selected for meat production and, comparatively, in an aboriginal Russian breed of Ushanka (USH). Reads were aligned to the reference genome bGalGal1.mat.broiler.GRCg7b and filtered to remove PCR duplicates and low-quality reads using BWA-MEM2 and bcftools software; 12,563,892 SNPs were produced for subsequent analyses. Compared to CRW and PRW, USH had a lower diversity and a higher genetic distinctiveness. Selective sweep regions and corresponding candidate genes were examined based on ZFST, hapFLK, and ROH assessment procedures. Twenty-seven prioritized chicken genes and the functional projection from human homologs suggest their importance for selection signals in the studied breeds. These genes have a functional relationship with such trait categories as body weight, muscles, fat metabolism and deposition, reproduction, etc., mainly aligned with the QTLs in the sweep regions. This information is pivotal for further executing genomic selection to enhance phenotypic traits.

Introgression and disruption of migration routes have shaped the genetic integrity of wildebeest populations

The blue wildebeest (Connochaetes taurinus) is a keystone species in savanna ecosystems from southern to eastern Africa, and is well known for its spectacular migrations and locally extreme abundance. In contrast, the black wildebeest (C. gnou) is endemic to southern Africa, barely escaped extinction in the 1900s and is feared to be in danger of genetic swamping from the blue wildebeest. Despite the ecological importance of the wildebeest, there is a lack of understanding of how its unique migratory ecology has affected its gene flow, genetic structure and phylogeography. Here, we analyze whole genomes from 121 blue and 22 black wildebeest across the genus' range. We find discrete genetic structure consistent with the morphologically defined subspecies. Unexpectedly, our analyses reveal no signs of recent interspecific admixture, but rather a late Pleistocene introgression of black wildebeest into the southern blue wildebeest populations. Finally, we find that migratory blue wildebeest populations exhibit a combination of long-range panmixia, higher genetic diversity and lower inbreeding levels compared to neighboring populations whose migration has recently been disrupted. These findings provide crucial insights into the evolutionary history of the wildebeest, and tangible genetic evidence for the negative effects of anthropogenic activities on highly migratory ungulates.

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

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

Insights for the Captive Management of South China Tigers Based on a Large-Scale Genetic Survey

There is an urgent need to find a way to improve the genetic diversity of captive South China tiger (SCT, Panthera tigris amoyensis), the most critically endangered taxon of living tigers, facing inbreeding depression. The genomes showed that 13 hybrid SCTs from Meihuashan were divided into two groups; one group included three individuals who had a closer relationship with pureblood SCTs than another group. The three individuals shared more that 40% of their genome with pureblood SCTs and might be potential individuals for genetic rescuing in SCTs. A large-scale genetic survey based on 319 pureblood SCTs showed that the mean microsatellite inbreeding coefficient of pureblood SCTs decreased significantly from 0.1789 to 0.0600 (p = 0.000009) and the ratio of heterozygous loci increased significantly from 38.5% to 43.2% (p = 0.02) after one individual of the Chongqing line joined the Suzhou line and began to breed in the mid-1980s, which is a reason why the current SCTs keep a moderate level of microsatellite heterozygosity and nucleotide diversity. However, it is important to establish a back-up population based on the three individuals through introducing one pureblood SCT into the back-up population every year. The back-up population should be an important reserve in case the pureblood SCTs are in danger in the future.

Exploring the Effects of Robertsonian Translocation 1/29 (Rob (1;29)) on Genetic Diversity in Minor Breeds of Spanish Berrenda Cattle via Genome-Wide Analysis

Most of the previous studies on the genetic variability in Spanish "Berrenda" breeds have been carried out using DNA microsatellites. The present work aimed to estimate the genetic diversity, population structure, and potential genetic differences among individuals of both Berrenda breeds and groups based on the presence of the Robertsonian chromosomal translocation, rob (1;29). A total of 373 samples from animals belonging to the two breeds, including 169 cases diagnosed as rob (1;29)-positive, were genotyped using an SNP50K chip. The genetic diversity at the breed level did not show significant differences, but it was significantly lower in those subpopulations containing the rob (1;29). Runs of homozygosity identified a region of homozygosity on chromosome 6, where the KIT (KIT proto-oncogene, receptor tyrosine kinase) gene, which determines the typical spotted coat pattern in both breeds, is located. The four subpopulations considered showed minor genetic differences. The regions of the genome that most determined the differences between the breeds were observed on chromosomes 4, 6, 18, and 22. The presence of this Robertsonian translocation did not result in sub-structuring within each of the breeds considered. To improve the reproductive performance of Berrenda breeds, it would be necessary to implement strategies considering the involvement of potential breeding stock carrying rob (1;29).

Differential patterns in runs of homozygosity in two mice lines under divergent selection for environmental variability for birth weight

Runs of homozygosity (ROH) are defined as long continuous homozygous stretches in the genome which are assumed to originate from a common ancestor. It has been demonstrated that divergent selection for variability in mice is possible and that low variability in birth weight is associated with robustness. To analyse ROH patterns and ROH-based genomic inbreeding, two mouse lines that were divergently selected for birth weight variability for 26 generations were used, with: 752 individuals for the high variability line (H-Line), 766 individuals for the low variability line (L-Line) and 74 individuals as a reference population. Individuals were genotyped using the high density Affymetrix Mouse Diversity Genotyping Array. ROH were identified using both the sliding windows (SW) and the consecutive runs (CR) methods. Inbreeding coefficients were calculated based on pedigree (FPED ) information, on ROH identified using the SW method (FROHSW ) and on ROH identified using the CR method (FROHCR ). Differences in genomic inbreeding were not consistent across generations and these parameters did not show clear differences between lines. Correlations between FPED and FROH were high, particularly for FROHSW . Moreover, correlations between FROHSW and FPED were even higher when ROH were identified with no restrictions in the number of heterozygotes per ROH. The comparison of FROH estimates between either of the selected lines were based on significant differences at the chromosome level, mainly in chromosomes 3, 4, 6, 8, 11, 15 and 19. ROH-based inbreeding estimates that were computed using longer homozygous segments had a higher relationship with FPED . Differences in robustness between lines were not attributable to a higher homozygosis in the L-Line, but maybe to the different distribution of ROH at the chromosome level between lines. The analysis identified a set of genomic regions for future research to establish the genomic basis of robustness.

Genome-Wide Detection for Runs of Homozygosity in Baoshan Pigs Using Whole Genome Resequencing

Baoshan pigs (BS) are a local breed in Yunnan Province that may face inbreeding owing to its limited population size. To accurately evaluate the inbreeding level of the BS pig population, we used whole-genome resequencing to identify runs of homozygosity (ROH) regions in BS pigs, calculated the inbreeding coefficient based on pedigree and ROH, and screened candidate genes with important economic traits from ROH islands. A total of 22,633,391 SNPS were obtained from the whole genome of BS pigs, and 201 ROHs were detected from 532,450 SNPS after quality control. The number of medium-length ROH (1-5 Mb) was the highest (98.43%), the number of long ROH (>5 Mb) was the lowest (1.57%), and the inbreeding of BS pigs mainly occurred in distant generations. The inbreeding coefficient FROH, calculated based on ROH, was 0.018 ± 0.016, and the FPED, calculated based on the pedigree, was 0.027 ± 0.028, which were positively correlated. Forty ROH islands were identified, containing 507 genes and 891 QTLs. Several genes were associated with growth and development (IGFALS, PTN, DLX5, DKK1, WNT2), meat quality traits (MC3R, ACSM3, ECI1, CD36, ROCK1, CACNA2D1), and reproductive traits (NPW, TSHR, BMP7). This study provides a reference for the protection and utilization of BS pigs.

Characterization of runs of Homozygosity revealed genomic inbreeding and patterns of selection in indigenous sahiwal cattle

Runs of homozygosity (ROH) are contiguous genomic regions, homozygous across all sites which arise in an individual due to the parents transmitting identical haplotypes to their offspring. The genetic improvement program of Sahiwal cattle after decades of selection needs re-assessment of breeding strategy and population phenomena. Hence, the present study was carried out to optimize input parameters in PLINK for ROH estimates, to explore ROH islands and assessment of pedigree and genome-based inbreeding in Sahiwal cattle. The sliding window approach with parameters standardized to define ROH for the specific population under study was used for the identification of runs. The optimum maximum gap, density, window-snp and window-threshold were 250 Kb, 120 Kb/SNP, 10, 0.05 respectively and ROH patterns were also characterized. ROH islands were defined as the short homozygous genomic regions shared by a large proportion of individuals in a population, containing significantly higher occurrences of ROH than the population specific threshold level. These were identified using the -homozyg-group function of the PLINK v1.9 program. Our results indicated that the Islands of ROH harbor a few candidate genes, ACAD11, RFX4, BANP, UBA5 that are associated with major economic traits. The average FPED (Pedigree based inbreeding coefficient), FROH (Genomic inbreeding coefficient), FHOM (Inbreeding estimated as the ratio of observed and expected homozygous genotypes), FGRM (Inbreeding estimated on genomic relationship method) and FGRM0.5 (Inbreeding estimated from the diagonal of a GRM with allele frequencies near to 0.5) were 0.009, 0.091, 0.035, -0.104 and -0.009, respectively. Our study revealed the optimum parameter setting in PLINK viz. maximal gaps between two SNPs, minimal density of SNPs in a segment (in kb/SNP) and scanning window size to identify ROH segments, which will enable ROH estimation more efficient and comparable across various SNP genotyping-based studies. The result further emphasized the significant role of genomics in unraveling population diversity, selection signatures and inbreeding in the ongoing Sahiwal breed improvement programs.

SNP heterozygosity, relatedness and inbreeding of whole genomes from the isolated population of the Faroe Islands

BACKGROUND

The population of the Faroe Islands is an isolated population but very little is known about it from whole genome sequencing. The population of about 50000 people has a high incidence of rare diseases e.g., 1:300 for Primary Carnitine Deficiency. A screening programme was implemented, and eleven persons were also whole genome sequenced at x37 coverage for diagnostic purposes of those cases that were not affected by the known mutations. The purpose of our study is to utilize the high coverage data to explore the genomic variation and the ancestral history of the population. We study the SNP heterozygosity, the pairwise relatedness from kinship, the inbreeding from runs of homozygosity ROH, and we find the minor allele frequency distribution. We estimate the population ancestry and the timing of the founding event by using the whole genomes from eight consenting individuals.

RESULTS

We find the number of SNPs and the heterozygosity for the eight individual samples, and for merged samples, for which we also study the relatedness. We find close relatedness between the supposedly unrelated individuals. From ROH, we interpret the high relatedness as an ancient property of the isolated population. A bottleneck event is estimated starting between years [Formula: see text] with a maximum consanguineous population in year [Formula: see text] and similarly consanguineous between years [Formula: see text]. The ancestry analysis shows the population descends from founders of [Formula: see text] European and [Formula: see text] Admixed American ancestry. A distinct clustering near the central European and British populations of the 1000 Genome Project is likely the result of the population isolation and genetic drift. The minor allele frequency distribution suggests many rare variants.

CONCLUSIONS

The ancestry is mainly European while the inbreeding is higher compared to European populations and population isolates. The Faroese population has inbreeding more like ancient Europeans. We discovered a bottlenecked and consanguineous population event and estimated it starting in the 1st-4th century as compared to the oldest archaeological findings from the 4th-6th century.

Whole-genome SNP markers reveal runs of homozygosity in indigenous cattle breeds of Pakistan

The runs of homozygosity (ROH) were identified in 14 Pakistani cattle breeds (n = 105) by genotyping with the Illumina 50 K SNP BeadChip. These breeds were categorized into Dairy, Dual, and Draft breeds based on their utility and production performance. We identified a total of 10,936 ROHs which mainly consisted of a high number of shorter segments (1-4 Mb). Dairy group exhibited the highest level of inbreeding (FROH: 0.078 ± 0.028) while the lowest (FROH: 0.002 ± 0.008) was observed in Dual group. In 48 genomic regions identified with a high frequency of ROH, 207 genes were detected in the three breed groups. A substantially higher number of ROH islands detected in dairy breeds indicated the impact of the positive selection pressure over the years. Important candidate genes and QTL were detected in the ROH islands associated with economic traits like milk production, reproduction, meat, carcass, and health traits in dairy cattle.

Identification of homozygosity-rich regions in the Holstein genome

In this study, 371 Holstein cows from six herds and 26 Holstein bulls, which were used in these herds, were genotyped by the Illumina BovineSNP50 array. For runs of homozygosity (ROH) identification, consecutive and sliding runs were performed by the detectRUNS and Plink software. The missing calls did not significantly affect the ROH data. The mean number of ROH identified by consecutive runs was 95.4 ± 2.7, and that by sliding runs was 86.0 ± 2.6 in cows, while this number for Holstein bulls was lower 58.9 ± 1.9. The length of the ROH segments varied from 1 Mb to over 16 Mb, with the largest number of ROH having a length of 1-2 Mb. Of the 29 chromosomes, BTA 14, BTA 16, and BTA 7 were the most covered by ROH. The mean coefficient of inbreeding across the herds was 0.111 ± 0.003 and 0.104 ± 0.004 based on consecutive and sliding runs, respectively, and 0.078 ± 0.005 for bulls based on consecutive runs. These values do not exceed those for Holstein cattle in North America. The results of this study confirmed the more accurate identification of ROH by consecutive runs, and also that the number of allowed heterozygous SNPs may have a significant effect on ROH data.

Genetic characterisation of the Connemara pony and the Warmblood horse using a within-breed clustering approach

BACKGROUND

The Connemara pony (CP) is an Irish breed that has experienced varied selection by breeders over the last fifty years, with objectives ranging from the traditional hardy pony to an agile athlete. We compared these ponies with well-studied Warmblood (WB) horses, which are also selectively bred for athletic performance but with a much larger census population. Using genome-wide single nucleotide polymorphism (SNP) and whole-genome sequencing data from 116 WB (94 UK WB and 22 European WB) and 36 CP (33 UK CP and 3 US CP), we studied the genomic diversity, inbreeding and population structure of these breeds.

RESULTS

The k-means clustering approach divided both the CP and WB populations into four genetic groups, among which the CP genetic group 1 (C1) associated with non-registered CP, C4 with US CP, WB genetic group 1 (W1) with Holsteiners, and W3 with Anglo European and British WB. Maximum and mean linkage disequilibrium (LD) varied significantly between the two breeds (mean from 0.077 to 0.130 for CP and from 0.016 to 0.370 for WB), but the rate of LD decay was generally slower in CP than WB. The LD block size distribution peaked at 225 kb for all genetic groups, with most of the LD blocks not exceeding 1 Mb. The top 0.5% harmonic mean pairwise fixation index (FST) values identified ontology terms related to cancer risk when the four CP genetic groups were compared. The four CP genetic groups were less inbred than the WB genetic groups, but C2, C3 and C4 had a lower proportion of shorter runs of homozygosity (ROH) (74 to 76% < 4 Mb) than the four WB genetic groups (80 to 85% < 4 Mb), indicating more recent inbreeding. The CP and WB genetic groups had a similar ratio of effective number of breeders (Neb) to effective population size (Ne).

CONCLUSIONS

Distinct genetic groups of individuals were revealed within each breed, and in WB these genetic groups reflected population substructure better than studbook or country of origin. Ontology terms associated with immune and inflammatory responses were identified from the signatures of selection between CP genetic groups, and while CP were less inbred than WB, the evidence pointed to a greater degree of recent inbreeding. The ratio of Neb to Ne was similar in CP and WB, indicating the influence of popular sires is similar in CP and WB.

Genome-Wide Single-Nucleotide Polymorphism-Based Genomic Diversity and Runs of Homozygosity for Selection Signatures in Equine Breeds

The horse, one of the most domesticated animals, has been used for several purposes, like transportation, hunting, in sport, or for agriculture-related works. Kathiawari, Marwari, Manipuri, Zanskari, Bhutia, Spiti, and Thoroughbred are the main breeds of horses, particularly due to their agroclimatic adaptation and role in any kind of strong physical activity, and these characteristics are majorly governed by genetic factors. The genetic diversity and phylogenetic relationship of these Indian equine breeds using microsatellite markers have been reported, but further studies exploring the SNP diversity and runs of homozygosity revealing the selection signature of breeds are still warranted. In our study, the identification of genes that play a vital role in muscle development is performed through SNP detection via the whole-genome sequencing approach. A total of 96 samples, categorized under seven breeds, and 620,721 SNPs were considered to ascertain the ROH patterns amongst all the seven breeds. Over 5444 ROH islands were mined, and the maximum number of ROHs was found to be present in Zanskari, while Thoroughbred was confined to the lowest number of ROHs. Gene enrichment of these ROH islands produced 6757 functional genes, with AGPAT1, CLEC4, and CFAP20 as important gene families. However, QTL annotation revealed that the maximum QTLs were associated with Wither's height trait ontology that falls under the growth trait in all seven breeds. An Equine SNP marker database (EqSNPDb) was developed to catalogue ROHs for all these equine breeds for the flexible and easy chromosome-wise retrieval of ROH along with the genotype details of all the SNPs. Such a study can reveal breed divergence in different climatic and ecological conditions.

Homozygosity mapping in the Kazakh national dog breed Tazy

The Tazy is a breed of sighthound common in Kazakhstan. The identification of runs of homozygosity (ROH) is an informative approach to assessing the history and possible patterns of directional selection pressure. To our knowledge, the present study is the first to provide an overview of the ROH pattern in the Tazy dogs from a genome-wide perspective. The ROH of the Tazy was found to be mainly composed of shorter segments (1-2 Mb), accounting for approximately 67% of the total ROH. The estimated ROH-based inbreeding coefficients (F_ROH) ranged from 0.028 to 0.058 with a mean of 0.057. Five genomic regions under positive selection were identified on chromosomes 18, 22, and 25. The regions on chromosomes 18 and 22 may be breed specific, while the region on chromosome 22 overlaps with regions of hunting traits in other hunting dog breeds. Among the 12 candidate genes located in these regions, the gene CAB39L may be a candidate that affects running speed and endurance of the Tazy dog. Eight genes could belong to an evolutionarily conserved complex as they were clustered in a large protein network with strong linkages. The results may enable effective interventions when incorporated into conservation planning and selection of the Tazy breed.

High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources

Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.

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 (F_ROH ) was 0.1 in BY_P, whereas the lowest F_ROH was 0.0743 in BY_C. In addition, the inbreeding coefficient of BY_R2 (F_ROH  = 0.0798) was higher than that of BY_R1 (F_ROH  = 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.

Pedigree and genome-based patterns of homozygosity in the South African Ayrshire, Holstein, and Jersey breeds

The erosion of genetic diversity limits long-term genetic gain and impedes the sustainability of livestock production. In the South African (SA) dairy industry, the major commercial dairy breeds have been applying estimated breeding values (EBVs) and/or have been participating in Multiple Across Country Evaluations (MACE). The transition to genomic estimated breeding values (GEBVs) in selection strategies requires monitoring of the genetic diversity and inbreeding of current genotyped animals, especially considering the comparatively small population sizes of global dairy breeds in SA. This study aimed to perform a homozygosity-based evaluation of the SA Ayrshire (AYR), Holstein (HST), and Jersey (JER) dairy cattle breeds. Three sources of information, namely 1) single nucleotide polymorphism (SNP) genotypes (3,199 animals genotyped for 35,572 SNPs) 2) pedigree records (7,885 AYR; 28,391 HST; 18,755 JER), and 3) identified runs of homozygosity (ROH) segments were used to quantify inbreeding related parameters. The lowest pedigree completeness was for the HST population reducing from a value of 0.990 to 0.186 for generation depths of one to six. Across all breeds, 46.7% of the detected ROH were between 4 megabase pairs (Mb) and 8 Mb in length. Two conserved homozygous haplotypes were identified in more than 70% of the JER population on Bos taurus autosome (BTA) 7. The JER breed displayed the highest level of inbreeding across all inbreeding coefficients. The mean (± standard deviation) pedigree-based inbreeding coefficient (FPED) ranged from 0.051 (±0.020) for AYR to 0.062 (±0.027) for JER, whereas SNP-based inbreeding coefficients (FSNP) ranged from 0.020 (HST) to 0.190 (JER) and ROH-based inbreeding coefficients, considering all ROH segment coverage (FROH), ranged from 0.053 (AYR) to 0.085 (JER). Within-breed Spearman correlations between pedigree-based and genome-based estimates ranged from weak (AYR: 0.132 between FPED and FROH calculated for ROH &lt;4Mb in size) to moderate (HST: 0.584 between FPED and FSNP). Correlations strengthened between FPED and FROH as the ROH length category was considered lengthened, suggesting a dependency on breed-specific pedigree depth. The genomic homozygosity-based parameters studied proved useful in investigating the current inbreeding status of reference populations genotyped to implement genomic selection in the three most prominent South African dairy cattle breeds.

Spatial genetic structure of European wild boar, with inferences on late-Pleistocene and Holocene demographic history

European wildlife has been subjected to intensifying levels of anthropogenic impact throughout the Holocene, yet the main genetic partitioning of many species is thought to still reflect the late-Pleistocene glacial refugia. We analyzed 26,342 nuclear SNPs of 464 wild boar (Sus scrofa) across the European continent to infer demographic history and reassess the genetic consequences of natural and anthropogenic forces. We found that population fragmentation, inbreeding and recent hybridization with domestic pigs have caused the spatial genetic structure to be heterogeneous at the local scale. Underlying local anthropogenic signatures, we found a deep genetic structure in the form of an arch-shaped cline extending from the Dinaric Alps, via Southeastern Europe and the Baltic states, to Western Europe and, finally, to the genetically diverged Iberian peninsula. These findings indicate that, despite considerable anthropogenic influence, the deeper, natural continental structure is still intact. Regarding the glacial refugia, our findings show a weaker signal than generally assumed, but are nevertheless suggestive of two main recolonization routes, with important roles for Southern France and the Balkans. Our results highlight the importance of applying genomic resources and framing genetic results within a species' demographic history and geographic distribution for a better understanding of the complex mixture of underlying processes.

Estimating inbreeding rate and effective population size in the Finnish Ayrshire population in the era of genomic selection

Genomic selection has been applied in dairy cattle breeding over the last decade. Using genomic information may speed up genetic gain as breeding values can be predicted reasonably accurately directly after birth. However, genetic diversity may decrease if the inbreeding rate per generation increases and the effective population size decreases. Despite many positive qualities of the Finnish Ayrshire, for example, high average protein yield and fertility, over time the breed has lost its place as the most common dairy breed in Finland. Thus, maintaining the genetic variability of the breed is becoming more important. The aim of our research was to estimate the impact of genomic selection on inbreeding rate and effective population size using both pedigree and genomic data. The genomic data included 46,914 imputed single nucleotide polymorphism (SNP) variants from 75,038 individuals, and the pedigree data included 2,770,025 individuals. All animals in the data were born between 2000 and 2020. Genomic inbreeding coefficients were estimated as the proportion of SNPs in runs of homozygosity (ROH) out of the total number of SNPs. The inbreeding rate was estimated by regressing the mean genomic inbreeding coefficients on birth years. Effective population size was then estimated based on the inbreeding rate. Additionally, effective population size was estimated from the mean increase in individual inbreeding using pedigree data. Introduction of genomic selection was assumed to have taken place gradually; years 2012-2014 were treated as a transition period from the traditional phenotype-based breeding value estimation to genomic-based estimation. The median length of the identified homozygous segments was 5.5 Mbp, and a slight increase in the proportion of segments over 10 Mbp was observed after 2010. The inbreeding rate decreased from 2000 to 2011 and subsequently increased slightly. The pedigree- and genomic-based estimates of inbreeding rate were similar to each other. The estimates of effective population size based on the regression method were very sensitive to the number of years considered; thus, the estimates were not very reliable. The effective population size estimated from the mean increase in individual inbreeding reached its highest value of 160 in 2011 and decreased to 150 after that. In addition, the generation interval in the sire path has decreased from 5.5 years to 3.5 years after genomic selection was implemented. Based on our results, after the implementation of genomic selection, the proportion of long ROH stretches has increased, the generation interval in the sire path has decreased, the inbreeding rate has increased and the effective population size has decreased. However, the effective population size is still at a good level, allowing for an efficient selection scheme in the Finnish Ayrshire breed.

The coefficients of inbreeding revealed by ROH study among inbred individuals belonging to each type of the first cousin marriage: A preliminary report from North India

BACKGROUND

Genome-wide runs of homozygosity (ROH) are appropriate to estimate genomic inbreeding, determine population history, unravel the genetic architecture of complex traits and disorders.

OBJECTIVE

The study sought to investigate and compare the actual proportion of homozygosity or autozygosity in the genomes of progeny of four subtypes of first cousin mating in humans, using both pedigree and genomic measures for autosomes and sex chromosomes.

METHODS

For this purpose, Illumina Global Screening Array-24 v1.0 BeadChip followed by cyto-ROH analysis through Illumina Genome Studio was used to characterise the homozygosity in five participants from North Indian state (Uttar Pradesh). PLINK v.1.9 software was used to estimate the genomic inbreeding coefficients viz. ROH-based inbreeding estimate (F_ROH) and homozygous loci-based inbreeding estimate (F_HOM).

RESULTS

A total of 133 ROH segments were detected with maximum number and genomic coverage in Matrilateral Parallel (MP) type and minimum in outbred individual. ROH pattern revealed that MP type has a higher degree of homozygosity than other subtypes. The comparison of F_ROH, F_HOM, and pedigree-based inbreeding estimate (F_PED) showed some difference in theoretical and realised proportion of homozygosity for sex-chromosomal loci but not for autosome for each type of consanguinity.

CONCLUSIONS

This is the very first study to compare and estimate the pattern of homozygosity among the kindreds of first cousin unions. However, a greater number of individuals from each type of marriage is required for statistical inference of no difference between theoretical and realized homozygosity among different degrees of inbreeding prevalent in humans worldwide.

Genome-wide copy number variant screening of Saudi schizophrenia patients reveals larger deletions in cases versus controls

Introduction

Genome-wide association studies have discovered common polymorphisms in regions associated with schizophrenia. No genome-wide analyses have been performed in Saudi schizophrenia subjects.

Methods

Genome-wide genotyping data from 136 Saudi schizophrenia cases and 97 Saudi controls in addition to 4,625 American were examined for copy number variants (CNVs). A hidden Markov model approach was used to call CNVs.

Results

CNVs in schizophrenia cases were twice as large on average than CNVs in controls (p = 0.04). The analyses focused on extremely large >250 kilobases CNVs or homozygous deletions of any size. One extremely large deletion was noted in a single case (16.5 megabases on chromosome 10). Two cases had an 814 kb duplication of chromosome 7 spanning a cluster of genes, including circadian-related loci, and two other cases had 277 kb deletions of chromosome 9 encompassing an olfactory receptors gene family. CNVs were also seen in loci previously associated with schizophrenia, namely a 16p11 proximal duplication and two 22q11.2 deletions.

Discussion

Runs of homozygosity (ROHs) were analyzed across the genome to investigate correlation with schizophrenia risk. While rates and sizes of these ROHs were similar in cases and controls, we identified 10 regions where multiple cases had ROHs and controls did not.

Pedigree reconstruction and population structure using SNP markers in Gir cattle

Our objective was to establish a SNPs panel for pedigree reconstruction using microarrays of different densities and evaluate the genomic relationship coefficient of the inferred pedigree, in addition to analyzing the population structure based on genomic analyses in Gir cattle. For parentage analysis and genomic relationship, 16,205 genotyped Gir animals (14,458 females and 1747 males) and 1810 common markers to the four SNP microarrays were used. For population structure analyses, including linkage disequilibrium, effective population size, and runs of homozygosity (ROH), genotypes from 21,656 animals were imputed. Likelihood ratio (LR) approach was used to reconstruct the pedigree, deepening the pedigree and showing it is well established in terms of recent information. Coefficients for each relationship category of the inferred pedigree were adequate. Linkage disequilibrium showed rapid decay. We detected a decrease in the effective population size over the last 50 generations, with the average generation interval around 9.08 years. Higher ROH-based inbreeding coefficient in a class of short ROH segments, with moderate to high values, was also detected, suggesting bottlenecks in the Gir genome. Breeding strategies to minimize inbreeding and avoid massive use of few proven sires with high genetic value are suggested to maintain genetic variability in future generations. In addition, we recommend reducing the generation interval to maximize genetic progress and increase effective population size.

Facial attractiveness is only weakly linked to genome-wide heterozygosity

Introduction

It has been frequently suggested that overall genomic heterozygosity and, particularly, heterozygosity of loci on the so-called major histocompatibility complex (MHC), which is responsible for the recognition of foreign substances/ pathogens and the recognition of self and non-self, is associated with better health and better resistance to infections and parasites. It has further been speculated that such a potentially beneficial heterozygosity can be detected through body odor and facial attractiveness.

Methods

On the basis of genome wide SNP data (713,014 SNPs) of participants from the Wisconsin Longitudinal Study, we therefore investigated whether homozygosity either on the MHC (measured as inbreeding coefficient) or genome-wide (measured as runs of homozygosity and inbreeding coefficient) is associated with rated facial attractiveness.

Results

Although we found that the genome-wide average length of homozygous segments and the genome-wide inbreeding coefficient are significantly negatively associated with some measures of facial attractiveness, if corrected for multiple testing, any significant association was no longer formally significant after correction. In addition, the variance in facial attractiveness explained by the genome wide homozygosity is very low (<0.15%). We did not find any significant association between the inbreeding coefficient on the MHC and facial attractiveness.

Discussion

We only find a weak association of genome- wide heterozygosity and facial attractiveness.

High-density SNP markers elucidate the genetic divergence and population structure of Noticiana sheep breed in the Mediterranean context

Among livestock species, sheep have played an early major role in the Mediterranean area. Italy has a long history of sheep breeding and, despite a dramatic contraction in numbers, still raise several local populations that may represent a unique source of genetic diversity. The Noticiana is a breed of the south-eastern part of Sicily appreciated both for its dairy products and for its resistance to harsh environment. In this study, the high-density Illumina Ovine SNP600K BeadChip array was used for the first genome-wide characterization of 48 individuals of Noticiana sheep to investigate its diversity, the genome structure and the relationship within the context of worldwide and Italian breeds. Moreover, the runs of homozygosity (ROH) pattern and the pairwise F_ST-outliers were examined. Noticiana reported moderate levels of genetic diversity. The high percentage of short and medium length ROH segments (93% under 4 Mb) is indicative of a within breed relatedness dating back to ancient times, despite the absence of management for the mating plans and the reduced population size. In the worldwide context, the Southern Italian, Spanish and Albanian breeds overlapped in a macro cluster which also included the Noticiana sheep. The results highlighted ancestral genetic components of Noticiana shared with Comisana breed, and showed the clear separation from the other Italian sheep. This is likely the consequence of the combined effects of genetic drift, small population size and reproductive isolation. ROH islands and F_ST-outliers approaches in Noticiana identified genes and QTLs involved in milk and meat production, as well as related to the local adaptation, and therefore are consistent with the phenotypic traits of the studied breed. Although a wider sampling could be useful to deepen the genomic survey on Noticiana, these results represent a crucial starting point for the characterization of an important local genetic resource, with a view of supporting the local economy and preserving the biodiversity of the sheep species.

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.

Fine-tuning genomic and pedigree inbreeding rates in equine population with a deep and reliable stud book: the case of the Pura Raza Española horse

BACKGROUND

Estimating inbreeding, which is omnipresent and inevitable in livestock populations, is a primary goal for management and animal breeding especially for those interested in mitigating the negative consequences of inbreeding. Inbreeding coefficients have been historically estimated by using pedigree information; however, over the last decade, genome-base inbreeding coefficients have come to the forefront in this field. The Pura Raza Española (PRE) horse is an autochthonous Spanish horse breed which has been recognised since 1912. The total PRE population (344,718 horses) was used to estimate Classical (F), Ballou's ancestral, Kalinowski's ancestral, Kalinowski's new and the ancestral history coefficient values. In addition, genotypic data from a selected population of 805 PRE individuals was used to determine the individual inbreeding coefficient using SNP-by-SNP-based techniques (methods of moments -F_HOM-, the diagonal elements of the genomic -F_G-, and hybrid matrixes -F_H-) and ROH measures (F_RZ). The analyse of both pedigree and genomic based inbreeding coefficients in a large and robust population such as the PRE horse, with proven parenteral information for the last 40 years and a high degree of completeness (over 90% for the last 70 years) will allow us to understand PRE genetic variability better and the correlations between the estimations will give the data greater reliability.

RESULTS

The mean values of the pedigree-based inbreeding coefficients ranged from 0.01 (F for the last 3 generations -F3-) to 0.44 (ancestral history coefficient) and the mean values of genomic-based inbreeding coefficients varied from 0.05 (F_RZ for three generations, F_H and F_HOM) to 0.11 (F_RZ for nine generations). Significant correlations were also found between pedigree and genomic inbreeding values, which ranged between 0.58 (F3 with F_HOM) and 0.79 (F with F_RZ). In addition, the correlations between F_RZ estimated for the last 20 generations and the pedigree-based inbreeding highlight the fact that fewer generations of genomic data are required when comparing total inbreeding values, and the opposite when ancient values are calculated.

CONCLUSIONS

Ultimately, our results show that it is still useful to work with a deep and reliable pedigree in pedigree-based genetic studies with very large effective population sizes. Obtaining a satisfactory parameter will always be desirable, but the approximation obtained with a robust pedigree will allow us to work more efficiently and economically than with massive genotyping.

Tracking footprints of artificial and natural selection signatures in breeding and non-breeding cats

Stray non-breeding cats (stray) represent the largest heterogeneous cat population subject to natural selection, while populations of the Siamese (SIAM) and Oriental Shorthair (OSH) breeds developed through intensive artificial selection for aesthetic traits. Runs of homozygosity (ROH) and demographic measures are useful tools to discover chromosomal regions of recent selection and to characterize genetic diversity in domestic cat populations. To achieve this, we genotyped 150 stray and 26 household non-breeding cats (household) on the Illumina feline 63 K SNP BeadChip and compared them to SIAM and OSH. The 50% decay value of squared correlation coefficients (r²) in stray (0.23), household (0.25), OSH (0.24) and SIAM (0.25) corresponded to a mean marker distance of 1.12 Kb, 4.55 Kb, 62.50 Kb and 175.07 Kb, respectively. The effective population size (N_e) decreased in the current generation to 55 in stray, 11 in household, 9 in OSH and 7 in SIAM. In the recent generation, the increase in inbreeding per generation (ΔF) reached its maximum values of 0.0090, 0.0443, 0.0561 and 0.0710 in stray, household, OSH and SIAM, respectively. The genomic inbreeding coefficient (F_ROH) based on ROH was calculated for three length categories. The F_ROH was between 0.014 (F_ROH60) and 0.020 (F_ROH5) for stray, between 0.018 (F_ROH60) and 0.024 (F_ROH5) for household, between 0.048 (F_ROH60) and 0.069 (F_ROH5) for OSH and between 0.053 (F_ROH60) and 0.073 (F_ROH5) for SIAM. We identified nine unique selective regions for stray through genome-wide analyses for regions with reduced heterozygosity based on F_ST statistics. Genes in these regions have previously been associated with reproduction (BUB1B), motor/neurological behavior (GPHN, GABRB3), cold-induced thermogenesis (DIO2, TSHR), immune system development (TSHR), viral carcinogenesis (GTF2A1), host immune response against bacteria, viruses, chemoattractant and cancer cells (PLCB2, BAHD1, TIGAR), and lifespan and aging (BUB1B, FGF23). In addition, we identified twelve unique selective regions for OSH containing candidate genes for a wide range of coat colors and patterns (ADAMTS20, KITLG, TYR, TYRO3-a MITF regulator, GPNMB, FGF7, RAB38) as well as congenital heart defects (PDE4D, PKP2) and gastrointestinal disorders (NLGN1, ALDH1B1). Genes in stray that represent unique selective events indicate, at least in part, natural selection for environmental adaptation and resistance to infectious disease, and should be the subject of future research. Stray cats represent an important genetic resource and have the potential to become a research model for disease resistance and longevity, which is why we recommend preserving semen before neutering.

Integrative QTL mapping and selection signatures in Groningen White Headed cattle inferred from whole-genome sequences

Here, we aimed to identify and characterize genomic regions that differ between Groningen White Headed (GWH) breed and other cattle, and in particular to identify candidate genes associated with coat color and/or eye-protective phenotypes. Firstly, whole genome sequences of 170 animals from eight breeds were used to evaluate the genetic structure of the GWH in relation to other cattle breeds by carrying out principal components and model-based clustering analyses. Secondly, the candidate genomic regions were identified by integrating the findings from: a) a genome-wide association study using GWH, other white headed breeds (Hereford and Simmental), and breeds with a non-white headed phenotype (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel, Dutch Belted, and Holstein Friesian); b) scans for specific signatures of selection in GWH cattle by comparison with four other Dutch traditional breeds (Dutch Friesian, Deep Red, Meuse-Rhine-Yssel and Dutch Belted) and the commercial Holstein Friesian; and c) detection of candidate genes identified via these approaches. The alignment of the filtered reads to the reference genome (ARS-UCD1.2) resulted in a mean depth of coverage of 8.7X. After variant calling, the lowest number of breed-specific variants was detected in Holstein Friesian (148,213), and the largest in Deep Red (558,909). By integrating the results, we identified five genomic regions under selection on BTA4 (70.2-71.3 Mb), BTA5 (10.0-19.7 Mb), BTA20 (10.0-19.9 and 20.0-22.7 Mb), and BTA25 (0.5-9.2 Mb). These regions contain positional and functional candidate genes associated with retinal degeneration (e.g., CWC27 and CLUAP1), ultraviolet protection (e.g., ERCC8), and pigmentation (e.g. PDE4D) which are probably associated with the GWH specific pigmentation and/or eye-protective phenotypes, e.g. Ambilateral Circumocular Pigmentation (ACOP). Our results will assist in characterizing the molecular basis of GWH phenotypes and the biological implications of its adaptation.

Genomic insights into the conservation and population genetics of two Chinese native goat breeds

Chinese goats are an important group of goats worldwide. However, there are few studies on the conservation priority, genetic relationship, and potential gene flow between Chinese and global goat breeds. Here, we genotyped 239 goats from conservation populations of the Chinese Guangfeng and Ganxi breeds using the GoatSNP50 BeadChip. The conservation priority, population structure, selection signatures and introgression of these goats were analyzed in the context of 36 global goat breeds. First, we showed that Guangfeng and Ganxi goats had the largest effective population sizes across the global breeds 13 generations ago. Nevertheless, Ganxi goats have recently experienced a high degree of inbreeding, resulting in their conservation priority based on total gene and allelic diversities being lower than that of most other Chinese breeds (including Guangfeng goats). Population structure and admixture analyses showed that an average of 18% of Guangfeng genomic components were introgressed from Boer goats approximately 18-yr ago. Next, we reconstructed the subfamily structure of the core populations of Guangfeng and Ganxi goats, and proposed reasonable conservation strategies for inbreeding management. Moreover, a list of candidate genes under selection for fertility, immunity, growth, and meat quality were detected in Guangfeng and Ganxi goats. Finally, we identified some genes related to body development and reproduction, which were introgressed from Boer goats and may be beneficial for improving performance and productivity of Guangfeng goats. In conclusion, this study not only provides new insights into the conservation and utilization of Guangfeng and Ganxi goats but also enriches our understanding of artificial introgression from exotic goats into Chinese local goats.

Challenges in quantifying genome erosion for conservation

Massive defaunation and high extinction rates have become characteristic of the Anthropocene. Genetic effects of population decline can lead populations into an extinction vortex, where declining populations show lower genetic fitness, in turn leading to lower populations still. The lower genetic fitness in a declining population due to a shrinking gene pool is known as genetic erosion. Three different types of genetic erosion are highlighted in this review: overall homozygosity, genetic load and runs of homozygosity (ROH), which are indicative of inbreeding. The ability to quantify genetic erosion could be a very helpful tool for conservationists, as it can provide them with an objective, quantifiable measure to use in the assessment of species at risk of extinction. The link between conservation status and genetic erosion should become more apparent. Currently, no clear correlation can be observed between the current conservation status and genetic erosion. However, the high quantities of genetic erosion in wild populations, especially in those species dealing with habitat fragmentation and habitat decline, may be early signs of deteriorating populations. Whole genome sequencing data is the way forward to quantify genetic erosion. Extra screening steps for genetic load and hybridization can be included, since they could potentially have great impact on population fitness. This way, the information yielded from genetic sequence data can provide conservationists with an objective genetic method in the assessment of species at risk of extinction. However, the great complexity of genome erosion quantification asks for consensus and bridging science and its applications, which remains challenging.

Genomic characterization of autozygosity and recent inbreeding trends in all major breeds of US dairy cattle

Maintaining a genetically diverse dairy cattle population is critical to preserving adaptability to future breeding goals and avoiding declines in fitness. This study characterized the genomic landscape of autozygosity and assessed trends in genetic diversity in 5 breeds of US dairy cattle. We analyzed a sizable genomic data set containing 4,173,679 pedigreed and genotyped animals of the Ayrshire, Brown Swiss, Guernsey, Holstein, and Jersey breeds. Runs of homozygosity (ROH) of 2 Mb or longer in length were identified in each animal. The within-breed means for number and the combined length of ROH were highest in Jerseys (62.66 ± 8.29 ROH and 426.24 ± 83.40 Mb, respectively; mean ± SD) and lowest in Ayrshires (37.24 ± 8.27 ROH and 265.05 ± 85.00 Mb, respectively). Short ROH were the most abundant, but moderate to large ROH made up the largest proportion of genome autozygosity in all breeds. In addition, we identified ROH islands in each breed. This revealed selection patterns for milk production, productive life, health, and reproduction in most breeds and evidence for parallel selective pressure for loci on chromosome 6 between Ayrshire and Brown Swiss and for loci on chromosome 20 between Holstein and Jersey. We calculated inbreeding coefficients using 3 different approaches, pedigree-based (F_PED), marker-based using a genomic relationship matrix (F_GRM), and segment-based using ROH (F_ROH). The average inbreeding coefficient ranged from 0.06 in Ayrshires and Brown Swiss to 0.08 in Jerseys and Holsteins using F_PED, from 0.22 in Holsteins to 0.29 in Guernsey and Jerseys using F_GRM, and from 0.11 in Ayrshires to 0.17 in Jerseys using F_ROH. In addition, the effective population size at past generations (5-100 generations ago), the yearly rate of inbreeding, and the effective population size in 3 recent periods (2000-2009, 2010-2014, and 2015-2018) were determined in each breed to ascertain current and historical trends of genetic diversity. We found a historical trend of decreasing effective population size in the last 100 generations in all breeds and breed differences in the effect of the recent implementation of genomic selection on inbreeding accumulation.

African wild dogs (Lycaon pictus) from the Kruger National Park, South Africa are currently not inbred but have low genomic diversity

African wild dogs (Lycaon pictus) have undergone severe population reductions and are listed as endangered on the International Union for Conservation of Nature Red List. Small, isolated populations have the potential to suffer from threats to their genetic diversity that may impact species viability and future survival. This study provides the first set of population-wide genomic data to address conservation concerns for this endangered species. Whole genome sequencing data were generated for 71 free-ranging African wild dogs from the Kruger National Park (KNP), South Africa, and used to estimate important population genomic parameters. Genomic diversity metrics revealed that variation levels were low; however, this African wild dog population showed low levels of inbreeding. Very few first- and second-order relationships were observed in this cohort, with most relationships falling into the third-order or distant category. Patterns of homozygosity could have resulted from historical inbreeding or a loss in genome variation due to a population bottleneck. Although the results suggest that this stronghold African wild dog population maintains low levels of inbreeding, likely due to their cooperative breeding system, it may lead to a continuous population decline when a reduced number of suitable mates are available. Consequently, the low genomic variation may influence species viability over time. This study highlights the importance of assessing population genomic parameters to set conservation priorities. Future studies should include the investigation of the potential of this endangered species to adapt to environmental changes considering the low genomic diversity in this population.

Conservation status and historical relatedness of South African communal indigenous goat populations using a genome-wide single-nucleotide polymorphism marker

Indigenous goats form the majority of populations in smallholder, low input, low output production systems and are considered an important genetic resource due to their adaptability to different production environments and support of communal farming. Effective population size (N_e), inbreeding levels, and the runs of homozygosity (ROHs) are effective tools for exploring the genetic diversity and understanding the demographic history in efforts to support breeding strategies to use and conserve genetic resources. Across populations, the current N_e of Gauteng was the lowest at 371 animals, while the historical N_e across populations suggests that the ancestor N_e has decreased by 53.86%, 44.58%, 42.16%, and 41.16% in Free State (FS), North West (NW), Limpopo (LP), and Gauteng (GP), respectively, over the last 971 generations. Genomic inbreeding levels related to ancient kinship (F_ROH > 5 Mb) were highest in FS (0.08 ± 0.09) and lowest in the Eastern Cape (EC) (0.02 ± 0.02). A total of 871 ROH island regions which include important environmental adaptation and hermo-tolerance genes such as IL10RB, IL23A, FGF9, IGF1, EGR1, MTOR, and MAPK3 were identified (occurring in over 20% of the samples) in FS (n = 37), GP (n = 42), and NW (n = 2) populations only. The mean length of ROH across populations was 7.76 Mb and ranged from 1.61 Mb in KwaZulu-Natal (KZN) to 98.05 Mb (GP and NW). The distribution of ROH according to their size showed that the majority (n = 1949) of the detected ROH were > 5 Mb in length compared to the other categories. Assuming two hypothetical ancestral populations, the populations from KZN and LP are revealed, supporting PC 1. The genomes of KZN and LP share a common origin but have substantial admixture from the EC and NW populations. The findings revealed that the occurrence of high N_e and autozygosity varied largely across breeds in communal indigenous goat populations at recent and ancient events when a genome-wide single-nucleotide polymorphism (SNP) marker was used. The use of Illumina goat SNP50K BeadChip shows that there was a migration route of communal indigenous goat populations from the northern part (LP) of South Africa to the eastern areas of the KZN that confirmed their historical relatedness and coincides with the migration periods of the Bantu nation.

Simultaneous testing of rule- and model-based approaches for runs of homozygosity detection opens up a window into genomic footprints of selection in pigs

BACKGROUND

Past selection events left footprints in the genome of domestic animals, which can be traced back by stretches of homozygous genotypes, designated as runs of homozygosity (ROHs). The analysis of common ROH regions within groups or populations displaying potential signatures of selection requires high-quality SNP data as well as carefully adjusted ROH-defining parameters. In this study, we used a simultaneous testing of rule- and model-based approaches to perform strategic ROH calling in genomic data from different pig populations to detect genomic regions under selection for specific phenotypes.

RESULTS

Our ROH analysis using a rule-based approach offered by PLINK, as well as a model-based approach run by RZooRoH demonstrated a high efficiency of both methods. It underlined the importance of providing a high-quality SNP set as input as well as adjusting parameters based on dataset and population for ROH calling. Particularly, ROHs ≤ 20 kb were called in a high frequency by both tools, but to some extent covered different gene sets in subsequent analysis of ROH regions common for investigated pig groups. Phenotype associated ROH analysis resulted in regions under potential selection characterizing heritage pig breeds, known to harbour a long-established breeding history. In particular, the selection focus on fitness-related traits was underlined by various ROHs harbouring disease resistance or tolerance-associated genes. Moreover, we identified potential selection signatures associated with ear morphology, which confirmed known candidate genes as well as uncovered a missense mutation in the ABCA6 gene potentially supporting ear cartilage formation.

CONCLUSIONS

The results of this study highlight the strengths and unique features of rule- and model-based approaches as well as demonstrate their potential for ROH analysis in animal populations. We provide a workflow for ROH detection, evaluating the major steps from filtering for high-quality SNP sets to intersecting ROH regions. Formula-based estimations defining ROHs for rule-based method show its limits, particularly for efficient detection of smaller ROHs. Moreover, we emphasize the role of ROH detection for the identification of potential footprints of selection in pigs, displaying their breed-specific characteristics or favourable phenotypes.

Genetic and Genomic Characterization of a New Beef Cattle Composite Breed (Purunã) Developed for Production in Pasture-Based Systems

Purunã is a composite beef cattle breed, developed in Southern Brazil by crossing the Angus, Charolais, Canchim, and Caracu breeds. The goal of this study was to perform the first genetic characterization of the Purunã breed, based on both pedigree and genomic information. For this, 100 randomly selected animals were genotyped, and 11,205 animals born from 1997 to 2019 had pedigree information. The genetic analyses performed were principal component analysis, admixture, phylogenic tree, pedigree and genomic inbreeding, linkage disequilibrium (LD), effective population size (Ne), consistency of the gametic phase, runs of homozygosity (ROH), heterozygosity-enriched regions (HERs), and functional analyses of the ROH and HER regions identified. Our findings indicate that Purunã is more genetically related to the Charolais, Canchim, and Angus breeds than Caracu or Nellore. The levels of inbreeding were shown to be small based on all the metrics evaluated and ranged from −0.009 to 0.029. A low (−0.12–0.31) correlation of the pedigree-based inbreeding compared to all the genomic inbreeding coefficients evaluated was observed. The LD average was 0.031 (±0.0517), and the consistency of the gametic phase was shown to be low for all the breed pairs, ranging from 0.42 to 0.27 to the distance of 20 Mb. The Ne values based on pedigree and genomic information were 158 and 115, respectively. A total of 1,839 ROHs were found, and the majority of them are of small length (&lt;4 Mb). An important homozygous region was identified on BTA5 with pathways related to behavioral traits (sensory perception, detection of stimulus, and others), as well as candidate genes related to heat tolerance (MY O 1A), feed conversion rate (RDH5), and reproduction (AMDHD1). A total of 1,799 HERs were identified in the Purunã breed with 92.3% of them classified within the 0.5–1 Mb length group, and 19 HER islands were identified in the autosomal genome. These HER islands harbor genes involved in growth pathways, carcass weight (SDCBP), meat and carcass quality (MT2A), and marbling deposition (CISH). Despite the genetic relationship between Purunã and the founder breeds, a multi-breed genomic evaluation is likely not feasible due to their population structure and low consistency of the gametic phase among them.

Genomic status of yellow-breasted bunting following recent rapid population decline

Global biodiversity is facing serious threats. However, knowledge of the genomic consequences of recent rapid population declines of wild organisms is limited. Do populations experiencing recent rapid population decline have the same genomic status as wild populations that experience long-term declines? Yellow-breasted Bunting (Emberiza aureola) is a critically endangered species that has been experiencing a recent rapid population decline. To answer the question, we assembled and annotated the whole genome of Yellow-breasted Bunting. Furthermore, we found high genetic diversity, low linkage disequilibrium, and low proportion of long runs of homozygosity in Yellow-breasted Bunting, suggesting that the populations following recent rapid declines have different genomic statuses from the population that experienced long-term population decline.

High-Density Genomic Characterization of Native Croatian Sheep Breeds

A recent comprehensive genomic analysis based on 50K SNP profiles has shown that the regional Balkan sheep populations have considerable genetic overlap but are distinctly different from surrounding breeds. All eight Croatian sheep breeds were represented by a small number of individuals per breed. Here, we genotyped 220 individuals representing the native Croatian sheep breeds (Istrian Sheep, Krk Island Sheep, Cres Island Sheep, Rab Island Sheep, Lika Pramenka, Pag Island Sheep, Dalmatian Pramenka, Dubrovnik Sheep) and mouflon using the Ovine Infinium® HD SNP BeadChip (606,006 SNPs). In addition, we included publicly available Balkan Pramenka and other Mediterranean sheep breeds. Our analyses revealed the complex population structure of Croatian sheep breeds and their origin and geographic barriers (island versus mainland). Migration patterns confirmed the historical establishment of breeds and the pathways of gene flow. Inbreeding coefficients (FROH&gt;2 Mb) between sheep populations ranged from 0.025 to 0.070, with lower inbreeding coefficients observed in Dalmatian Pramenka and Pag Island Sheep and higher inbreeding in Dubrovnik sheep. The estimated effective population size ranged from 61 to 1039 for Krk Island Sheep and Dalmatian Pramenka, respectively. Higher inbreeding levels and lower effective population size indicate the need for improved conservation management to maintain genetic diversity in some breeds. Our results will contribute to breeding and conservation strategies of native Croatian sheep breeds.

Identification of Runs of Homozygosity Islands and Genomic Estimated Inbreeding Values in Caqueteño Creole Cattle (Colombia)

The Caqueteño Creole (CAQ) is a native breed of cattle from the Caquetá department (Colombia), adapted to tropical conditions, which is extremely important to production systems in those regions. However, CAQ is poorly studied. In this sense, population structure studies associated with runs of homozygosity (ROH) analysis would allow for a better understanding of CAQ. Through ROH analysis, it is possible to reveal genetic relationships between individuals, measure genome inbreeding levels, and identify regions associated with traits of economic interest. Samples from a CAQ population (n = 127) were genotyped with the Bovine HD BeadChip (777,000 SNPs) and analyzed with the PLINK 1.9 program to estimate FROH and ROH islands. We highlighted a decrease in inbreeding frequency for FROH 4−8 Mb, 8−16 Mb, and >16 Mb classes, indicating inbreeding control in recent matings. We also found genomic hotspot regions on chromosomes 3, 5, 6, 8, 16, 20, and 22, where chromosome 20 harbored four hotspots. Genes in those regions were associated with fertility and immunity traits, muscle development, and environmental resistance, which may be present in the CAQ breed due to natural selection. This indicates potential for production systems in tropical regions. However, further studies are necessary to elucidate the CAQ production objective.

Demographic history differences between Hispanics and Brazilians imprint haplotype features

Admixture is known to greatly impact the genetic landscape of a population and, while genetic variation underlying human phenotypes has been shown to differ among populations, studies on admixed subjects are still scarce. Latin American populations are the result of complex demographic history, such as 2 or 3-way admixing events, bottlenecks and/or expansions, and adaptive events unique to the American continent. To explore the impact of these events on the genetic structure of Latino populations, we evaluated the following haplotype features: linkage disequilibrium, shared identity by descent segments, runs of homozygosity, and extended haplotype homozygosity (integrated haplotype score) in Latinos represented in the 1000 Genome Project along with array data from 171 Brazilians sampled in the South and Southeast regions of Brazil. We found that linkage disequilibrium decay relates to the amount of American and African ancestry. The extent of identity by descent sharing positively correlates with historical effective population sizes, which we found to be steady or growing, except for Puerto Ricans and Colombians. Long runs of homozygosity, a particular instance of autozygosity, was only enriched in Peruvians and Native Americans. We used simulations to account for random sampling and linkage disequilibrium to filter positive selection indexes and found 244 unique markers under selection, 26 of which are common to 2 or more populations. Some markers exhibiting positive selection signals had estimated time to the most recent common ancestor consistent with human adaptation to the American continent. In conclusion, Latino populations present highly divergent haplotype characteristics that impact genetic architecture and underlie complex phenotypes.

Genome-wide detection of genetic structure and runs of homozygosity analysis in Anhui indigenous and Western commercial pig breeds using PorcineSNP80k data

Background

Runs of homozygosity (ROH) are continuous homozygous regions typically located in the DNA sequence of diploid organisms. Identifications of ROH that lead to reduced performance can provide valuable insight into the genetic architecture of complex traits. Here, we systematically investigated the population genetic structure of five Anhui indigenous pig breeds (AHIPs), and compared them to those of five Western commercial pig breeds (WECPs). Furthermore, we examined the occurrence and distribution of ROHs in the five AHIPs and estimated the inbreeding coefficients based on the ROHs (F_ROH) and homozygosity (F_HOM). Finally, we identified genomic regions with high frequencies of ROHs and annotated candidate genes contained therein.

Results

The WECPs and AHIPs were clearly differentiated into two separate clades consistent with their geographical origins, as revealed by the population structure and principal component analysis. We identified 13,530 ROHs across all individuals, of which 4,555 and 8,975 ROHs were unique to AHIPs and WECPs, respectively. Most ROHs identified in our study were short (< 10 Mb) or medium (10–20 Mb) in length. WECPs had significantly higher numbers of short ROHs, and AHIPs generally had longer ROHs. F_ROH values were significantly lower in AHIPs than in WECPs, indicating that breed improvement and conservation programmes were successful in AHIPs. On average, F_ROH and F_HOM values were highly correlated (0.952–0.991) in AHIPs and WECPs. A total of 27 regions had a high frequency of ROHs and contained 17 key candidate genes associated with economically important traits in pigs. Among these, nine candidate genes (CCNT2, EGR2, MYL3, CDH13, PROX1, FLVCR1, SETD2, FGF18, and FGF20) found in WECPs were related to muscular and skeletal development, whereas eight candidate genes (CSN1S1, SULT1E1, TJP1, ZNF366, LIPC, MCEE, STAP1, and DUSP) found in AHIPs were associated with health, reproduction, and fatness traits.

Conclusion

Our findings provide a useful reference for the selection and assortative mating of pig breeds, laying the groundwork for future research on the population genetic structures of AHIPs, ultimately helping protect these local varieties.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08583-9.

Genomic integration to identify molecular biomarkers associated with indicator traits of gastrointestinal nematode resistance in sheep

This study aimed to integrate GWAS and structural variants to propose possible molecular biomarkers related to gastrointestinal nematode resistance traits in Santa Inês sheep. The phenotypic records FAMACHA, haematocrit, white blood cell count, red blood cell count, haemoglobin, platelets and egg counts per gram of faeces were collected from 700 naturally infected animals, belonging to four Brazilian flocks. A total of 576 animals were genotyped using the Ovine SNP12k BeadChip and were imputed using a reference population with Ovine SNP50 BeadChip. The GWAS approaches were based on SNPs, haplotypes, CNVs and ROH. The overlapping between the significant genomic regions detected from all approaches was investigated, and the results were integrated using a network analysis. Genes related to the immune system were found, such as ABCB1, IL6, WNT5A and IRF5. Genomic regions containing candidate genes and metabolic pathways involved in immune responses, inflammatory processes and immune cells affecting parasite resistance traits were identified. The genomic regions, biological processes and candidate genes uncovered could lead to biomarkers for selecting more resilient sheep and improving herd welfare and productivity. The results obtained are the start point to identify molecular biomarkers related to indicator traits of gastrointestinal nematode resistance in Santa Inês sheep.