Conservation genomic analysis of domestic and wild pig populations from the Iberian Peninsula

Hybridization in an isolated population of blesbok and red hartebeest

Hybridization in antelope species has been widely reported in South African national parks and provincial reserves as well as on private land due to anthropogenic effects. In a closed management setting, hybridization may occur due to the crossbreeding of closely related species with unequal sex ratios, resulting in either sterile or fertile offspring. In this study, we used molecular techniques to evaluate the risk of anthropogenic hybridization between blesbok (Damaliscus pygargus phillipsi) and red hartebeest (Alcelaphus buselaphus caama) in an isolated group that purposely included the two species with unequal sex ratios (one red hartebeest male and 19 male and female blesbok). Genetic analysis based on microsatellites confirmed the presence of seven hybrid individuals. Mitochondrial analysis verified that hybridization occurred between blesbok females and the red hartebeest male. STRUCTURE and NEWHYBRIDS classified the hybrids as F1. It is suspected that the hybrid individuals were sterile as the males had undeveloped testes and only F1 hybrids were detected. Thus, the risk of hybridization between these two species may be limited in the wild. In captive settings, genetic monitoring should be included in management plans for blesbok and red hartebeest to ensure that the long-term consequences of wasted reproductive effort are limited.

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.

Functional and evolutionary analysis of host Synaptogyrin-2 in porcine circovirus type 2 susceptibility

Mammalian evolution has been influenced by viruses for millions of years, leaving signatures of adaptive evolution within genes encoding for viral interacting proteins. Synaptogyrin-2 (SYNGR2) is a transmembrane protein implicated in promoting bacterial and viral infections. A genome-wide association study of pigs experimentally infected with porcine circovirus type 2b (PCV2b) uncovered a missense mutation (SYNGR2 p.Arg63Cys) associated with viral load. In this study, CRISPR/Cas9-mediated gene editing of the porcine kidney 15 (PK15, wtSYNGR2+p.63Arg) cell line generated clones homozygous for the favorable SYNGR2 p.63Cys allele (emSYNGR2+p.63Cys). Infection of edited clones resulted in decreased PCV2 replication compared to wildtype PK15 (P<0.05), with consistent effects across genetically distinct PCV2b and PCV2d isolates. Sequence analyses of wild and domestic pigs (n>700) revealed the favorable SYNGR2 p.63Cys allele is unique to domestic pigs and more predominant in European than Asian breeds. A haplotype defined by the SYNGR2 p.63Cys allele was likely derived from an ancestral haplotype nearly fixed within European (0.977) but absent from Asian wild boar. We hypothesize that the SYNGR2 p.63Cys allele arose post-domestication in ancestral European swine. Decreased genetic diversity in homozygotes for the SYNGR2 p.63Cys allele compared to SYNGR2 p.63Arg, corroborates a rapid increase in frequency of SYGNR2 p.63Cys via positive selection. Signatures of adaptive evolution across mammalian species were also identified within SYNGR2 intraluminal loop domains, coinciding with the location of SYNGR2 p.Arg63Cys. Therefore, SYNGR2 may reflect a novel component of the host-virus evolutionary arms race across mammals with SYNGR2 p.Arg63Cys representing a species-specific example of putative adaptive evolution.

Characterization analysis of Rongchang pig population based on the Zhongxin-1 Porcine Breeding Array PLUS

OBJECTIVE

To carry out a comprehensive production planning of the existing Rongchang pig population from both environmental and genetic aspects, and to establish a closed population with stable genetic diversity and strict pathogen control, it is necessary to fully understand the genetic background of the population.

METHODS

We genotyped 54 specific pathogen free (SPF) Rongchang pigs using the Zhongxin-1 Porcine Breeding Array PLUS, calculated their genetic diversity parameters and constructed their families. In addition, we also counted the runs of homozygosity (ROH) of each individual and calculated the value of inbreeding coefficient based on ROH for each individual.

RESULTS

Firstly, the results of genetic diversity analysis showed that the effective population size (Ne) of this population was 3.2, proportion of polymorphic markers (PN) was 0.515, desired heterozygosity (He) and observed heterozygosity (Ho) were 0.315 and 0.335. Ho was higher than He, indicating that the heterozygosity of all the selected loci was high. Secondly, combining the results of genomic relatedness analysis and cluster analysis, it was found that the existing Rongchang pig population could be divided into four families. Finally, we also counted the ROH of each individual and calculated the inbreeding coefficient value accordingly, whose mean value was 0.09.

CONCLUSION

Due to the limitation of population size and other factors, the genetic diversity of this Rongchang pig population is low. The results of this study can provide basic data to support the development of Rongchang pig breeding program, the establishment of SPF Rongchang pig closed herd and its experimental utilization.

Identification of breed-specific genomic variants in Colombian Creole pig breeds by whole-genome sequencing

Dissecting genetic variation of local breeds is important for the success of conservation. In this research, we investigated the genomic variation of Colombian Creole (CR) pigs, with a focus on the breed-specific variants in the exonic region of 34 genes with reported effects on adaptive and economic traits. Seven individuals of each of the three CR breeds (CM, Casco de Mula; SP, San Pedreño; and ZU, Zungo) were whole-genome sequenced along with 7 Iberian (IB) pigs and 7 pigs of each of the four most used cosmopolitan (CP) breeds (Duroc, Landrace × Large White, and Pietrain). Molecular variability in CR (6,451,218 variants; from 3,919,242, in SP, to 4,648,069, in CM) was comparable to that in CP, but higher than in IB. For the investigated genes, SP pigs displayed less exonic variants (178) than ZU (254), CM (263), IB (200), and the individual CP genetic types (201 to 335). Sequence variation in these genes confirmed the resemblance of CR to IB and indicates that CR pigs, particularly ZU and CM, are not exempt from selective introgression of other breeds. A total of 50 exonic variants were identified as being potentially specific to CR, including a high-impact deletion in the intron between exons 15 and 16 of the leptin receptor gene, which was only found in CM and ZU. The identification of breed-specific variants in genes related to adaptive and economical traits can bolster the understanding of the role of gene-environment interactions on local adaptation and points the way for effective breeding and conservation of CR pigs.

Application potential of chicken DNA chip in domestic pigeon species - Preliminary results

Introducing the SNP technology to pigeon breeding will enhance the competitiveness of a sector that produces one of the healthiest and best quality meats. The present study aimed to test the applicability of the Illumina Chicken_50K_CobbCons array on 24 domestic pigeon individuals from the Mirthys hybrids and Racing pigeon breeds. A total of 53,313 SNPs were genotyped. Principal component analysis shows a significant overlap between the two groups. The chip performed poorly in this data set, with a call rate per sample of 0.474 (49%). The low call rate was likely due to an increase in the evolutionary distance. A total of 356 SNPs were retained after a relatively strict quality control. We have demonstrated that it is technically feasible to use a chicken microarray chip on pigeon samples. Presumably, with a larger sample size and by assigning phenotypic data, efficiency would be improved, allowing more thorough analyses, such as genome-wide association studies.

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.

Expression Profile of Housekeeping Genes and Tissue-Specific Genes in Multiple Tissues of Pigs

Pigs have become an ideal model system for human disease research and development and an important farm animal that provides a valuable source of nutrition. To profile the all-sided gene expression and their biological functions across multiple tissues, we conducted a comprehensive analysis of gene expression on a large scale around the side of housekeeping genes (HKGs), tissue specific genes (TSGs), and the co-expressed genes in 14 various tissues. In this study, we identified 2351 HKGs and 3018 TSGs across tissues, among which 4 HKGs (COX1, UBB, OAZ1/NPFF) exhibited low variation and high expression levels, and 31 particular TSGs (e.g., PDC, FKBP6, STAT2, and COL1A1) were exclusively expressed in several tissues, including endocrine brain, ovaries, livers, backfat, jejunum, kidneys, lungs, and longissimus dorsi muscles. We also obtained 17 modules with 230 hub genes (HUBGs) by weighted gene co-expression network analysis. On the other hand, HKGs functions were enriched in the signaling pathways of the ribosome, spliceosome, thermogenesis, oxidative phosphorylation, and nucleocytoplasmic transport, which have been highly suggested to involve in the basic biological tissue activities. While TSGs were highly enriched in the signaling pathways that were involved in specific physiological processes, such as the ovarian steroidogenesis pathway in ovaries and the renin-angiotensin system pathway in kidneys. Collectively, these stable, specifical, and co-expressed genes provided useful information for the investigation of the molecular mechanism for an understanding of the genetic and biological processes of complex traits in pigs.

Lost pigs of Angola: Whole genome sequencing reveals unique regions of selection with emphasis on metabolism and feed efficiency

Angola, in the western coast of Africa, has been through dramatic social events that have led to the near-disappearance of native swine populations, and the recent introduction of European exotic breeds has also contributed to the erosion of this native swine repertoire. In an effort to investigate the genetic basis of native pigs in Angola (ANG) we have generated whole genomes from animals of a remote local pig population in Huambo province, which we have compared with 78 genomes of European and Asian pig breeds as well as European and Asian wild boars that are currently in public domain. Analyses of population structure showed that ANG pigs grouped within the European cluster and were clearly separated from Asian pig breeds. Pairwise F_ST ranged from 0.14 to 0.26, ANG pigs display lower levels of genetic differentiation towards European breeds. Finally, we have identified candidate regions for selection using a complementary approach based on various methods. All results suggest that selection towards feed efficiency and metabolism has occurred. Moreover, all analysis identified CDKAL1 gene, which is related with insulin and cholesterol metabolism, as a candidate gene overlapping signatures of selection unique to ANG pigs. This study presents the first assessment of the genetic relationship between ANG pigs and other world breeds and uncovers selection signatures that may indicate adaptation features unique to this important genetic resource.

Optimization of the Thawing Protocol for Iberian Boar Sperm

Simple Summary

Limited attention is paid to sperm thawing protocols, and their study could be relevant to endangered species or breeds, especially for cryopreserved material present in the existing gene banks. The aim of this study was to determine the most optimal thawing protocol for Iberian boar sperm through testing different thawing rates and different modifications of the thawing extender. Based on the findings, the most optimal results were obtained thawing at 70 °C for 8 s with the inclusion of cyclodextrins loaded with cholesterol (CLC) in the extender, revealing the importance of adapting the thawing protocols.

Abstract

Thawing protocols have been barely studied, and their modifications may lead to a substantial improvement in post-thawing sperm quality, which could be of great relevance to existing sperm banks, such as those for Iberian pig breeds with varieties in danger of extinction. For that, the study aimed to evaluate different thawing rates and to evaluate modifications in the composition of the thawing extender (basic pH to 8–8.2, incorporation of cyclodextrins loaded with cholesterol [CLC] and the incorporation of ion chelators [EDTA and EGTA]). After thawing, overall sperm motility and kinematic parameters, acrosome status and sperm membrane integrity were evaluated. The most optimal results were obtained with the thawing rate reaching 70 °C for 8 s with the inclusion of 12.5 mg of CLC/500 × 10⁶ spermatozoa in the thawing extender, which showed an improvement compared to the control at 70 °C. In conclusion, to adapt the thawing conditions may be relevant, especially for endangered species or breeds such as some varieties of Iberian pig, since this process could also be used in samples cryopreserved in gene banks.

Genomic inbreeding and runs of homozygosity analysis of indigenous cattle populations in southern China

Runs of homozygosity (ROH) are continuous homozygous segments from the common ancestor of parents. Evaluating ROH pattern can help to understand inbreeding level and genetic basis of important traits. In this study, three representative cattle populations including Leiqiong cattle (LQC), Lufeng cattle (LFC) and Hainan cattle (HNC) were genotyped using the Illumina BovineHD SNPs array (770K) to assess ROH pattern at genome wide level. Totally, we identified 26,537 ROH with an average of 153 ROH per individual. The sizes of ROH ranged from 0.5 to 53.26Mb, and the average length was 1.03Mb. The average of F_ROH ranged from 0.10 (LQC) to 0.15 (HNC). Moreover, we identified 34 ROH islands (with frequency > 0.5) across genome. Based on these regions, we observed several breed-specific candidate genes related to adaptive traits. Several common genes related to immunity (TMEM173, MZB1 and SIL1), and heat stress (DNAJC18) were identified in all three populations. Three genes related to immunity (UGP2), development (PURA) and reproduction (VPS54) were detected in both HNC and LQC. Notably, we identified several breed-specific genes related to sperm development (BRDT and SPAG6) and heat stress (TAF7) in HNC, and immunity (CDC23 and NME5) and development (WNT87) in LFC. Our findings provided valuable insights into understanding the genomic homozygosity pattern and promoting the conservation of genetic resources of Chinese indigenous cattle.

Runs of homozygosity in Sable Island feral horses reveal the genomic consequences of inbreeding and divergence from domestic breeds

BACKGROUND

Understanding inbreeding and its impact on fitness and evolutionary potential is fundamental to species conservation and agriculture. Long stretches of homozygous genotypes, known as runs of homozygosity (ROH), result from inbreeding and their number and length can provide useful population-level information on inbreeding characteristics and locations of signatures of selection. However, the utility of ROH for conservation is limited for natural populations where baseline data and genomic tools are lacking. Comparing ROH metrics in recently feral vs. domestic populations of well understood species like the horse could provide information on the genetic health of those populations and offer insight into how such metrics compare between managed and unmanaged populations. Here we characterized ROH, inbreeding coefficients, and ROH islands in a feral horse population from Sable Island, Canada, using ~41 000 SNPs and contrasted results with those from 33 domestic breeds to assess the impacts of isolation on ROH abundance, length, distribution, and ROH islands.

RESULTS

ROH number, length, and ROH-based inbreeding coefficients (F_ROH) in Sable Island horses were generally greater than in domestic breeds. Short runs, which typically coalesce many generations prior, were more abundant than long runs in all populations, but run length distributions indicated more recent population bottlenecks in Sable Island horses. Nine ROH islands were detected in Sable Island horses, exhibiting very little overlap with those found in domestic breeds. Gene ontology (GO) enrichment analysis for Sable Island ROH islands revealed enrichment for genes associated with 3 clusters of biological pathways largely associated with metabolism and immune function.

CONCLUSIONS

This study indicates that Sable Island horses tend to be more inbred than their domestic counterparts and that most of this inbreeding is due to historical bottlenecks and founder effects rather than recent mating between close relatives. Unique ROH islands in the Sable Island population suggest adaptation to local selective pressures and/or strong genetic drift and highlight the value of this population as a reservoir of equine genetic variation. This research illustrates how ROH analyses can be applied to gain insights into the population history, genetic health, and divergence of wild or feral populations of conservation concern.

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.

Identification of Homozygous Regions With Adverse Effects on the Five Economic Traits of Duroc Pigs

Runs of homozygosity (ROH) are widely used to estimate genomic inbreeding, which is linked to inbreeding depression on phenotypes. However, the adverse effects of specific homozygous regions on phenotypic characteristics are rarely studied in livestock. In this study, the 50 K SNP data of 3,770 S21 Duroc (American origin) and 2,096 S22 Duroc (Canadian origin) pigs were used to investigate the harmful ROH regions on five economic traits. The results showed that the two Duroc lines had different numbers and distributions of unfavorable ROHs, which may be related to the different selection directions and intensities between the two lines. A total of 114 and 58 ROH segments were found with significant adverse effects on the economic traits of S21 and S22 pigs, respectively. Serval pleiotropic ROHs were detected to reduce two or multiple phenotypic performances in two Duroc populations. Candidate genes in these shared regions were mainly related to growth, fertility, immunity, and fat deposition. We also observed that some ROH genotypes may cause opposite effects on different traits. This study not only enhances our understanding of the adverse effects of ROH on phenotypes, but also indicates that ROH information could be incorporated into breeding programs to estimate and control the detrimental effects of homozygous regions.

Genetic Differentiation among Livestock Breeds—Values for Fst

Simple Summary

The degree of relationship among livestock breeds can be quantified by the F_st statistic, which measures the extent of genetic differentiation between them. An F_st value of 0.1 has often been taken as indicating that two breeds are indeed genetically distinct, but this concept has not been evaluated critically. Here, F_st values have been collated for the six major livestock species: cattle, sheep, goats, pigs, horses, and chickens. These values are remarkably variable both within and between species, demonstrating that F_st > 0.1 is not a reliable criterion for breed distinctiveness. However, the large body of F_st data accumulated in the last 20–30 years represents an untapped database that could contribute to the development of interdisciplinary research involving livestock breeds.

Abstract

(1) Background: The F_st statistic is widely used to characterize between-breed relationships. F_st = 0.1 has frequently been taken as indicating genetic distinctiveness between breeds. This study investigates whether this is justified. (2) Methods: A database was created of 35,080 breed pairs and their corresponding F_st values, deduced from microsatellite and SNP studies covering cattle, sheep, goats, pigs, horses, and chickens. Overall, 6560 (19%) of breed pairs were between breeds located in the same country, 7395 (21%) between breeds of different countries within the same region, 20,563 (59%) between breeds located far apart, and 562 (1%) between a breed and the supposed wild ancestor of the species. (3) Results: General values for between-breed F_st were as follows, cattle: microsatellite 0.06–0.12, SNP 0.08–0.15; sheep: microsatellite 0.06–0.10, SNP 0.06–0.17; horses: microsatellite 0.04–0.11, SNP 0.08–0.12; goats: microsatellite 0.04–0.14, SNP 0.08–0.16; pigs: microsatellite 0.06–0.27, SNP 0.15–0.22; chickens: microsatellite 0.05–0.28, SNP 0.08–0.26. (4) Conclusions: (1) Large amounts of F_st data are available for a substantial proportion of the world’s livestock breeds, (2) the value for between-breed F_st of 0.1 is not appropriate owing to its considerable variability, and (3) accumulated F_st data may have value for interdisciplinary research.

Dental formula variations in wild and domestic Sus scrofa: is the first premolar agenesis an evolutionary trend?

We investigated the occurrence of congenital dental abnormalities in 94 wild boars Sus scrofa scrofa, and 20 Iberian pigs Sus scrofa domesticus from Doñana (DNP) and Sierra de Andújar Natural Park (ANP) southern Spain. The only dental variation found was agenesis of the lower first premolar, either deciduous (Dp1) or permanent (P1). We analyzed prevalence variations using: odds ratios, Fisher’s exact test and logistic-regression models. The lower first premolar has often been reported to be absent in modern domestic breeds of pigs, but it is usually reported as present in wild boar. However, we found a similar occurrence of agenesis in wild boar as in Iberian pigs. A common genotypic background between wild boar and Iberian pig populations of ANP sharing the same ecosystems may explain this result. When considering only wild boar, unexpected differences between populations were detected. The wild boar from ANP showed higher estimated probabilities of agenesis than those from DNP. Environmental conditions can also represent influencing factors for changes within subspecies. Our findings reinforce the idea that in suids the congenital absence of premolars probably has an evolutionary relevance, which could represent a generalized trend towards the reduction of the dental formula from the plesiomorphic placental dentition.

Genetic diversity of the Central European wild boar (Sus scrofa scrofa) population and domestic pig (Sus scrofa domesticus) breeds based on a microsatellite DNA locus

The results of studies of the genetic structure of the Central European wild boar (Sus scrofa scrofa) population and four breeds of domestic pigs (Duroc, Yorkshire, Large White and Landrace) bred in the Central Black Earth region of Russia are presented in this work. Based on 12 microsatellite loci, a significant ( p <0.05) decrease in the level of genetic variability in bred breeds was shown. The expected heterozygosity and Shannon index were as follows: in the wild boar, Ho = 0.763 ± 0.026, I = 1.717 ± 0.091; in the Duroc breed, Ho = 0.569 ± 0.068, I = 1.191 ± 0.157; in the Landrace, Ho = 0.618 ± 0.062, I = 1.201 ± 0.147; in the Large White, Ho = 0.680 ± 0.029, I = 1.362 ± 0.074; and in the Yorkshire, Ho = 0.642 ± 0.065, I = 1.287 ± 0.156. The results of checking genotypic Hardy–Weinberg equilibrium based on the G-test of maximum likelihood demonstrated that the overwhelming majority of loci in the wild boar population were in the state of said equilibrium. By contrast, in pig breed populations, some loci demonstrated a significant deviation from the indicated equilibrium. In addition, the Yorkshire, Large White, and Landrace populations had loci, for which the hypothesis of neutrality was reliably rejected based on the results of the Ewens–Watterson test. The revealed private alleles, characteristic of the wild boar and breeds, can later be used to identify them. The ordination of the centroids of different herds in the space of the first two principal coordinates based on the matrix of pairwise estimates of Nei’s genetic distances showed that the most distant populations are the Duroc and Boar breeds, and the most genetically close are the Yorkshire and Landrace breeds. The closest to the wild boar population was the Large White breed. The assessment of the effective size, carried out using the method based on the linkage disequilibrium and the molecular coancestry method, showed that in all studied groups, including the wild boar population, the effective size was less than 100 individuals. The low effective size of the wild boar population (Ne = 21.8, Neb = 4.0) is probably caused by the death and shooting of animals due to Pestis africana suum.

A continent-wide high genetic load in African buffalo revealed by clines in the frequency of deleterious alleles, genetic hitchhiking and linkage disequilibrium

A high genetic load can negatively affect population viability and increase susceptibility to diseases and other environmental stressors. Prior microsatellite studies of two African buffalo (Syncerus caffer) populations in South Africa indicated substantial genome-wide genetic load due to high-frequency occurrence of deleterious alleles. The occurrence of these alleles, which negatively affect male body condition and bovine tuberculosis resistance, throughout most of the buffalo's range were evaluated in this study. Using available microsatellite data (2-17 microsatellite loci) for 1676 animals from 34 localities (from 25°S to 5°N), we uncovered continent-wide frequency clines of microsatellite alleles associated with the aforementioned male traits. Frequencies decreased over a south-to-north latitude range (average per-locus Pearson r = -0.22). The frequency clines coincided with a multilocus-heterozygosity cline (adjusted R2 = 0.84), showing up to a 16% decrease in southern Africa compared to East Africa. Furthermore, continent-wide linkage disequilibrium (LD) at five linked locus pairs was detected, characterized by a high fraction of positive interlocus associations (0.66, 95% CI: 0.53, 0.77) between male-deleterious-trait-associated alleles. Our findings suggest continent-wide and genome-wide selection of male-deleterious alleles driven by an earlier observed sex-chromosomal meiotic drive system, resulting in frequency clines, reduced heterozygosity due to hitchhiking effects and extensive LD due to male-deleterious alleles co-occurring in haplotypes. The selection pressures involved must be high to prevent destruction of allele-frequency clines and haplotypes by LD decay. Since most buffalo populations are stable, these results indicate that natural mammal populations, depending on their genetic background, can withstand a high genetic load.

Assessment of the Genetic Diversity of a Local Pig Breed Using Pedigree and SNP Data

Herein, the genetic diversity of the local Přeštice Black-Pied pig breed was assessed by the simultaneous analysis of the pedigree and single nucleotide polymorphism (SNP) data. The information about sire line, dam, date of birth, sex, breeding line, and herd for 1971 individuals was considered in the pedigree analysis. The SNP analysis (n = 181) was performed using the Illumina PorcineSNP60 BeadChip kit. The quality of pedigree and SNPs and the inbreeding coefficients (F) and effective population size (Ne) were evaluated. The correlations between inbreeding based on the runs of homozygosity (FROH) and pedigree (FPED) were also calculated. The average FPED for all animals was 3.44%, while the FROH varied from 10.81% for a minimum size of 1 Mbp to 3.98% for a minimum size of 16 Mbp. The average minor allele frequency was 0.28 ± 0.11. The observed and expected within breed heterozygosities were 0.38 ± 0.13 and 0.37 ± 0.12, respectively. The Ne, obtained using both the data sources, reached values around 50 animals. Moderate correlation coefficients (0.49-0.54) were observed between FPED and FROH. It is necessary to make decisions that stabilize the inbreeding rate in the long-term using optimal contribution selection based on the available SNP data.

Further host-genomic characterization of total antibody response to PRRSV vaccination and its relationship with reproductive performance in commercial sows: genome-wide haplotype and zygosity analyses

Background

The possibility of using antibody response (S/P ratio) to PRRSV vaccination measured in crossbred commercial gilts as a genetic indicator for reproductive performance in vaccinated crossbred sows has motivated further studies of the genomic basis of this trait. In this study, we investigated the association of haplotypes and runs of homozygosity (ROH) and heterozygosity (ROHet) with S/P ratio and their impact on reproductive performance.

Results

There was no association (P-value ≥ 0.18) of S/P ratio with the percentage of ROH or ROHet, or with the percentage of heterozygosity across the whole genome or in the major histocompatibility complex (MHC) region. However, specific ROH and ROHet regions were significantly associated (P-value ≤ 0.01) with S/P ratio on chromosomes 1, 4, 5, 7, 10, 11, 13, and 17 but not (P-value ≥ 0.10) with reproductive performance. With the haplotype-based genome-wide association study (GWAS), additional genomic regions associated with S/P ratio were identified on chromosomes 4, 7, and 9. These regions harbor immune-related genes, such as SLA-DOB, TAP2, TAPBP, TMIGD3, and ADORA. Four haplotypes at the identified region on chromosome 7 were also associated with multiple reproductive traits. A haplotype significantly associated with S/P ratio that is located in the MHC region may be in stronger linkage disequilibrium (LD) with the quantitative trait loci (QTL) than the previously identified single nucleotide polymorphism (SNP) (H3GA0020505) given the larger estimate of genetic variance explained by the haplotype than by the SNP.

Conclusions

Specific ROH and ROHet regions were significantly associated with S/P ratio. The haplotype-based GWAS identified novel QTL for S/P ratio on chromosomes 4, 7, and 9 and confirmed the presence of at least one QTL in the MHC region. The chromosome 7 region was also associated with reproductive performance. These results narrow the search for causal genes in this region and suggest SLA-DOB and TAP2 as potential candidate genes associated with S/P ratio on chromosome 7. These results provide additional opportunities for marker-assisted selection and genomic selection for S/P ratio as genetic indicator for litter size in commercial pig populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00676-5.

Detection and Visualization of Heterozygosity-Rich Regions and Runs of Homozygosity in Worldwide Sheep Populations

In this study, we chose 17 worldwide sheep populations of eight breeds, which were intensively selected for different purposes (meat, milk, or wool), or locally-adapted breeds, in order to identify and characterize factors impacting the detection of runs of homozygosity (ROH) and heterozygosity-rich regions (HRRs) in sheep. We also applied a business intelligence (BI) tool to integrate and visualize outputs from complementary analyses. We observed a prevalence of short ROH, and a clear distinction between the ROH profiles across populations. The visualizations showed a fragmentation of medium and long ROH segments. Furthermore, we tested different scenarios for the detection of HRR and evaluated the impact of the detection parameters used. Our findings suggest that HRRs are small and frequent in the sheep genome; however, further studies with higher density SNP chips and different detection methods are suggested for future research. We also defined ROH and HRR islands and identified common regions across the populations, where genes related to a variety of traits were reported, such as body size, muscle development, and brain functions. These results indicate that such regions are associated with many traits, and thus were under selective pressure in sheep breeds raised for different purposes. Interestingly, many candidate genes detected within the HRR islands were associated with brain integrity. We also observed a strong association of high linkage disequilibrium pattern with ROH compared with HRR, despite the fact that many regions in linkage disequilibrium were not located in ROH regions.

Estimates of genomic inbreeding and identification of candidate regions that differ between Chinese indigenous sheep breeds

Background

A run of homozygosity (ROH) is a consecutive tract of homozygous genotypes in an individual that indicates it has inherited the same ancestral haplotype from both parents. Genomic inbreeding can be quantified based on ROH. Genomic regions enriched with ROH may be indicative of selection sweeps and are known as ROH islands. We carried out ROH analyses in five Chinese indigenous sheep breeds; Altay sheep (n = 50 individuals), Large-tailed Han sheep (n = 50), Hulun Buir sheep (n = 150), Short-tailed grassland sheep (n = 150), and Tibetan sheep (n = 50), using genotypes from an Ovine Infinium HD SNP BeadChip.

Results

A total of 18,288 ROH were identified. The average number of ROH per individual across the five sheep breeds ranged from 39 (Hulun Buir sheep) to 78 (Large-tailed Han sheep) and the average length of ROH ranged from 0.929 Mb (Hulun Buir sheep) to 2.544 Mb (Large-tailed Han sheep). The effective population size (Ne) of Altay sheep, Large-tailed Han sheep, Hulun Buir sheep, Short-tailed grassland sheep and Tibetan sheep were estimated to be 81, 78, 253, 238 and 70 five generations ago. The highest ROH-based inbreeding estimate (F_ROH) was 0.0808 in Large-tailed Han sheep, whereas the lowest F_ROH was 0.0148 in Hulun Buir sheep. Furthermore, the highest proportion of long ROH fragments (> 5 Mb) was observed in the Large-tailed Han sheep breed which indicated recent inbreeding. In total, 49 ROH islands (the top 0.1% of the SNPs most commonly observed in ROH) were identified in the five sheep breeds. Three ROH islands were common to all the five sheep breeds, and were located on OAR2: 12.2–12.3 Mb, OAR12: 78.4–79.1 Mb and OAR13: 53.0–53.6 Mb. Three breed-specific ROH islands were observed in Altay sheep (OAR15: 3.4–3.8 Mb), Large-tailed Han sheep (ORA17: 53.5–53.8 Mb) and Tibetan sheep (ORA5:19.8–20.2 Mb). Collectively, the ROH islands harbored 78 unique genes, including 19 genes that have been documented as having associations with tail types, adaptation, growth, body size, reproduction or immune response.

Conclusion

Different ROH patterns were observed in five Chinese indigenous sheep breeds, which reflected their different population histories. Large-tailed Han sheep had the highest genomic inbreeding coefficients and the highest proportion of long ROH fragments indicating recent inbreeding. Candidate genes in ROH islands could be used to illustrate the genetic characteristics of these five sheep breeds. Our findings contribute to the understanding of genetic diversity and population demography, and help design and implement breeding and conservation strategies for Chinese sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-021-00608-9.

An Overview of the Use of Genotyping Techniques for Assessing Genetic Diversity in Local Farm Animal Breeds

Simple Summary

The number of local farm animal breeds is declining worldwide. However, these breeds have different degrees of genetic diversity. Measuring genetic diversity is important for the development of conservation strategies and, therefore, various genomic analysis techniques are available. The aim of the present work was to shed light on the use of these techniques in diversity studies of local breeds. In summary, a total of 133 worldwide studies that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that over time, almost all available genomic techniques were used and various diversity parameters were calculated. Therefore, the present results provide a comprehensive overview of the application of these techniques in the field of local breeds. This can provide helpful insights into the advancement of the conservation of breeds with high genetic diversity.

Abstract

Globally, many local farm animal breeds are threatened with extinction. However, these breeds contribute to the high amount of genetic diversity required to combat unforeseen future challenges of livestock production systems. To assess genetic diversity, various genotyping techniques have been developed. Based on the respective genomic information, different parameters, e.g., heterozygosity, allele frequencies and inbreeding coefficient, can be measured in order to reveal genetic diversity between and within breeds. The aim of the present work was to shed light on the use of genotyping techniques in the field of local farm animal breeds. Therefore, a total of 133 studies across the world that examined genetic diversity in local cattle, sheep, goat, chicken and pig breeds were reviewed. The results show that diversity of cattle was most often investigated with microsatellite use as the main technique. Furthermore, a large variety of diversity parameters that were calculated with different programs were identified. For 15% of the included studies, the used genotypes are publicly available, and, in 6%, phenotypes were recorded. In conclusion, the present results provide a comprehensive overview of the application of genotyping techniques in the field of local breeds. This can provide helpful insights to advance the conservation of breeds.

Host Genetic Diversity and Infectious Diseases. Focus on Wild Boar, Red Deer and Tuberculosis

Host genetic diversity tends to limit disease spread in nature and buffers populations against epidemics. Genetic diversity in wildlife is expected to receive increasing attention in contexts related to disease transmission and human health. Ungulates such as wild boar (Sus scrofa) and red deer (Cervus elaphus) are important zoonotic hosts that can be precursors to disease emergence and spread in humans. Tuberculosis is a zoonotic disease with relevant consequences and can present high prevalence in wild boar and red deer populations. Here, we review studies on the genetic diversity of ungulates and determine to what extent these studies consider its importance on the spread of disease. This assessment also focused on wild boar, red deer, and tuberculosis. We found a disconnection between studies treating genetic diversity and those dealing with infectious diseases. Contrarily, genetic diversity studies in ungulates are mainly concerned with conservation. Despite the existing disconnection between studies on genetic diversity and studies on disease emergence and spread, the knowledge gathered in each discipline can be applied to the other. The bidirectional applications are illustrated in wild boar and red deer populations from Spain, where TB is an important threat for wildlife, livestock, and humans.

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

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

Intra- and interbreed genetic heterogeneity and divergence in four commercial pig breeds based on microsatellite markers

In recent years, there has been an increasing amount of attention paid to the genetic health of domesticated animals and its relationship with the level of inbreeding and genetic diversity. At the same time, insufficient attention is still paid to the study of intrabreed genetic diversity and intrabreed stratification. The main goal of our work was to analyze the intra- and interbreed genetic diversity of commercial pig breeds on the basis of DNA microsatellite (MS-DNA) polymorphism. In total, the work used data for 3,308 pigs, which represented 11 herds. The animals belonged to four commercial pig breeds – Duroc (DR), Yorkshire (YR), Landrace (LN) and Large White (LW). 12 microsatellite loci recommended by ISAG-FAO and arranged in one multiplex panel (S0101, S0155, S0228, S0355, S0386, SW24, SW240, SW72, SW857, SW911, SW936, SW951) were used as DNA markers. When analyzing the intra- and interbreed variability of 11 herds, we found that all studied breeds significantly differed in terms of the proportion of both rare and the most common alleles. At the same time, the noted differences were determined, first of all, by the variability between individual herds within their breed. The location of herd centroids is random and is not consistent with their breed affiliation at all. When individuals belonging to the same breed are combined, the centroids of pig breeds in the space of first two axes from a Principal Coordinate Analysis form two clusters. The first one contains the only red pig breed (DR) used in the analysis, while the second one contains white pig breeds. In six pig herds the Ne estimates were below 50 inds., in two herds they were in the range of 50–100 inds., and finally in three herds the Ne estimates exceeded 100 inds. The analysis of the genetic variability of pigs of four commercial breeds showed that the high level of interbreed differences is caused, first of all, by the high variability among pig herds within each studied breed. Such intrabreed stratification can be formed due to the manifestation of many causes: different genetic basis of the founders of intrabreed genealogical groups, geographical isolation, different directions of selection within individual herds, exchange of animals between separate herds, the use of inbreeding in the practice of selection together with isolation, etc. Important consequences of intrabreed stratification are an increase in the level of interherd diversity (which is not lower than the level of interbreed diversity) against the background of a decrease in variability within individual herds, as well as a significant deficit of heterozygotes and an increase in the role of negative genetic and demographic processes. Thus, the existence of genetic heterogeneity within commercial pig breeds should be considered as an essential element in the history of their formation and breeding.

Whole-genome SNP markers reveal conservation status, signatures of selection, and introgression in Chinese Laiwu pigs

Laiwu pigs are a Chinese indigenous breed that is renowned for its exceptionally high intramuscular fat content (average greater than 6%), providing an excellent genetic resource for the genetic improvement of meat quality of modern commercial pigs. To uncover genetic diversity, population structure, signature of selection, and potential exotic introgression in this breed, we sampled 238 Laiwu pigs from a state-supported conservation population and genotyped these individuals using GeneSeek 80K SNP BeadChip. We then conducted in-depth population genetics analyses for the Laiwu pig in a context of 1,116 pigs from 42 Eurasian diverse breeds. First, we show that the current Laiwu population has more abundant genetic diversity than the population of 18 years ago likely due to gene flow from European commercial breeds. Both neighbor-joining (NJ) and principal component analyses indicate the introgression of European haplotypes into Laiwu pigs. The admixture analysis reveals that an average 26.66% of Laiwu genetic components are of European origin. Then, we assigned the tested individuals to different families according to their clustering patterns in the NJ tree and proposed a family-based conservation strategy to reduce the risk of inbreeding depression in Laiwu pigs. Next, we explored three statistics (ROH and iHS and EigenGWAS) to identify a list of candidate genes for fat deposition, reproduction, and growth in Laiwu pigs. Last, we detected a strong signature of introgression from European pigs into Laiwu pigs at the GPC6 locus that regulates the growth of developing long bones. Further association analyses indicate that the introgressed GPC6 haplotype likely contributed to the improvement of growth performance in Laiwu pigs. Altogether, this study not only benefits the better conservation of the Laiwu pig, but also advances our knowledge of the poorly understood effect of human-mediated introgression on phenotypic traits in Chinese indigenous pigs.

Genome-Wide Analysis Revealed Homozygosity and Demographic History of Five Chinese Sheep Breeds Adapted to Different Environments

Homozygosity of long sequence genotypes are a result of parents transmitting identical haplotypes, which can be used to estimate their auto-zygosity. Therefore, we used high-density SNP Chip data to characterize the auto-zygosity of each breed according to the occurrence and distribution of runs of homozygosity (ROH). Subsequently, we identified the genomic regions with high runs of homozygosity frequencies within individuals of each breed. We selected 96 sheep samples from five local Chinese sheep breeds belonging to different geographical locations. We identified 3046 ROHs within the study breed individuals, among which the longer segments (>1–5 Mb) were dominant. On average, ROH segments covered about 12% of the genomes; the coverage rate of OAR20 was the lowest and that of OAR2 was the highest. The distribution analysis of runs of homozygosity showed that the detected ROH mainly distributed between >26 and 28 Mb. The Hetian and Hu sheep showed the lowest ROH distribution. The estimation of homozygosity level reflects the history of modern and ancient inbreeding, which may affect the genomes of Chinese indigenous sheep breeds and indicate that some animals have experienced recent self-pollination events (Yabuyi, Karakul and Wadi). In these sheep breeds, the genomic regions were assumed to be under selection signatures frequently in line with long ROH. These regions included candidate genes associated with disease resistance traits (5S_rRNA), the innate and adaptive immune response (HERC2 and CYFIP1), digestion and metabolism (CENPJ), growth (SPP1), body size and developments (GJB2 and GJA3). This study highlighted new insights into the ROH patterns and provides a basis for future breeding and conservation strategies of Chinese sheep breeds.

High-density single nucleotide polymorphism chip-based conservation genetic analysis of indigenous pig breeds from Shandong Province, China

Objective

Shandong indigenous pig breeds are important Chinese pig resources. Their progressive population decline in recent decades has attracted attention towards their conservation. Conservation genetics of these indigenous breeds are essential for developing a conservation and utilization scheme.

Methods

A high-density single nucleotide polymorphism (HD-SNP) chip-based comparative analysis of genetic characteristics was performed for seven Shandong indigenous pig breeds in the context of five Western commercial breeds.

Results

The results showed that Shandong indigenous pig breeds varied greatly in genetic diversity, effective population size, inbreeding level, and genetic distance with the Western commercial breeds. Specifically, Laiwu and Dapulian displayed low genetic diversity, and had a genetically distant relationship with the Western commercial breeds (average F statistics [F_ST] value of 0.3226 and 0.2666, respectively). Contrastingly, the other five breeds (Yantai, Licha, Yimeng, Wulain, and Heigai) displayed high genetic diversity within breed and had some extent of mixture pattern with the Western commercial breeds, especially Duroc and Landrace (F_ST values from 0.1043 to 0.2536). Furthermore, intensive gene flow was discovered among the seven Shandong indigenous breeds, particularly Wulian, Licha, and Heigai, as indicated by the large cluster formed in the principal component analysis scatterplot and small population differentiation (average of 0.1253) among them.

Conclusion

Our study advances the understanding of genetic characteristics of Shandong indigenous breeds and provides essential information for developing an appropriate conservation and utilization scheme for these breeds.

Genetic Variability in Local and Imported Germplasm Chicken Populations as Revealed by Analyzing Runs of Homozygosity

Simple Summary

To maintain the uniqueness of conserved chicken populations of local and imported breeds is of great importance. In this study, we genotyped small populations belonging to 14 breeds and 7 crossbreds using an Illumina Chicken 60K SNP (Single Nucleotide Polymorphisms) BeadChip and looked for appropriate methods to characterize their purity/variability. It was not straightforward to identify crossbred individuals, and the best approach was based on calculating the length and number of homozygous regions, or runs of homozygosity (ROH), in the populations studied. The latter enabled most accurate identification of crossbreds and can be served as an effective tool in testing genome-wide purity of chicken breeds.

Abstract

Preserving breed uniqueness and purity is vitally important in developing conservation/breeding programs for a germplasm collection of rare and endangered chicken breeds. The present study was aimed at analyzing SNP genetic variability of 21 small local and imported purebred and F₁ crossbred populations and identifying crossbreeding events via whole-genome evaluation of runs of homozygosity (ROH). The admixture models more efficiently reflected population structure, pinpointing crossbreeding events in the presence of ancestral populations but not in their absence. Multidimensional scaling and F_ST-based analyses did not discriminate properly between purebred populations and F₁ crossbreds, especially when comparing related breeds. When applying the ROH-based approach, more and longer ROHs were revealed in purebred individuals/populations, suggesting this as an effective implement in genome-wide analysis of germplasm breed purity.

Genome-Wide Assessment of Runs of Homozygosity in Chinese Wagyu Beef Cattle

Runs of homozygosity (ROH) are continuous homozygous regions that generally exist in the DNA sequence of diploid organisms. Identifications of ROH leading to reduction in performance can provide valuable insight into the genetic architecture of complex traits. Here, we evaluated genome-wide patterns of homozygosity and their association with important traits in Chinese Wagyu beef cattle. We identified a total of 29,271 ROH segments from 462 animals. Within each animal, an average number of ROH was 63.36 while an average length was 62.19 Mb. To evaluate the enrichment of ROH across genomes, we initially identified 280 ROH regions by merging ROH events across all individuals. Of these, nine regions containing 154 candidate genes, were significantly associated with six traits (body height, chest circumference, fat coverage, backfat thickness, ribeye area, and carcass length; p < 0.01). Moreover, we found 26 consensus ROH regions with frequencies exceeding 10%, and several regions overlapped with QTLs, which are associated with body weight, calving ease, and stillbirth. Among them, we observed 41 candidate genes, including BCKDHB, MAB21L1, SLC2A13, FGFR3, FGFRL1, CPLX1, CTNNA1, CORT, CTNNBIP1, and NMNAT1, which have been previously reported to be related to body conformation, meat quality, susceptibility, and reproductive traits. In summary, we assessed genome-wide autozygosity patterns and inbreeding levels in Chinese Wagyu beef cattle. Our study identified many candidate regions and genes overlapped with ROH for several important traits, which could be unitized to assist the design of a selection mating strategy in beef cattle.

How to study runs of homozygosity using PLINK? A guide for analyzing medium density SNP data in livestock and pet species

Background

PLINK is probably the most used program for analyzing SNP genotypes and runs of homozygosity (ROH), both in human and in animal populations. The last decade, ROH analyses have become the state-of-the-art method for inbreeding assessment. In PLINK, the --homozyg function is used to perform ROH analyses and relies on several input settings. These settings can have a large impact on the outcome and default values are not always appropriate for medium density SNP array data. Guidelines for a robust and uniform ROH analysis in PLINK using medium density data are lacking, albeit these guidelines are vital for comparing different ROH studies. In this study, 8 populations of different livestock and pet species are used to demonstrate the importance of PLINK input settings. Moreover, the effects of pruning SNPs for low minor allele frequencies and linkage disequilibrium on ROH detection are shown.

Results

We introduce the genome coverage parameter to appropriately estimate F_ROH and to check the validity of ROH analyses. The effect of pruning for linkage disequilibrium and low minor allele frequencies on ROH analyses is highly population dependent and such pruning may result in missed ROH. PLINK’s minimal density requirement is crucial for medium density genotypes and if set too low, genome coverage of the ROH analysis is limited. Finally, we provide recommendations for the maximal gap, scanning window length and threshold settings.

Conclusions

In this study, we present guidelines for an adequate and robust ROH analysis in PLINK on medium density SNP data. Furthermore, we advise to report parameter settings in publications, and to validate them prior to analysis. Moreover, we encourage authors to report genome coverage to reflect the ROH analysis’ validity. Implementing these guidelines will substantially improve the overall quality and uniformity of ROH analyses.

Genomic diversity, linkage disequilibrium and selection signatures in European local pig breeds assessed with a high density SNP chip

Genetic characterization of local breeds is essential to preserve their genomic variability, to advance conservation policies and to contribute to their promotion and sustainability. Genomic diversity of twenty European local pig breeds and a small sample of Spanish wild pigs was assessed using high density SNP chips. A total of 992 DNA samples were analyzed with the GeneSeek Genomic Profiler (GGP) 70 K HD porcine genotyping chip. Genotype data was employed to compute genetic diversity, population differentiation and structure, genetic distances, linkage disequilibrium and effective population size. Our results point out several breeds, such as Turopolje, Apulo Calabrese, Casertana, Mora Romagnola and Lithuanian indigenous wattle, having the lowest genetic diversity, supported by low heterozygosity and very small effective population size, demonstrating the need of enhanced conservation strategies. Principal components analysis showed the clustering of the individuals of the same breed, with few breeds being clearly isolated from the rest. Several breeds were partially overlapped, suggesting genetic closeness, which was particularly marked in the case of Iberian and Alentejana breeds. Spanish wild boar was also narrowly related to other western populations, in agreement with recurrent admixture between wild and domestic animals. We also searched across the genome for loci under diversifying selection based on F_ST outlier tests. Candidate genes that may underlie differences in adaptation to specific environments and productive systems and phenotypic traits were detected in potentially selected genomic regions.

Genome-Wide Scan for Runs of Homozygosity Identifies Candidate Genes in Three Pig Breeds

Runs of homozygosity (ROH) are contiguous homozygous genotype segments in the genome that are present in an individual since the identical haplotypes are inherited from each parent. The aim of this study was to investigate the frequency and distribution of ROH in the genomes of Landrace, Songliao black and Yorkshire pigs. We calculated two types of genome inbreeding coefficients and their correlation, including the inbreeding coefficient based on ROH (F_ROH) and the inbreeding coefficient based on the difference between the observed and expected number of homozygous genotypes (F_HOM). Furthermore, we identified candidate genes in the genomic region most associated with ROH. We identified 21,312 ROH in total. The average number of ROH per individual was 32.99 ± 0.38 and the average length of ROH was 6.40 ± 0.070 Mb in the three breeds. The F_ROH results showed that Yorkshire pigs exhibited the highest level of inbreeding (0.092 ± 0.0015) and that Landrace pigs exhibited the lowest level of inbreeding (0.073 ± 0.0047). The average correlation between F_ROH and F_HOM was high (0.94) within three breeds. The length of ROH provides insight into the inbreeding history of these three pig breeds. In this study, Songliao black pigs presented a higher frequency and average length of long ROH (>40 Mb) compared with those of Landrace and Yorkshire pigs, which indicated greater inbreeding in recent times. Genes related to reproductive traits (GATM, SPATA46, HSD17B7, VANGL2, DAXX, CPEB1), meat quality traits (NR1I3, APOA2, USF1) and energy conversion (NDUFS2) were identified within genomic regions with a high frequency of ROH. These genes could be used as target genes for further marker-assisted selection and genome selection.

Genomic Analysis Reveals Specific Patterns of Homozygosity and Heterozygosity in Inbred Pigs

The inbred strain of miniature pig is an ideal model for biomedical research due to its high level of homozygosity. In this study, we investigated genetic diversity, relatedness, homozygosity, and heterozygosity using the Porcine SNP60K BeadChip in both inbred and non-inbred Wuzhishan pigs (WZSPs). Our results from multidimensional scaling, admixture, and phylogenetic analyses indicated that the inbred WZSP, with its unique genetic properties, can be utilized as a novel genetic resource for pig genome studies. Inbreeding depression and run of homozygosity (ROH) analyses revealed an average of 61 and 12 ROH regions in the inbred and non-inbred genomes of WZSPs, respectively. By investigating ROH number, length, and distribution across generations, we further briefly studied the impacts of recombination and demography on ROH in these WZSPs. Finally, we explored the SNPs with higher heterozygosity across generations and their potential functional implications in the inbred WZSP. We detected 56 SNPs showing constant heterozygosity with He = 1 across six generations in inbred pigs, while only one was found in the non-inbred population. Among these SNPs, we observed nine SNPs located in swine RefSeq genes, which were found to be involved in signaling and immune processes. Together, our findings indicate that the inbred-specific pattern of homozygosity and heterozygosity in inbred pigs can offer valuable insights for elucidating the mechanisms of inbreeding in farm animals.

A guinea fowl genome assembly provides new evidence on evolution following domestication and selection in galliformes

The helmeted guinea fowl Numida meleagris belongs to the order Galliformes. Its natural range includes a large part of sub‐Saharan Africa, from Senegal to Eritrea and from Chad to South Africa. Archaeozoological and artistic evidence suggest domestication of this species may have occurred about 2,000 years BP in Mali and Sudan primarily as a food resource, although villagers also benefit from its capacity to give loud alarm calls in case of danger, of its ability to consume parasites such as ticks and to hunt snakes, thus suggesting its domestication may have resulted from a commensal association process. Today, it is still farmed in Africa, mainly as a traditional village poultry, and is also bred more intensively in other countries, mainly France and Italy. The lack of available molecular genetic markers has limited the genetic studies conducted to date on guinea fowl. We present here a first‐generation whole‐genome sequence draft assembly used as a reference for a study by a Pool‐seq approach of wild and domestic populations from Europe and Africa. We show that the domestic populations share a higher genetic similarity between each other than they do to wild populations living in the same geographical area. Several genomic regions showing selection signatures putatively related to domestication or importation to Europe were detected, containing candidate genes, most notably EDNRB2, possibly explaining losses in plumage coloration phenotypes in domesticated populations.

Inbreeding and effective population size in French dairy sheep: Comparison between genomic and pedigree estimates

Before availability of dense SNP data, genetic diversity was characterized and managed with pedigree-based information. Besides this classical approach, 2 methodologies have been proposed in recent years to characterize and manage diversity from dense SNP data: the SNP-by-SNP approach and the alternative based on runs of homozygosity (ROH). The establishment of criteria to identify ROH is a current constraint in the literature dealing with ROH. The objective of this study was, using a medium-density SNP chip, to quantify by 3 methods (pedigree, SNP-by-SNP, and ROH) the genetic diversity on 5 selected French dairy sheep subpopulations and breeds and to assess the effect of the definition of ROH on these estimates. The data set available included individuals from the breeds Basco-Béarnaise, Manech Tête Noire, Manech Tête Rousse, and 2 subpopulations of Lacaune: Lacaune Confederation and Lacaune Ovitest. Animals were genotyped with the Illumina OvineSNP50 BeadChip (Illumina Inc., San Diego, CA). After filtering, the genomic data included 38,287 autosomal SNP and 8,700 individuals, which comprised 72,803 animals in the pedigree. The results indicated that no significant differences were observed in effective population size estimates obtained from pedigree or genomic (SNP-by-SNP or ROH) information. In general, estimates of effective population size were above 200 in Lacaune Confederation and Lacaune Ovitest subpopulations and below 200 in Basco-Béarnaise, Manech Tête Noire, and Manech Tête Rousse breeds. The minimum length that constituted a ROH, the minimum number of SNP that constituted a ROH, as well as the minimum density and the maximum distance allowed between 2 homozygous SNP are ROH-defining factors with important implications in the estimation of the rate of inbreeding. The ROH-based rates of inbreeding in concordance with those obtained from pedigree information require a specific set of values. This particular set of values is different from that identified to obtain ROH-based rates of inbreeding similar to those obtained on a SNP-by-SNP basis. Factors to define ROH do not change the results much unless extreme values are considered, although further research on ROH-based inbreeding is still required.

Genomic structure of a crossbred Landrace pig population

Single nucleotide polymorphism (SNP) markers are used to study population structure and conservation genetics, which permits assessing similarities regarding the linkage disequilibrium and information about the relationship among individuals. To investigate the population genomic structure of 300 females and 25 males from a commercial maternal pig line we analyzed linkage disequilibrium extent, inbreeding coefficients using genomic and conventional pedigree data, and population stratification. The average linkage disequilibrium (r²) was 0.291 ± 0.312 for all adjacent SNPs, distancing less than 100 Kb (kilobase) between markers. The average inbreeding coefficients obtained from runs of homozygosity (ROH) and pedigree analyses were 0.119 and 0.0001, respectively. Low correlation was observed between the inbreeding coefficients possibly as a result of genetic recombination effect accounted for the ROH estimates or caused by pedigree identification errors. A large number of long ROHs might indicate recent inbreeding events in the studied population. A total of 36 homozygous segments were found in more than 30% of the population and these ROH harbor genes associated with reproductive traits. The population stratification analysis indicated that this population was possibly originated from two distinct populations, which is a result from crossings between the eastern and western breeds used in the formation of the line. Our findings provide support to understand the genetic structure of swine populations and may assist breeding companies to avoid a high level of inbreeding coefficients to maintain genetic diversity, showing the effectiveness of using genome-wide SNP information for quantifying inbreeding when the pedigree was incomplete or incorrect.

Diversity across major and candidate genes in European local pig breeds

The aim of this work was to analyse the distribution of causal and candidate mutations associated to relevant productive traits in twenty local European pig breeds. Also, the potential of the SNP panel employed for elucidating the genetic structure and relationships among breeds was evaluated. Most relevant genes and mutations associated with pig morphological, productive, meat quality, reproductive and disease resistance traits were prioritized and analyzed in a maximum of 47 blood samples from each of the breeds (Alentejana, Apulo-Calabrese, Basque, Bísara, Majorcan Black, Black Slavonian (Crna slavonska), Casertana, Cinta Senese, Gascon, Iberian, Krškopolje (Krškopoljski), Lithuanian indigenous wattle, Lithuanian White Old Type, Mora Romagnola, Moravka, Nero Siciliano, Sarda, Schwäbisch-Hällisches Schwein (Swabian Hall pig), Swallow-Bellied Mangalitsa and Turopolje). We successfully analyzed allelic variation in 39 polymorphisms, located in 33 candidate genes. Results provide relevant information regarding genetic diversity and segregation of SNPs associated to production and quality traits. Coat color and morphological trait-genes that show low level of segregation, and fixed SNPs may be useful for traceability. On the other hand, we detected SNPs which may be useful for association studies as well as breeding programs. For instance, we observed predominance of alleles that might be unfavorable for disease resistance and boar taint in most breeds and segregation of many alleles involved in meat quality, fatness and growth traits. Overall, these findings provide a detailed catalogue of segregating candidate SNPs in 20 European local pig breeds that may be useful for traceability purposes, for association studies and for breeding schemes. Population genetic analyses based on these candidate genes are able to uncover some clues regarding the hidden genetic substructure of these populations, as the extreme genetic closeness between Iberian and Alentejana breeds and an uneven admixture of the breeds studied. The results are in agreement with available knowledge regarding breed history and management, although largest panels of neutral markers should be employed to get a deeper understanding of the population's structure and relationships.

Genetic diversity, population structure and phylogenetic relationships of three indigenous pig breeds from Jiangxi Province, China, in a worldwide panel of pigs

Dongxiang Spotted, Pingxiang Two-End-Black and Yushan Black pigs are three indigenous breeds in Jiangxi Province, China, that have been listed in the national conservation program for Chinese indigenous livestock germplasm. Here, we investigated the genetic diversity and population structure of the nucleus populations of these three breeds in a worldwide context of European and Chinese pigs using Illumina Porcine 60K chips. Our data indicate that Dongxiang Spotted and Yushan Black pigs have recently experienced severe inbreeding. The two breeds show large runs-of-homozygosity values, long-range extents of linkage disequilibrium and reduced observed heterozygosity. In contrast, Pingxiang Two-End-Black pigs are rich in genetic diversity and have few inbred individuals. Both phylogenetic and admixture analyses illustrate that Dongxiang Spotted and Yushan Black pigs are genetically close to their geographical neighbors in East China and that Pingxiang Two-End-Black pigs have a close relationship with three other Two-End-Black breeds in central China. Finally, we reconstructed the family structures of the three breeds and propose a reliable breeding strategy to better conserve these breeds.

Conservation status and historical relatedness of Italian cattle breeds

Background

In the last 50 years, the diversity of cattle breeds has experienced a severe contraction. However, in spite of the growing diffusion of cosmopolite specialized breeds, several local cattle breeds are still farmed in Italy. Genetic characterization of breeds represents an essential step to guide decisions in the management of farm animal genetic resources. The aim of this work was to provide a high-resolution representation of the genome-wide diversity and population structure of Italian local cattle breeds using a medium-density single nucleotide polymorphism (SNP) array.

Results

After quality control filtering, the dataset included 31,013 SNPs for 800 samples from 32 breeds. Our results on the genetic diversity of these breeds agree largely with their recorded history. We observed a low level of genetic diversity, which together with the small size of the effective populations, confirmed that several breeds are threatened with extinction. According to the analysis of runs of homozygosity, evidence of recent inbreeding was strong in some local breeds, such as Garfagnina, Mucca Pisana and Pontremolese. Patterns of genetic differentiation, shared ancestry, admixture events, and the phylogenetic tree, all suggest the presence of gene flow, in particular among breeds that originate from the same geographical area, such as the Sicilian breeds. In spite of the complex admixture events that most Italian cattle breeds have experienced, they have preserved distinctive characteristics and can be clearly discriminated, which is probably due to differences in genetic origin, environment, genetic isolation and inbreeding.

Conclusions

This study is the first exhaustive genome-wide analysis of the diversity of Italian cattle breeds. The results are of significant importance because they will help design and implement conservation strategies. Indeed, efforts to maintain genetic diversity in these breeds are needed. Improvement of systems to record and monitor inbreeding in these breeds may contribute to their in situ conservation and, in view of this, the availability of genomic data is a fundamental resource.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0406-x) contains supplementary material, which is available to authorized users.

The effects of recent changes in breeding preferences on maintaining traditional Dutch chicken genomic diversity

Traditional Dutch chicken breeds are marginalised breeds of ornamental and cultural-historical importance. In the last decades, miniaturising of existing breeds (so called neo-bantam) has become popular and resulted in alternatives to original large breeds. However, while backcrossing is increasing the neo-bantams homozygosity, genetic exchange between breeders may increase their genetic diversity. We use the 60 K SNP array to characterise the genetic diversity, demographic history, and level of inbreeding of Dutch heritage breeds, and particularly of neo-bantams. Commercial white layers are used to contrast the impact of management strategy on genetic diversity and demography. A high proportion of alleles was found to be shared between large fowls and neo-bantams, suggesting gene flow during neo-bantams development. Population admixture analysis supports these findings, in addition to revealing introgression from neo-bantams of the same breed and of phenotypically similar breeds. The prevalence of long runs of homozygosity (ROH) confirms the importance of recent inbreeding. A high diversity in management, carried out in small breeding units explains the high heterogeneity in diversity and ROH profile displayed by traditional breeds compared to commercial lines. Population bottlenecks may explain the long ROHs in large fowls, while repetitive backcrossing for phenotype selection may account for them in neo-bantams. Our results highlight the importance of using markers to inform breeding programmes on potentially harmful homozygosity to prevent loss of genetic diversity. We conclude that bantamisation has generated unique and identifiable genetic diversity. However, this diversity can only be preserved in the near future through structured breeding programmes.

Genetic diversity, breed composition and admixture of Kenyan domestic pigs

The genetic diversity of African pigs, whether domestic or wild has not been widely studied and there is very limited published information available. Available data suggests that African domestic pigs originate from different domestication centers as opposed to international commercial breeds. We evaluated two domestic pig populations in Western Kenya, in order to characterize the genetic diversity, breed composition and admixture of the pigs in an area known to be endemic for African swine fever (ASF). One of the reasons for characterizing these specific populations is the fact that a proportion of indigenous pigs have tested ASF virus (ASFv) positive but do not present with clinical symptoms of disease indicating some form of tolerance to infection. Pigs were genotyped using either the porcine SNP60 or SNP80 chip. Village pigs were sourced from Busia and Homabay counties in Kenya. Because bush pigs (Potamochoerus larvatus) and warthogs (Phacochoerus spp.) are known to be tolerant to ASFv infection (exhibiting no clinical symptoms despite infection), they were included in the study to assess whether domestic pigs have similar genomic signatures. Additionally, samples representing European wild boar and international commercial breeds were included as references, given their potential contribution to the genetic make-up of the target domestic populations. The data indicate that village pigs in Busia are a non-homogenous admixed population with significant introgression of genes from international commercial breeds. Pigs from Homabay by contrast, represent a homogenous population with a "local indigenous' composition that is distinct from the international breeds, and clusters more closely with the European wild boar than African wild pigs. Interestingly, village pigs from Busia that tested negative by PCR for ASFv genotype IX, had significantly higher local ancestry (>54%) compared to those testing positive, which contained more commercial breed gene introgression. This may have implication for breed selection and utilization in ASF endemic areas. A genome wide scan detected several regions under preferential selection with signatures for pigs from Busia and Homabay being very distinct. Additionally, there was no similarity in specific genes under selection between the wild pigs and domestic pigs despite having some broad areas under similar selection signatures. These results provide a basis to explore possible genetic determinants underlying tolerance to infection by ASFv genotypes and suggests multiple pathways for genetically mediated ASFv tolerance given the diversity of selection signatures observed among the populations studied.

Assessment of runs of homozygosity islands and estimates of genomic inbreeding in Gyr (Bos indicus) dairy cattle

Background

Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (F_ROH), genomic relationship matrix approach (F_GRM) and based on the observed versus expected number of homozygous genotypes (F_HOM), and from pedigree-based coefficient (F_PED).

Results

ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH_{1–2 Mb}) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of F_ROH, F_GRM, F_HOM, and F_PED approaches. F_PED estimates ranged from 0.00 to 0.327 and F_ROH from 0.001 to 0.201. Low to moderate correlations were observed between F_PED-F_ROH and F_GRM-F_ROH, with values ranging from −0.11 to 0.51. Low to high correlations were observed between F_ROH-F_HOM and moderate between F_PED-F_HOM and F_GRM-F_HOM. Correlations between F_ROH from different lengths and F_PED gradually increased with ROH length.

Conclusions

Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low F_PED-F_ROH correlations for small segments indicate that F_PED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that F_ROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4365-3) contains supplementary material, which is available to authorized users.

Mapping of genome-wide copy number variations in the Iranian indigenous cattle using a dense SNP data set

The objective of this study was to present the first map of the copy number variations (CNVs) in Iranian indigenous cattle based on a high-density single nucleotide polymorphism (SNP) dataset. A total of 90 individuals were genotyped using the Illumina BovineHD BeadChip containing 777 962 SNPs. The QuantiSNP algorithm was used to perform a genome-wide CNV detection across autosomal genome. After merging the overlapping CNV, a total of 221 CNV regions were identified encompassing 36.4 Mb or 1.44% of the bovine autosomal genome. The length of the CNV regions ranged from 3.5 to 2252.8 Kb with an average of 163.8 Kb. These regions included 147 loss (66.52%) and 74 gain (33.48%) events containing a total of 637 annotated Ensembl genes. Gene ontology analysis revealed that most of genes in the CNV regions were involved in environmental responses, disease susceptibility and immune system functions. Furthermore, 543 of these genes corresponded to the human orthologous genes, which involved in a wide range of biological functions. Altogether, 73% of the 221 CNV regions overlapped either completely or partially with those previously reported in other cattle studies. Moreover, novel CNV regions involved several quantitative trait loci (QTL)-related to adaptative traits of Iranian indigenous cattle. These results provided a basis to conduct future studies on association between CNV regions and phenotypic variations in the Iranian indigenous cattle.

Weighted likelihood inference of genomic autozygosity patterns in dense genotype data

Background

Genomic regions of autozygosity (ROA) arise when an individual is homozygous for haplotypes inherited identical-by-descent from ancestors shared by both parents. Over the past decade, they have gained importance for understanding evolutionary history and the genetic basis of complex diseases and traits. However, methods to infer ROA in dense genotype data have not evolved in step with advances in genome technology that now enable us to rapidly create large high-resolution genotype datasets, limiting our ability to investigate their constituent ROA patterns.

Methods

We report a weighted likelihood approach for inferring ROA in dense genotype data that accounts for autocorrelation among genotyped positions and the possibilities of unobserved mutation and recombination events, and variability in the confidence of individual genotype calls in whole genome sequence (WGS) data.

Results

Forward-time genetic simulations under two demographic scenarios that reflect situations where inbreeding and its effect on fitness are of interest suggest this approach is better powered than existing state-of-the-art methods to infer ROA at marker densities consistent with WGS and popular microarray genotyping platforms used in human and non-human studies. Moreover, we present evidence that suggests this approach is able to distinguish ROA arising via consanguinity from ROA arising via endogamy. Using subsets of The 1000 Genomes Project Phase 3 data we show that, relative to WGS, intermediate and long ROA are captured robustly with popular microarray platforms, while detection of short ROA is more variable and improves with marker density. Worldwide ROA patterns inferred from WGS data are found to accord well with those previously reported on the basis of microarray genotype data. Finally, we highlight the potential of this approach to detect genomic regions enriched for autozygosity signals in one group relative to another based upon comparisons of per-individual autozygosity likelihoods instead of inferred ROA frequencies.

Conclusions

This weighted likelihood ROA inference approach can assist population- and disease-geneticists working with a wide variety of data types and species to explore ROA patterns and to identify genomic regions with differential ROA signals among groups, thereby advancing our understanding of evolutionary history and the role of recessive variation in phenotypic variation and disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4312-3) contains supplementary material, which is available to authorized users.