Genome-wide scan for runs of homozygosity identifies candidate genes in Wannan Black pigs

Article information

Anim Biosci. 2021;34(12):1895-1902
Publication date (electronic) : 2021 March 10
doi : https://doi.org/10.5713/ab.20.0679
1College of Animal Science and Technology, Anhui Agricultural University, Hefei, Anhui 230036, China
2Anhui Province Key Laboratory of Local Livestock and Poultry Genetic Resource Conservation and Bio-Breeding, Anhui Agricultural University, Hefei, 230036, China
3Key Laboratory of Pig Molecular Quantitative Genetics of Anhui Academy of Agricultural Sciences, Anhui Provincial Key Laboratory of Livestock and Poultry Product Safety Engineering, Institute of Animal Husbandry and Veterinary Medicine, Anhui Academy of Agricultural Sciences, Hefei, Anhui 230031, China
4Key Laboratory of Animal (Poultry) Genetics Breeding and Reproduction, Ministry of Agriculture and Rural Affairs, Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100193, China
*Corresponding Authors: Zongjun Yin, Tel: +86-551-65786303, Fax: +86-551-65786303, E-mail: yinzongjun@ahau.edu.cn. Yueyun Ding, Tel: +86-0551-65785662, Fax: +86-0551-65785662, E-mail: dingyueyun@ahau.edu.cn
a

These authors contributed equally to this work.

Received 2020 September 28; Revised 2020 December 5; Accepted 2021 February 7.

Abstract

Objective

Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that can reveal inbreeding levels, selection pressure, and mating schemes. In this study, ROHs were evaluated in Wannan Black pigs to assess the inbreeding levels and the genome regions with high ROH frequency.

Methods

In a previous study, we obtained 501.52 GB of raw data from resequencing (10×) of the genome and identified 21,316,754 single-nucleotide variants in 20 Wannan Black pig samples. We investigated the number, length, and frequency of ROH using resequencing data to characterize the homozygosity in Wannan Black pigs and identified genomic regions with high ROH frequencies.

Results

In this work, 1,813 ROHs (837 ROHs in 100 to 500 kb, 449 ROHs in 500 to 1,000 kb, 527 ROHs in >1,000 kb) were identified in all samples, and the average genomic inbreeding coefficient (FROH) in Wannan Black pigs was 0.5234. Sixty-one regions on chromosomes 2, 3, 7, 8, 13, 15, and 16 harbored ROH islands. In total, 105 genes were identified in 42 ROH islands, among which some genes were related to production traits.

Conclusion

This is the first study to identify ROH across the genome of Wannan Black pigs, the Chinese native breed of the Anhui province. Overall, Wannan Black pigs have high levels of inbreeding due to the influence of ancient and recent inbreeding due to the genome. These findings are a reliable resource for future studies and contribute to save and use the germplasm resources of Wannan Black pigs.

INTRODUCTION

Pigs were among the first species that were domesticated. They provide an important meat resource and serve as a biomedical animal model for human [1]. Domestication of wild boars began approximately 10,000 years ago [2], and hundreds of pig breeds evolved worldwide due to natural selection and artificial selection pressure. For a long time, pig breeding was the main method of increasing the economic performance of pig production [3]. The fertility, growth speed, and meat quality of pigs is significantly improved over that of the wild boar.

Molecular genetics is key to understanding the biodiversity and evolutionary relationships, and breeding programs [4]. DNA technologies are improving continuously and being increasingly applied to investigate the genotype of animals subjected to breeding [5]. Compared with traditional selection, molecular-assisted selection has advantages of seed selection and shortened breeding period in pigs. However, increased inbreeding is an inevitable consequence of genetic selection in livestock populations [6]. Inbreeding induces impaired performance traits (inbreeding depression) and reduced genetic variation and therefore, is an important factor to consider in animal breeding practices [79]. Inbreeding is usually estimated from the pedigree information of the animal populations. Incomplete records and neglectful historical inbreeding leads to the high error rate estimation and underestimation of the true inbreeding levels [10]. With the reducing cost of resequencing, it is possible to use genome-wide single-nucleotide polymorphism (SNP) information to assess true population genomic inbreeding levels [11].

Runs of homozygosity (ROH) are contiguous lengths of homozygous genotypes that are present in an individual due to parents transmitting identical haplotypes to their offspring [12]. Longer haplotypes are inherited from recent common ancestors and shorter haplotypes are inherited from distant ancestors. ROHs are not randomly distributed across the genome [13], and the distribution of ROH islands across the genome shows that the genome was subjected to selective pressure [14]. Furthermore, ROH has been widely used as a predictor of whole genome inbreeding levels. Past studies have shown that genomic inbreeding coefficient (FROH) is a more accurate surrogate of inbreeding than the inbreeding coefficient estimated from the pedigree [15].

To the best of our knowledge, numerous reports have revealed that ROH has been used to explore the selection signatures in sheep [16], horses [17], and Western pig breeds [18]. However, few studies have focused on ROH analysis in Chinese pig breeds. Wannan Black pig has been bred for more than 100 years in mountainous area of southern Anhui province. Wannan Black pig were historically bred in dark and damp environment and two important factors were chosen for their breeding: reproductive performance and adaptive capacity [19]. As a native Chinese breed endemic to the local agriculture which exhibits disease resistance, high prolificacy, and high fat deposition [20], Wannan Black pig makes an important contribution to the local commercial pork production. In our previous work, we identified several candidate genes that have mutations in Wannan Black pigs that are significantly associated with production traits [2123]. In recent decades, Wannan Black pigs have been in various unfavorable predicaments, leading to their low population density and diversity, greatly influenced by commercial, lean pig genotypes. From 1982 to 2019, the numbers of Wannan black pig sows changed from 6,688 to 360. Importantly, the risk of high level of inbreeding increased in small populations where the choice of breeding mates is limited [24]. On the other hand, some inbreeding depression in fitness-related traits of Wannan Black pigs has already been reported (Ministry of Agriculture and Rural Affairs of the People’s Republic of China, http://www.moa.gov.cn/). In order to provide reference for rational and efficient utilization of Wannan black pig germplasm resources, the frequency, number, and FROH of ROHs in Wannan black pig population were detected.

MATERIALS AND METHODS

Sequencing, genetic diversity and genetic relationship analysis

Our research was approved by the Anhui Agricultural University Animal Ethics Committee under permission No. AHAU20140215. In this study, 20 ear tissues of Wannan Black pigs (10 females and 10 males) were collected from the Wannan Black pig conservation farm, Huangshan City, Anhui province, People’s Republic of China). In a previous work, these samples were re-sequenced on the Illumina HiSeq X Ten platform (Illumina, San Diego, CA, USA) at Novogene Biotech Co., Ltd. (Beijing, China) with an average depth of ten generations. The analytical procedures related to the resequencing were adapted from Zhang et al [25]. Based on sequencing data, 501.52 G of raw data of the Wannan Black pig genome were obtained, which have been submitted to the National Center for Biotechnology Information (NCBI) database under the accession number PRJNA524263. A total of 21,316,754 SNPs were identified. We used the PLINK’-genome’ command to analyze the genetic relationship between the 20 pigs, as described by Purrcell et al [26].

Runs of homozygosity detection criteria

ROH were identified for each animal using PLINK v1.07 software [26]. The following criteria were chosen for ROH estimation [2729]: i) the minimum length of the filter input regions was set to 1 Mb; ii) one heterozygous and five missing calls were allowed per window to account for genotyping error; iii) the minimum number of SNPs was set to 100; and iv) the minimum quality of bases in the filter input regions was set to 10.

Inbreeding coefficient estimation

Based on ROH data, genomic inbreeding for each animal was estimated from ROH (FROH) as the ratio of the total length of the genome covered by ROH to the total length of the genome covered by SNPs or sequences, which is given as follows:

FROH=LROHLgenome

Where LROH is the total length of an individual’s ROH in the autosomal genome, and Lgenome is the length of the autosomal genome covered by the SNPs, which was 2,435 G in this study.

Detection of common runs of homozygosity

To identify genomic regions with a high frequency of ROH, we take the population as a unit, in which the ROH ratio of SNP sites in ROH was counted for each SNP site. Manhattan map was drawn according to ROH ratio of each SNP site. The ROH ratio where the top 1% is located is taken as the threshold line of high frequency SNP, and ROH island is obtained according to the distribution of SNP sites exceeding the threshold in the genome [30].

The gene content of the ROH islands was annotated using the annotation database provided by NCBI (https://www.ncbi.nlm.nih.gov). To further analyze the functions of the identified genes, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) analyses were performed using Database for Annotation, Visualization and Integrated Discovery (DAVID, v. 6.8) [31]. An extensive accurate literature search was then performed.

RESULTS

The results of the genetic relationship analysis of the 20 Wannan Black pigs (Supplementary Table S1) showed that the animals were not related, which was useful for elucidating the ROH information.

Runs of homozygosity detection and inbreeding coefficient

The descriptive statistics of ROH number and length by class are given in Table 1 and Figure 1. In total, 1,813 ROHs were identified in 20 Wannan Black pigs. Among the identified ROHs, 837 ranged between 100 to 500 kb, 449 between 500 to 1,000 kb, and 527 ROHs were >1,000 kb, with a mean ROH length of 0.69 MB. The longest segment was found in Sus Scrofa chromosome (SSC) 2, which was 3.42 Mb, and the shortest segment was 0.01 Mb. The total ROH number of Wannan Black pigs was composed mostly of a high number of shorter segments (100 to 500 kb and 500 to 1,000 kb) which accounted for 70.94% of all ROH detected. Although the number of ROHs of 100 to 500 kb was the largest, the proportion of the genome covered by them was relatively small compared to that of ROH segments longer than 1,000 kb.

Descriptive statistics of ROH number and length (in kb) by ROH length class

Figure 1

Distribution of the runs of homozygosity.

The numbers of ROH on each chromosome and the percentage of chromosomes covered by ROH in Wannan Black pigs are shown in Figure 2. The results indicated that ROHs were distributed unevenly, with SSC13 and SSC5 showing the highest (n = 267) and lowest (n = 48) number of ROHs, respectively. The highest genome coverage by ROH was observed on SSC6 (38.84%), whereas the lowest coverage was on SSC3 (9.40%).

Figure 2

Number of runs of homozygosity (ROH) longer per chromosome (bars) and average percentage of each chromosome covered by ROH (red line).

From the sequencing data, we inferred that 2,435,262,063 bp long autosomes were covered by SNPs and 1,274,547,767 bp long ROH were present on autosomes. Based on this information, the average FROH of Wannan Black pigs was 0.523.

Identification of runs of homozygosity islands

By determining the frequency of SNPs in ROH, we identified the genomic regions that were most commonly associated with ROHs in Wannan Black pigs. The results were plotted against the positions of the SNPs along the autosome. A total of 61 ROH islands were identified on SSC2, 3, 7, 8, 13, 15, and 16 (Figure 3; Supplementary Table S2). The longest ROH island was observed on SSC3, whereas the shortest one was observed on SSC13.

Figure 3

Manhattan plot of incidence of each single-nucleotide polymorphism (SNP) in the runs of homozygosity (ROH) across individuals. The red line represents the 30% threshold.

Functional annotation of genes in the runs of homozygosity islands

Furthermore, the 61 ROH islands harbored 105 candidate genes, whereas 19 ROH islands did not harbor any genes (Supplementary Table S2). The function of the genes in the identified regions was analyzed (Supplementary Table S3, S4). Eight genes including eukaryotic translation initiation factor 1 (EIF1), nucleoporin 54 (Nup54), mitochondrial ribosomal protein S9 (MRPS9), and ribosomal protein S7 (Rps7) were related to signal translation processes, and four genes (four and a half LIM domains 2 [FHL2], leucine rich repeat containing 4C [LRRC4C], cytoplasmic polyadenylation element binding protein 4 ([PEB4], and EPH receptor A3 [EPHA3]) were involved in development. The remaining genes were involved in the immune system, nucleotide metabolism, amino acid metabolism, cellular community, and other important biological processes. Based on the results of the literature query, five genes were associated with economic traits in animals (Table 2).

Candidate genes located in genomic regions with a high frequency of runs of homozygosity associated with economic traits

DISCUSSION

ROH are continuous homozygous segments of the DNA sequence in diploid genomes and are widely distributed in humans and livestock populations [32]. Inbreeding level has a significant influence on abundance, length, and number of ROH. In addition, genetic drift, natural and artificial selection, population bottlenecks among several other factors promote ROH buildup in a population [33]. However, ROH segments are not randomly distributed across the genome. In Chinese Merino sheep, Ovies aries chromosome (OAR) 21 and OAR3 exhibited the highest and lowest coverage across the chromosomes by ROH, respectively [34]. Xu et al [28] reported that in Jinhua pigs, the highest number of ROH per chromosome was on SSC6 with 1,672 segments, and the lowest was on SSC16 with 485 segments. In agreement with these findings, the SSC13 of Wannan Black pigs harbored the highest number of ROHs (n = 267), whereas SSC5 had the lowest number (n = 48).

The information regarding the abundance, length, and number of ROHs on the chromosome is valuable in exploring the demographic history of livestock species. With the availability of genomics tools, numerous polymorphisms can be synchronously genotyped across the genome and help to identify long stretches of homozygous genotypes [35]. In this work, the distribution and frequency of ROH in Wannan Black pig population were analyzed by resequencing (10×) data. Several studies have reported that long ROHs are associated with recent events of inbreeding, whereas short ROHs have been broken by repeated recombination and represent inbreeding that took place several generations ago [36,37]. A total of 1,813 ROHs were identified and short segments (lengths <1,000 kb) account for 70.93%. Wannan Black pig may therefore have undergone inbreeding in a distant generation, which concurs with what we know of the history of the Wannan Black pig.

Conventionally, coefficients of inbreeding were estimated from pedigree records. However, native breed conservation may suffer imperfections, such as the inaccurate pedigree records and chaotic consanguinity. Hence, ROH can be used as a predictor of whole genome inbreeding levels. FROH is defined as the proportion of the autosomal genome in ROH exceeding a specified length. McQuillan et al [38] first reported that FROH strongly correlates with the inbreeding coefficient of pedigrees in European populations. Our results showed that the average FROH in Wannan Black pigs was higher than that of other pig breeds (Złotnicka Spotted, Polish Landrace, 0.287, 0.171, respectively) [39]. Li et al [40] analyzed the genetic diversity of the same Wannan Black pig population (from the same conservation farm), and reported that these animals have more allele frequencies of microsatellite loci deviated from that expected under Hardy-Weinberg Equilibrium than introduced pig breeds (Duroc, Berkshire, Pietrain, Landrance, and Yorkshire) and other native pig breeds of Anhui province (Dinyuan pig, Wei pig, and Anqing six-end-white pig), revealing that Wannan Black pig has a higher inbreeding level than other pigs. With a small population size, inbreeding increases at a higher rate, leading to loss of alleles, which is impossible to counterbalance without migration [41]. Introducing new consanguinity, and increasing the effective population size of Wannan Black pig is necessary in the future.

Several candidate genes associated with production traits were identified in the ROH islands. Calmodulin-lysine N-methyltransferase (CAMKMT) is located on SSC3 of Wannan Black pig, which is significantly correlated with the number of embryonic cells and bone surface area [42]. The polymorphisms SNP07, SNP28, and SNP31 of CAMKMT polymorphisms are associated with sheep growth traits [43]. CPEB4 is involved in bone, muscle, fat, and lung development, which is related to growth traits in Large White pig populations by genome wide association study [44]. The expression of collagen type III alpha 1 chain (Col3a1) was reported to affect intramuscular collagen, which is associated with meat quality in pigs [45]. Co-analysis of ROH and the signature selection identified GULP PTB domain containing engulfment adaptor 1 (GULP1) and kelch repeat and BTB domain containing 8 (KBTBD8) genes. GULP1 is located on SSC15 and plays key roles in cancer suppression and immune responses [46], whereas KBTBD8 is involved in cytoskeleton arrangement, regulation of cell morphology, and idiopathic short stature [47,48]. These two genes might be associated with the characteristics of Wannan Black pigs that adapt to natural and artificial selection. Therefore, elucidating the potential involvement of GULP1 and KBTBD8 in pig productive performance warrants further research.

IMPLICATIONS

We identified runs of homozygosity s across the genome of the Wannan Black pig, the Chinese native breeds of the Anhui province, for the first time. The frequency, numbers, and FROH were obtained and reflected the inbreeding history in the Wannan Black pigs. Our results showed that ancient and recent inbreeding had an influence on the genome and revealed a high level of inbreeding in the existing population. Furthermore, several candidate genes associated with important biological processes were identified. In conclusion, our work will help to save and utilize the germplasm resource of Wannan Black pig.

Supplementary Information

Notes

CONFLICT OF INTEREST

We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.

ACKNOWLEDGMENTS

This work was supported by the National Natural Science Foundation of China (31972531), Special Projects of the Central Government Guiding Local Science and Technology Development of Anhui Provincial (202007d06020016), Anhui Provincial Science and Technology Major Project (17030 701008, 17030701061) and the Anhui Provincial Modern Industrial Technology System of Swine.

References

1. Chen H, Huang T, Zhang Z, et al. Genome-wide association studies and meta-analysis reveal novel quantitative trait loci and pleiotropic loci for swine head-related traits. J Anim Sci 2017;95:2354–66. https://doi.org/10.2527/jas.2016.1137 .
2. Maselli V, Polese G, Larson G, et al. A dysfunctional sense of smell: the irreversibility of olfactory evolution in free-living pigs. Evol Biol 2014;41:229–39. https://doi.org/10.1007/s11692-013-9262-3 .
3. Novoselova TA, Meuwissen MPM, Jongbloed AW, Huirne RBM. Expected economic performance of genetic modification in pork production NJAS-Wageningen. J Life Sci 2013. 64–659–15. https://doi.org/10.1016/j.njas.2012.03.002 .
4. Muhammad H, Iqbal Z, Bashir Q. Genetic polymorphism among five species of Indus River carps on the basis of random amplified polymorphic DNA (RAPD) analysis. Rend Lincei Sci Fis Nat 2017;28:73–80. https://doi.org/10.1007/s12210-016-0576-9 .
5. Gorlov IF, Kolosov YA, Shirokova NV, et al. GDF9 gene polymorphism and its association with litter size in two Russian sheep breeds. Rend Lincei Sci Fis Nat 2018;29:61–6. https://doi.org/10.1007/s12210-017-0659-2 .
6. Forutan M, Mahyari SA, Baes C, et al. Inbreeding and runs of homozygosity before and after genomic selection in North American Holstein cattle. BMC Genomics 2018;19:98. https://doi.org/10.1186/s12864-018-4453-z .
7. Christensen OF, Nielsen B, Su G, et al. A bivariate genomic model with additive, dominance and inbreeding depression effects for sire line and three-way crossbred pigs. Genet Sel Evol 2019;51:45. https://doi.org/10.1186/s12711-019-0486-2 .
8. Do C, Yang C, Choi J, et al. The outcomes of selection in a closed herd on a farm in operation. Asian-Australas J Anim Sci 2015;28:1244–51. https://doi.org/10.5713/ajas.14.0962 .
9. Ouborg NJ, Pertoldi C, Loeschcke V, Bijlsma R, Hedrick PW. Conservation genetics in transition to conservation genomics. Trends Genet 2010;26:177–87. https://doi.org/10.1016/j.tig.2010.01.001 .
10. Sell-Kubiak E, Czarniecki A, Strabel T. Challenges in inbreeding estimation of large populations based on Polish Holstein-Friesian cattle pedigree. J Appl Genet 2018;59:313–23. https://doi.org/10.1007/s13353-018-0441-x .
11. Simino J, Shi G, Arnett D, Broeckel U, Hunt SC, Rao DC. Variants on chromosome 6p22.3 associated with blood pressure in the HyperGEN study: follow-up of FBPP quantitative trait loci. Am J Hypertens 2011;24:1227–33. https://doi.org/10.1038/ajh.2011.140 .
12. Purfield DC, Berry DP, McParland S, Bradley DG. Runs of homozygosity and population history in cattle. BMC Genet 2012;13:70. https://doi.org/10.1186/1471-2156-13-70 .
13. Zavarez LB, Utsunomiya YT, Carmo AS, et al. Assessment of autozygosity in Nellore cows (Bos indicus) through high-density SNP genotypes. Front Genet 2015;6:5. https://doi.org/10.3389/fgene.2015.00005 .
14. Kim ES, Cole JB, Huson H, et al. Effect of artificial selection on runs of homozygosity in US Holstein cattle. PLoS One 2013;8:e80813. https://doi.org/10.1371/journal.pone.0080813 .
15. Ceballos FC, Joshi PK, Clark DW, Ramsay M, Wilson JF. Runs of homozygosity: windows into population history and trait architecture. Nat Rev Genet 2018;19:220–34. https://doi.org/10.1038/nrg.2017.109 .
16. Purfield DC, McParland S, Wall E, Berry DP. The distribution of runs of homozygosity and selection signatures in six commercial meat sheep breeds. PLoS One 2017;12:e0176780. https://doi.org/10.1371/journal.pone.0176780 .
17. Metzger J, Karwath M, Tonda R, et al. Runs of homozygosity reveal signatures of positive selection for reproduction traits in breed and non-breed horses. BMC Genomics 2015;16:764. https://doi.org/10.1186/s12864-015-1977-3 .
18. Shi L, Wang L, Liu J, et al. Estimation of inbreeding and identification of regions under heavy selection based on runs of homozygosity in a Large White pig population. J Anim Sci Biotechnol 2020;11:46. https://doi.org/10.1186/s40104-020-00447-0 .
19. China National Commission of Animal Genetic Resources. Animal genetic resources in China pigs Beijing, China: China Agricultural Press; 2011.
20. Ding YY, Zhang W, Zhang MQ, et al. Functional and association studies of the cholesteryl ester transfer protein (CETP) gene in a Wannan Black pig model. Anim Genet 2015;46:702–6. https://doi.org/10.1111/age.12370 .
21. Ding X, Ding Y, Zhang C, Huang L, He X, Yin Z, et al. Polymorphisms and sequencing analysis in 5′ upstream and intron 2 region of H-FABP gene in Anhui Wannan Black pig population. J Anim Vet Adv 2011;10:1568–72.
22. Zhang XD, Zhu HY, Zhou J, et al. Relationship between polymorphisms in exon 10 of FSHR gene and litter size in swine. Genet Mol Res 2015;14:8252–61. https://doi.org/10.4238/2015.July.27.13 .
23. Ding X, Zhang X, Yang Y, et al. Polymorphism, expression of natural resistance-associated macrophage protein 1 encoding gene (NRAMP1) and Its association with immune traits in pigs. Asian-Australas J Anim Sci 2014;27:1189–95. https://doi.org/10.5713/ajas.2014.14017 .
24. Tang GQ, Xuea J, Lian MJ, et al. Inbreeding and genetic diversity in three imported swine breeds in china using pedigree data. Asian-Australas J Anim Sci 2013;26:755–65. https://doi.org/10.5713/ajas.2012.12645 .
25. Zhang W, Yang M, Wang YL, et al. Genomic analysis reveals selection signatures of the Wannan Black pig during domestication and breeding. Asian-Australas J Anim Sci 2020;33:712–21. https://doi.org/10.5713/ajas.19.0289 .
26. Purcell S, Neale B, Todd-Brown K, et al. PLINK: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 2007;81:559–75. https://doi.org/10.1086/519795 .
27. Peripolli E, Metzger J, de Lemos MVA, et al. Autozygosity islands and ROH patterns in Nellore lineages: evidence of selection for functionally important traits. BMC Genomics 2018;19:680. https://doi.org/10.1186/s12864-018-5060-8 .
28. Xu Z, Sun H, Zhang Z, et al. Assessment of autozygosity derived from runs of homozygosity in Jinhua pigs disclosed by sequencing data. Front Genet 2019;10:274. https://doi.org/10.3389/fgene.2019.00274 .
29. Zhao G, Zhang T, Liu Y, et al. Genome-wide assessment of runs of homozygosity in Chinese Wagyu beef cattle. Animals 2020;10:1425. https://doi.org/10.3390/ani10081425 .
30. Mastrangelo S, Ciani E, Sardina MT, et al. Runs of homozygosity reveal genome-wide autozygosity in Italian sheep breeds. Anim Genet 2018;49:71–81. https://doi.org/10.1111/age.12634 .
31. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009;4:44–57. https://doi.org/10.1038/nprot.2008.211 .
32. Gibson J, Morton NE, Collins A. Extended tracts of homozygosity in outbred human populations. Hum Mol Genet 2006;15:789–95. https://doi.org/10.1093/hmg/ddi493 .
33. Peripolli E, Munari DP, Silva MVGB, Lima ALF, Irgang R, Baldi F. Runs of homozygosity: current knowledge and applications in livestock. Anim Genet 2017;48:255–71. https://doi.org/10.1111/age.12526 .
34. He SG, Di J, Han B, et al. Genome-wide scan for runs of homozygosity identifies candidate genes related to economically important traits in Chinese merino. Animals 2020;10:524. https://doi.org/10.3390/ani10030524 .
35. Bosse M, Megens HJ, Madsen O, et al. Regions of homozygosity in the porcine genome: consequence of demography and the recombination landscape. PLoS Genet 2012;8:e1003100. https://doi.org/10.1371/journal.pgen.1003100 .
36. Marras G, Gaspa G, Sorbolini S, et al. Analysis of runs of homozygosity and their relationship with inbreeding in five cattle breeds farmed in Italy. Anim Genet 2015;46:110–21. https://doi.org/10.1111/age.12259 .
37. Gomez-Raya L, Rauw WM, Dunkeiberger JR, Dekkers JCM. Autozygosity and genetic differentiation of landrace and large white pigs as revealed by the genetic analyses of crossbreds. Front Genet 2019;10 https://doi.org/10.3389/fgene.2019.00739 .
38. McQuillan R, Leutenegger A-L, Abdel-Rahman R, et al. Runs of homozygosity in european populations. Am J Hum Genet 2008;83:658. https://doi.org/10.1016/j.ajhg.2008.10.009 .
39. Szmatola T, Jasielczuk I, Semik-Gurgul E, et al. Detection of runs of homozygosity in conserved and commercial pig breeds in Poland. J Anim Breed Genet 2020; https://doi.org/10.1111/jbg.12482 .
40. Li Q, Wang C, Yang J, et al. Genetic diversity analysis based on microsatellite markers in 5 Anhui local and 5 introduced pig breeds. Swine Prod 2018;3:67–72.
41. Woolliams J. Genetic contributions and inbreeding. p. 147–165. In : Oldenbroek K, ed. Utilization and conservation of farm animal genetic resources Wageningen, the Netherlands: Wageningen Academic Publishers; 2007. p. 147–65.
42. Abdelmohsen K, Srikantan S, Yang X, et al. Ubiquitin-mediated proteolysis of HuR by heat shock. EMBO J 2009;28:1271–82. https://doi.org/10.1038/emboj.2009.67 .
43. Ma X, Guan L, Xuan J, et al. Effect of polymorphisms in the CAMKMT gene on growth traits in Ujumqin sheep. Anim Genet 2016;47:618–22. https://doi.org/10.1111/age.12455 .
44. Meng Q, Wang K, Liu X, et al. Identification of growth trait related genes in a Yorkshire purebred pig population by genome-wide association studies. Asian-Australas J Anim Sci 2017;30:462–9. https://doi.org/10.5713/ajas.16.0548 .
45. Bao X, Zeng Y, Wei S, et al. Developmental changes of Col3a1 mRNA expression in muscle and their association with intramuscular collagen in pigs. J Genet Genomics 2007;34:223–8. https://doi.org/10.1016/s1673-8527(07)60023-x .
46. Maldonado L, Brait M, Begum S, et al. GULP1, a potential tumor suppressor gene in ovarian tumors and its utility as a biomarker. Cancer Res 2010;70:4891. https://doi.org/10.1158/1538-7445.Am10-4891 .
47. Luhrig S, Kolb S, Mellies N, Nolte J. The novel BTB-kelch protein, KBTBD8, is located in the Golgi apparatus and translocates to the spindle apparatus during mitosis. Cell Div 2013;8:3. https://doi.org/10.1186/1747-1028-8-3 .
48. Avila-Arcos MC, McManus KF, Sandoval K, et al. Population history and gene divergence in native Mexicans inferred from 76 human exomes. Mol Biol Evol 2020;37:994–1006. https://doi.org/10.1093/molbev/msz282 .

Article information Continued

Figure 1

Distribution of the runs of homozygosity.

Figure 2

Number of runs of homozygosity (ROH) longer per chromosome (bars) and average percentage of each chromosome covered by ROH (red line).

Figure 3

Manhattan plot of incidence of each single-nucleotide polymorphism (SNP) in the runs of homozygosity (ROH) across individuals. The red line represents the 30% threshold.

Table 1

Descriptive statistics of ROH number and length (in kb) by ROH length class

ROH length (kb) ROH number Percentage (%) Mean length (Mb) (mean±SD) Genome coverage (%)
100–500 837 46.17 0.27±0.11 9.16
500–1,000 449 24.77 0.72±0.14 13.35
>1,000 527 29.07 1.34±0.88 29.05
Total >100 1,813 100 0.69±8.81 51.56

ROH, runs of homozygosity.

Table 2

Candidate genes located in genomic regions with a high frequency of runs of homozygosity associated with economic traits

Chromosome Location (bp) Size (bp) Gene Gene function
3 95,854,107 – 95,918,246 64,139 calmodulin-lysine N-methyltransferase Growth
15 92,482,944 – 93,357,026 874,082 GULP PTB domain containing engulfment adaptor 1 Immune responses
15 93,556,914 – 93,595,678 38,764 collagen type III alpha 1 chain Meat quality
16 50,185,113 – 505,01,620 316,507 cytoplasmic polyadenylation element binding protein 4 Growth
13 48,486,919 – 48,499,992 13,073 kelch repeat and BTB domain containing 8 Growth