Genetic variability, structure and assignment of Spanish and French pig populations based on a large sampling

Admixture and breed traceability in European indigenous pig breeds and wild boar using genome-wide SNP data

Preserving diversity of indigenous pig (Sus scrofa) breeds is a key factor to (i) sustain the pork chain (both at local and global scales) including the production of high-quality branded products, (ii) enrich the animal biobanking and (iii) progress conservation policies. Single nucleotide polymorphism (SNP) chips offer the opportunity for whole-genome comparisons among individuals and breeds. Animals from twenty European local pigs breeds, reared in nine countries (Croatia: Black Slavonian, Turopolje; France: Basque, Gascon; Germany: Schwabisch-Hällisches Schwein; Italy: Apulo Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano, Sarda; Lithuania: Indigenous Wattle, White Old Type; Portugal: Alentejana, Bísara; Serbia: Moravka, Swallow-Bellied Mangalitsa; Slovenia: Krškopolje pig; Spain: Iberian, Majorcan Black), and three commercial breeds (Duroc, Landrace and Large White) were sampled and genotyped with the GeneSeek Genomic Profiler (GGP) 70 K HD porcine genotyping chip. A dataset of 51 Wild Boars from nine countries was also added, summing up to 1186 pigs (~ 49 pigs/breed). The aim was to: (i) investigate individual admixture ancestries and (ii) assess breed traceability via discriminant analysis on principal components (DAPC). Albeit the mosaic of shared ancestries found for Nero Siciliano, Sarda and Moravka, admixture analysis indicated independent evolvement for the rest of the breeds. High prediction accuracy of DAPC mark SNP data as a reliable solution for the traceability of breed-specific pig products.

Individual identification and population assignment with microsatellite markers: an evaluation in two Indian pig populations

Ability of a set of 24 microsatellite markers for individual identification and their potential for breed assignment of individuals was evaluated in 2 Indian pig populations. The cumulative probabilities of identity of 2 random individuals within a population, even with selected set of 5 loci (CGA, S0026, S0228, S0355, SW936) were 2.87×10–8 (Assamese) and 9.66×10–8 (Anakamali) and from 2 different population was 1.13×10–12. However, the population assignment precision even with all the 24 loci was only 80 (Assamese) and 88% (Ankamali). These results suggested that although this set of markers can be safely employed for identification of individuals but their utility for breed allocation in Indian pigs needs further authentication before they can be practically used for such purposes.

SINE jumping contributes to large-scale polymorphisms in the pig genomes

BACKGROUND

Molecular markers based on retrotransposon insertion polymorphisms (RIPs) have been developed and are widely used in plants and animals. Short interspersed nuclear elements (SINEs) exert wide impacts on gene activity and even on phenotypes. However, SINE RIP profiles in livestock remain largely unknown, and not be revealed in pigs.

RESULTS

Our data revealed that SINEA1 displayed the most polymorphic insertions (22.5 % intragenic and 26.5 % intergenic), followed by SINEA2 (10.5 % intragenic and 9 % intergenic) and SINEA3 (12.5 % intragenic and 5.0 % intergenic). We developed a genome-wide SINE RIP mining protocol and obtained a large number of SINE RIPs (36,284), with over 80 % accuracy and an even distribution in chromosomes (14.5/Mb), and 74.34 % of SINE RIPs generated by SINEA1 element. Over 65 % of pig SINE RIPs overlap with genes, most of them (> 95 %) are in introns. Overall, about one forth (23.09 %) of the total genes contain SINE RIPs. Significant biases of SINE RIPs in the transcripts of protein coding genes were observed. Nearly half of the RIPs are common in these pig breeds. Sixteen SINE RIPs were applied for population genetic analysis in 23 pig breeds, the phylogeny tree and cluster analysis were generally consistent with the geographical distributions of native pig breeds in China.

CONCLUSIONS

Our analysis revealed that SINEA1-3 elements, particularly SINEA1, are high polymorphic across different pig breeds, and generate large-scale structural variations in the pig genomes. And over 35,000 SINE RIP markers were obtained. These data indicate that young SINE elements play important roles in creating new genetic variations and shaping the evolution of pig genome, and also provide strong evidences to support the great potential of SINE RIPs as genetic markers, which can be used for population genetic analysis and quantitative trait locus (QTL) mapping in pig.

Evaluation of genetic richness among Vietnamese native pig breeds using microsatellite markers

We have elucidated genetic relationships of Vietnamese native pigs (VNP) using preliminarily collected samples by a single-nucleotide polymorphism (SNP) array. In order to confirm our previous results and compare with the results of a previous study using microsatellite (MS) markers, we aimed to characterize genetic diversity and population structure in wider varieties (24 breeds from 21 Provinces) of VNP across the country using 20 polymorphic MS markers recommended by ISAG/FAO (International Society for Animal Genetics/Food and Agriculture Organization) for diversity study. In this study, we collected 1,136 DNA samples of the VNPs and three exotic breeds. Our results revealed that the average number of alleles and allelic richness across the loci in VNPs were 10.0 and 7.6, which were higher than those of exotic breeds. Genomic components among VNPs were subjected to the sampling locations. Interestingly, Co Binh Thuan showed remarkable genetic feature compared to the other VNPs, because the habitation of Co Binh Thuan was relatively far from the other breeds. The results of this study provided useful information for exploitation, conservation, and development trends of the VNP breeds. More recently, African swine fever caused significant damage to most of the VNP populations. Therefore, our findings will help a reconstruction scheme of the VNP genetic resources.

Extent of third-order linkage disequilibrium in a composite line of Iberian pigs

BACKGROUND

Previous studies on linkage disequilibrium have investigated second order linkage disequilibrium in animal and plant populations. The objective of this paper was to investigate the genome-wide levels of third order linkage disequilibrium in a composite line founded by admixture of four Iberian pig strains. A model for the generation of third order linkage disequilibrium by population admixture is proposed. A computer Expectation-Maximization algorithm is developed and applied to the estimation of third order linkage disequilibrium at inter- and intra-chromosomal level using 26,347 SNPs typed in 306 sows. The relationship of third order linkage disequilibrium with physical distance was investigated over 35 million triplets in SSC12. Basic and normalized estimates of inter and intra-chromosomal third order linkage disequilibrium are reported.

RESULTS

Genome-wide analyses revealed that third order linkage disequilibrium is rather common among linked loci in this Iberian pig line. It is shown that population admixture of multiple populations may explain the observed levels of third order linkage disequilibrium although it could be generated by genetic drift. Third order linkage disequilibrium decreases rapidly up to 4 Mb and then declines slowly. The short distances between consecutive markers explain the maintenance of the observed third order linkage disequilibria levels when using a model incorporating the break-up of disequilibrium by recombination. Genome-wide testing also revealed that only 3.6% of the normalized estimates were different from 1, - 1, 0, or from a not well-defined situation in which there is only one possible value for the third order linkage disequilibrium parameter, given allele frequencies and pairwise linkage disequilibria parameters.

CONCLUSIONS

Third order linkage disequilibrium is common among linked markers in the analyzed pig line and may have been generated by population admixture of multiple populations or by genetic drift. As with second order linkage disequilibrium, the absolute value of the third order linkage disequilibrium decreases with physical distance. Normalization of third order linkage disequilibrium should be avoided for closely linked bi-allelic loci.

Landscape genomic approach to detect selection signatures in locally adapted Brazilian swine genetic groups

Samples of 191 animals from 18 different Brazilian locally adapted swine genetic groups were genotyped using Illumina Porcine SNP60 BeadChip in order to identify selection signatures related to the monthly variation of Brazilian environmental variables. Using BayeScan software, 71 SNP markers were identified as FST outliers and 60 genotypes (58 markers) were found by Samβada software in 371 logistic models correlated with 112 environmental variables. Five markers were identified in both methods, with a Kappa value of 0.073 (95% CI: 0.011-0.134). The frequency of these markers indicated a clear north-south country division that reflects Brazilian environmental differences in temperature, solar radiation, and precipitation. Global spatial territory correlation for environmental variables corroborates this finding (average Moran's I = 0.89, range from 0.55 to 0.97). The distribution of alleles over the territory was not strongly correlated with the breed/genetic groups. These results are congruent with previous mtDNA studies and should be used to direct germplasm collection for the National gene bank.

Genetic characterization of local Criollo pig breeds from the Americas using microsatellite markers

Little is known about local Criollo pig genetic resources and relationships among the various populations. In this paper, genetic diversity and relationships among 17 Criollo pig populations from 11 American countries were assessed with 24 microsatellite markers. Heterozygosities, F-statistics, and genetic distances were estimated, and multivariate, genetic structure and admixture analyses were performed. The overall means for genetic variability parameters based on the 24 microsatellite markers were the following: mean number of alleles per locus of 6.25 ± 2.3; effective number of alleles per locus of 3.33 ± 1.56; allelic richness per locus of 4.61 ± 1.37; expected and observed heterozygosity of 0.62 ± 0.04 and 0.57 ± 0.02, respectively; within-population inbreeding coefficient of 0.089; and proportion of genetic variability accounted for by differences among breeds of 0.11 ± 0.01. Genetic differences were not significantly associated with the geographical location to which breeds were assigned or their country of origin. Still, the NeighborNet dendrogram depicted the clustering by geographic origin of several South American breeds (Criollo Boliviano, Criollo of northeastern Argentina wet, and Criollo of northeastern Argentina dry), but some unexpected results were also observed, such as the grouping of breeds from countries as distant as El Salvador, Mexico, Ecuador, and Cuba. The results of genetic structure and admixture analyses indicated that the most likely number of ancestral populations was 11, and most breeds clustered separately when this was the number of predefined populations, with the exception of some closely related breeds that shared the same cluster and others that were admixed. These results indicate that Criollo pigs represent important reservoirs of pig genetic diversity useful for local development as well as for the pig industry.

Genetic characterization of four native Italian shepherd dog breeds and analysis of their relationship to cosmopolitan dog breeds using microsatellite markers

Very little research into genetic diversity of Italian native dog breeds has been carried out so far. In this study we aimed to estimate and compare the genetic diversity of four native Italian shepherd dog breeds: the Maremma, Bergamasco, Lupino del Gigante and Oropa shepherds. Therefore, some cosmopolitan dog breeds, which have been widely raised in Italy for a long time past, have also been considered to check possible influence of these dog populations on the Italian autochthonous breeds considered here. A total of 212 individuals, belonging to 10 different dog breeds, were sampled and genotyped using 18 autosomal microsatellite loci. We analyzed the genetic diversity of these breeds, within breed diversity, breed relationship and population structure. The 10 breeds considered in this study were clearly genetically differentiated from each other, regardless of current population sizes and the onset of separate breeding history. The level of genetic diversity explained 20% of the total genetic variation. The level of H E found here is in agreement with that found by other studies. The native Italian breeds showed generally higher genetic diversity compared with the long established, well-defined cosmopolitan dog breeds. As the Border Collie seems closer to the Italian breeds than the other cosmopolitan shepherd dogs considered here, a possible utilization of this breed to improve working performance in Italian traditional working shepherd dogs cannot be ignored. The data and information found here can be utilized in the organization of conservation programs planned to reduce inbreeding and to minimize loss of genetic variability.

Microsatellite based phylogeny and bottleneck studies of Iranian indigenous goat populations

Genetic analyses, structure, and bottlenecks were examined in six populations of Iranian indigenous goat using 13 microsatellite loci. The overall heterozygosity, polymorphism information content (PIC), and Shannon index values were 0.80, 0.74, and 2.14, respectively, indicating high genetic diversity. Both a phylogenetic tree and factorial correspondence analysis grouped the populations into two major clusters. Signatures for bottleneck events in the populations were examined by two methods, which suggested that bottlenecks had occurred in two Tali and Markhoz populations, whereas other populations (Raeini, Korki jonobe Khorasan, Lori, and Najdi) showed no signature of a genetic bottleneck in the recent past. The results showed that Iranian goats have high genetic diversity and may be of value to alternative breeding and conservation programs.

Empirical Selection of Informative Microsatellite Markers within Co-ancestry Pig Populations Is Required for Improving the Individual Assignment Efficiency

The Lanyu is a miniature pig breed indigenous to Lanyu Island, Taiwan. It is distantly related to Asian and European pig breeds. It has been inbred to generate two breeds and crossed with Landrace and Duroc to produce two hybrids for laboratory use. Selecting sets of informative genetic markers to track the genetic qualities of laboratory animals and stud stock is an important function of genetic databases. For more than two decades, Lanyu derived breeds of common ancestry and crossbreeds have been used to examine the effectiveness of genetic marker selection and optimal approaches for individual assignment. In this paper, these pigs and the following breeds: Berkshire, Duroc, Landrace and Yorkshire, Meishan and Taoyuan, TLRI Black Pig No. 1, and Kaohsiung Animal Propagation Station Black pig are studied to build a genetic reference database. Nineteen microsatellite markers (loci) provide information on genetic variation and differentiation among studied breeds. High differentiation index (FST) and Cavalli-Sforza chord distances give genetic differentiation among breeds, including Lanyu's inbred populations. Inbreeding values (FIS) show that Lanyu and its derived inbred breeds have significant loss of heterozygosity. Individual assignment testing of 352 animals was done with different numbers of microsatellite markers in this study. The testing assigned 99% of the animals successfully into their correct reference populations based on 9 to 14 markers ranking D-scores, allelic number, expected heterozygosity (HE) or FST, respectively. All miss-assigned individuals came from close lineage Lanyu breeds. To improve individual assignment among close lineage breeds, microsatellite markers selected from Lanyu populations with high polymorphic, heterozygosity, FST and D-scores were used. Only 6 to 8 markers ranking HE, FST or allelic number were required to obtain 99% assignment accuracy. This result suggests empirical examination of assignment-error rates is required if discernible levels of co-ancestry exist. In the reference group, optimum assignment accuracy was achievable achieved through a combination of different markers by ranking the heterozygosity, FST and allelic number of close lineage populations.

Variation of cats under domestication: genetic assignment of domestic cats to breeds and worldwide random-bred populations

Both cat breeders and the lay public have interests in the origins of their pets, not only in the genetic identity of the purebred individuals, but also in the historical origins of common household cats. The cat fancy is a relatively new institution with over 85% of its 40-50 breeds arising only in the past 75 years, primarily through selection on single-gene aesthetic traits. The short, yet intense cat breed history poses a significant challenge to the development of a genetic marker-based breed identification strategy. Using different breed assignment strategies and methods, 477 cats representing 29 fancy breeds were analysed with 38 short tandem repeats, 148 intergenic and five phenotypic single nucleotide polymorphisms. Results suggest the frequentist method of Paetkau (single nucleotide polymorphisms = 0.78, short tandem repeats = 0.88) surpasses the Bayesian method of Rannala and Mountain (single nucleotide polymorphisms = 0.56, short tandem repeats = 0.83) for accurate assignment of individuals to the correct breed. Additionally, a post-assignment verification step with the five phenotypic single nucleotide polymorphisms accurately identified between 0.31 and 0.58 of the misassigned individuals raising the sensitivity of assignment with the frequentist method to 0.89 and 0.92 for single nucleotide polymorphisms and short tandem repeats respectively. This study provides a novel multistep assignment strategy and suggests that, despite their short breed history and breed family groupings, a majority of cats can be assigned to their proper breed or population of origin, that is, race.

Genetic diversity, population structure and subdivision of local Balkan pig breeds in Austria, Croatia, Serbia and Bosnia-Herzegovina and its practical value in conservation programs

BACKGROUND

At present the Croatian Turopolje pig population comprises about 157 breeding animals. In Austria, 324 Turopolje pigs originating from six Croatian founder animals are registered. Multiple bottlenecks have occurred in this population, one major one rather recently and several more older and moderate ones. In addition, it has been subdivided into three subpopulations, one in Austria and two in Croatia, with restricted gene flow. These specificities explain the delicate situation of this endangered Croatian lard-type pig breed.

METHODS

In order to identify candidate breeding animals or gene pools for future conservation breeding programs, we studied the genetic diversity and population structure of this breed using microsatellite data from 197 individuals belonging to five different breeds.

RESULTS

The genetic diversity of the Turopolje pig is dramatically low with observed heterozygosities values ranging from 0.38 to 0.57. Split into three populations since 1994, two genetic clusters could be identified: one highly conserved Croatian gene pool in Turopoljski Lug and the"Posavina" gene pool mainly present in the Austrian population. The second Croatian subpopulation in Lonjsko Polje in the Posavina region shows a constant gene flow from the Turopoljski Lug animals.

CONCLUSIONS

One practical conclusion is that it is necessary to develop a "Posavina" boar line to preserve the "Posavina" gene pool and constitute a corresponding population in Croatia. Animals of the highly inbred herd in Turopoljski Lug should not be crossed with animals of other populations since they represent a specific phenotype-genotype combination. However to increase the genetic diversity of this herd, a program to optimize its sex ratio should be carried out, as was done in the Austrian population where the level of heterozygosity has remained moderate despite its heavy bottleneck in 1994.

Combining US and Brazilian microsatellite data for a meta-analysis of sheep (Ovis aries) breed diversity: facilitating the FAO Global Plan of Action for Conserving Animal Genetic Resources

Microsatellites are commonly used to understand genetic diversity among livestock populations. Nevertheless, most studies have involved the processing of samples in one laboratory or with common standards across laboratories. Our objective was to identify an approach to facilitate the merger of microsatellite data for cross-country comparison of genetic resources when samples were not evaluated in a single laboratory. Eleven microsatellites were included in the analysis of 13 US and 9 Brazilian sheep breeds (N = 706). A Bayesian approach was selected and evaluated with and without a shared set of samples analyzed by each country. All markers had a posterior probability of greater than 0.5, which was higher than predicted as reasonable by the software used. Sensitivity analysis indicated no difference between results with or without shared samples. Cluster analysis showed breeds to be partitioned by functional groups of hair, meat, or wool types (K = 7 and 12 of STRUCTURE). Cross-country comparison of hair breeds indicated substantial genetic distances and within breed variability. The selected approach can facilitate the merger and analysis of microsatellite data for cross-country comparison and extend the utility of previously collected molecular markers. In addition, the result of this type of analysis can be used in new and existing conservation programs.

Genetic variability and individual assignment of Chinese indigenous sheep populations (Ovis aries) using microsatellites

The purpose of this study was to assess the genetic characteristics of six breeds of Chinese local sheep using 19 microsatellite loci and to effectively validate statistical methods for individual assignment based on informative microsatellites. All the six breeds deviated from Hardy-Weinberg equilibrium expectations, while the majority of markers complied. The polymorphism information content (PIC) of overall loci for the six populations ranged from 0.283 (SRCRSP5) to 0.852 (OarVH72). Tibetan sheep were the most diverse population with the highest mean allelic richness (6.895), while Ujmuqin (UQ) harboured the lowest allelic richness (6.000). The F-statistics for the six populations were F(IS)  = -0.172, F(IT)  = -0.082 and F(ST)  = 0.077, respectively. Furthermore, the pair-wise F(IS) revealed a moderate genetic differentiation among populations (P < 0.01), indicating that all breeds can be considered genetically independent entities. The lowest genetic differentiation was between Tengchong (TC) and UQ (F(ST)  = 0.041), and the highest one was between TC and Fat-tailed Han (F(ST)  = 0.111). In comparing the three statistical models, we note that the seven microsatellite loci (MAF65, OarJMP58, SRCRSP9, MCM140, OarAE129, BM8125 and SRCRSP5) commonly used for individual assignment will ensure a powerful detection of individual origin, with accuracy up to 91.87%, when the likelihood-based method is used. Overall, these findings shed light onto the genetic characteristics of Chinese indigenous sheep and offer a set of microsatellite loci that is simple, economic and highly informative for individual assignment of Chinese sheep.