Assessing genetic diversity in Italian goat populations using AFLP® markers

Amplified fragment length polymorphism (AFLP) markers were used to investigate the genetic variation in a sample of seven goat (Capra hircus) populations. A total of 210 individuals (30 per population) were analysed using seven selected AFLP primer combinations that produced 219 clear polymorphisms. Four autochthonous goat breeds (Bionda dell'Adamello, Frisa, Orobica and Verzaschese), two primary populations, one from the Lombardy Alps (Val di Livo) and the other from Sardinia island (Sarda) and a reference cosmopolitan breed (Saanen) were included in the analysis. The expected heterozygosity (Het) did not differ significantly among breeds (range 0.21-0.24). No breed specific markers were identified. The variability at AFLP loci was largely maintained within breeds, as indicated by the coefficient of genetic differentiation (Gst) value (0.11). Dice similarities calculated between pairs of individuals belonging to the same or to different breeds largely overlapped. Bootstrapping on markers indicated that the coefficient of variation (CV) of the genetic indexes tested decreases only marginally by adding markers over 100 AFLPs. Cluster analysis based on standard genetic distance between breeds indicates that Sarda is the most distant population, while Bionda, Frisa, Verzaschese and Val di Livo seem to be highly related populations. Interestingly, Saanen is closer than Orobica to the other four goat populations of the Lombardy Alps. Principal co-ordinates analysis based on Dice similarities confirms these observations. Genetic diversity of the goat populations investigated confirms what is expected on the basis of their geographical location. Results from Orobica are not correlated with geographical distances and may reflect undocumented migrations and gene flows and identify an original genetic resource.     

Signatures of selection in five Italian cattle breeds detected by a 54K SNP panel

In this study we used a medium density panel of SNP markers to perform population genetic analysis in five Italian cattle breeds. The BovineSNP50 BeadChip was used to genotype a total of 2,935 bulls of Piedmontese, Marchigiana, Italian Holstein, Italian Brown and Italian Pezzata Rossa breeds. To determine a genome-wide pattern of positive selection we mapped the F _st values against genome location. The highest F _st peaks were obtained on BTA6 and BTA13 where some candidate genes are located. We identified selection signatures peculiar of each breed which suggest selection for genes involved in milk or meat traits. The genetic structure was investigated by using a multidimensional scaling of the genetic distance matrix and a Bayesian approach implemented in the STRUCTURE software. The genotyping data showed a clear partitioning of the cattle genetic diversity into distinct breeds if a number of clusters equal to the number of populations were given. Assuming a lower number of clusters beef breeds group together. Both methods showed all five breeds separated in well defined clusters and the Bayesian approach assigned individuals to the breed of origin. The work is of interest not only because it enriches the knowledge on the process of evolution but also because the results generated could have implications for selective breeding programs.     

A machine learning approach for the identification of population-informative markers from high-throughput genotyping data: application to several pig breeds

Single nucleotide polymorphisms (SNPs) able to describe population differences can be used for important applications in livestock, including breed assignment of individual animals, authentication of mono-breed products and parentage verification among several other applications. To identify the most discriminating SNPs among thousands of markers in the available commercial SNP chip tools, several methods have been used. Random forest (RF) is a machine learning technique that has been proposed for this purpose. In this study, we used RF to analyse PorcineSNP60 BeadChip array genotyping data obtained from a total of 2737 pigs of 7 Italian pig breeds (3 cosmopolitan-derived breeds: Italian Large White, Italian Duroc and Italian Landrace, and 4 autochthonous breeds: Apulo-Calabrese, Casertana, Cinta Senese and Nero Siciliano) to identify breed informative and reduced SNP panels using the mean decrease in the Gini Index and the Mean Decrease in Accuracy parameters with stability evaluation. Other reduced informative SNP panels were obtained using Delta, Fixation index and principal component analysis statistics, and their performances were compared with those obtained using the RF-defined panels using the RF classification method and its derived Out Of Bag rates and correct prediction proportions. Therefore, the performances of a total of six reduced panels were evaluated. The correct assignment of the animals to its breed was close to 100% for all tested approaches. Porcine chromosome 8 harboured the largest number of selected SNPs across all panels. Many SNPs were included in genomic regions in which previous studies identified signatures of selection or genes (e.g. ESR1, KITL and LCORL) that could contribute to explain, at least in part, phenotypically or economically relevant traits that might differentiate cosmopolitan and autochthonous pig breeds. Random forest used as preselection statistics highlighted informative SNPs that were not the same as those identified by other methods. This might be due to specific features of this machine learning methodology. It will be interesting to explore if the adaptation of RF methods for the identification of selection signature regions could be able to describe population-specific features that are not captured by other approaches.     

Assessing genetic diversity in indigenous Veneto chicken breeds using AFLP markers

Genetic variation in four indigenous chicken breeds from the Veneto region of Italy was assessed using amplified fragment length polymorphism (AFLP) markers. A total of 99 individuals were analysed using three AFLP primer combinations that produced 70 polymorphisms. Four indigenous Veneto chicken breeds (Ermellinata, Padovana, Pépoi and Robusta) and a reference broiler line were included in the analysis. Breed-specific markers were identified in each breed. The expected heterozygosity did not differ significantly among the indigenous Veneto chicken breeds and the broiler line. The coefficient of gene variation (Gst) value across loci indicated that almost half of the total variability was observed among breeds. Nei's standard genetic distance between pairs of breeds showed that the distance between the broiler line and the Pépoi breed was greater than the distances between the broiler line and the other three chicken breeds. Cluster analysis based on standard genetic distances between breeds indicated that the Padovana and Pépoi breeds were closely related. Factorial analysis based on a binary matrix of the AFLP data showed a clear distinction of all breeds.     

Using 2-DE for the differentiation of local chicken breeds

The aim of this study was to apply a proteomic approach for the characterisation of local chicken breeds. The experiment involved a total of 29 males of three local Italian chicken breeds: Pépoi, Padovana and Ermellinata di Rovigo. Sarcoplasmic protein fractions of breast muscle were analysed by 2-DE. Image analysis followed by statistical analysis enabled to differentiate groups of individuals based on the similarities of protein expression. Individuals were distinguished into clusters and groups, corresponding to the breed of origin. Distances among individuals, calculated using data on spot volumes, were used to draw a neighbour-joining tree, showing clear individual and breed grouping. The most relevant spots regarding breed differentiation were detected; 11 were identified by MS revealing preliminary evidences on the mechanisms of the breed differentiation process. The results evidenced the ability of proteomic analyses for the characterisation of chicken breeds.     

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30–35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.     

Genetic characterization and breed assignment in five Italian sheep breeds using microsatellite markers

The knowledge of the genetic relationship and admixture among neighbouring populations is crucial for conservation efforts. The aim of this study was to analyse the genetic diversity of five Italian sheep breeds (Appenninica, Garfagnina Bianca, Massese, Pomarancina and Zerasca) using a panel of 24 microsatellite markers. Blood samples from 226 individuals belonging to the aforementioned populations were obtained and genotyped. All the investigated breeds showed a significant heterozygote deficiency caused by the high level of inbreeding indicated also by the high level of F_IS (0.146). Genetic differentiation between breeds was moderate (F_ST = 0.05) but significant and the individuals could be assigned to their breeds with an high success rate even if the inter-individual distances showed that few animals clustered separately from the other individuals of the same breed, especially for Pomarancina breed. The genetic distances reflect the historical knowledge of these breeds and some patterns of ancestral and recent gene flow between neighbour populations arise. The clustering analysis detects the presence of six clusters. Massese and Zerasca breeds were grouped together as well as Appenninica and Pomarancina with the latter forming two distinct clusters equally represented. The formation of this last breed is occurred with the absorption of individuals of the Appenninica breed and the gene flow probably continued in these recent years allowing the presence of a population substructure for Pomarancina breed. Such substructure supports the high level of heterozygote deficiency found for this breed despite the relatively high population size. The five populations analysed presented some genetic similarities but a clear uniqueness of the populations has been showed for almost all of them. Special attention to monitor genetic variability and to organize mating plans should be given especially for the three endangered breeds.     

Synaptonemal complex analysis of the Holstein-Friesian, Piemontese and Simmental breeds of Bos taurus taurus

The synaptonemal complex (SC) of specimens of Bos taurus taurus from the Holstein-Friesian, Piemontese, and Simmental breeds, was analysed. The analysis included quantification of the frequency of various types of abnormalities in the SC, and the frequency of cells with SC abnormalities. All animals had 29 autosomal bivalents and one sexual bivalent and the most frequently recorded abnormality was pairing failure. The number of cells with abnormalities in the Holstein-Friesian breed was 29.41%, in the Piemontese breed was 30.00% and in the Simmental breed it was 29.54%. The subspecies Bos taurus taurus had 29.63% of cells showing abnormalities with 57.33% of these abnormalities occurring in zygotene and 42.67% occurring in pachytene. Statistical analyses showed that there were no significant differences in the number of cells with SC abnormalities among the breeds studied. The frequency of cells with abnormalities, and the efect on the fertility of the Holstein-Friesian, Piemontese and Simmental breeds are discussed.     

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.     

Genomic variation and population structure detected by single nucleotide polymorphism arrays in Corriedale,Merino and Creole sheep

The aim of this study was to investigate the genetic diversity within and among three breeds of sheep: Corriedale, Merino and Creole. Sheep from the three breeds (Merino n = 110, Corriedale n = 108 and Creole n = 10) were genotyped using the Illumina Ovine SNP50 beadchip^®. Genetic diversity was evaluated by comparing the minor allele frequency (MAF) among breeds. Population structure and genetic differentiation were assessed using STRUCTURE software, principal component analysis (PCA) and fixation index (F_ST). Fixed markers (MAF = 0) that were different among breeds were identified as specific breed markers. Using a subset of 18,181 single nucleotide polymorphisms (SNPs), PCA and STUCTURE analysis were able to explain population stratification within breeds. Merino and Corriedale divergent lines showed high levels of polymorphism (89.4% and 86% of polymorphic SNPs, respectively) and moderate genetic differentiation (F_ST = 0.08) between them. In contrast, Creole had only 69% polymorphic SNPs and showed greater genetic differentiation from the other two breeds (F_ST = 0.17 for both breeds). Hence, a subset of molecular markers present in the OvineSNP50 is informative enough for breed assignment and population structure analysis of commercial and Creole breeds.     

Genetic Polymorphism of Russian,European, and Asian Chicken Breeds as Revealed with DNA and Protein Markers

The variation in polymorphic DNA (RAPD and minisatellite) and protein markers was compared for nine Russian chicken breeds differing in morphological and productivity types and in origin, three European egg breeds, and three meat breeds of the Asian origin. Genetic diversity indices were calculated for each breed group and each marker type and were used to construct dendrograms of genetic similarity. In all breed groups, minisatellites and RAPD markers revealed higher genetic diversity as compared with protein markers. With any type of markers, genetic diversity of the Russian and Asian meat breeds proved to be significantly higher than that of the European egg breeds. The differentiating potentialities of molecular and genetic biochemical markers at the breed level and the origin of the Russian chicken breeds are discussed.     

Genetic polymorphism of Russian,European, and Asian chicken breeds as revealed with DNA and protein markers

The variation in polymorphic DNA (RAPD and minisatellite) and protein markers was compared for nine Russian chicken breeds differing in morphological and productivity types and in origin, three European egg breeds, and three broiler breeds of the Asian origin. Genetic diversity indices were calculated for each breed group and each marker type and were used to construct dendrograms of genetic similarity. In all breed groups, minisatellites and RAPD markers revealed higher genetic diversity as compared with protein markers. With any type of markers, genetic diversity of the Russian and Asian broiler breeds proved to be significantly higher than that of the European egg breeds. The differentiating potentialities of molecular and genetic biochemical markers at the breed level and the origin of the Russian chicken breeds are discussed.     

Genome‐wide genetic variation discovery in Chinese Taihu pig breeds using next generation sequencing

The Chinese Taihu pig breeds are an invaluable component of the world's pig genetic resources, and they are the most prolific breeds of swine in the world. In this study, the genomes of 252 pigs of the six indigenous breeds in the Taihu Lake region were sequenced using the genotyping by genome reducing and sequencing approach. A total of 950 million good reads were obtained using an Illumina Hiseq2000 at an average depth of 13× (for SNP calling) and an average coverage of 2.3%. In total, 122 632 indels, 31 444 insertions, 44 056 deletions and 455 CNVs (copy number variants) were identified in the genomes of the pigs. Approximately 2.3% of these genetic markers were mapped to gene exon regions, and 25% were in QTL regions related to economically important traits. The KEGG pathway or GO enrichment analyses revealed that genetic variants assumed to be large‐effect mutations were significantly overrepresented in 22 SNP, 56 indel, 26 insertion, 28 deletion and three CNV gene sets. A total of 343 breed‐specific SNPs were also identified in the six Chinese indigenous pigs. The findings from this study can contribute to future investigations of the genetic diversity, population structure, positive selection signals and molecular evolutionary history of these pigs at the genome level and can serve as a valuable reference for improving the breeding and cultivation of these pigs.     

Identification and functional characterization of copy number variations in diverse chicken breeds

Background

The detection and functional characterization of genomic structural variations are important for understanding the landscape of genetic variation in the chicken. A recently recognized aspect of genomic structural variation, called copy number variation (CNV), is gaining interest in chicken genomic studies. The aim of the present study was to investigate the pattern and functional characterization of CNVs in five characteristic chicken breeds, which will be important for future studies associating phenotype with chicken genome architecture.

Results

Using a commercial 385 K array-based comparative genomic hybridization (aCGH) genome array, we performed CNV discovery using 10 chicken samples from four local Chinese breeds and the French breed Houdan chicken. The female Anka broiler was used as a reference. A total of 281 copy number variation regions (CNVR) were identified, covering 12.8 Mb of polymorphic sequences or 1.07% of the entire chicken genome. The functional annotation of CNVRs indicated that these regions completely or partially overlapped with 231 genes and 1032 quantitative traits loci, suggesting these CNVs have important functions and might be promising resources for exploring differences among various breeds. In addition, we employed quantitative PCR (qPCR) to further validate several copy number variable genes, such as prolactin receptor, endothelin 3 (EDN3), suppressor of cytokine signaling 2, CD8a molecule, with important functions, and the results suggested that EDN3 might be a molecular marker for the selection of dark skin color in poultry production. Moreover, we also identified a new CNVR (chr24: 3484617–3512275), encoding the sortilin-related receptor gene, with copy number changes in only black-bone chicken.

Conclusions

Here, we report a genome-wide analysis of the CNVs in five chicken breeds using aCGH. The association between EDN3 and melanoblast proliferation was further confirmed using qPCR. These results provide additional information for understanding genomic variation and related phenotypic characteristics.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-934) contains supplementary material, which is available to authorized users.     

Insights into the genetic history of Green‐legged Partridgelike fowl: mtDNA and genome‐wide SNP analysis

The Green‐legged Partridgelike (GP) fowl, an old native Polish breed, is characterised by reseda green‐coloured shanks rather than yellow, white, slate or black commonly observed across most domestic breeds of chicken. Here, we investigate the origin, genetic relationships and structure of the GP fowl using mtDNA D‐loop sequencing and genome‐wide SNP analysis. Genome‐wide association analysis between breeds enables us to verify the genetic control of the reseda green shank phenotype, a defining trait for the breed. Two mtDNA D‐loop haplogroups and three autosomal genetic backgrounds are revealed. Significant associations of SNPs on chromosomes GGA24 and GGAZ indicate that the reseda green leg phenotype is associated with recessive alleles linked to the W and Id loci. Our results provide new insights into the genetic history of European chicken, indicating an admixd origin of East European traditional breeds of chicken on the continent, as supported by the presence of the reseda green phenotype and the knowledge that the GP fowl as a breed was developed before the advent of commercial stocks.     

Tracking the emergence of a new breed using 49,034 SNP in sheep

Domestic animals are unique in that they have been organised into managed populations called breeds. The strength of genetic divergence between breeds may vary dependent on the age of the breed, the scenario under which it emerged and the strength of reproductive isolation it has from other breeds. In this study, we investigated the Gulf Coast Native breed of sheep to determine if it contains lines of animals that are sufficiently divergent to be considered separate breeds. Allele sharing and principal component analysis (PCA) using nearly 50,000 SNP loci revealed a clear genetic division that corresponded with membership of either the Florida or Louisiana Native lines. Subsequent analysis aimed to determine if the strength of the divergence exceeded that found between recognised breed pairs. Genotypes from 14 breeds sampled from Europe and Asia were used to obtain estimates of pair-wise population divergence measured as F _ST. The divergence separating the Florida and Louisiana Native (F _ST = 6.2%) was approximately 50% higher than the average divergence separating breeds developed within the same region of Europe (F _ST = 4.2%). This strongly indicated that the two Gulf Coast Native lines are sufficiently different to be considered separate breeds. PCA using small SNP sets successfully distinguished between the Florida and Louisiana Native animals, suggesting that allele frequency differences have accumulated across the genome. This is consistent with a population history involving geographic separation and genetic drift. Suggestive evidence was detected for divergence at the poll locus on sheep chromosome 10; however drift at neutral markers has been the largest contributor to the genetic separation observed. These results document the emergence of populations that can be considered separate breeds, with practical consequences for bio-conservation priorities, animal registration and the establishment of separate breed societies.     

Genetic Diversity Among Turkish Native Chickens, Denizli and Gerze, Estimated by Microsatellite Markers

The genetic diversity of the Turkish native chicken breeds Denizli and Gerze was evaluated with 10 microsatellite markers. We genotyped a total of 125 individuals from five subpopulations. Among loci, the mean number of alleles was 7.5, expected heterozygosity (H (e)) was 0.665, PIC value was 0.610, and Wright's fixation index was 0.301. H (e) was higher in the Denizli breed (0.656) than in the Gerze breed (0.475). The PIC values were 0.599 and 0.426 for Denizli and Gerze, respectively. A phylogenetic tree was constructed using genetic distance and the neighbor-joining method. Its topology reflects the general pattern of genetic differentiation among the Denizli and Gerze breeds. The present study suggests that Denizli and Gerze subpopulations have a rich genetic diversity. The information about Denizli and Gerze breeds estimated by microsatellite analysis may also be useful as an initial guide in defining objectives for designing future investigations of genetic variation and developing conservation strategies.     

Mitochondrial D-loop sequence variation among Italian horse breeds

The genetic variability of the mitochondrial D-loop DNA sequence in seven horse breeds bred in Italy (Giara, Haflinger, Italian trotter, Lipizzan, Maremmano, Thoroughbred and Sarcidano) was analysed. Five unrelated horses were chosen in each breed and twenty-two haplotypes were identified. The sequences obtained were aligned and compared with a reference sequence and with 27 mtDNA D-loop sequences selected in the GenBank database, representing Spanish, Portuguese, North African, wild horses and an Equus asinus sequence as the outgroup. Kimura two-parameter distances were calculated and a cluster analysis using the Neighbour-joining method was performed to obtain phylogenetic trees among breeds bred in Italy and among Italian and foreign breeds. The cluster analysis indicates that all the breeds but Giara are divided in the two trees, and no clear relationships were revealed between Italian populations and the other breeds. These results could be interpreted as showing the mixed origin of breeds bred in Italy and probably indicate the presence of many ancient maternal lineages with high diversity in mtDNA sequences.     

Genetic diversity among chicken breeds estimated through randomly amplified polymorphic DNA.

The randomly amplified polymorphic DNA (RAPD) markers were used to detect polymorphism among five breeds of chicken i.e. White Leghorn and Rhodes Island Red (selected for part period egg production and egg mass respectively), Red Cornish and White Plymouth Rock (selected for early body weights) and Kadaknath (native breed). Twelve of the fifty random primers screened yielded distinct polymorphic RAPD profiles. Of the total 96 fragments amplified, about 25% showed polymorphism. Using the RAPD data matrix, the within population and between population genetic similarity was estimated. The selected improved breeds showed higher within population genetic similarity in comparison to the native breed. The two meat type breeds showed a high level of genetic similarity between themselves. The White Leghorn breed showed a low genetic similarity with other breeds. The native breed showed highest similarity with Rhodes Island Red. The dendogram was constructed to show phylogenetic relationship among these breeds. As expected, the genetic distances were lowest within similar type breeds and were highest between dissimilar type breeds. The results indicated the effectiveness of RAPD in detecting polymorphism between chicken populations and their applicability in population studies and establishing genetic relationships among the chicken populations.