Next Article in Journal
An Integrated Linkage Map of Three Recombinant Inbred Populations of Pea (Pisum sativum L.)
Previous Article in Journal
The Difference in the Proportions of Deleterious Variations within and between Populations Influences the Estimation of FST
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of the Genetic Architecture of Pigs Subjected to Breeding Intensification

1
Don State Agrarian University, 346493 Persianovski, Russia
2
Centre for Strategic Planning and Management of Biomedical Health Risks Disabled, 119435 Moscow, Russia
3
Federal Research Center for Animal Husbandry Named after Academy Member LK. Ernst, 142132 Dubrovitsy, Russia
4
Mendeleev University of Chemical Technology, 125480 Moscow, Russia
*
Author to whom correspondence should be addressed.
Genes 2022, 13(2), 197; https://doi.org/10.3390/genes13020197
Submission received: 26 December 2021 / Revised: 18 January 2022 / Accepted: 21 January 2022 / Published: 22 January 2022

Abstract

:
Pigs are strategically important animals for the agricultural industry. An assessment of genetic differentiation between pigs, undergone and not undergone to selection intensification, is of particular interest. Our research was conducted on two groups of Large White pigs grown on the same farm but in different years. A total of 165 samples were selected with 78 LW_А (n = 78, the Russian selection) and LW_B (n = 87, a commercial livestock). For genotyping, we used GeneSeek® GGP Porcine HD Genomic Profiler v1 (Illumina Inc, San Diego, CA, USA). To define breeding characteristics of selection, we used smoothing FST and segment identification of HBD (Homozygous-by-Descent). The results of smoothing FST showed 20 areas of a genome with strong ejection regions of the genome located on all chromosomes except SSC2, SSC3, and SSC8. The average realized autozygosity in Large White pigs of native selection was in (LW_A)—0.21, in LW_В—0.29. LW_А showed 13,338 HBD segments, 171 per one animal, and LW_B—15,747 HBD segments, 181 per one animal. The ejections found by the smoothing FST method were partially localized in the HBD regions. In these areas, the genes ((NCBP1, PLPPR1, GRIN3A, NBEA, TRPC4, HS6ST3, NALCN, SMG6, TTC3, KCNJ6, IKZF2, OBSL1, CARD10, ETV6, VWF, CCND2, TSPAN9, CDH13, CEP128, SERPINA11, PIK3CG, COG5, BCAP29, SLC26A4) were defined. The revealed genes can be of special interest for further studying their influence on an organism of an animal since they can act as candidate genes for selection-significant traits.

1. Introduction

A human noticed long ago that some traits are displayed in animals of one species differently, and he selected the individuals with traits interesting to him and raised them artificially on his farm. Thus, there appeared a selection traditionally based on breeding and selection aimed at fixing desirable traits in the population. The development of genetics and molecular biology made it possible to considerably raise the efficiency of selection and breeding work, which also influenced the rate of improvement of animals. A long-term, targeted selection by the same breeding important traits leads to the appearance of the so-called selection signatures in the genome of farm animals and is associated with specific traits [1].
Selection signatures were found for the genes-candidates linked to such industrial characteristics as growth rate and development [2,3], reproduction [4], meat quality [5], biological regulation and metabolism [6], level of immunity, adaptability [7], processes providing development of a brain [8], and morphological variations, the form of ears [9], and body length [10].
When searching for selection signatures, the FST method is extremely interesting. Value FST is a differentiation measure between populations [11]. The FST locus value is calculated as the ratio of the variance of allele frequencies between populations and the sum of variances within and between populations. The locus with more significant values of FST in comparison with other loci may indicate positive selection [12,13]. To date, various modifications of this method have been presented, but undoubtedly, it remains the most widespread and reliable one to identify genome traits of selection signatures between observable populations [14,15].
Method FST is used to identify adjacent regions of the genome during selection and is useful for analyzing distantly related populations as it reveals subtle differences between them [16]. Smoothed method FST is based on the model of Nicholson pure drift [17], according to which separate SNP clusters in genome windows are calculated as average values.
Selection signatures can also be localized in homozygous areas of various lengths [1]. When creating breeds of farm animals, the accumulation of homozygosity allows pure-bred animals to possess certain qualities and steadily pass them on to their offspring [18]. The use of homozygosity patterns of ROH (Runs Of Homozygosity) makes it possible to reveal long homozygous areas in a genome [19]. In the research, it is reasonable to use the method proposed by Drouet and Gaultier based on the model of HBD multiple classes (homozygous-by-Descent). This method allows estimating autozygosity according to the age of the ancestors.
Pigs are strategically important animals for agriculture. In this connection, over the past decades, targeted work has been carried out to increase selection-significant indicators [20]. This enabled us to considerably increase the pigs’ reproductive, growth, and meat indexes. However, alongside high efficiency in pigs, there appeared various anomalies, congenital defects, problems with limbs, and susceptibility to various diseases [21]. Thereby an assessment of genetical differentiation between the pigs that have undergone and have not undergone selection intensification is of particular interest. In this aspect, we conducted research using pigs raised on the territory of the Russian Federation at different periods of time and identified the signatures of selection in pigs due to trends in different socio-economic conditions [22]. However, studying modern livestock of Large White pigs has shown that formations of selection signatures are influenced by selection centers themselves since each of them implements its own breeding strategy. In this connection, we focused our work on pigs raised on one farm of the Russian Federation, but at different periods of time. In addition to the FST method, we investigated homozygosity areas using a model of plural HBD classes, defined genome signatures identified by two methods, determined QTL enrichment, and positioned genes in these areas.

2. Materials and Methods

2.1. Animals

Anesthesia, euthanasia, or any animal sacrifice was not used to conduct this study. This study did not involve any endangered or protected species. According to standard monitoring procedures and guidelines, the participating holdings specialists collected tissue samples, following the ethical protocols outlined in the Directive 2010/63/EU (2010). The pig ear samples (ear pluck) were obtained as a general breeding monitoring procedure. The collection of ear samples is a standard practice in pig breeding [23].

2.2. Sampling and Genotyping

For our work, we chose Large White pigs, which were kept on the same farm but in different years. Pigs of the LW_A group belonged to the Russian selection, which is based on pigs of the Large White breed, imported from England in 1923–1931. Long-term breeding work, taking into account the local climate, changed the English type of Large White pigs, and, in fact, a new domestic Large White breed was created, which at that time surpassed the English in many respects. At the beginning of the 21st century, these pigs almost completely disappeared, and imported pigs began to be imported into the Russian Federation [24]. LW_A were distinguished by good adaptation and resistance to various diseases were less whimsical to the conditions of keeping and feeding, but the pigs of imported selection were significantly superior in growth rate, reproductive performance, and thinner fat. Group LW_A (date of assembly 2008–2010) and LW_B (date of assembly 2014–2016). Pigs of the LW_A group belonged to the Russian selection, and pigs LW_B belonged to commercial livestock, which was delivered to the farm from Europe in 2013. For work, 165 samples were selected, 78 LW_А and 87 LW_B. Genomic DNA was extracted from tissue (ear pinch) using a set of DNA-Extran-2 reagents (OOO NPF Sintol, Russia). For genotyping, we used GeneSeek® GGP Porcine HD Genomic Profiler v1, which included 68,516 SNPs evenly distributed with an average spacing of 25 kb. (Illumina Inc, San Diego, CA, USA). The total genotyping rate was 0.99.

2.3. Data Analysis

To make relations between populations visual we conducted SVD (singular value decomposition) by means of basic svd function in R. The Heatmap graph was plotted on the basis of the GRM matrix. To define selection characteristics, we used methods of smoothing FST. For smoothing FST a filtration of the data hwe 1 × 10 −7 maf 0.01—geno 0.2—mind 0.2—indep-pairwise 50 5 0.2 42,442 variants passed the QC filters and were retained for further analysis. To filter the noise obtained as a result of the FST calculation, the lokern smoothing algorithm was applied: Kernel Regression Smoothing with Local or Global Plug-in Bandwidth of the lokern package in R [25] with the n.out = 424 parameter, which approximately corresponds to one point for every 100 SNPs and allows smoothing SNP data set. The value of the x.out parameter was used to match the smoothed values against the SNP reference map and its position. The smoothest FST values, corresponding to 0.999%, were identified and translated into genomic positions of Sus scrofa 11.1, and the gene content of each region was analyzed.
To identify HBD segments and to assess autozygosity (or the coefficient of inbreeding), we used the multiple HBD classes model presented in RZooRoH package [26,27]. The method was insensitive to MAF filtration and rather resistant to the structure of LD. In this connection, the data filtration by MAF and LD was not conducted. The Rk coefficients were set from 2 to 516 (2, 4, 8, 16, 32, 64, 128, 256, 512). The inbreeding coefficient was calculated as the sum of autozygosity for all HBD classes. The total number of HBD segments, the average number per individual, the average length of the segment per individual, and the distribution of segments (and their average length) on the chromosomes were assessed for each group. We defined SNP frequency (%) in the found HBD segments and, for each group, chose top HBD provided that HBD frequency was not less than 60% and included at least 10 SNPs. Based on the results of the 2 methods (FST and HBD), regions of the genome and genes were identified, probably associated with the intensification of the selection process in commercial pigs.

2.4. Search and the Analysis of QTL Enrichment

The search of QTL, genes, and the QTL enrichment analysis performed in Genomic Annotation in Livestock for positional candidate LOci (GALLO) was an R package Ensembl genome browser [28], and also a literature search was also carried out manually for the presence of data on the associations of genes with any traits in humans and animals.

3. Results

To assess the genetic structure of the studied populations of Large White pigs, we used SVD and Heatmap. Figure 1A,B show that pigs of groups LW_A and LW_B have their own individual cluster.
The results of smoothing FST showed 20 regions of the genome with strong outliers located on all chromosomes, with the exception of SSC2, SSC3, and SSC8 (additional Table S1). These areas overlap with quantitative trait loci (QTLs), of which Meat and Carcass traits were most represented (Figure 2A). Based on the analysis of QTL enrichment with the most significant enrichment, the signs of pH 24 hr post-mortem (lion), meat color a *, and Cortisol level were identified (Figure 2B).
On average, the realized autozygosity was 0.21 for LW_A and 0.29 for LW_B (Figure 3). The maximum values (0.34) were recorded in the LW_B group and the minimum (0.15) in the LW_A group.
The class Rk = 128 contributed greatly to autozygosity in LW_B (proportion of the genome about 0.1) (Figure 4A,B; additional Table S2). Herewith, the contribution of the classes Rk = 128 (0.05) and the class Rk = 256 (0.05) can be traced in LW_A. Variations in individual levels of autozygosity in both groups did not deviate much from the midline (Figure 4C,D). HBD class Rk = 512 was absent in all pigs in the study groups.
In general, LW_A had 13,338 HBD segments, an average of 171 per animal; LW_B has 15,747 HBD segments, an average of 181 per animal. The largest length of HBD segments was determined for LW_A (138.52 Mb, 1759 Number SNP, SSC1). The average length of HBD segments for LW_A was about 2.47 Mb (54 Number SNP), for LW_B 3.41 Mb (77 Number SNP). The average length of HDB segments on chromosomes (taking into account different classes) is shown in Figure 5. Segments of the Rk = 2 class were determined only in pigs from the LW_A group. These segments were located in SSC1, and their average length was 108.67 Mb. Segments of class Rk = 4 for LW_A were defined at SSC1, SSC2, SSC4, SSC5, SSC6, SSC9, SSC13, SSC17; LW_B has SSC1, SSC2, SSC8, SSC9, SSC13, and SSC15. HBD segments of class Rk = 8 were absent in LW_A on SSC10 and in LW_B on SSC10 and SSC16. Starting from Rk = 16 and further up to Rk = 256, HBD segments were relatively evenly distributed on all chromosomes in LW_A and LW_B.
The SNP frequencies (%) in the detected HBD were estimated for each group of pigs and plotted against the position of the SNP in the autosomes (Figure 6 and Figure 7).
For each group, the top HBD were selected, provided that the HBD frequency was at least 60% and at least 10 SNPs were included. As a result, LW_A has 4 regions located in SSC1 (Table 1). LW_B has 10 regions, of which 5 were in SSC1, 2 in SSC6, and one each in SSC10, SSC14, and SSC15. The topHBD regions did not overlap between groups.
In both groups, the top HBD areas overlapped with QTLs, among which the most represented were the signs of Meat and Carcass traits (Figure 8). In pigs LW_A, relative to LW_B, the QTL type of Exterior, Health and Production was more represented. In LW_B pigs, QTL Reproduction was more represented. Using the analysis of QTL enrichment in the top HBD, LW_A had the most represented characteristics: Ph Logissmus Dorsi, carcass weight (hot), body weight (weaning), backfat at last rib, average daily gain. LW_B has shoulder subcutaneous fat thickness, shear force, loin muscle area, fat area percentage in the carcass, estimated carcass lean content, and dressing percentage.
In the top HBD, genes encoding proteins and having names (according to Ensembl genome browser 104) were identified, but genes encoding lncRNA, snoRNA, snRNA were also represented in these regions with a high frequency (additional Table S1).
The ejections detected by the smoothing FST method were partially localized in the HBD regions. Table 2 shows the areas in which high outliers were found in FST smoothing and their frequencies in the HBD areas.

4. Discussion

The average length of HBD, chromosome distribution, and the genome proportion covered by HBD can be used as indicators of the origin and history of a population, as well as reflecting events of artificial selection. The Large White breed was created in the 1870s–1880s and officially received its name in 1885 [24]. According to the results of our research, the shortest HBD was inherited from ancestors about 128 years ago, which in general was exactly the period of the formation of the Large White breed. Subsequently, the Large White breed took part in the creation and improvement of most modern European breeds, as well as local breeds created on the territory of the USSR. Our studies showed that the total autozygosity of LW_B pigs was 0.29 proportion of the genome, while the contribution of the Rk = 128 class segments was about 0.1. On this basis, we may assume that ancestors contributed greatly to the autozygosity of this group about 64 years ago. This period was superimposed on the period of growth of intensification processes in pig breeding and the formation of commercial livestock, characterized by high production indicators [29]. In its turn, the total autozygosity in pigs of the LW_A group was 0.21 proportion of the genome, but the dominant contribution of any of the classes was not observed because the contribution of ancestors about 64 (Rk = 128) and 128 (Rk = 256) years ago amounted to 0.05 shares of autozygosity. It can be assumed that the period of intensification was reflected in the LW_A population, but in a much smaller volume, which made it possible to preserve the signatures inherited from more distant ancestors.
Common autozygosity areas in a population identify selection hotspots [30]. In the groups under study, the top HBD was determined, provided that their frequency was less than 60% and included at least 10 SNPs. In general, in both groups, QTLs for Meat and Carcass were most represented, but in group LW_A they stood out to a greater extent with traits associated with carcass weight, and in LW_B—traits associated with obtaining lean pork.
Ejections identified by FST smoothing were partially localized in HBD segments. For example, SNP rs81349176 (SSC1) occured with a frequency of 0.73 in HBD in LW_A. The functional significance of SNP rs81349176 was difficult to interpret since it was localized in the intragenic region. However, it is interesting to mention that the adjacent genes CDH19 and CDH7 belong to the cadherin family, which plays a key role in the regulation of adhesion. Dysregulation of adhesion molecules often causes various diseases, including inflammation and tumors [31,32,33]. Earlier, we also suggested a hypothesis about the connection of cadregins with Capped Hock in Pig [34]. In the future, it is vital to study cadregins in more detail in terms of limb tumors since this is a source of significant economic losses in pig production [35,36].
Genes were determined in the regions identified by smoothing FST (NCBP1, PLPPR1, GRIN3A, NBEA, TRPC4, HS6ST3, NALCN, SMG6, TTC3, KCNJ6, IKZF2, OBSL1, CARD10, ETV6, VWF, CCND2, TSPAN9128, CDH13, PIK3CG, COG5, BCAP29, SLC26A4), which may be selection signals related to growth, conformation, health, reproductive performance, and meat quality. The identified genes were associated with the structural and functional work of the cell (CARD10, TSPAN9, PIK3CG, COG5, BCAP29). These genes are involved in the transmission of apoptosis signals (CARD10, BCAP29), play a certain role in the regulation of development, activation, growth, and motility of cells (TSPAN9), as well as in maintaining the structural and functional integrity of the epithelium and regulation of cytotoxicity in NK cells (PIK3CG), morphology and functions of the Golgi complex (COG5) [37].
It is interesting to note the variants of the genes influencing the function of a hemopoiesis and hereditary diseases of the circulatory, cardiovascular system, and other pathologies (HS6ST3, IKZF2, ETV6, SMG6, SLC26A4, VWF). The genes participate in proliferation and differentiation, adhesion, migration, inflammation, fibrillation, and other various processes (HS6ST3), in the regulation of development of lymphocytes (IKZF2), play a role in the hemopoiesis and malignant transformation (ETV6), and are associated with an increased risk of ischemic heart disease (SMG6), inherited hearing loss in domestic animals (SLC26A4).
The genes regulating important components of the nervous system have also been identified (KCNJ6, TRPC4, NALCN, GRIN3A, NBEA). The gene mutations are associated with severe developmental delay, facial dysmorphism and mental retardation, reduced cognitive ability (KCNJ6, NALCN). They play an important role in dopamine-related processes, including addiction and attention (TRPC4), physiological and pathological processes in the central nervous system (GRIN3A), autism (NBEA), etc.
The revealed genes can be of special interest for further studying of their impact on an organism of an animal as they can represent themselves as genes-candidates bound to the physiological features of an organism providing high industrial indexes, but also probably associated with changes of the nervous system, various pathologies, including the number of circulatory and cardiovascular systems.

5. Conclusions

Investigations aimed at studying the presence and localization of selection signatures (FST), as well as the identification of areas of homozygosity (HBD) in two groups of pigs bred at different times on the same farm, enabled us to identify differences between populations. The presence of QTLs located in areas of homozygosity and associated with traits, the improvement of which was aimed at selection breeding work, makes such areas the most promising for the search for potential candidate genes associated with the level of productivity and the presence of diseases.
In the genome regions determined by using the FST and HBD methods, we identified genes that may have contributed to the changes associated with the intensification of the selection process in pigs. In general, the results presented in our work show promising prospects for genome scanning using FST and HBD methods for studying population history, as well as for identifying genomic regions and genes associated with important economic traits and various pathologies.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes13020197/s1, Table S1: Gene ID Fst, HBD_LW_A, HBD_LW_B; Table S2: Summary statistics of realized autozygosity in each HBD class (as a percentage of the genome).

Author Contributions

A.K.: data curation, L.G.: data curation, project administration, writing—original draft, writing—review and editing, M.K., T.R.: writing—original draft, N.B.: investigation, E.R.: visualization, I.B.: investigation, O.K.: validation, Y.P.: Investigation, O.T.: investigation, S.B.: data curation, investigation, software, writing—original draft, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Russian Scienсe Foundation (RSF) Project No. 19-76-10012.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhan, H.; Zhang, S.; Zhang, K.; Peng, X.; Xie, S.; Li, X.; Zhao, S.; Ma, Y. Genome-Wide Patterns of Homozygosity and Relevant Characterizations on the Population Structure in Piétrain Pigs. Genes 2020, 11, 577. [Google Scholar] [CrossRef]
  2. Wang, Y.; Ding, X.; Tan, Z.; Xing, K.; Yang, T.; Pan, Y.; Mi, S.; Sun, D.; Wang, C. Genome-wide association study for reproductive traits in a Large White pig population. Anim. Genet. 2018, 49, 127–131. [Google Scholar] [CrossRef] [Green Version]
  3. Zhu, Y.; Li, W.; Yang, B.; Zhang, Z.; Ai, H.; Ren, J.; Huang, L. Signatures of Selection and Interspecies Introgression in the Genome of Chinese Domestic Pigs. Genome Biol. Evol. 2017, 9, 2592–2603. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Chen, M.; Wang, J.; Wang, Y.; Wu, Y.; Fu, J.; Liu, J.-F. Genome-wide detection of selection signatures in Chinese indigenous Laiwu pigs revealed candidate genes regulating fat deposition in muscle. BMC Genet. 2018, 19, 31. [Google Scholar] [CrossRef] [PubMed]
  5. Jeong, H.; Song, K.-D.; Seo, M.; Caetano-Anollés, K.; Kim, J.; Kwak, W.; Oh, J.-D.; Kim, E.; Jeong, D.K.; Cho, S.; et al. Exploring evidence of positive selection reveals genetic basis of meat quality traits in Berkshire pigs through whole genome sequencing. BMC Genet. 2015, 16, 104. [Google Scholar] [CrossRef] [Green Version]
  6. Szmatoła, T.; Jasielczuk, I.; Semik-Gurgul, E.; Szyndler-Nędza, M.; Blicharski, T.; Szulc, K.; Skrzypczak, E.; Gurgul, A. Detection of runs of homozygosity in conserved and commercial pig breeds in Poland. J. Anim. Breed. Genet. 2020, 137, 571–580. [Google Scholar] [CrossRef] [PubMed]
  7. Zhang, Z.; Zhang, Q.; Xiao, Q.; Sun, H.; Gao, H.; Yang, Y.; Chen, J.; Li, Z.; Xue, M.; Ma, P.; et al. Distribution of runs of homozygosity in Chinese and Western pig breeds evaluated by reduced-representation sequencing data. Anim. Genet. 2018, 49, 579–591. [Google Scholar] [CrossRef]
  8. Chen, L.; Tian, S.; Jin, L.; Guo, Z.; Zhu, D.; Jing, L.; Che, T.; Tang, Q.; Chen, S.; Zhang, L.; et al. Genome-wide analysis reveals selection for Chinese Rongchang pigs. Front. Agric. Sci. Eng. 2017, 4, 319. [Google Scholar] [CrossRef] [Green Version]
  9. Wilkinson, S.; Lu, Z.H.; Megens, H.-J.; Archibald, A.; Haley, C.; Jackson, I.; Groenen, M.; Crooijmans, R.P.M.A.; Ogden, R.; Wiener, P. Signatures of Diversifying Selection in European Pig Breeds. PLoS Genet. 2013, 9, e1003453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Li, X.; Yang, S.; Dong, K.; Tang, Z.; Li, K.; Fan, B.; Wang, Z.; Liu, B. Identification of positive selection signatures in pigs by comparing linkage disequilibrium variances. Anim. Genet. 2017, 48, 600–605. [Google Scholar] [CrossRef]
  11. Moravčíková, N.; Kasarda, R.; Vostrý, L.; Krupová, Z.; Krupa, E.; Lehocká, K.; Olšanská, B.; Trakovická, A.; Nádaský, R.; Židek, R.; et al. Analysis of selection signatures in the beef cattle genome. Czech J. Anim. Sci. 2019, 64, 491–503. [Google Scholar] [CrossRef] [Green Version]
  12. Nielsen, R.; Williamson, S.; Kim, Y.; Hubisz, M.J.; Clark, A.G.; Bustamante, C. Genomic scans for selective sweeps using SNP data. Genome Res. 2005, 15, 1566–1575. [Google Scholar] [CrossRef] [Green Version]
  13. Weigand, H.; Leese, F. Detecting signatures of positive selection in non-model species using genomic data. Zool. J. Linn. Soc. 2018, 184, 528–583. [Google Scholar] [CrossRef]
  14. Orlov, Y.L.; Anashkina, A.A. Life: Computational Genomics Applications in Life Sciences. Life 2021, 11, 1211. [Google Scholar] [CrossRef]
  15. Wang, K.; Wu, P.-X.; Chen, D.-J.; Zhou, J.; Yang, X.-D.; Jiang, A.-A.; Ma, J.-D.; Tang, Q.-Z.; Xiao, W.-H.; Jiang, Y.-Z.; et al. Genome-wide scan for selection signatures based on whole-genome re-sequencing in Landrace and Yorkshire pigs. J. Integr. Agric. 2021, 20, 1898–1906. [Google Scholar] [CrossRef]
  16. Porto-Neto, L.R.; Sonstegard, T.S.; Liu, G.E.; Bickhart, D.M.; Da Silva, M.V.; Machado, M.A.; Utsunomiya, Y.T.; Garcia, J.F.; Gondro, C.; Van Tassell, C.P. Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping. BMC Genom. 2013, 14, 876. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Nicholson, G.; Smith, A.V.; Jónsson, F.; Gústafsson, Ó.; Stefánsson, K.; Donnelly, P. Assessing population differentiation and isolation from single-nucleotide polymorphism data. J. R. Stat. Soc. Ser. B (Stat. Methodol.) 2002, 647, 695–715. [Google Scholar] [CrossRef]
  18. Bakoev, S.; Kolosov, A.; Bakoev, F.; Kostyunina, O.; Bakoev, N.; Romanets, T.; Koshkina, O.; Getmantseva, L. Analysis of Homozygous-by-Descent (HBD) Segments for Purebred and Crossbred Pigs in Russia. Life 2021, 11, 861. [Google Scholar] [CrossRef] [PubMed]
  19. Бакoев, С.; Гетманцева, Л. Метoды oценки инбридинга и пoдписей селекции сельскoхoзяйственных живoтных на oснoве прoтяженных гoмoзигoтных oбластей. Дoстижения науки и техники АПК 2019, 33, 11. [Google Scholar]
  20. Rauw, W.M.; Rydhmer, L.; Kyriazakis, I.; Øverland, M.; Gilbert, H.; Dekkers, J.C.; Hermesch, S.; Bouquet, A.; Izquierdo, E.G.; Louveau, I.; et al. Prospects for sustainability of pig production in relation to climate change and novel feed resources. J. Sci. Food Agric. 2020, 100, 3575–3586. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Derks, M.F.L.; Megens, H.-J.; Bosse, M.; Lopes, M.S.; Harlizius, B.; Groenen, M.A.M. A systematic survey to identify lethal recessive variation in highly managed pig populations. BMC Genom. 2017, 18, 858. [Google Scholar] [CrossRef] [Green Version]
  22. Bakoev, S.; Getmantseva, L.; Kostyunina, O.; Bakoev, N.; Prytkov, Y.; Usatov, A.; Tatarinova, T.V. Genome-wide analysis of genetic diversity and artificial selection in Large White pigs in Russia. PeerJ 2021, 9, e11595. [Google Scholar] [CrossRef] [PubMed]
  23. Kunhareang, S.; Zhou, H.; Hickford, J. Rapid DNA extraction of pig ear tissues. Meat Sci. 2010, 85, 589–590. [Google Scholar] [CrossRef]
  24. Getmantseva, L.; Bakoev, S.; Bakoev, N.; Karpushkina, T.; Kostyunina, O. Mitochondrial DNA Diversity in Large White Pigs in Russia. Animals 2020, 10, 1365. [Google Scholar] [CrossRef]
  25. Muller, H.-G. Smooth optimum kernel estimators of densities, regression curves and modes. Ann. Stat. 1984, 12, 766–774. [Google Scholar] [CrossRef]
  26. Druet, T.; Gautier, M. A model-based approach to characterize individual inbreeding at both global and local genomic scales. Mol. Ecol. 2017, 26, 5820–5841. [Google Scholar] [CrossRef] [Green Version]
  27. Bertrand, A.R.; Kadri, N.K.; Flori, L.; Gautier, M.; Druet, T. RZooRoH: An R package to characterize individual genomic autozygosity and identify homozygous-by-descent segments. Methods Ecol. Evol. 2019, 10, 860–866. [Google Scholar] [CrossRef]
  28. Fonseca, P.A.S.; Suárez-Vega, A.; Marras, G.; Cánovas, Á. GALLO: An R package for genomic annotation and integration of multiple data sources in livestock for positional candidate loci. GigaScience 2020, 9, giaa149. [Google Scholar] [CrossRef] [PubMed]
  29. Amills, M.; Clop, A.; Ramírez, O.; Pérez-Enciso, M. Origin and Genetic Diversity of Pig Breeds. eLS 2010. [Google Scholar] [CrossRef] [Green Version]
  30. Purfield, D.C.; McParland, S.; Wall, E.; Berry, D. The distribution of runs of homozygosity and selection signatures in six commercial meat sheep breeds. PLoS ONE 2017, 12, e0176780. [Google Scholar] [CrossRef] [Green Version]
  31. Becker, S.F.; Langhe, R.; Huang, C.; Wedlich, D.; Kashef, J. Giving the right tug for migration: Cadherins in tissue movements. Arch. Biochem. Biophys. 2012, 524, 30–42. [Google Scholar] [CrossRef] [PubMed]
  32. Ding, X.; Zhang, Y.; Huang, Y.; Liu, S.; Lu, H.; Sun, T. Cadherin-11 involves in synovitis and increases the migratory and invasive capacity of fibroblast-like synoviocytes of osteoarthritis. Int. Immunopharmacol. 2015, 26, 153–161. [Google Scholar] [CrossRef] [PubMed]
  33. Langhe, R.P.; Gudzenko, T.; Bachmann, M.; Becker, S.F.; Gonnermann, C.; Winter, C.; Abbruzzese, G.; Alfandari, D.; Kratzer, M.-C.; Franz, C.M.; et al. Cadherin-11 localizes to focal adhesions and promotes cell–substrate adhesion. Nat. Commun. 2016, 7, 10909. [Google Scholar] [CrossRef] [Green Version]
  34. Getmantseva, L.; Kolosova, M.; Bakoev, F.; Zimina, A.; Bakoev, S. Genomic Regions and Candidate Genes Linked to Capped Hock in Pig. Life 2021, 11, 510. [Google Scholar] [CrossRef]
  35. Le, T.H.; Christensen, O.F.; Nielsen, B.; Sahana, G. Genome-wide association study for conformation traits in three Danish pig breeds. Genet. Sel. Evol. 2017, 49, 12. [Google Scholar] [CrossRef] [PubMed]
  36. Teixeira, D.L.; Salazar, L.C.; Enriquez-Hidalgo, D.; Boyle, L.A. Assessment of Animal-Based Pig Welfare Outcomes on Farm and at the Abattoir: A Case Study. Front. Veter. Sci. 2020, 7, 576942. [Google Scholar] [CrossRef]
  37. Wang, L.; Guo, Y.; Huang, W.-J.; Ke, X.; Poyet, J.-L.; Manji, G.A.; Merriam, S.; Glucksmann, M.; DiStefano, P.S.; Alnemri, E.S.; et al. CARD10 Is a Novel Caspase Recruitment Domain/Membrane-associated Guanylate Kinase Family Member That Interacts with BCL10 and Activates NF-κB. J. Biol. Chem. 2001, 276, 21405–21409. [Google Scholar] [CrossRef] [Green Version]
Figure 1. SVD (A) and Heatmap (B) for pigs LW_A and LW_B. (LW_A—Large White Russian selection, LW-B—a commercial Large White).
Figure 1. SVD (A) and Heatmap (B) for pigs LW_A and LW_B. (LW_A—Large White Russian selection, LW-B—a commercial Large White).
Genes 13 00197 g001
Figure 2. QTLs in genomic regions with strong outliers according to smoothing FST. ((A)—percent QTL type. (B)—QTL enrichment analysis (the more intensive the red shade, the more significant is enrichment; the area of circles is proportional to quantity QTL; richness factor—the attitude of QTL quantity, annotated in the study areas, to the total number of each QTL in the reference database).
Figure 2. QTLs in genomic regions with strong outliers according to smoothing FST. ((A)—percent QTL type. (B)—QTL enrichment analysis (the more intensive the red shade, the more significant is enrichment; the area of circles is proportional to quantity QTL; richness factor—the attitude of QTL quantity, annotated in the study areas, to the total number of each QTL in the reference database).
Genes 13 00197 g002
Figure 3. Implemented autozygosity in pigs. ((A)—distribution of the values of the inbreeding coefficient for LW_A; (B)—distribution of the values of the inbreeding coefficient for LW_B).
Figure 3. Implemented autozygosity in pigs. ((A)—distribution of the values of the inbreeding coefficient for LW_A; (B)—distribution of the values of the inbreeding coefficient for LW_B).
Genes 13 00197 g003
Figure 4. The proportion of autozygosity in different HBD classes. ((A)—the proportion of autozygosity in different classes of HBD in LW_A; (B)—the proportion of autozygosity in different classes of HBD in LW_B; (C)—variations in individual levels of autozygosity in LW_A; (D)—variations in individual levels of autozygosity in LW_B).
Figure 4. The proportion of autozygosity in different HBD classes. ((A)—the proportion of autozygosity in different classes of HBD in LW_A; (B)—the proportion of autozygosity in different classes of HBD in LW_B; (C)—variations in individual levels of autozygosity in LW_A; (D)—variations in individual levels of autozygosity in LW_B).
Genes 13 00197 g004
Figure 5. The average length of HDB segments on chromosomes in different classes in pigs.
Figure 5. The average length of HDB segments on chromosomes in different classes in pigs.
Genes 13 00197 g005
Figure 6. Manhattan Plot-% SNP in HBD for LW_A.
Figure 6. Manhattan Plot-% SNP in HBD for LW_A.
Genes 13 00197 g006
Figure 7. Manhattan Plot-% SNP in HBD for LW_B.
Figure 7. Manhattan Plot-% SNP in HBD for LW_B.
Genes 13 00197 g007
Figure 8. QTLs type in HBD areas.
Figure 8. QTLs type in HBD areas.
Genes 13 00197 g008
Table 1. Top HBD pig LW_A and LW_B.
Table 1. Top HBD pig LW_A and LW_B.
GroupChromnSNPFromToGene
LW_A1136370473364301394FHL5, GPR63, NDUFAF4, KLHL32
LW_A114116300444117138985PYGO1, DNAAF4, C15orf65, CCPG1, PIGB, RAB27, RSL24D1
LW_A117149261646150160033ZNF407, CNDP1, CNDP2, DIPK1C, C18orf63, CYB5A, FBX015, TIMM21
LW_A146154796451158940596DSEL, CDH19, CDH7, SERPINB10, SERPINB7, SERPINB12, SERPINB8, SERPINB5, SERPINB2, SERPINB13, SERPINB11, VDS4B, PHLPP1, KDSR
LW_B12454194039153513159MAN1A1, FAM184A, ASF1A, CEP85L, SLC35F1, NUS1, DCBLD1, ROS1, VGLL2, RFX6, GPRC6A, FAM162B, KPNA5, ZUP1, RSPH4A, PTP4A1, PHF3, ADGRB3, LMBRD1, COL19A1, COL9A1, FAM135A, SMAP1, B3GAT2, OGFRL1, RIMS1, KCNQ5, DPPA5, OOEP, CYB5R4, MRAP2, CEP162
LW_B1105656395657069421RNGTT
LW_B1116198623262371709-
LW_B114225017887225614369CEMIP2, ABHD17B, C9orf85, GDA
LW_B118227781084229301982TRPM6, C9orf40, NMRK1, CARNMT1, OSTF1, PCSK5
LW_B6206126167163135016PEG3, AURKC, ZNF304, ZNF772, ZNF773, ZNF550, ZNF606, ZNF135, ZNF329, ZNF274, ZNF8, RPS5, ZNF584, ZNF446, SLC27A5, ZBTB45, TRIM28, CHMP2A, UBE2M, MZF1
LW_B684107551091113309905RBBP8, CABLES1, TMEM241, RIOK3, RMC1, NPC1, ANKRD29, LAMA3, TTC39C, CABYR, OSBPL1A, HRH4, ZNF521, SS18, PSMA8, TAF4B, KCTD1, AQP4, CHST9, CDH2
LW_B10142837398128880540-
LW_B1419104827372105533891MYOF, CEP55, FFAR4, RBP4, PDE6C, FRA10AC1, LGI1, SLC35G1, PLCE1
LW_B15222606975426699994-
Table 2. Intersection of the smoothing FST and HBD areas.
Table 2. Intersection of the smoothing FST and HBD areas.
LW_ALW_BGene
1SSC1: 239389749–2433677270.380.51NCBP1, PLPPR1, GRIN3A
2SSC1: 1559862860.730.31-
3SSC4: 79218630.080.07-
4SSC4:40772747–414349850.270.45-
5SSC4: 1193385110.230.31lncRNA
6SSC5: 60515812–670770930.190.48ETV6, VWF, CCND2, TSPAN9
7SSC5: 102941380.120.31CARD10
8SSC6: 51546940.120.38CDH13
9SSC7: 103380850–1156664930.240.07CEP128, SERPINA11
10SSC9: 106642564–1197574340.080.52PIK3CG, COG5, BCAP29, SLC26A4
11SSC10: 49802418–502773300.230.31-
12SSC11: 113356140.230.17NBEA
SSC11: 13345472–135274220.230.31TRPC4
SSC11: 65881573–698728680.300.66HS6ST3, NALCN
13SSC12: 481810490.120.38SMG6
14SSC13: 200778735–2012458270.040.48TTC3, KCNJ6
15SSC14: 138917940.120.21-
16SSC15: 115593713–1215615030.230.38IKZF2, OBSL1
17SSC17: 81626130.190.31-
18SSC18: 458019970.080.24-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kolosov, A.; Getmantseva, L.; Kolosova, M.; Romanets, T.; Bakoev, N.; Romanets, E.; Bakoeva, I.; Kostyunina, O.; Prytkov, Y.; Tretiakova, O.; et al. Investigation of the Genetic Architecture of Pigs Subjected to Breeding Intensification. Genes 2022, 13, 197. https://doi.org/10.3390/genes13020197

AMA Style

Kolosov A, Getmantseva L, Kolosova M, Romanets T, Bakoev N, Romanets E, Bakoeva I, Kostyunina O, Prytkov Y, Tretiakova O, et al. Investigation of the Genetic Architecture of Pigs Subjected to Breeding Intensification. Genes. 2022; 13(2):197. https://doi.org/10.3390/genes13020197

Chicago/Turabian Style

Kolosov, Anatoly, Lyubov Getmantseva, Maria Kolosova, Timofey Romanets, Nekruz Bakoev, Elena Romanets, Ilona Bakoeva, Olga Kostyunina, Yuri Prytkov, Olga Tretiakova, and et al. 2022. "Investigation of the Genetic Architecture of Pigs Subjected to Breeding Intensification" Genes 13, no. 2: 197. https://doi.org/10.3390/genes13020197

APA Style

Kolosov, A., Getmantseva, L., Kolosova, M., Romanets, T., Bakoev, N., Romanets, E., Bakoeva, I., Kostyunina, O., Prytkov, Y., Tretiakova, O., & Bakoev, S. (2022). Investigation of the Genetic Architecture of Pigs Subjected to Breeding Intensification. Genes, 13(2), 197. https://doi.org/10.3390/genes13020197

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop