Skip to main content

Population dynamics of potentially harmful haplotypes: a pedigree analysis

Abstract

Background

The identification of low-frequency haplotypes, never observed in homozygous state in a population, is considered informative on the presence of potentially harmful alleles (candidate alleles), putatively involved in inbreeding depression. Although identification of candidate alleles is challenging, studies analyzing the dynamics of potentially harmful alleles are lacking. A pedigree of the highly endangered Gochu Asturcelta pig breed, including 471 individuals belonging to 51 different families with at least 5 offspring each, was genotyped using the Axiom PigHDv1 Array (658,692 SNPs). Analyses were carried out on four different cohorts defined according to pedigree depth and at the whole population (WP) level.

Results

The 4,470 Linkage Blocks (LB) identified in the Base Population (10 individuals), gathered a total of 16,981 alleles in the WP. Up to 5,466 (32%) haplotypes were statistically considered candidate alleles, 3,995 of them (73%) having one copy only. The number of alleles and candidate alleles varied across cohorts according to sample size. Up to 4,610 of the alleles identified in the WP (27% of the total) were present in one cohort only. Parentage analysis identified a total of 67,742 parent-offspring incompatibilities. The number of mismatches varied according to family size. Parent-offspring inconsistencies were identified in 98.2% of the candidate alleles and 100% of the LB in which they were located. Segregation analyses informed that most potential candidate alleles appeared de novo in the pedigree. Only 17 candidate alleles were identified in the boar, sow, and paternal and maternal grandparents and were considered segregants.

Conclusions

Our results suggest that neither mutation nor recombination are the major forces causing the apparition of candidate alleles. Their occurrence is more likely caused by Allele-Drop-In events due to SNP calling errors. New alleles appear when wrongly called SNPs are used to construct haplotypes. The presence of candidate alleles in either parents or grandparents of the carrier individuals does not ensure that they are true alleles. Minimum Allele Frequency thresholds may remove informative alleles. Only fully segregant candidate alleles should be considered potentially harmful alleles. A set of 16 candidate genes, potentially involved in inbreeding depression, is described.

Peer Review reports

Background

Linkage disequilibrium (LD) may shape the genome causing its organization in block-like structures [1, 2], varying in size, usually called Linkage disequilibrium Blocks (LB). LB were defined as strong LD stretches flanked by genomic areas with higher recombination rates [3]. Although there exist different approaches to identify LB in non-related individuals [4, 5], LB show large variation both among populations and between individuals within population. LD patterns may vary due to genetic factors such as mutation and recombination but also to population events including selection, demographic changes, or drift [6, 7]. This is why, for practical applications, the ascertainment of patterns of haplotype conservation within LB is of primary interest.

Although LB are assumed to provide insights for disease association mapping and population genetics, there is an increasing interest in using them to explain the apparition of inbreeding depression in highly selected livestock populations [8, 9]. Since deleterious alleles are assumed to be fatal, low-frequency haplotypes, that are never observed in homozygous state in a population, are considered potentially lethal recessive genetic variants [8, 10,11,12]. However, the identification of low-frequency haplotypes requires a large number of genotyped individuals to ensure that they do not appear in homozygous state by chance [13, 14].

Detection of such genetic variants in large populations is difficult and subject to the occurrence of both false positives and false negatives [9]. Therefore, identification of potentially lethal haplotypes is usually carried out in two steps: first, their frequencies and their probabilities of being observed as homozygotes in a population are computed; and second, further restrictions such as their presence in one or both parents or grand-parents are applied [8, 9, 11]. The use of pedigree information may be even more important in populations with small to moderate sizes because lethal recessive variants have the potential to drift rapidly to higher frequencies due to inbreeding [15].

Although the potential role of mutation, demography, and selection on the frequency of deleterious alleles has been discussed before [7], little is known about the dynamics of potentially lethal haplotypes. Assuming that genetic variants causing inherited defects should expand in a population from founder events involving carriers [16], this research aims at the ascertainment of segregation patterns of haplotypes with low probabilities of being observed at a homozygous state using a small pedigree of the highly inbred Gochu Asturcelta pig breed.

The breeding history of Gochu Asturcelta pig can make the available pedigree informative for uncovering the dynamics of potentially lethal alleles. The highly endangered Gochu Asturcelta pig breed derives from four founders only [17, 18]. Wrong management practices, including full-sib matings, caused a sudden increase in inbreeding early after the start of the conservation programme [19] and the occurrence of inbreeding-related events such as stillborn parturitions that have been reported before in other inbred pig populations [20,21,22]. This scenario advised the implementation of a strict minimum coancestry mating policy [17] allowing, at present, a normal production and reproduction performance and, theoretically, the maintenance of lethal alleles in low frequency in the population and even purging.

In this research, carriers of potentially harmful haplotypes will be identified and the dynamics of such haplotypes, including Mendelian inheritance and parents-offspring inconsistencies, analyzed. Causes of deviation from the rules of Mendelian inheritance and the occurrence of alleles that are additional to the parental genotypes will be discussed. Segregating LB, carrying potentially lethal recessive haplotypes, will be subject to gene enrichment analyses to contribute to the knowledge of genomic areas potentially involved in inbreeding depression in pig.

Methods

Pedigree, cohorts and genotyping

A Gochu Asturcelta pig pedigree including 471 individuals (fully described in Supplementary Table S1), belonging to 51 different families (descendants of the same parental couple), was analyzed. This pedigree is a sample of that previously analyzed in Arias et al. (including 534 individuals, and 76 families) which was edited removing families with less than five offspring typed [23]. However, to gain consistency, pedigrees of two small families giving two boars (168 and 658) with good reproduction success were kept. Individuals were genotyped using the Axiom Porcine Genotyping Array (Axiom_PigHDv1; 658,692 SNPs). The typed individuals derived from 17 genotyped boars and 35 genotyped sows. Family size (number of offspring per parental couple) varied from 5 to 34, with 9 offspring on average. When necessary, the pedigree was visualized using the library visPedigree [24] of R environment.

For descriptive purposes, the individuals in the pedigree (whole population; WP) were split according to pedigree depth (equivalents to complete generations; t; [25]):

  • Base population (BP) formed by 10 different individuals including: (i) two founders of the Gochu Asturcelta pig breed; (ii) four descendants of two untyped founders of the breed; and (iii) four direct descendants of the typed founders.

  • G12 cohort, formed by 52 individuals with t ≤ 2 equivalents to complete generations and not included in the BP.

  • G23 cohort, formed by 281 individuals with pedigree depth varying from t > 2 to t ≤ 3.

  • G3 cohort, formed by 128 individuals with t > 3 in their pedigree.

The software Axiom Analysis Suite v4.0.3 (Thermo Fisher Scientific, Waltham, MA) was used to create standard .ped, .map, and .vcf files. SNPs were mapped using the Sscrofa genome build 11.1 [26]. Only autosomal chromosomes with known positions were considered. Following previous approaches, genotypes were only filtered on Mendelian errors to avoid the presence of null and false alleles [23, 27]. A total of 503,043 SNPs (372,402 of them polymorphic) with a minimum call rate of 0.97 were retained for further analysis. Missing genotypes were imputed using BEAGLE v5.4 [28, 29] using default parameters.

Identification of LB and candidate haplotypes

Haplotype blocks (LB), defined following Gabriel et al. [3], were identified at a base population level using PLINK v1.9 [30] following the default procedure implemented in HAPLOVIEW v4.1 [1]. Since Veroneze et al. reported that the average extent of LD in pig was about 400 kb [31], the option --blocks-max-kb was fitted to 400. Considering that our dataset was previously filtered for Mendelian errors, no Minimum Allele Frequency thresholds were applied for the LB identification. Alleles (arbitrarily coded using three digits) and genotypes were identified per each LB, using our own R code and allelic frequency, q, computed.

Following Howard et al., potentially lethal recessive haplotypes were identified by computing haplotype frequency at both the population-wide and the within-cohort level [8]. The expected number of homozygotes (E[H]) will be calculated as E[H] = q2N, where q is the frequency of the haplotype and N is the total number of haplotyped individuals in each cohort (10, 52, 281, and 128, respectively) and the whole population (471). Assuming Poisson distribution, the probability of observing no homozygous haplotypes (O[H]) given E[H] as P(O[H] = 0|E[H]) = e− E[H].

Haplotypes with P(O[H] = 0|E[H]) < 0.05 were considered candidate harmful haplotypes. Candidate haplotypes were identified at both the whole population level and the cohort level.

When necessary, relationships among candidate haplotypes were summarized with Venn diagrams constructed using the web-based software InteractiVenn [32].

Parents-offspring inconsistencies

Haplotypes within LB were arbitrarily coded as codominant (microsatellite-like) marker alleles at the whole population level. The program COLONY 2.0.7.0 [33] was used to identify parent-offspring inconsistencies (i.e. departures of Mendelian inheritance). Analyses would inform on the consistency of the inheritance patterns across generations and families. The inconsistent haplotypes were identified using the genotypes inferred by the full-likelihood method [34] implemented in the COLONY and further assigned to the Offspring, the Father, or the Mother for each parent-offspring trio. The full-likelihood method [34] implemented in the COLONY pedigree uses the available pedigree, including family structures (known and excluded paternal and maternal sibships), and genotypes of related individuals. COLONY was run using the following settings: (a) polygamy with inbreeding as mating system; and (b) medium run lengths to search for the best assignment in the simulated annealing algorithm. Since parentage was previously verified using COLONY and SNP array data, the option Known paternity was fitted for all individuals in dataset.

Following Arias et al., parent-offspring incompatibility rates were quantified as: (a) Mean error rate per locus (el) as el = ml/nt ; and (b) Mean error rate per allele (ea) computed as ea = ma/2nt, were ml is the number of single-locus genotypes including at least one allelic error, ma, the number of allelic mismatches, and nt, the number of replicated single-locus genotypes [27].

Complementary analyses

Recombination rate (ρ) was estimated on phased genotypes using a machine-learning approach implemented in the R package FastEPRR [35]. Estimates were carried out on each porcine autosome with overlapping sliding windows under default settings for diploid organisms and window size fitted to 50 kb and step length fitted to 25 kb. To ascertain the relationships between LB and ρ, LB and ρ windows were overlapped using the intersectBed function of the BedTools software [36].

Furthermore, since genetic features such as Copy Number Variations (CNV) have been shown to cause calling errors giving wrong SNP alleles [27], the LB identified were overlapped, using BedTools as well, with 344 CNV regions previously identified in the Gochu Asturcelta pig population (see Supplementary Table S5 of Arias et al. [37]).

Enrichment and functional annotation analyses

Following previous analyses [38], candidate haplotypes were subject to enrichment analyses using the BioMart tool [39]. Protein-coding genes found within these candidate haplotypes identified were retrieved from the Ensembl Genes 91 database, based on the Sscrofa v11.1 porcine reference genome. All the genes identified in the genomic areas spanned by candidate haplotypes were processed using the functional annotation tool implemented in DAVID Bioinformatics resources 6.8 [40] to determine enriched functional terms. An enrichment score of 1.3, which is equivalent to the Fisher exact test P-value of 0.05 [40], was used as a threshold to define the significantly enriched functional terms in comparison to the whole porcine reference genome background. Selection of significant composite annotation terms (clusters) using enrichment score as a criterion for selection rather than single annotation terms as independent statistically significant entities supports the identification of biological functions. In other words, this strategy allows to consider relationships between Gene Ontology annotation terms, by moving the analysis of biological function from the level of single genes to that of biological processes [40, 41].

Results

General overview

Up to 4,470 LB were identified in the BP, gathering a total of 16,981 alleles in the whole typed population (3.8 ± 2.9 alleles per LB, on average) and covering 236.2 Mb of the porcine genome (52.8 ± 93.9 kb on average; Table 1; Supplementary Table S2). Most LB (2,003; 45%) had two alleles only, 1,475 LB (33%) had more than four alleles, 232 (5%) had ten alleles or more, and 16 (0.4%) had 20 alleles or more. Furthermore, the number of alleles identified in each of the cohorts defined showed a large variation (Fig. 1). In the BP, the 4,470 LB had a total of 11,900 alleles. In the other cohorts the number of alleles increased with the number of individuals typed: 12,441 different alleles were identified in the 52 individuals forming cohort G12; 14,639 alleles were identified in the 281 individuals assigned to cohort G23; and 12,966 alleles were identified in the 128 individuals belonging to the cohort G3. The observed homozygosity in the whole population was 0.464.

Table 1 Number of Linkage Blocks (LB), mean length (in kb) and mean number of alleles per LB identified in the pedigree of Gochu Asturcelta analyzed per porcine chromosome. Minimum and maximum values are in brackets. Furthermore, mean recombination rate (ρ ± s.d.; in cM/Mb) and maximum Rho value (in brackets) are also given
Fig. 1
figure 1

Allelic frequencies per cohort defined according to pedigree depth (BP, G12, G23 and G3), and in the whole typed population. Plot illustrates the total number of alleles identified (in blue), the number of candidate alleles (in orange), and the number of alleles that are identified in one cohort only (in grey). The number of individuals typed in each of the cohorts and in the whole population are in brackets

The number of candidate alleles tended to vary across cohorts with the number of individuals typed as well (Fig. 1). Up to 1,516 (on 1,075 LB), 958 (on 694 LB), 3,230 (on 1,524 LB), and 2,193 (on 1,259 LB) alleles were considered candidate in the BP and cohorts G12, G23 and G3, respectively. Finally, 5,466 candidate alleles (32% of the total; gathering 8,587 copies) identified in the whole typed population, on 2,031 different LB (45% of the total), were considered candidates. Most candidate alleles identified at the whole population level (3,995; 73% of the total) had one copy only, 794 (15%) had two copies, 445 (8%) had three or four copies, and 232 (4%) from five (63) to nine (29) copies. Only three candidate alleles identified in the WP (two of them identified in the G23 cohort as well) had 1 homozygous genotype each (Supplementary Table S3).

Up to 4,610 of the alleles identified at the whole population level (27% of the total) were present in one cohort only (Fig. 1): 260 alleles were identified in the BP only and 528, 2,416, and 1,406 alleles were only identified within cohorts G12, G23, and G3, respectively. Of them, 181, 525, 2,410, and 1,399 were considered candidates within each cohort. Notably, six out of these 4,610 alleles (104 on LB196, 111 on LB645, 109 on LB1656, 112 on LB2925, 106 on LB3194, and 109 on LB3315) were not considered candidates at the WP level. Figure 2 illustrates the pedigrees of the carriers of these six non-candidate alleles. In five cases segregation involved one parent (sow 343 for alleles 104 on LB196 and 109 on LB1656; boar 486 for alleles 112 on LB2925 and 109 LB3315; and sow 332 for allele 111 on LB645) in which the alleles appear de novo and is transmitted to the offspring (from 9 to 11 descendants). Furthermore, one of these candidate alleles (allele 106 on LB3194; Fig. 2F) was identified de novo in 11 offspring of sow 167 (mated with two different boars) in which that allele was not present.

Fig. 2
figure 2

Pedigrees showing the segregation of six non-candidate alleles in the Gochu Asturcelta population analyzed. Boars are in squares and sows are in circles. Black squares and circles identify the carrier individuals. Plots A and B show the pedigree of sow 343 for alleles 104 on LB196 and 109 on LB1656, respectively. Plots C and D show the pedigree of boar 486 for alleles 112 on LB2925 and 109 on LB3315, respectively. Plot E shows the pedigree of sow 332 for allele 111 on LB645; and Plot F shows the pedigree of sow 167 (mated with two different boars) for allele 106 on LB3194. Only carrier offspring are included

Parent-offspring inconsistencies

A total of 67,742 parent-offspring incompatibilities (ea = 0.016), involving 9,624 different alleles, were identified on 3,405 LB (el = 0.761) across the 471 individuals of the analyzed pedigree (Supplementary Table S4). Up to 5,229 parent-offspring incompatibilities involved the two alleles at a locus. A total of 63,016 parent-offspring incompatibilities were assigned to the offspring, 2,556 to the father, and 2,170 to the mother. The 3,405 LB on which paternal inconsistencies were identified were, on average, bigger (64.9 ± 104.4 kb vs. 14.1 ± 13.9 kb), gathered a higher number of alleles (4.3 ± 3.1 vs. 2.0 ± 0.3), and were less homozygous (observed homozygosity of 0.226 vs. 0.250) than the remnant 1,065 LB (Supplementary Table S2).

Parents-offspring allelic incompatibilities varied across both cohorts and families according to sample size. In this respect, the BP gathered 1,009 errors (1.5% of the total), and cohorts G12, G23, and G3 gathered 6,609 (9.8%), 52,675 (77.8%) and 7,450 (11.0%) errors, respectively (Supplementary Table S4). Figure 3 illustrates how both total and mean number of mismatches per family varied according to family size. The two more-sized families (family 24 with 31 offspring and family 12 with 34 offspring) gathered a total of 11,326 and 14,742 mismatches, respectively.

Fig. 3
figure 3

Dispersion plots constructed according to family size (on the Y-axis) and the total number of parents-offspring inconsistencies per family (on the X-axis; Plot A) and the mean number of parents-offspring inconsistencies per family (on the X-axis; Plot B)

Figure 4 illustrates the relationships between the 5,466 candidate alleles identified in the WP and their corresponding LB, and those in which parent-offspring inconsistencies were identified. Parent-offspring inconsistencies were identified in 98.2% of the candidate alleles and 100% of the corresponding LB.

Fig. 4
figure 4

Venn diagrams summarizing the relationships between the 5,466 candidate alleles identified in the whole population and the 3,405 alleles involved in parents-offspring incompatibilities (Plot A). The relationships between the 3,405 LB on which parent-offspring allelic incompatibilities were identified and the 2,031 LB carrying candidate (potentially harmful) haplotypes at the whole population level are also summarized (Plot B)

Segregation patterns of candidate alleles

The segregation patterns of the candidate alleles showed considerable variation. Figure 5 contrasts a fully segregant candidate allele which is transmitted from the founder boar 20 to the daughter 67 and five grandsons (Plot A) with three pedigrees in which the identification of candidate alleles could be considered appearance de novo except for its identification, in some cases in grandparents or great-grandparents.

Fig. 5
figure 5

Pedigrees showing the segregation of four candidate alleles in the Gochu Asturcelta population analyzed. Boars are in squares and sows are in circles. Black squares and circles identify the carrier individuals. Plot A shows the fully segregant allele 101 on LB3529 which is transmitted from the founder boar 20 to the daughter 67 and five grandsons. Plots B (allele 102 on LB746), C (allele 105 on LB2481), and D (allele 103 on LB706) show pedigrees in which candidate alleles could be considered as identified de novo except for their identification in grandparents or great-grandparents

Figure 6 summarizes the frequency in which candidate alleles in the WP were identified in the fathers, the mothers, the paternal grandfathers, or the maternal grandfathers of the carrier individuals (see Supplementary Tables S2 and S3 as well). Up to 5,169 candidate alleles identified in the WP (94.6% of the total), gathering 6,774 copies (79% of the total copies of the candidate alleles), could not be identified in any ancestor of the carrier individuals.

Fig. 6
figure 6

Venn diagram illustrating the frequency in which the 5,466 candidate alleles identified in the whole population could be identified in the fathers, the mothers, the paternal grandfathers, or the maternal grandfathers of the carrier

A total of 188 alleles (3.4%), gathering 834 copies (1% of the total), were present in either the Father or the Mother of the carrier individuals but not in other ancestors. Moreover, 43 alleles (0.8%), gathering 138 copies (0.2%), were identified in at least one maternal or paternal Grand Parent of the carrier individuals but not in the parents (either sow or boar) of the carrier individuals.

Finally, a total of 66 candidate alleles (1.2%), gathering 470 copies (5.5%), were jointly identified in at least one parent (either sow or boar) and at least one maternal or paternal Grand Parent of the carrier individuals. Only 17 candidate alleles (0.3%), gathering 111 copies (1.3%), located on SSC4, SSC6, and SSC13, were identified in the boar, sow, and paternal and maternal grandparents and were, therefore, classified as segregating (Table 2). The LB on which these 17 segregating candidate alleles were identified lengthened 102.8 kb on average (min. 8.9 kb; max. 386.3 kb) and had 5.9 alleles on average (min. 3; max. 12). Figure 7 illustrates the variation of these 17 segregating candidate haplotypes. Interestingly, all of them were carried by founder individuals and were transmitted to progeny with pedigree depth up to 2.5 equivalents to complete generations. These 17 candidate alleles were selected for enrichment analyses.

Table 2 List of the 17 Linkage Blocks (LB) identified in the Gochu Asturcelta pig pedigree analyzed and selected for enrichment analyses. The porcine chromosome (SSC), position (in bp) of the start and end, the total Length (in bp), number of alleles (haplotypes), and the number of homozygous genotypes observed are given for each LB
Fig. 7
figure 7

Pedigree illustrating the segregation of 17 segregant candidate alleles in the Gochu Asturcelta population analyzed. Boars are in squares and sows are in circles. Open squares and circles identify the non-carrier individuals. Letters beside the identification of the individuals mean that the individual carries segregant candidate haplotypes in the following LB (allele code in brackets): (A) LB1260 (101), LB1271 (101), LB3444 (103), LB3445 (106), LB3446 (106), LB3448 (104), LB3449 (104), LB3450 (103), LB3452 (104), LB3453 (103), LB3495 (101), LB3515 (101), LB3525 (101), LB3526 (101), LB3529 (101); (B) LB3444 (103), LB3445 (106), LB3446 (106), LB3448 (104), LB3449 (104), LB3450 (103), LB3452 (104), LB3453 (103), LB3495 (101), LB3515 (101), LB3525 (101), LB3526 (101), LB3529 (101); (C) LB3444 (103), LB3446 (106), LB3448 (104), LB3449 (104), LB3450 (103), LB3452 (104), LB3453 (103), LB3495 (101), LB3515 (101), LB3525 (101), LB3526 (101), LB3529 (101); (D) LB1260 (101), LB1271 (101), LB3495 (101), LB3515 (101), LB3525 (101), LB3526 (101), LB3529 (101); (E) LB1260 (101), LB1271 (101), LB1781 (103), LB1788 (103); (F) LB3525 (101), LB3526 (101), LB3529 (101); (G) LB1260 (101), LB1271 (101); (H) LB1781 (103), LB1788 (103); (I) LB1788 (103); and (J) LB1271 (101). Individuals carrying the same haplotypic combination are in the same color

Relationships between LB, recombination rate and CNV regions

Mean ρ estimated for the whole typed population was 2.2 ± 3.9 4Ner (Table 1). To ascertain the possible influence of recombination in the occurrence of new alleles within LB, a total of 24,752 windows, with ρ higher than the mean, overlapped with the 4,470 LB identified. A total of 212 LB, gathering a total of 967 alleles (5.7% of the total number of alleles identified), overlapped with at least one window with ρ 2.5 s.d. above the mean, considered recombination hotspots (Supplementary Table S2). These 212 LB lengthened 53.8 ± 89.1 kb and gathered 4.6 ± 3.6 alleles, on average. The mean observed homozygosity in these 212 LB (0.453) was similar to that observed at the whole population level as well (Supplementary Table S2).

Furthermore, a total of 177 LB, gathering 804 different alleles (4.7% of the total number of alleles identified; 4.5 alleles per LB) overlapped with 72 different CNV regions (Supplementary Table S2) previously identified in Gochu Asturcelta pig [38]. These 177 LB lengthened 78.3 ± 115.2 kb on average and their mean number of alleles per locus was 4.7 ± 3.6. The observed homozygosity in these 177 LB was 0.457.

Only six different LB on SSC2, SSC6 and SSC11 (LB686, LB688, LB2023, LB2024, LB3086 and LB3087), gathering a total of 45 different alleles (7.5 alleles per LB), overlapped with both recombination hotspots and CNV regions. The observed homozygosity in these 177 LB was 0.459.

Annotation and Enrichment analyses

A total of 16 protein coding-genes were located in the bounds of the 17 LB gathering the segregating candidate alleles, such as the syntaxin binding protein 5 L (STXBP5L) gene, the myosin heavy chain 15 (MYH15), the intraflagellar transport 57 (IFT57) gene, and the lipase I (LIPI) gene. A full description of these candidate genes, including their identification, description and location, retrieved from the Ensembl Genes 91 database, is given in Supplementary Table S5.

The genes identified formed two functional annotation clusters (AC; Supplementary Table S6), only one of them significant, putatively involved in immunity. The annotation cluster AC1 (enrichment score = 2.86; IPR007110:Immunoglobulin-like domain), included five candidate genes: CD86 molecule (CD86) gene, immunoglobulin like domain containing receptor 1 (ILDR1) gene, HERV-H LTR-associating 2 (HHLA2) gene, cell adhesion molecule 2 (CADM2) gene, and the roundabout guidance receptor 1 (ROBO1) gene. The non-significant annotation cluster AC2 (enrichment score = 0.32; AC2 GO:0016021  integral component of membrane) includes up to six candidate genes, four of them included in AC1 as well (CD86, ILDR1, HHLA2 and CADM2), but also, the GABA type A receptor associated protein like 2 (GABARAPL2) gene and the ectonucleotide pyrophosphatase/phosphodiesterase 2 (ENPP2) gene.

Discussion

The identification of potentially harmful alleles, by definition in very low frequency, is highly dependent on sample size. This is important because projects aiming at the identification of potentially harmful alleles require large sample sizes [8, 13, 42, 43]. Differences in sample size may partially explain the large differences in the identification of candidate haplotypes across cohorts and the WP in the current research. However, the close relationship between sample size and the number of potential candidate alleles in our data (Fig. 1), together with the high proportion of candidate alleles to the number of alleles identified (32% at the WP level), suggests that most candidate haplotypes identified are likely to be spurious.

Here, LB are identified at the BP, therefore making it possible to assume that most, if not all, allelic variants to be identified at the WP level should be present in the BP. However, our results substantially depart from that expectation. The number of alleles identified in the BP was 70% of those identified at the WP level only. Furthermore, it is intuitive to assume that only potentially harmful (candidate) allelic variants identified in the BP could be identified as candidate haplotypes in the WP. This is important because the identification of the same alleles in different generations is expected to give confidence on the allelic frequencies computed [11]. However, the number of alleles in low frequency probably being observed in heterozygous state only in the WP was 3.6-fold higher than those in the BP.

These facts illustrate a scenario in which the identification of haplotypes in very low frequency is subject to a high degree of uncertainty.

It is worth noting that our approach tried to avoid uncertainty as much as possible. Their identification at the BP gave a relatively low number of LB to be analyzed. When the presence of founders is not considered, the number of haplotypes substantially increases, virtually covering all analyzed genome adding uncertainty to the identification of potentially harmful haplotypes. When no pedigree was included in the analyses, PLINK identified roughly four-fold LB (16,793) in the whole typed population. Furthermore, the GHap R package [5] gave the identification of 20,681 different LB carrying a total of 308,497 alleles.

Allele-Drop-in events

Events such as mutation or recombination may influence haplotypic frequencies [7]. In our case, however, since a significant number of alleles appear de novo in cohorts defined by their pedigree depth (Figs. 1, 2 and 5), we could discard that mutation causes the excess of variation identified. Although the role of recombination in the occurrence of new haplotypes within LB cannot be neglected, in our data the overlap of the LB identified with genomic areas potentially considered recombination hotspots is low and, therefore, recombination cannot explain by itself the excess of new alleles identified.

Although the alleles identified within LB fit well with the Mendelian expectations of inheritance, a non-negligible number of parent-offspring inconsistencies were identified. Both the candidate alleles and their LB coincided substantially with those in which paternity errors were identified (Figs. 4 and 6). Therefore, it is likely that causes of paternal incompatibility are related to those causing the occurrence of new alleles in low frequency at the WP level.

The paternal incompatibilities found are likely dependent on the occurrence of Allele-Drop-In (ADI; i.e. alleles that are additional to the parental genotypes) events. Figures 2 and 6 illustrate this. Validation of potentially harmful haplotypes usually relies upon their inheritance from genotyped ancestors [11]. However, in our data the identification of the candidate alleles in the parents and grandparents of the carrier individuals is extremely low (Fig. 6).

The analyzed dataset was edited to remove those SNPs in which Mendelian Errors were identified [23, 27]. In SNP arrays data, a relatively low number of Mendelian Errors identified in a very high number of loci are considered ADI events caused by the influence of genomic alterations, namely Copy Number Variations, SINEs, or LINEs [27]. Interestingly, consistently with the parental incompatibilities identified in the current research, Arias et al. reported that such Mendelian Errors increase with family size [27].

It is known that not all SNP genotyping errors depart from Mendelian rules of inheritance and, therefore, would not be identifiable. We suggest that such SNP calling errors causing ADI events are likely to be the main cause of the occurrence of ADI haplotypes. When constructing haplotypes in a LB, the presence of “hidden” calling errors at a locus, probably affecting a high number of loci in low frequency, being consistent with parental genotypes, causes the apparition of wrong alleles [27]. When used to construct haplotypes, the wrongly called SNP genotypes can cause the apparition of new variants. Since the wrongly called SNPs are in low frequency [27], ADI haplotypes are likely to appear in frequencies consistent with putatively harmful alleles, potentially darkening this kind of studies. Furthermore, since calling errors can be caused by genomic alterations such as Copy Number Variations, which can be inherited from parents to offspring [38], the calling errors could also be “inherited” in a relatively high frequency (Fig. 2). This may darkness the identification of potentially lethal alleles; their frequency may be relatively high if they are associated with heterozygous advantage due to positive pleiotropic effects.

The dense Gochu Asturcelta pedigree analyzed is relatively small. When the target is to identify haplotypes in very small frequency using populations including thousands of individuals, the number of “true” potentially harmful haplotypes will be extremely low if compared with those arising from ADI events, and their identification must be considered with caution.

This particularly applies to the criteria applied to remove false candidate alleles from data. In this respect, the separated assessment of the presence of the candidate haplotypes in either the parents or the (maternal and paternal) grandparents of the carrier individuals [8, 9, 11, 43, 44], although no doubt adding confidence to the correct identification of potentially harmful alleles, may not be enough. Some studies apply Minimum Allele Frequency (MAF) thresholds (e.g. MAF > 2%; [9] to remove “false” potentially harmful alleles from datasets. In our data, no candidate allele identified in the WP reached MAF = 1%, and MAF values for the 17 segregating candidate alleles identified varied from MAF = 0.53% (alleles 103 on LB1781 and 105 on LB3445) to MAF = 0.58% (alleles 101 on LB1271, LB3525, LB3526, and LB3529; Supplementary Table S3). In previous analyses using SNP arrays data, it has been suggested that the use of MAF thresholds may not be advisable if the goal is to keep in arrays data truly informative alleles [27]. In the case of haplotype construction, ADI alleles may easily exceed MAF thresholds, particularly if the analyzed population is large and includes big offspring size. ADI haplotypes may exceed MAF thresholds even if, as in our case, alleles considered candidates should never be observed as in homozygous state.

The assessment of the segregation of the potential candidate haplotypes identified at least two generations back from the carrier individuals can be a conservative approach to dealing with the identification of potentially lethal haplotypes.

Putative role of candidate genes in inbreeding depression

Enrichment analyses were performed on the genomic areas of the LB in which segregating candidate alleles were identified to assess the possible role of the candidate genes identified in inbreeding depression. This is a difficult task because inbreeding depression effects are known to be heterogeneous along the genome, and the architecture of inbreeding depression is expected to differ between populations [45]. Causes of such differences may rely upon incomplete linkage disequilibrium between the causative lethal allele and the haplotypes inferred using marker information [9, 13, 14, 46, 47].

In any case, the 17 segregating candidate alleles identified fit well with expectations and the breeding history of the Gochu Asturcelta pig breed: they expanded in the population from founders [16], boar 20 and sow 71, and could have been purged from the population after the implementation of strict mating policies in the conservation programme [17].

A significant part of the candidate genes identified, such as the T-cell antigen receptor CD86 gene, involved in nutrient absorption and gut barrier function in pigs [48], ILDR2, involved in the formation of epithelial cells barrier [49], HHLA2, belonging to the B7 gene family can substantially regulate T-cell activity [50], GABARAPL2 involved in macroautophagy which plays essential roles during development and immunological phenomena [51], the ABCC13 gene, involved in xenobiotic defense [52], or the IFT57 gene, involved in inflammatory responses often associated with elevated cytokine concentrations [53], are expected to participate defense and immune response. However, it has been suggested that genes associated to immunity can have a pleiotropic effect [54]. Furthermore, it is worth mentioning that the GABARAPL2 gene can be involved in wound healing [55, 56], the HHLA2 gene is involved in sow fecundity and litter size by interacting with peri-implantation endometrium in Yorkshire pigs [57] and that the IFT57 gene has been associated with the regulation of the embryonic development causing short stature and brachymesophalangia in humans [58]. Although pleiotropy could be the rule in the mammal genome [44, 59] and, therefore, it is always possible to find consistency between any list of candidate genes obtained at random and a given trait under study, the candidate genes located in the bounds of the haplotypes considered candidate in this study would fit well with expectations of possible effects of potentially harmful alleles on inbreeding depression.

Conclusions

The population dynamics of potentially harmful haplotypes suggest that most alleles identified in very low frequencies in a pedigree can be considered Allele-Drop-In events. The number of candidate alleles probably never identified in homozygous state in the pedigree analyzed substantially exceeded the expectations and it cannot be explained by mutation or recombination. The de novo occurrence of haplotypes in a LB is likely to be related to SNP calling errors. In cases of dense pedigrees, such calling errors can cause that these de novo haplotypes identified in a parent or grandparent being erroneously identified in the offspring or grand-offspring as well despite its lack of biological meaning. The number of segregating haplotypes (being identifiable at least in both the parents and the grandparents of the carrier individuals) is very low, and it would be advisable to restrict the identification of potentially harmful haplotypes to that scenario. Furthermore, the use of haplotypic segregation can be a promising strategy to identify calling errors in SNP genotyping.

Data availability

All genotypic data (4,470 LB and 471 individuals) necessary to replicate the results presented are in a plain text file (16.5 MB) that can be downloaded from https://drive.google.com/file/d/1NYKNsuk6IqutnpCCLpUjDl6O7CQL6Gx6/. SNP array data are under currently analysis within the “AutoGenome” project framework. However, the data set used is available from the corresponding author upon reasonable request.

Abbreviations

BP:

Base population

4N e r :

Recombination units, where Ne is the effective population size and r the recombination rate for the window

G12:

Cohort formed by individuals with t ≤ 2

G23:

Cohort formed by individuals with pedigree depth varying from t > 2 to t ≤ 3

G3:

Cohort formed by individuals with t > 3

LB:

Linkage disequilibrium Blocks

LD:

Linkage disequilibrium

t :

Equivalents to complete generations

WP:

Whole population

References

  1. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5.

    Article  CAS  PubMed  Google Scholar 

  2. Stephens M, Scheet P. Accounting for decay of linkage disequilibrium in haplotype inference and missing-data imputation. Am J Hum Genet. 2005;76(3):449–62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The structure of Haplotype blocks in the Human Genome. Science. 2002;296(5576):2225–9.

    Article  CAS  PubMed  Google Scholar 

  4. Kim SA, Cho C-S, Kim S-R, Bull SB, Yoo YJ. A new haplotype block detection method for dense genome sequencing data based on interval graph modeling of clusters of highly correlated SNPs. Bioinformatics. 2018;34(3):388–97.

    Article  CAS  PubMed  Google Scholar 

  5. Utsunomiya YT, Milanesi M, Utsunomiya ATH, Ajmone-Marsan P, Garcia JF. GHap: an R package for genome-wide haplotyping. Bioinformatics. 2016;32(18):2861–2.

    Article  CAS  PubMed  Google Scholar 

  6. Pritchard JK, Przeworski M. Linkage disequilibrium in humans: models and data. Am J Hum Genet. 2001;69(1):1–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Bosse M, Megens H-J, Derks MFL, de Cara ÁMR, Groenen MAM. Deleterious alleles in the context of domestication, inbreeding, and selection. Evol Appl. 2019;12(1):6–17.

    Article  PubMed  Google Scholar 

  8. Howard DM, Pong-Wong R, Knap PW, Woolliams JA. Use of haplotypes to identify regions harbouring lethal recessive variants in pigs. Genet Sel Evol. 2017;49(1):57.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Hoff JL, Decker JE, Schnabel RD, Taylor JF. Candidate lethal haplotypes and causal mutations in Angus cattle. BMC Genomics. 2017;18(1):799.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Jenko J, McClure MC, Matthews D, McClure J, Johnsson M, Gorjanc G, et al. Analysis of a large dataset reveals haplotypes carrying putatively recessive lethal and semi-lethal alleles with pleiotropic effects on economically important traits in beef cattle. Genet Sel Evol. 2019;51(1):9.

    Article  PubMed  PubMed Central  Google Scholar 

  11. VanRaden PM, Olson KM, Null DJ, Hutchison JL. Harmful recessive effects on fertility detected by absence of homozygous haplotypes. J Dairy Sci. 2011;94(12):6153–61.

    Article  CAS  PubMed  Google Scholar 

  12. Zhang C, MacNeil MD, Kemp RA, Dyck MK, Plastow GS. Putative loci causing early embryonic mortality in Duroc Swine. Front Genet. 2018;9:655.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Pausch H, Schwarzenbacher H, Burgstaller J, Flisikowski K, Wurmser C, Jansen S, et al. Homozygous haplotype deficiency reveals deleterious mutations compromising reproductive and rearing success in cattle. BMC Genomics. 2015;16(1):312.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Taylor JF, Schnabel RD, Sutovsky P, Review. Genomics of bull fertility. Animal. 2018;12:s172–83.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Leroy G. Inbreeding depression in livestock species: review and meta-analysis. Anim Genet. 2014;45(5):618–28.

    Article  CAS  PubMed  Google Scholar 

  16. Yin T, Wensch-Dorendorf M, Simianer H, Swalve HH, König S. Assessing the impact of natural service bulls and genotype by environment interactions on genetic gain and inbreeding in organic dairy cattle genomic breeding programs. Animal. 2014;8(6):877–86.

    Article  CAS  PubMed  Google Scholar 

  17. Menéndez J, Álvarez I, Fernández I, Goyache F. Genealogical analysis of the Gochu Asturcelta pig breed: insights for conservation. Czech J Anim Sci. 2016;61:140–3.

    Article  Google Scholar 

  18. Menéndez J, Álvarez I, Fernández I, Menéndez-Arias NA, Goyache F. Assessing performance of single‐sample molecular genetic methods to estimate effective population size: empirical evidence from the endangered Gochu Asturcelta pig breed. Ecol Evol. 2016;6:4971–80.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Menéndez J, Álvarez I, Fernández I, de la Roza B, Goyache F. Multiple paternity in domestic pigs under equally probable natural matings – a case study in the endangered Gochu Asturcelta pig breed. Arch Anim Breed. 2015;58:217–20.

    Article  Google Scholar 

  20. Saura M, Fernández A, Varona L, Fernández AI, de Cara MÁR, Barragán C, et al. Detecting inbreeding depression for reproductive traits in Iberian pigs using genome-wide data. Genet Sel Evol. 2015;47:1.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Zhang Y, Zhuo Y, Ning C, Zhou L, Liu J-F. Estimate of inbreeding depression on growth and reproductive traits in a large White pig population. G3. 2022;12:jkac118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Rodrigáñez J, Toro MA, Rodriguez MC, Silió L. Effect of founder allele survival and inbreeding depression on litter size in a closed line of large White pigs. Anim Sci. 1998;67:573–82.

    Article  Google Scholar 

  23. Arias KD, Gutiérrez JP, Fernández I, Álvarez I, Goyache F. Approaching autozygosity in a small pedigree of Gochu Asturcelta pigs. Genet Sel Evol. 2023;55:74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Luan S. visPedigree. visPedigree: A package for tidying and drawing animal pedigree. 2018. https://github.com/luansheng/visPedigree. 2018.

  25. Maignel L, Boichard D, Verrier E. Genetic variability of French dairy breeds estimated from pedigree information. Interbull Bull. 1996;49–54.

  26. Groenen MAM, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491(7424):393–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Arias KD, Álvarez I, Gutiérrez JP, Fernandez I, Menéndez J, Menéndez-Arias NA, et al. Understanding mendelian errors in SNP arrays data using a Gochu Asturcelta pig pedigree: genomic alterations, family size and calling errors. Sci Rep. 2022;12:19686.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Browning BL, Zhou Y, Browning SR. A one-penny Imputed Genome from Next-Generation reference panels. Am J Hum Genet. 2018;103(3):338–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Browning BL, Tian X, Zhou Y, Browning SR. Fast two-stage phasing of large-scale sequence data. Am J Hum Genet. 2021;108(10):1880–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Chang CC, Chow CC, Tellier LCAM, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience. 2015;4:7.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Veroneze R, Lopes PS, Guimarães SEF, Silva FF, Lopes MS, Harlizius B, et al. Linkage disequilibrium and haplotype block structure in six commercial pig lines. J Anim Sci. 2013;91:3493–501.

    Article  CAS  PubMed  Google Scholar 

  32. Heberle H, Meirelles GV, da Silva FR, Telles GP, Minghim R. InteractiVenn: a web-based tool for the analysis of sets through Venn diagrams. BMC Bioinformatics. 2015;16:169.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Wang J. Pedigree reconstruction from poor quality genotype data. Heredity. 2019;122:719–28.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wang J. Computationally efficient sibship and parentage assignment from multilocus marker data. Genetics. 2012;191(1):183–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Gao F, Ming C, Hu W, Li H. New Software for the fast estimation of Population Recombination Rates (FastEPRR) in the genomic era. G3 Genes|Genom|Genet. 2016;6:1563–71.

    Article  CAS  PubMed Central  Google Scholar 

  36. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Arias KD, Gutiérrez JP, Fernández I, Álvarez I, Goyache F. Copy number variation regions differing in segregation patterns span different sets of genes. Animals. 2023;13:2351.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Arias KD, Gutiérrez JP, Fernández I, Menéndez-Arias NA, Álvarez I, Goyache F. Segregation patterns and inheritance rate of copy number variations regions assessed in a Gochu Asturcelta pig pedigree. Gene. 2023;854:147111.

    Article  CAS  PubMed  Google Scholar 

  39. Kinsella RJ, Kahari A, Haider S, Zamora J, Proctor G, Spudich G et al. Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database. 2011;bar030–030.

  40. Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    Article  CAS  PubMed  Google Scholar 

  41. Tipney H, Hunter L. An introduction to effective use of enrichment analysis software. Hum Genomics. 2010;4:202.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Derks MFL, Megens H-J, Bosse M, Visscher J, Peeters K, Bink MCAM, et al. A survey of functional genomic variation in domesticated chickens. Genet Sel Evol. 2018;50(1):17.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Ben Braiek M, Fabre S, Hozé C, Astruc J-M, Moreno-Romieux C. Identification of homozygous haplotypes carrying putative recessive lethal mutations that compromise fertility traits in French lacaune dairy sheep. Genet Sel Evol. 2021;53(1):41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Kadri NK, Sahana G, Charlier C, Iso-Touru T, Guldbrandtsen B, Karim L et al. T Leeb editor 2014 A 660-Kb deletion with antagonistic effects on fertility and milk production segregates at high frequency in nordic red cattle: additional evidence for the common occurrence of balancing selection in Livestock. PLoS Genet 10 1 e1004049.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Hervás-Rivero C, Srihi H, López-Carbonell D, Casellas J, Ibáñez-Escriche N, Negro S, et al. Genomic scanning of Inbreeding Depression for Litter size in two varieties of Iberian pigs. Genes. 2023;14(10):1941.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Cole JB, Null DJ, VanRaden PM. Phenotypic and genetic effects of recessive haplotypes on yield, longevity, and fertility. J Dairy Sci. 2016;99(9):7274–88.

    Article  CAS  PubMed  Google Scholar 

  47. Fritz S, Capitan A, Djari A, Rodriguez SC, Barbat A, Baur A et al. Detection of Haplotypes Associated with Prenatal Death in Dairy Cattle and Identification of Deleterious Mutations in GART, SHBG and SLC37A2. Veitia RA, editor. PLoS ONE. 2013;8:e65550.

  48. Zhao X, Schindell B, Li W, Ni L, Liu S, Wijerathne CUB, et al. Distribution and localization of porcine calcium sensing receptor in different tissues of weaned piglets1. J Anim Sci. 2019;97(6):2402–13.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Higashi T, Katsuno T, Kitajiri S, Furuse M. Deficiency of Angulin-2/ILDR1, a tricellular tight Junction-Associated membrane protein, causes deafness with cochlear hair cell degeneration in mice. PLoS ONE. 2015;10(3):e0120674.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Ahangar NK, Khalaj-Kondori M, Alizadeh N, Mokhtarzadeh A, Baghbanzadeh A, Shadbad MA, et al. Silencing tumor-intrinsic HHLA2 potentiates the anti-tumoral effect of paclitaxel on MG63 cells: another side of immune checkpoint. Gene. 2023;855:147086.

    Article  CAS  PubMed  Google Scholar 

  51. Mizushima N, Yoshimori T, Levine B. Methods in mammalian autophagy research. Cell. 2010;140(3):313–26.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Silverton L, Dean M, Moitra K. Variation and evolution of the ABC transporter genes ABCB1, ABCC1, ABCG2, ABCG5 and ABCG8: implication for pharmacogenetics and disease. Drug Metabol Drug Interact. 2011;26(4):169–79.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Wann A, Knight M. A role for IFT88/the primary cilium in the inflammatory response to interleukin-1. Cilia. 2012;1:P60.

    Article  PubMed Central  Google Scholar 

  54. Horodyska J, Wimmers K, Reyer H, Trakooljul N, Mullen AM, Lawlor PG, et al. RNA-seq of muscle from pigs divergent in feed efficiency and product quality identifies differences in immune response, growth, and macronutrient and connective tissue metabolism. BMC Genomics. 2018;19:791.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Brigger D, Torbett BE, Chen J, Fey MF, Tschan MP. Inhibition of GATE-16 attenuates ATRA-induced neutrophil differentiation of APL cells and interferes with autophagosome formation. Biochem Biophys Res Commun. 2013;438:283–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. Cho HS, Park SY, Kim SM, Kim WJ, Jung JY. Autophagy-related protein MAP1LC3C plays a crucial role in Odontogenic differentiation of Human Dental Pulp cells. Tissue Eng Regen Med. 2020;18:265–77.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Zhou C, Cheng X, Meng F, Wang Y, Luo W, Zheng E, et al. Identification and characterization of circRNAs in peri-implantation endometrium between Yorkshire and Erhualian pigs. BMC Genomics. 2023;24:412.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Thevenon J, Duplomb L, Phadke S, Eguether T, Saunier A, Avila M, et al. Autosomal recessive IFT57 hypomorphic mutation cause ciliary transport defect in unclassified oral-facial-digital syndrome with short stature and brachymesophalangia. Clin Genet. 2016;90(6):509–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Brown SDM, Lad HV. The dark genome and pleiotropy: challenges for precision medicine. Mamm Genome. 2019;30(7):212–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The Authors kindly thank the members of the GochuAsturcelta Breeders Association (ACGA; https://www.gochuasturcelta.org/) and, particularly Dr. Juan Menéndez, for their support and kind collaboration. The genotyping service was carried out at CEGEN-PRB3-ISCIII; it is supported by grant PT17/0019, of the PE I+D+i 2013-2016, funded by ISCIII and ERDF. Authors are indebted to Dr. Gonçalo Carvalho (Thermo Fisher Scientific) for his willingness and kind contribution to the execution of this project.

Funding

This research was partially funded by AEI-FEDER, grant no. PID2019-103951RB/AEI/10.13039/501100011033. AutoGenome K.D.A. was funded by AEI-ESF, grant no. PRE2020-092905.

Author information

Authors and Affiliations

Authors

Contributions

F.G. and I.A conceived and planned the project; K.D.A., F.G., and J.P.G. wrote the paper; K.D.A., F.G., and I.F. did the data analyses; K.D.A., F.G., and J.P.G. interpreted models; I.A and I.F. undertook sampling and discussed and interpreted genetic data in light of the statistical and breeding evidence; I.A. did the laboratory work. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Félix Goyache.

Ethics declarations

Ethics approval and consent to participate

SERIDA is adhered to the Ethical Committee in Research of the University of Oviedo (Spain), which ensures that all research with biological agents follows Good Laboratory Practices and European and Spanish regulations on biosecurity under the Regulation of February 13th, 2014 (BOPA no. 47, February 26th, 2014). The use of samples in the current research followed the “three-R principles”: R1: replacement, as Gochu Asturcelta pig breed has not been characterized thus far at the haplotype level, and the current study is thus not redundant with other data; R2: reduction, as samples were originally collected for previous projects being shared to maximize their contribution to research and knowledge production; and R3: refinement, as samples were obtained using standard veterinary procedures in full respect with animal welfare and minimizing stress. Tissue and hair root samples used in this project were collected by veterinary practitioners working for the Gochu Asturcelta pig Breeders’ Association (ACGA), with the permission and in presence of the owners. For this reason, permission from the Ethical Committee in Research of the University of Oviedo was not required. In all instances, ACGA veterinarians followed standard procedures and relevant national guidelines to ensure appropriate animal care. No animals were sacrificed during the execution of the current research.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Arias, K.D., Fernández, I., Gutiérrez, J.P. et al. Population dynamics of potentially harmful haplotypes: a pedigree analysis. BMC Genomics 25, 487 (2024). https://doi.org/10.1186/s12864-024-10407-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10407-x

Keywords