Medium-density Single Nucleotide Polymorphism Research Articles

The Use of Genomic Screening for the Detection of Chromosomal Abnormalities in the Domestic Horse: Five New Cases of 65,XXY Syndrome in the Pura Raza Español Breed.

Sex chromosomal abnormalities are a well-established cause of reproductive failure in domestic horses. Because of its difficult diagnosis, the Pura Raza Español breeding program established a routine screening for chromosomal abnormalities in all the horses prior to enrolling in the studbook. This genomic procedure combines an initial assessment based on the results from Short Tandem Repeat (STR) parentage testing followed by a Single-Nucleotide Polymorphism (SNP) based copy number aberration (CNA) confirmative analysis in positive cases. Using this methodology, we identified five new individuals carrying a 65,XXY chromosomal number aberration (CNA) among 27,330 foals enrolled over the past two reproductive seasons. The animals were initially flagged as CNA candidates due to abnormal results in STR testing. Subsequent analysis genotyping using an STR sex-linked dedicated panel and a medium-density SNP array in ECAX and ECAY confirmed the diagnosis as 65,XXY carriers. Four cases (upon sample availability) underwent further analysis using in situ fluorescent hybridization with ECAX and ECAY probes, showing identical results. Phenotypic analysis revealed abnormal gonad development in one of the cases, showing that the remaining four had a normal reproductive morphology. To our knowledge, this represents the largest number of horses exhibiting the equine form of Klinefelter syndrome (65,XXY) reported to date. Our study highlights the importance of genomic screening in the accurate detection of chromosomal abnormalities in horses.

Evaluation of DNA yield from various tissue and sampling sources for use in single nucleotide polymorphism panels

Genetics studies are used by wildlife managers and researchers to gain inference into a population of a species of interest. To gain these insights, microsatellites have been the primary method; however, there currently is a shift from microsatellites to single nucleotide polymorphisms (SNPs). With the different DNA requirements between microsatellites and SNPs, an investigation into which samples can provide adequate DNA yield is warranted. Using samples that were collected from previous genetic projects from regions in the USA from 2014 to 2021, we investigated the DNA yield of eight sample categories to gain insights into which provided adequate DNA to be used in ddRADseq or already developed high- or medium-density SNP panels. We found seven sample categories that met the DNA requirements for use in all three panels, and one sample category that did not meet any of the three panels requirements; however, DNA integrity was highly variable and not all sample categories that met panel DNA requirements could be considered high quality DNA. Additionally, we used linear random-effects models to determine which covariates would have the greatest influence on DNA yield. We determined that all covariates (tissue type, storage method, preservative, DNA quality, time until DNA extraction and time after DNA extraction) could influence DNA yield.

Comprehensive analysis of runs of homozygosity and heterozygosity in Holstein cattle on the basis of medium and high density SNP panels and large population sample

Abstract This study reports runs of homozygosity (ROH) and heterozygosity (ROHet) distributed in a large population of Holstein cattle on the basis of two microarrays of medium (50k; 2163 animals; 54 609 SNPs) and high single nucleotide polymorphism (SNP) density (HD; 600 animals; 777 692 SNPs). To assess the inbreeding values of Holstein cattle, the ROH-based genomic inbreeding coefficient (FROH) was calculated. The comparison of SNP panels suggested that FROH values above 4 Mb should be considered for panels of medium densities as a relatively reliable measure of inbreeding. Moreover, ROH hotspots and coldspots were identified and compared between the HD and 50k SNP panels and were carefully examined for association with production and functional traits. The obtained results pinpointed genomic regions presumably under selection pressure in Holstein cattle. The regions overlapped with a large number of genes, including GHR, GBF1, SUMF1, CCL28, NIM1K, U6, BTRC and FABP1, many of which are involved in important Holstein cattle characteristics. We also found that some ROH hotspots and coldspots identified with the HD panel were not detected with the 50k panel, mainly because of insufficient SNP density in certain genomic regions. This suggests that using medium-density panels might not be the best choice when precise identification of ROH patterns is the main goal. In summary, in this work, we confirmed that a high-density SNP panel compared to a medium-density SNP panel allows for more precise identification of ROH patterns, especially in the case of short ROH that could be associated with ancestral inbreeding.

PSXII-10 Development of a 70k Snp Genotyping Array for American Mink (Neogale vison)

Abstract Mink is an important fur animal species; however, little is known about the genomics of complex traits in mink. The development of single nucleotide polymorphism (SNP) array has accelerated genomics application in animal breeding by facilitating genotyping and associated analyses. The current project aimed to design a medium-density SNP array in American mink for genomic studies and selection in this species. For this purpose, whole genome sequence with 30x coverage was generated for 100 individuals from the Canadian Centre for Fur Animal Research (CCFAR) at the Dalhousie Faculty of Agriculture (Truro, NS, Canada) and Millbank Fur Farm (Rockwood, ON, Canada). Variant calling was done using Samtools and GATK pipelines. A total of 8,373,854 bi-allelic SNPs identified by both pipelines were used for an initial selection of 91K SNPs for submission to the Affymetrix scoring system. After filtering, 62,375 SNPs were selected for probe design and production of an Axiom Mink Genotyping Array 70K. This 70K SNP array was used to genotype 2,973 mink for validation. After quality control of genotypes, 47,159 SNPs were classified as high-resolution SNPs, and among them 24,237 were polymorphic. A preliminary genome-wide association study of 1,356 mink for body length at harvest identified a SNP (AX-647623051) significantly associated (P < 1.26 × 10-5) with the trait. In conclusion, this 70K mink SNP array could be used for genomic studies and potentially as a starting point for genomic selection programs in American mink.

Genetic structure and genome-wide association study of the traditional Kazakh horses.

Horses are traditionally used in Kazakhstan as a source of food and as working and saddle animals as well. Here, for the first time, microarray-based medium-density single nucleotide polymorphism (SNP) genotyping of six traditionally defined types and breeds of indigenous Kazakh horses was conducted to reveal their genetic structure and find markers associated with animal size and weight. The results showed that the predefined separation between breeds and sampled populations was not supported by the molecular data. The lack of genetic variation between breeds and populations was revealed by the principal component analysis, ADMIXTURE, and distance-based analyses, as well as the general population parameters expected and observed heterozygosity (He and Ho) and between-group fixation index (Fst). The analysis revealed that the studied types and breeds should be considered as a single breed, namely the 'Kazakh horse'. The comparison with previously published data on global horse breed diversity revealed the relatively high level of individual diversity of Kazakh horses in comparison with the well-known foreign breeds. The Mongolian and Tuva breeds were identified as the closest horse landraces, demonstrating similar patterns of internal variability. The genome-wide association analysis was performed for animal size and weight as the traits directly related with the meat productivity of horses. The analysis identified a set of 60 SNPs linked with horse genes involved in the regulation of processes of development of connective tissues and the bone system, neural system, immune system regulation, and other processes. The present study is novel and introduces Kazakh horses as a promising genetic source for horse breeding and selection both on the domestic and international levels.

Population genomics of American mink using genotype data.

Understanding the genetic structure of the target population is critically important to develop an efficient genomic selection program in domestic animals. In this study, 2,973 American mink of six color types from two farms (Canadian Centre for Fur Animal Research (CCFAR), Truro, NS and Millbank Fur Farm (MFF), Rockwood, ON) were genotyped with the Affymetrix Mink 70K panel to compute their linkage disequilibrium (LD) patterns, effective population size (Ne), genetic diversity, genetic distances, and population differentiation and structure. The LD pattern represented by average r 2, decreased to <0.2 when the inter-marker interval reached larger than 350kb and 650kb for CCFAR and MFF, respectively, and suggested at least 7,700 and 4,200 single nucleotide polymorphisms (SNPs) be used to obtain adequate accuracy for genomic selection programs in CCFAR and MFF respectively. The Ne for five generations ago was estimated to be 76 and 91 respectively. Our results from genetic distance and diversity analyses showed that American mink of the various color types had a close genetic relationship and low genetic diversity, with most of the genetic variation occurring within rather than between color types. Three ancestral genetic groups was considered the most appropriate number to delineate the genetic structure of these populations. Black (in both CCFAR and MFF) and pastel color types had their own ancestral clusters, while demi, mahogany, and stardust color types were admixed with the three ancestral genetic groups. This study provided essential information to utilize the first medium-density SNP panel for American mink in their genomic studies.

Assessment of Genetic Diversity, Runs of Homozygosity, and Signatures of Selection in Tropical Milking Criollo Cattle Using Pedigree and Genomic Data.

The objective of this study was to evaluate the genetic diversity of the Tropical Milking Criollo cattle (TMC) breed in Mexico through parameters derived from pedigree and genomic information assessment. The pedigree file consisted of 3780 animals. Seventy-nine bovines were genotyped with the medium-density single nucleotide polymorphism chip and considered a reference population for pedigree analysis. The effective population size and the probability of gene origin used to assess the evolution of genetic diversity were calculated with pedigree information. Inbreeding coefficients were evaluated based on pedigree (), the genomic relationship matrix (), and runs of homozygosity () of different length classes. The average inbreeding was 2.82 ± 2.66%, −0.7 ± 3.8%, and 10.9 ± 3.0% for , , and , respectively. Correlation between and was significant only for runs of homozygosity > 4 Mb, indicating the of a population with an average equivalent complete generation of five only recovers the most recent inbreeding. The parameters of the probability of gene origin indicated the existence of genetic bottlenecks and the loss of genetic diversity in the history of the TMC cattle population; however, pedigree and genomic information revealed the existence of current sufficient genetic diversity to design a sustainable breeding program.

Replacement of microsatellite markers by imputed medium-density SNP arrays for parentage control in German warmblood horses

In horses, parentage control is currently performed based on an internationally standardized panel of 17 microsatellite (MS) markers comprising 12 mandatory and five optional markers. Unlike MS, single nucleotide polymorphism (SNP) profiles support a wider portfolio of genomic applications, including parentage control. A transition to SNP-based parentage control is favorable, but requires additional efforts for ensuring generation-overlapping availability of marker genotypes of the same type. To avoid double genotyping of either parents or offspring for changing to SNP technology and enable efficient transition, we tested whether MS genotypes used for parentage control could be reliably imputed from a medium-density SNP panel in German warmblood horses. Imputation accuracy was tested in a tenfold cross-validation with two approaches: within breed (option A) and across breeds (option B). Average imputation accuracies of 97.98% (A) and 96.17% (B) were achieved, respectively. Due to interbreed differences in genotyping rates, five MS markers of low genotyping rate (GTR; < 90%) could be imputed with higher accuracy within breed (98.18%) than across breeds (90.73%). MS markers with high GTR performed homogeneously well in option B (98.44%) and showed slightly lower accuracy in option A (97.90%). Among these markers, AHT5 proved to be problematic for imputation regardless of the approach, revealing accuracies of 86.40% (A) and 88.70% (B). Better results for MS markers with high GTR and savings in computational processing justified the choice of option B for routine implementation. To date, more than 9500 horses have undergone the new parentage control based on imputed MS genotypes.

Impact of linkage disequilibrium heterogeneity along the genome on genomic prediction and heritability estimation

BackgroundCompared to medium-density single nucleotide polymorphism (SNP) data, high-density SNP data contain abundant genetic variants and provide more information for the genetic evaluation of livestock, but it has been shown that they do not confer any advantage for genomic prediction and heritability estimation. One possible reason is the uneven distribution of the linkage disequilibrium (LD) along the genome, i.e., LD heterogeneity among regions. The aim of this study was to effectively use genome-wide SNP data for genomic prediction and heritability estimation by using models that control LD heterogeneity among regions.MethodsThe LD-adjusted kinship (LDAK) and LD-stratified multicomponent (LDS) models were used to control LD heterogeneity among regions and were compared with the classical model that has no such control. Simulated and real traits of 2000 dairy cattle individuals with imputed high-density (770K) SNP data were used. Five types of phenotypes were simulated, which were controlled by very strongly, strongly, moderately, weakly and very weakly tagged causal variants, respectively. The performances of the models with high- and medium-density (50K) panels were compared to verify that the models that controlled LD heterogeneity among regions were more effective with high-density data.ResultsCompared to the medium-density panel, the use of the high-density panel did not improve and even decreased prediction accuracies and heritability estimates from the classical model for both simulated and real traits. Compared to the classical model, LDS effectively improved the accuracy of genomic predictions and unbiasedness of heritability estimates, regardless of the genetic architecture of the trait. LDAK applies only to traits that are mainly controlled by weakly tagged causal variants, but is still less effective than LDS for this type of trait. Compared with the classical model, LDS improved prediction accuracy by about 13% for simulated phenotypes and by 0.3 to ~ 10.7% for real traits with the high-density panel, and by ~ 1% for simulated phenotypes and by − 0.1 to ~ 6.9% for real traits with the medium-density panel.ConclusionsGrouping SNPs based on regional LD to construct the LD-stratified multicomponent model can effectively eliminate the adverse effects of LD heterogeneity among regions, and greatly improve the efficiency of high-density SNP data for genomic prediction and heritability estimation.

Analysis of the Genetic Diversity and Population Structure of Four Senegalese Sheep Breeds Using Medium-Density Single-Nucleotide Polymorphisms.

Simple SummaryThis paper reported genetic parameters of four Senegalese sheep breeds, in relation to inbreeding, diversity and genetic proximity. The results provide informations on genetic conservation and adaptability of the breeds in the Senegalese context.In Senegal, sheep breeds have adapted to their environment and play a key socio-economic role. This study aimed to explore the genetic diversity and structure of four Senegalese sheep breeds (Peul-peul, Djallonke, Touabire, and Ladoum) and their relationships with global sheep breeds. To that end, forty-seven sheep were genotyped using the OvineSNP50 BeadChip, and these genotypic data were analysed with those of 73 sheep breeds representative of worldwide ovine diversity (2729 animals). The average observed heterozygosity (Ho) ranged from 0.293 in Djallonke sheep to 0.339 in Touabire sheep. The estimated Fis values were low, ranging from 0.019 for Ladoum to 0.034 for Peul-peul sheep. The estimated Fst values were low (0.003–0.044) among the trypanosusceptible breeds (Peul-peul, Touabire, and Ladoum) but high between the previous breeds and the trypanotolerant Djallonke breed (0.075–0.116), indicating better genetic conservation of the Djallonke sheep. A principal component analysis revealed clustering of the Senegalese sheep breeds according to their geographic distribution. However, owing to genetic improvement practices, the introgression of Touabire sheep blood seems to have reshaped the genetic landscape of the trypanosusceptible sheep breeds in Senegal. The Senegalese sheep breeds showed lower genetic diversity than their presumed ancestral sheep breeds of the Middle East. They also presented some relatedness with Caribbean sheep breeds, which reveals their contribution to the global genetic diversity and to the development of Caribbean sheep breeds.

Phylogenetic relationships of Argentinean Creole with other Latin American Creole Cattle as revealed by a medium density Single Nucleotide polymorphism microarray

The demographic history of Creole cattle in Argentina and in Latin America dates back to the time of the Spanishcolonization. This study aimed to investigate the potential use of a medium­density SNP array to describe cattlefromthe most representative andoldest herds of the Argentine Creole cattle breedregisteredat the time of theconstitution of the Argentinean Creole Cattle Breeders Association and to explore the phylogenetic relationship withCreole cattle fromother LatinAmericancountries. Toachieve this goal, genotypes from51 animals on34,008autosomal SNP were used to generate genetic distance matrices based on the proportion of shared identical­by­statealleles among individual animals and animals clustered according to their origin, analyzed by the PLINKprogram.Aneighbor­joining phylogenetic tree based on pairwise genetic distance was constructed using PHYLIP and wasprepared for visualization using FigTree. Amultidimensional scaling analysis was performed to evaluate the levelof relationship in terms of genetic distance among the different animal clusters. Genetic distances between animalsvaried from0.186 to 0.357 when considering all pairs of animals, and from0.186 to 0.338 when considering Creolepairs. The dendrogramobtained showed three major clusters. Cluster 1 included Latin American Creole cattle fromColombia, Guadalupe, Paraguay, and Uruguay, and the reference groups of Holstein and Jersey cattle. Cluster 2contained exclusively Patagonian Creole cattle, while the third cluster included the remaining Argentinean Creoles.The genetic relationship patterns obtainedvia multidimensional scaling showeda close relationship among fourgroups of Creole animals fromArgentina. The closeness between clusters can be explained in part on the basis ofearly migration of animals that gave rise to founders herds at some Argentinean locations. The outcomes of thisstudycontribute toabetter understandingof the compositionof the earlyfounder herds of Creole cattle inArgentina and the relationship with other Latin America Creole cattle populations.

Genome-Wide Analysis Reveals Selection Signatures Involved in Meat Traits and Local Adaptation in Semi-Feral Maremmana Cattle.

The Maremmana cattle is an ancient Podolian-derived Italian breed raised in semi-wild conditions with distinctive morphological and adaptive traits. The aim of this study was to detect potential selection signatures in Maremmana using medium-density single nucleotide polymorphism array. Putative selection signatures were investigated combining three statistical approaches designed to quantify the excess of haplotype homozygosity either within (integrated haplotype score, iHS) or among pairs of populations (Rsb and XP-EHH), and contrasting the Maremmana with a single reference population composed of a pool of seven Podolian-derived Italian breeds. Overall, the three haplotype-based analyses revealed selection signatures distributed over 19 genomic regions. Of these, six relevant candidate regions were identified by at least two approaches. We found genomic signatures of selective sweeps spanning genes related to mitochondrial function, muscle development, growth, and meat traits (SCIN, THSD7A, ETV1, UCHL1, and MYOD1), which reflects the different breeding schemes between Maremmana (semi-wild conditions) and the other Podolian-derived Italian breeds (semi-extensive). We also identified several genes linked to Maremmana adaptation to the environment of the western-central part of Italy, known to be hyperendemic for malaria and other tick-borne diseases. These include several chemokine (C-C motif) ligand genes crucially involved in both innate and adaptive immune responses to intracellular parasite infections and other genes playing key roles in pulmonary disease (HEATR9, MMP28, and ASIC2) or strongly associated with malaria resistance/susceptibility (AP2B1). Our results provide a glimpse into diverse selection signatures in Maremmana cattle and can be used to enhance our understanding of the genomic basis of environmental adaptation in cattle.

Leveraging the 3000 Rice Genome Data for Computational Design of Polymorphic Markers in a Local Rice Variety Lacking Sequence Data

DNA markers can detect DNA sequence variations in the genome, and they are useful for genetic studies, DNA fingerprinting, and genotype-based selection in breeding programs. Rice, as one of the model plants for genetic and genomic studies, has abundant DNA markers stored in various online databases. Selecting markers in rice is not limited by marker availability but rather by their polymorphism in the target population. We developed a computational method to screen millions of single nucleotide polymorphism (SNP) markers listed in IRRI 3000 rice genome database in order to find a subset of markers that are polymorphic in an F2 mapping population created from a cross between a parental line with a known genome sequence and a local Indonesian variety with no genome sequence data. The parental lines were genotyped using an affordable medium-density SNP array. The genotype data was cross-referenced with the rice genome database to perform phylogenetic analysis and identify accessions clusters with the highest genetic similarities to each parental line. The cluster data was then used to identify monomorphic SNP candidates within the cluster but exhibit consistent polymorphism between the two clusters. Using this method, we obtained a SNP marker set for a segment in rice chromosome 8 with 76.19% polymorphism rate, which is much higher than the expected 1.06% polymorphism rate if the SNP markers were chosen randomly. The improved polymorphism rate was also observed when the method was applied to other random chromosome segments and randomly chosen parental candidates.

Inference of Ancestries and Heterozygosity Proportion and Genotype Imputation in West African Cattle Populations.

Several studies have evaluated computational methods that infer the haplotypes from population genotype data in European cattle populations. However, little is known about how well they perform in African indigenous and crossbred populations. This study investigates: (1) global and local ancestry inference; (2) heterozygosity proportion estimation; and (3) genotype imputation in West African indigenous and crossbred cattle populations. Principal component analysis (PCA), ADMIXTURE, and LAMP-LD were used to analyse a medium-density single nucleotide polymorphism (SNP) dataset from Senegalese crossbred cattle. Reference SNP data of East and West African indigenous and crossbred cattle populations were used to investigate the accuracy of imputation from low to medium-density and from medium to high-density SNP datasets using Minimac v3. The first two principal components differentiated Bos indicus from European Bos taurus and African Bos taurus from other breeds. Irrespective of assuming two or three ancestral breeds for the Senegalese crossbreds, breed proportion estimates from ADMIXTURE and LAMP-LD showed a high correlation (r ≥ 0.981). The observed ancestral origin heterozygosity proportion in putative F1 crosses was close to the expected value of 1.0, and clearly differentiated F1 from all other crosses. The imputation accuracies (estimated as correlation) between imputed and the real data in crossbred animals ranged from 0.142 to 0.717 when imputing from low to medium-density, and from 0.478 to 0.899 for imputation from medium to high-density. The imputation accuracy was generally higher when the reference data came from the same geographical region as the target population, and when crossbred reference data was used to impute crossbred genotypes. The lowest imputation accuracies were observed for indigenous breed genotypes. This study shows that ancestral origin heterozygosity can be estimated with high accuracy and will be far superior to the use of observed individual heterozygosity for estimating heterosis in African crossbred populations. It was not possible to achieve high imputation accuracy in West African crossbred or indigenous populations based on reference data sets from East Africa, and population-specific genotyping with high-density SNP assays is required to improve imputation.

A publicly available repository of ROH islands reveals signatures of selection in different livestock and pet species

BackgroundRuns of homozygosity (ROH) have become the state-of-the-art method for analysis of inbreeding in animal populations. Moreover, ROH are suited to detect signatures of selection via ROH islands and are used in other applications, such as genomic prediction and genome-wide association studies (GWAS). Currently, a vast amount of single nucleotide polymorphism (SNP) data is available online, but most of these data have never been used for ROH analysis. Therefore, we performed a ROH analysis on large medium-density SNP datasets in eight animal species (cat, cattle, dog, goat, horse, pig, sheep and water buffalo; 442 different populations) and make these results publicly available.ResultsThe results include an overview of ROH islands per population and a comparison of the incidence of these ROH islands among populations from the same species, which can assist researchers when studying other (livestock) populations or when looking for similar signatures of selection. We were able to confirm many known ROH islands, for example signatures of selection for the myostatin (MSTN) gene in sheep and horses. However, our results also included multiple other ROH islands, which are common to many populations and not identified to date (e.g. on chromosomes D4 and E2 in cats and on chromosome 6 in sheep).ConclusionsWe are confident that our repository of ROH islands is a valuable reference for future studies. The discovered ROH island regions represent a unique starting point for new studies or can be used as a reference for future studies. Furthermore, we encourage authors to add their population-specific ROH findings to our repository.

Genomic selection in American mink (Neovison vison) using a single-step genomic best linear unbiased prediction model for size and quality traits graded on live mink.

Genomic selection relies on single-nucleotide polymorphisms (SNPs), which are often collected using medium-density SNP arrays. In mink, no such array is available; instead, genotyping by sequencing (GBS) can be used to generate marker information. Here, we evaluated the effect of genomic selection for mink using GBS. We compared the estimated breeding values (EBVs) from single-step genomic best linear unbiased prediction (SSGBLUP) models to the EBV from ordinary pedigree-based BLUP models. We analyzed seven size and quality traits from the live grading of brown mink. The phenotype data consisted of ~20,600 records for the seven traits from the mink born between 2013 and 2016. Genotype data included 2,103 mink born between 2010 and 2014, mostly breeding animals. In total, 28,336 SNP markers from 391 scaffolds were available for genomic prediction. The pedigree file included 29,212 mink. The predictive ability was assessed by the correlation (r) between progeny trait deviation (PTD) and EBV, and the regression of PTD on EBV, using 5-fold cross-validation. For each fold, one-fifth of animals born in 2014 formed the validation set. For all traits, the SSGBLUP model resulted in higher accuracies than the BLUP model. The average increase in accuracy was 15% (between 3% for fur clarity and 28% for body weight). For three traits (body weight, silky appearance of the under wool, and guard hair thickness), the difference in r between the two models was significant (P < 0.05). For all traits, the regression slopes of PTD on EBV from SSGBLUP models were closer to 1 than regression slopes from BLUP models, indicating SSGBLUP models resulted in less bias of EBV for selection candidates than the BLUP models. However, the regression coefficients did not differ significantly. In conclusion, the SSGBLUP model is superior to conventional BLUP model in the accurate selection of superior animals, and, thus, it would increase genetic gain in a selective breeding program. In addition, this study shows that GBS data work well in genomic prediction in mink, demonstrating the potential of GBS for genomic selection in livestock species.

Homozygosity indicators in canine MHC region of the Standard Poodle and Leonberger populations

Homozygosity indicators in canine MHC region of the Standard Poodle and Leonberger populations

The genetic heritage of Alpine local cattle breeds using genomic SNP data

BackgroundAssessment of genetic diversity and population structure provides important control metrics to avoid genetic erosion, inbreeding depression and crossbreeding between exotic and locally-adapted cattle breeds since these events can have deleterious consequences and eventually lead to extinction. Historically, the Alpine Arc represents an important pocket of cattle biodiversity with a large number of autochthonous breeds that provide a fundamental source of income for the entire regional economy. By using genotype data from medium-density single nucleotide polymorphism (SNP) arrays, we performed a genome-wide comparative study of 23 cattle populations from the Alpine Arc and three cosmopolitan breeds.ResultsAfter filtering, we obtained a final genotyping dataset consisting of 30,176 SNPs for 711 individuals. The local breeds showed high or intermediate values of genetic diversity compared to the highly selected cosmopolitan breeds. Patterns of genetic differentiation, multidimensional scaling, admixture analysis and the constructed phylogenetic tree showed convergence, which indicates the presence of gene flow among the breeds according to both geographic origin and historical background. Among the most differentiated breeds, we identified the modern Brown cattle. In spite of admixture events, several local breeds have preserved distinctive characteristics, which is probably due to differences in genetic origin and geographic location.ConclusionsThis study represents one of the most comprehensive genome-wide analysis of the Alpine cattle breeds to date. Using such a large dataset that includes the majority of the local breeds found in this region, allowed us to expand knowledge on the evaluation and status of Alpine cattle biodiversity. Our results indicate that although many of the analyzed local breeds are listed as endangered, they still harbor a large amount of genetic diversity, even when compared to some cosmopolitan breeds. This finding, together with the reconstruction of the phylogeny and the relationships between these Alpine Arc cattle breeds, provide crucial insights not only into the improvement of genetic stocks but also into the implementation of future conservation strategies.

Identification of Ancestry Informative Marker (AIM) Panels to Assess Hybridisation between Feral and Domestic Sheep.

Simple SummaryOnce present in the entirety of Europe, mouflon (wild sheep) became extinct due to intense hunting, but remnant populations survived and became feral on the Mediterranean islands of Corsica and Sardinia. Although now protected by regional laws, Sardinian mouflon is threatened by crossbreeding with domestic sheep causing genetic hybridisation. The spread of domestic genes can be detrimental for wild populations as it dilutes the genetic features that characterise them. This work aimed to identify diagnostic tools that could be applied to monitor the level of hybridisation between mouflon and domestic sheep. Tens of thousands of genetic markers known as single nucleotide polymorphisms (SNPs) were screened and we identified the smallest number of SNPs necessary to discriminate between pure mouflon and sheep. We produced four SNP panels of different sizes which were able to assess the hybridisation level of a mouflon and we verified that the SNP panels efficacy is independent of the domestic sheep breed involved in the hybrid. The implementation of these results into actual diagnostic tools will help the conservation of this unique and irreplaceable mouflon population, and the methodology applied can easily be transferred to other case studies of interest. Hybridisation of wild populations with their domestic counterparts can lead to the loss of wildtype genetic integrity, outbreeding depression, and loss of adaptive features. The Mediterranean island of Sardinia hosts one of the last extant autochthonous European mouflon (Ovis aries musimon) populations. Although conservation policies, including reintroduction plans, have been enforced to preserve Sardinian mouflon, crossbreeding with domestic sheep has been documented. We identified panels of single nucleotide polymorphisms (SNPs) that could act as ancestry informative markers able to assess admixture in feral x domestic sheep hybrids. The medium-density SNP array genotyping data of Sardinian mouflon and domestic sheep (O. aries aries) showing pure ancestry were used as references. We applied a two-step selection algorithm to this data consisting of preselection via Principal Component Analysis followed by a supervised machine learning classification method based on random forest to develop SNP panels of various sizes. We generated ancestry informative marker (AIM) panels and tested their ability to assess admixture in mouflon x domestic sheep hybrids both in simulated and real populations of known ancestry proportions. All the AIM panels recorded high correlations with the ancestry proportion computed using the full medium-density SNP array. The AIM panels proposed here may be used by conservation practitioners as diagnostic tools to exclude hybrids from reintroduction plans and improve conservation strategies for mouflon populations.

Concordance rate between copy number variants detected using either high- or medium-density single nucleotide polymorphism genotype panels and the potential of imputing copy number variants from flanking high density single nucleotide polymorphism haplotypes in cattle

BackgroundThe trading of individual animal genotype information often involves only the exchange of the called genotypes and not necessarily the additional information required to effectively call structural variants. The main aim here was to determine if it is possible to impute copy number variants (CNVs) using the flanking single nucleotide polymorphism (SNP) haplotype structure in cattle. While this objective was achieved using high-density genotype panels (i.e., 713,162 SNPs), a secondary objective investigated the concordance of CNVs called with this high-density genotype panel compared to CNVs called from a medium-density panel (i.e., 45,677 SNPs in the present study). This is the first study to compare CNVs called from high-density and medium-density SNP genotypes from the same animals. High (and medium-density) genotypes were available on 991 Holstein-Friesian, 1015 Charolais, and 1394 Limousin bulls. The concordance between CNVs called from the medium-density and high-density genotypes were calculated separately for each animal. A subset of CNVs which were called from the high-density genotypes was selected for imputation. Imputation was carried out separately for each breed using a set of high-density SNPs flanking the midpoint of each CNV. A CNV was deemed to be imputed correctly when the called copy number matched the imputed copy number.ResultsFor 97.0% of CNVs called from the high-density genotypes, the corresponding genomic position on the medium-density of the animal did not contain a called CNV. The average accuracy of imputation for CNV deletions was 0.281, with a standard deviation of 0.286. The average accuracy of imputation of the CNV normal state, i.e. the absence of a CNV, was 0.982 with a standard deviation of 0.022. Two CNV duplications were imputed in the Charolais, a single CNV duplication in the Limousins, and a single CNV duplication in the Holstein-Friesians; in all cases the CNV duplications were incorrectly imputed.ConclusionThe vast majority of CNVs called from the high-density genotypes were not detected using the medium-density genotypes. Furthermore, CNVs cannot be accurately predicted from flanking SNP haplotypes, at least based on the imputation algorithms routinely used in cattle, and using the SNPs currently available on the high-density genotype panel.