Mitochondrial D-loop analysis reveals low diversity in Mangalica pigs and their relationship to historical specimens

Population Subdivision and Migration Assessment of Mangalica Pig Breeds Based on Pedigree Analysis

In conserving the genetic diversity of domestic animal breeds, strategies that emphasise between-breed diversity may not be optimal, as they neglect within-breed variation. The aim of the present study was to assess the extent of population subdivision in three Mangalica pig breeds and the contribution of migration to their substructure. Wright's FST coefficient was calculated based on genealogical data with breeding animals born between 1981 and 2023, with three colour variants (Blonde, Swallow-Belly and Red). These Wright's FST coefficients were analysed using multidimensional scaling to reveal the population substructure. The average FST coefficient was 0.04 for the Blonde breed and 0.047 for the Swallow-Belly and Red Mangalica breeds, while these parameters were lower in the active herds at 0.03 and 0.04, respectively. The migration of individuals between herds was 61.63% for the Blonde breed and 75.53% and 63.64% for the Swallow-Belly and Red Magalica breeds, respectively. No population substructure was observed in any of the Mangalica breeds, which can be explained by the extensive migration between herds.

An insight into the runs of homozygosity distribution and breed differentiation in Mangalitsa pigs

Mangalitsa pigs exhibit three distinct coat color patterns based on which they are described as Red, Blond, and Swallow-bellied. The current study investigated genome-wide diversity and selection signatures in the three breeds using fixation index, runs of homozygosity and population structure analyses. The analyses were originally based on quality-controlled data on 77 Mangalitsa animals from Germany, including 23 Blond, 30 Swallow-bellied and 24 Red Mangalitsa genotyped with a customized version of the ProcineSNP60 v2 Genotyping Bead Chip. Also, 20 Hungarian Mangalitsa genotypes were included as outgroup data for comparison. Estimates of observed heterozygosity were 0.27, 0.28, and 0.29, and inbreeding coefficients estimated based on runs of homozygosity were 24.11%, 20.82%, and 16.34% for Blond, Swallow-bellied and Red Mangalitsa, respectively. ROH islands were detected in all breeds, however, none of these were shared amongst them. The KIF16B gene previously reported to play a role in synaptic signaling was found in a ROH island (SSC17: 16–26) in Swallow-bellied Mangalitsa. The same gene was found to harbor a significantly differentiated SNP (MARC0032380) while contrasting either Blond or Red to Swallow-belied Mangalitsa. In the Red Mangalitsa, some ROH islands were associated with genes that play a role in meat quality traits, i.e., ABCA12, VIL1, PLSCR5, and USP37. Our population structure analysis highlighted a separation of the three breeds, but also showed the closest relatedness between Red and Blond Mangalitsa pigs. Findings of this study improve our understanding of the diversity in the three breeds of Mangalitsa pigs.

Genetic diversity and population structure of six autochthonous pig breeds from Croatia, Serbia, and Slovenia

BACKGROUND

The importance of local breeds as genetic reservoirs of valuable genetic variation is well established. Pig breeding in Central and South-Eastern Europe has a long tradition that led to the formation of several local pig breeds. In the present study, genetic diversity parameters were analysed in six autochthonous pig breeds from Slovenia, Croatia and Serbia (Banija spotted, Black Slavonian, Turopolje pig, Swallow-bellied Mangalitsa, Moravka and Krskopolje pig). Animals from each of these breeds were genotyped using microsatellites and single nucleotide polymorphisms (SNPs). The results obtained with these two marker systems and those based on pedigree data were compared. In addition, we estimated inbreeding levels based on the distribution of runs of homozygosity (ROH) and identified genomic regions under selection pressure using ROH islands and the integrated haplotype score (iHS).

RESULTS

The lowest heterozygosity values calculated from microsatellite and SNP data were observed in the Turopolje pig. The observed heterozygosity was higher than the expected heterozygosity in the Black Slavonian, Moravka and Turopolje pig. Both types of markers allowed us to distinguish clusters of individuals belonging to each breed. The analysis of admixture between breeds revealed potential gene flow between the Mangalitsa and Moravka, and between the Mangalitsa and Black Slavonian, but no introgression events were detected in the Banija spotted and Turopolje pig. The distribution of ROH across the genome was not uniform. Analysis of the ROH islands identified genomic regions with an extremely high frequency of shared ROH within the Swallow-bellied Mangalitsa, which harboured genes associated with cholesterol biosynthesis, fatty acid metabolism and daily weight gain. The iHS approach to detect signatures of selection revealed candidate regions containing genes with potential roles in reproduction traits and disease resistance.

CONCLUSIONS

Based on the estimation of population parameters obtained from three data sets, we showed the existence of relationships among the six pig breeds analysed here. Analysis of the distribution of ROH allowed us to estimate the level of inbreeding and the extent of homozygous regions in these breeds. The iHS analysis revealed genomic regions potentially associated with phenotypic traits and allowed the detection of genomic regions under selection pressure.

The full mitochondrial genomes of Mangalica pig breeds and their possible origin

The mitogenomes of one animal of each of the three Mangalica breeds, Blonde, Red, and Swallow-belly were assembled from reads obtained by Next Generation Sequencing of the three genomes. Features of the mitogenomes were identical in the three breeds, apart from a second tRNA-Val gene on the L strand in Swallow-belly. Phylogenetic comparison of the three mitogenomes with 112 full mtDNA sequences clearly put Mangalicas into the European clade. Comparing the mitogenome of eight Mangalica animals revealed particular differences between them. The mitogenome of some Mangalicas was closely related to the Croatian Turopolje breed and this indicates either the common origin of their maternal lineages or admixture of some populations of the breeds. However, the origin of the mitogenome of certain purebred Mangalicas kept in the Hungarian Mangalica Gene Reserve still remains unknown.

The 'heritability' of domestication and its functional partitioning in the pig

We propose to estimate the proportion of variance explained by regression on genome-wide markers (or genomic heritability) when wild/domestic status is considered the phenotype of interest. This approach differs from the standard Fst in that it can accommodate genetic similarity between individuals in a general form. We apply this strategy to complete genome data from 47 wild and domestic pigs from Asia and Europe. When we partitioned the total genomic variance into components associated to subsets of single nucleotide polymorphisms (SNPs) defined in terms of their annotation, we found that potentially deleterious non-synonymous mutations (9566 SNPs) explained as much genetic variance as the whole set of 25 million SNPs. This suggests that domestication may have affected protein sequence to a larger extent than regulatory or other kinds of mutations. A pathway-guided analysis revealed ovarian steroidogenesis and leptin signaling as highly relevant in domestication. The genomic regression approach proposed in this study revealed molecular processes not apparent through typical differentiation statistics. We propose that at least some of these processes are likely new discoveries because domestication is a dynamic process of genetic selection, which may not be completely characterized by a static metric like Fst. Nevertheless, and despite some particularly influential mutation types or pathways, our analyses tend to rule out a simplistic genetic basis for the domestication process: neither a single pathway nor a unique set of SNPs can explain the process as a whole.

Romanian wild boars and Mangalitza pigs have a European ancestry and harbour genetic signatures compatible with past population bottlenecks

We aimed to analyse the genetic diversity of Romanian wild boars and to compare it with that from other wild boar and pig populations from Europe and Asia. Partial sequencing of the mitochondrial encoded cytochrome b (MT-CYB) gene from 36 Romanian wild boars and 36 domestic pigs (Mangalitza, Bazna and Vietnamese breeds) showed that the diversity of Romanian wild boars and Mangalitza pigs is fairly reduced, and that most of the members of these two populations share a common MT-CYB haplotype. Besides, in strong contrast with the Bazna animals, Romanian wild boars and Mangalitza swine did not carry Asian variants at the MT-CYB locus. The autosomal genotyping of 18 Romanian wild boars with the Illumina Porcine SNP60 BeadChip revealed that their genetic background is fundamentally European, even though signs of a potential Near Eastern ancestry (~25%) were detectable at K = 4 (the most significant number of clusters), but not at higher K-values. Admixture analysis also showed that two wild boars are of a hybrid origin, which could be explained by the mating of feral animals with domestic pigs. Finally, a number of Romanian wild boars displayed long runs of homozygosity, an observation that is consistent with the occurrence of past population bottlenecks and the raise of inbreeding possibly due to overhunting or to the outbreak of infectious diseases.

Genome sequencing and analysis of Mangalica, a fatty local pig of Hungary

BACKGROUND

Mangalicas are fatty type local/rare pig breeds with an increasing presence in the niche pork market in Hungary and in other countries. To explore their genetic resources, we have analysed data from next-generation sequencing of an individual male from each of three Mangalica breeds along with a local male Duroc pig. Structural variations, such as SNPs, INDELs and CNVs, were identified and particular genes with SNP variations were analysed with special emphasis on functions related to fat metabolism in pigs.

RESULTS

More than 60 Gb of sequence data were generated for each of the sequenced individuals, resulting in 11× to 19× autosomal median coverage. After stringent filtering, around six million SNPs, of which approximately 10% are novel compared to the dbSNP138 database, were identified in each animal. Several hundred thousands of INDELs and about 1,000 CNV gains were also identified. The functional annotation of genes with exonic, non-synonymous SNPs, which are common in all three Mangalicas but are absent in either the reference genome or the sequenced Duroc of this study, highlighted 52 genes in lipid metabolism processes. Further analysis revealed that 41 of these genes are associated with lipid metabolic or regulatory pathways, 49 are in fat-metabolism and fatness-phenotype QTLs and, with the exception of ACACA, ANKRD23, GM2A, KIT, MOGAT2, MTTP, FASN, SGMS1, SLC27A6 and RETSAT, have not previously been associated with fat-related phenotypes.

CONCLUSIONS

Genome analysis of Mangalica breeds revealed that local/rare breeds could be a rich source of sequence variations not present in cosmopolitan/industrial breeds. The identified Mangalica variations may, therefore, be a very useful resource for future studies of agronomically important traits in pigs.

Introgression and isolation contributed to the development of Hungarian Mangalica pigs from a particular European ancient bloodline

Background

Mangalica breeds are indigenous to Hungary and their breeding history dates back to about 200–250 years ago. They are fat-type pigs and have a rare curly hair phenotype. The aim of our study was to establish the relationships between these unique breeds and other European breeds.

Results

Based on a core sequence of 382 bp present in 2713 mitochondrial D-loop sequences from pigs belonging to 38 local breeds from nine countries, five cosmopolitan breeds and wild boars from 14 countries, we identified 164 haplotypes. More than half of the 2713 sequences belonged to either four haplotypes characteristic of continental European breeds or two haplotypes characteristic of British/cosmopolitan breeds; each haplotype is present in more than 100 individuals. Most Mangalica individuals belonged either to one of these common continental European haplotypes or to two Mangalica-specific haplotypes that were absent in all other breeds. In addition, we identified the ancestral mitochondrial D-loop signature present in these 2713 sequences and found that ~ 80% carried the European ancient signatures, ANC-Aside and ANC-Cside or their closely related signatures, while most of the remaining sequences carried a modern Asian signature, ANC-Easia. Mangalica individuals carried the ANC-Aside signature, but not the ANC-Cside or ANC-Easia signatures.

Conclusions

In all the Mangalica individuals, a unique ancient European signature was found in the mitochondrial DNA D-loop region, but they belonged almost exclusively to either certain very abundant European or two Mangalica-specific D-loop haplotypes. This indicates that the present-day Mangalica population in Hungary evolved either by introgression of other European breeds and wild boars or via total isolation after the divergence of European ancient porcine bloodlines.