The present study aimed at analyzing the genetic variability of indigenous goat breeds (Capra hircus) using the MHC-associated microsatellite markers BF1, BM1818, BM1258, DYMS1, and SMHCC1. The following breeds were included: Chinese Xuhuai, Indian Changthangi and Pashmina, Kenyan Small East African (SEA) and Galla, and Albanian Vendi. To examine genetic variability, the levels of heterozigosity, degrees of inbreeding, and genetic differences among the breeds were analyzed. The mean number of alleles ranged from nine in the Galla to 14.5 in the Vendi breed. The mean observed heterozygosity and mean expected heterozygosity varied from 0.483 in the Vendi to 0.577 in the Galla breed, and from 0.767 in the SEA to 0.879 in the Vendi breed, respectively. Significant loss of heterozygosity (p < 0.01) indicated that these loci were not in Hardy-Weinberg equilibrium. The mean FIS values ranged from 0.3299 in the SEA to 0.4605 in the Vendi breed with a mean value of 0.3623 in all breeds (p < 0.001). Analysis of molecular variance indicated that 7.14% and 4.74% genetic variation existed among the different breeds and geographic groups, whereas 92.86% and 95.26% existed in the breeds and the geographic groups, respectively (p < 0.001). The microsatellite marker analysis disclosed a high degree of genetic polymorphism. Loss of heterozygosity could be due to genetic drift and endogamy. The genetic variation among populations and geographic groups does not indicate a correlation of genetic differences with geographic distance.
genetic variability; indigenous goat breeds; major histocompatibility complex; microsatellites
Genetic variation in Barbari goats, a highly prolific breed distributed widely in the northern part of India, known for better milk and meat quality, was studied as a part of genetic characterization and conservation. The genomic DNA from 50 unrelated Barbari goats were amplified via PCR with a panel of 21 microsatellite markers, and resolved through 6 per cent denaturing polyacrylamide gel electrophoresis followed by silver staining. The number of alleles ranged from 4 to 11, with allele sizes ranging from 88 to 220 bp. The distribution of allele frequencies was between 0.0104 and 0.5208. Polymorphism information content varied from 0.5563 to 0.8348. The population was not in Hardy-Weinberg equilibrium for all except two microsatellite loci (ILSTS044 and ILSTS060). The observed heterozygosity ranged from 0.8478 to 1.0000 while the expected heterozygosity ranged from 0.6208 to 0.8509. Based on the results of the present study, there is a good scope for exploiting the genetic variability in the Barbari goats for further improvement of performance.
Barbari goats; genetic characterization; heterozygosity; microsatellite markers; polymorphism
Ovine rinderpest or goat plague is an economically important and contagious viral disease of sheep and goats, caused by the Peste des petits ruminants virus (PPRV). Differences in susceptibility to goat plague among different breeds and water buffalo exist. The host innate immune system discriminates between pathogen associated molecular patterns and self antigens through surveillance receptors known as Toll like receptors (TLR). We investigated the role of TLR and cytokines in differential susceptibility of goat breeds and water buffalo to PPRV. We examined the replication of PPRV in peripheral blood mononuclear cells (PBMC) of Indian domestic goats and water buffalo and demonstrated that the levels of TLR3 and TLR7 and downstream signalling molecules correlation with susceptibility vs resistance. Naturally susceptible goat breeds, Barbari and Tellichery, had dampened innate immune responses to PPRV and increased viral loads with lower basal expression levels of TLR 3/7. Upon stimulation of PBMC with synthetic TLR3 and TLR7 agonists or PPRV, the levels of proinflammatory cytokines were found to be significantly higher while immunosuppressive interleukin (IL) 10 levels were lower in PPRV resistant Kanni and Salem Black breeds and water buffalo at transcriptional level, correlating with reduced viralloads in infected PBMC. Water buffalo produced higher levels of interferon (IFN) α in comparison with goats at transcriptional and translational levels. Pre-treatment of Vero cells with human IFNα resulted in reduction of PPRV replication, confirming the role of IFNα in limiting PPRV replication. Treatment with IRS66, a TLR7 antagonist, resulted in the reduction of IFNα levels, with increased PPRV replication confirming the role of TLR7. Single nucleotide polymorphism analysis of TLR7 of these goat breeds did not show any marked nucleotide differences that might account for susceptibility vs resistance to PPRV. Analyzing other host genetic factors might provide further insights on susceptibility to PPRV and genetic polymorphisms in the host.
The aim of this study was to investigate allele frequencies at the CSN1S2 locus in two Chinese dairy goat breeds and the effects of its variation on dairy goat economic traits. Seven hundred and eight goats from Xinong Saanen (XS, n = 268) and Guanzhong (GZ, N = 440) breeds were selected. The milk samples of 268 XS goats were collected during the middle of lactation, body size parameters (708 goats) and daily milk yield (202 goats) were registered. The RFLP (restriction fragment length polymorphism) and SSCP (single strand conformation polymorphism) were used to detect the polymorphisms in CSN1S2. The Hardy-Weinberg (HW) equilibrium and the associations between body size, milk yield and composition and the genotypes were calculated. The results revealed that only A and F CSN1S2 alleles were found in the two Chinese dairy goat breeds. Allelic frequencies of A and F were 0.795, 0.205 and 0.739, 0.261 in Xinong Saanen and Guanzhong population respectively. Xinong Saanen breed was in Hardy-Weinberg equilibrium, while Guanzhong breed deviated from Hardy-Weinberg equilibrium (p<0.05). The association of polymorphism with economic traits indicated that the goats with FF genotype have higher milk fat and total solid concentration than those with AA and AF genotypes (p<0.05).
CSN1S2; SNP; Dairy Goats; Milk; A and F Alleles
Twelve Chinese indigenous goat populations were genotyped for twenty-six microsatellite markers recommended by the EU Sheep and Goat Biodiversity Project. A total of 452 goats were tested. Seventeen of the 26 microsatellite markers used in this analysis had four or more alleles. The mean expected heterozygosity and the mean observed heterozygosity for the population varied from 0.611 to 0.784 and 0.602 to 0.783 respectively. The mean FST (0.105) demonstrated that about 89.5% of the total genetic variation was due to the genetic differentiation within each population. A phylogenetic tree based on the Nei (1978) standard genetic distance displayed a remarkable degree of consistency with their different geographical origins and their presumed migration throughout China. The correspondence analysis did not only distinguish population groups, but also confirmed the above results, classifying the important populations contributing to diversity. Additionally, some specific alleles were shown to be important in the construction of the population structure. The study analyzed the recent origins of these populations and contributed to the knowledge and genetic characterization of Chinese indigenous goat populations. In addition, the seventeen microsatellites recommended by the EU Sheep and Goat Biodiversity Project proved to be useful for the biodiversity studies in goat breeds.
genetic relationship; microsatellite; goat; Chinese indigenous population
China has numerous native domestic goat breeds, however, extensive studies are focused on the genetic diversity within the fewer breeds and limited regions, the population demograogic history and origin of Chinese goats are still unclear. The roles of geographical structure have not been analyzed in Chinese goat domestic process. In this study, the genetic relationships of Chinese indigenous goat populations were evaluated using 30 microsatellite markers.
Forty Chinese indigenous populations containing 2078 goats were sampled from different geographic regions of China. Moderate genetic diversity at the population level (HS of 0.644) and high population diversity at the species level (HT value of 0.737) were estimated. Significant moderate population differentiation was detected (FST value of 0.129). Significant excess homozygosity (FIS of 0.105) and recent population bottlenecks were detected in thirty-six populations. Neighbour-joining tree, principal components analysis and Bayesian clusters all revealed that Chinese goat populations could be subdivided into at least four genetic clusters: Southwest China, South China, Northwest China and East China. It was observed that the genetic diversity of Northern China goats was highest among these clusters. The results here suggested that the goat populations in Southwest China might be the earliest domestic goats in China.
Our results suggested that the current genetic structure of Chinese goats were resulted from the special geographical structure, especially in the Western China, and the Western goat populations had been separated by the geographic structure (Hengduan Mountains and Qinling Mountains-Huaihe River Line) into two clusters: the Southwest and Northwest. It also indicated that the current genetic structure was caused by the geographical origin mainly, in close accordance with the human’s migration history throughout China. This study provides a fundamental genetic profile for the conservation of these populations and better to understand the domestication process and origin of Chinese goats.
Iran is an area of particular interest for investigating goat diversity. Archaeological remains indicate early goat domestication (about 10 000 years ago) in the Iranian Zagros Mountains as well as in the high Euphrates valley and southeastern Anatolia. In addition, mitochondrial DNA data of domestic goats and wild ancestors (C. aegagrusor bezoar) suggest a pre-domestication management of wild populations in southern Zagros and central Iranian Plateau. In this study genetic diversity was assessed in seven Iranian native goat breeds, namely Markhoz, Najdi, Taleshi, Khalkhali, Naini, native Abadeh and Turki-Ghashghaei. A total of 317 animals were characterized using 14 microsatellite loci. Two Pakistani goat populations, Pahari and Teddy, were genotyped for comparison.
Iranian goats possess a remarkable genetic diversity (average expected heterozygosity of 0.671 across loci, 10.7 alleles per locus) mainly accounted for by the within-breed component (GST = 5.9%). Positive and highly significant FIS values in the Naini, Turki-Ghashghaei, Abadeh and Markhoz breeds indicate some level of inbreeding in these populations. Multivariate analyses cluster Iranian goats into northern, central and western groups, with the western breeds relatively distinct from the others. Pakistani breeds show some relationship with Iranian populations, even if their position is not consistent across analyses. Gene flow was higher within regions (west, north, central) compared to between regions but particularly low between the western and the other two regions, probably due to the isolating topography of the Zagros mountain range. The Turki-Ghashghaei, Najdi and Abadeh breeds are reared in geographic areas where mtDNA provided evidence of early domestication. These breeds are highly variable, located on basal short branches in the neighbor-joining tree, close to the origin of the principal component analysis plot and, although highly admixed, they are quite distinct from those reared on the western side of the Zagros mountain range.
These observations call for further investigation of the nuclear DNA diversity of these breeds within a much wider geographic context to confirm or re-discuss the current hypothesis (based on maternal lineage data) of an almost exclusive contribution of the eastern Anatolian bezoar to the domestic goat gene pool.
Kisspeptins are the peptide products of KISS1 gene, which operate via the G - protein-coupled receptor GPR54. These peptides have emerged as essential upstream regulators of neurons secreting gonadotropin-releasing hormone (GnRH), the major hypothalamic node for the stimulatory control of the hypothalamic–pituitary– gonadal (HPG) axis. The present study detected the polymorphisms of caprine KISS1 gene in three goat breeds and investigated the associations between these genetic markers and litter size.
Three goat breeds (n = 680) were used to detect single nucleotide polymorphisms (SNPs) in the coding regions with their intron–exon boundaries and the proximal flanking regions of KISS1 gene by DNA sequencing and PCR–RFLP. Eleven novel SNPs (g.384G>A, g.1147T>C, g.1417G>A, g.1428_1429delG, g.2124C>T, g.2270C>T, g.2489T>C, g.2510G>A, g.2540C>T, g.3864_3865delCA and g.3885_3886insACCCC) were identified. It was shown that Xinong Saanen and Guanzhong goat breeds were in Hardy-Weinberg disequilibrium at g.384G>A locus (P < 0.05). Both g.2510G>A and g.2540C>T loci were closely linked in Xinong Saanen (SN), Guanzhong (GZ) and Boer (BG) goat breeds (r2 > 0.33). The g.384G>A, g.2489T>C, g.2510G>A and g.2540C>T SNPs were associated with litter size (P<0.05). Individuals with AATTAATT combinative genotype of SN breed (SC) and TTAATT combinative genotype of BG breed (BC) had higher litter size than those with other combinative genotypes in average parity. The results extend the spectrum of genetic variation of the caprine KISS1 gene, which might contribute to goat genetic resources and breeding.
This study explored the genetic polymorphism of KISS1 gene, and indicated that four SNPs may play an important role in litter size. Their genetic mechanism of reproduction in goat breeds should be further investigated. The female goats with SC1 (AATTAATT) and BC7 (TTAATT) had higher litter size than those with other combinative genotypes in average parity and could be used for the development of new breeds of prolific goats. Further research on a large number of animals is required to confirm the link with increased prolificacy in goats.
Combinative genotype; SNP; PCR-RFLP; Candidate gene
To investigate the genetic diversity of seven Chinese indigenous meat goat breeds (Tibet goat, Guizhou white goat, Shannan white goat, Yichang white goat, Matou goat, Changjiangsanjiaozhou white goat and Anhui white goat), explain their genetic relationship and assess their integrity and degree of admixture, 302 individuals from these breeds and 42 Boer goats introduced from Africa as reference samples were genotyped for 11 microsatellite markers. Results indicated that the genetic diversity of Chinese indigenous meat goats was rich. The mean heterozygosity and the mean allelic richness (AR) for the 8 goat breeds varied from 0.697 to 0.738 and 6.21 to 7.35, respectively. Structure analysis showed that Tibet goat breed was genetically distinct and was the first to separate and the other Chinese goats were then divided into two sub-clusters: Shannan white goat and Yichang white goat in one cluster; and Guizhou white goat, Matou goat, Changjiangsanjiaozhou white goat and Anhui white goat in the other cluster. This grouping pattern was further supported by clustering analysis and Principal component analysis. These results may provide a scientific basis for the characteristization, conservation and utilization of Chinese meat goats.
Meat Goat; Genetic Diversity; Genetic Differentiation; Microsatellite
Introgression of Sahel livestock genes southwards in West Africa may be favoured by human activity and the increase of the duration of the dry seasons since the 1970’s. The aim of this study is to assess the gene flow patterns in Burkina Faso goat and to ascertain the most likely factors influencing geographic patterns of genetic variation in the Burkina Faso goat population.
A total of 520 goat were sampled in 23 different locations of Burkina Faso and genotyped for a set of 19 microsatellites. Data deposited in the Dryad repository: http://dx.doi.org/10.5061/dryad.41h46j37. Although overall differentiation is poor (FST = 0.067 ± 0.003), the goat population of Burkina Faso is far from being homogeneous. Barrier analysis pointed out the existence of: a) genetic discontinuities in the Central and Southeast Burkina Faso; and b) genetic differences within the goat sampled in the Sahel or the Sudan areas of Burkina Faso. Principal component analysis and admixture proportion scores were computed for each population sampled and used to construct interpolation maps. Furthermore, Population Graph analysis revealed that the Sahel and the Sudan environmental areas of Burkina Faso were connected through a significant number of extended edges, which would be consistent with the hypothesis of long-distance dispersal. Genetic variation of Burkina Faso goat followed a geographic-related pattern. This pattern of variation is likely to be related to the presence of vectors of African animal trypanosomosis. Partial Mantel test identified the present Northern limit of trypanosome vectors as the most significant landscape boundary influencing the genetic variability of Burkina Faso goat (p = 0.008). The contribution of Sahel goat genes to the goat populations in the Northern and Eastern parts of the Sudan-Sahel area of Burkina Faso was substantial. The presence of perennial streams explains the existence of trypanosome vectors. The South half of the Nakambé river (Southern Ouagadougou) and the Mouhoun river loop determined, respectively, the Eastern and Northern limits for the expansion of Sahelian goat genes. Furthermore, results from partial Mantel test suggest that the introgression of Sahelian goat genes into Djallonké goat using human-influenced genetic corridors has a limited influence when compared to the biological boundary defined by the northern limits for the distribution of the tsetse fly. However, the genetic differences found between the goat sampled in Bobo Dioulasso and the other populations located in the Sudan area of Burkina Faso may be explained by the broad goat trade favoured by the main road of the country.
The current analysis clearly suggests that genetic variation in Burkina Faso goat: a) follows a North to South clinal; and b) is affected by the distribution of the tsetse fly that imposes a limit to the Sahelian goat expansion due to their trypanosusceptibility. Here we show how extensive surveys on livestock populations can be useful to indirectly assess the consequences of climate change and human action in developing countries.
The major histocompatibility complex (MHC) is a crucial component of the vertebrate immune system and shows extremely high levels of genetic polymorphism. The extraordinary genetic variation is thought to be ancient polymorphisms maintained by balancing selection. However, introgression from related species was recently proposed as an additional mechanism. Here we provide evidence for introgression at the MHC in Alpine ibex (Capra ibex ibex). At a usually very polymorphic MHC exon involved in pathogen recognition (DRB exon 2), Alpine ibex carried only two alleles. We found that one of these DRB alleles is identical to a DRB allele of domestic goats (Capra aegagrus hircus). We sequenced 2489 bp of the coding and non-coding regions of the DRB gene and found that Alpine ibex homozygous for the goat-type DRB exon 2 allele showed nearly identical sequences (99.8%) to a breed of domestic goats. Using Sanger and RAD sequencing, microsatellite and SNP chip data, we show that the chromosomal region containing the goat-type DRB allele has a signature of recent introgression in Alpine ibex. A region of approximately 750 kb including the DRB locus showed high rates of heterozygosity in individuals carrying one copy of the goat-type DRB allele. These individuals shared SNP alleles both with domestic goats and other Alpine ibex. In a survey of four Alpine ibex populations, we found that the region surrounding the DRB allele shows strong linkage disequilibria, strong sequence clustering and low diversity among haplotypes carrying the goat-type allele. Introgression at the MHC is likely adaptive and introgression critically increased MHC DRB diversity in the genetically impoverished Alpine ibex. Our finding contradicts the long-standing view that genetic variability at the MHC is solely a consequence of ancient trans-species polymorphism. Introgression is likely an underappreciated source of genetic diversity at the MHC and other loci under balancing selection.
The major histocompatibility complex (MHC), a crucial component of the defense against pathogens, contains the most polymorphic functional genes in vertebrate genomes. The extraordinary genetic variation is generally considered to be ancient. We investigated whether a previously neglected mechanism, introgression from related species, provides an additional source of MHC variation. We show that introgression from domestic goat dramatically increased genetic variation at the MHC of Alpine ibex, a species that had nearly gone extinct during the 18th century, but has been restored to large numbers since. We show that Alpine ibex share one of only two alleles at a generally highly polymorphic MHC locus with domestic goats and that the chromosomal region containing the goat-type allele has a signature of recent introgression. Our finding contradicts the long-standing view that ancient trans-species polymorphism is the sole source of the extraordinary genetic variability at the MHC. Instead, we show that in Alpine ibex introgression generated genetic diversity at a MHC locus. Our study supports the view that loci favoring genetic polymorphism may be susceptible to adaptive introgression from related species and will encourage future research to identify unexpected signatures of introgression.
Assessing genetic biodiversity and population structure of minor breeds through the information provided by neutral molecular markers, allows determination of their extinction risk and to design strategies for their management and conservation. Analysis of microsatellite loci is known to be highly informative in the reconstruction of the historical processes underlying the evolution and differentiation of animal populations. Guadarrama goat is a threatened Spanish breed which actual census (2008) consists of 3057 females and 203 males distributed in 22 populations more or less isolated. The aim of this work is to study the genetic status of this breed through the analysis of molecular data from 10 microsatellites typed in historic and actual live animals.
The mean expected heterozygosity across loci within populations ranged from 0.62 to 0.77. Genetic differentiation measures were moderate, with a mean FST of 0.074, GST of 0.081 and RST of 0.085. Percentages of variation among and within populations were 7.5 and 92.5, respectively. Bayesian clustering analyses pointed out a population subdivision in 16 clusters, however, no correlation between geographical distances and genetic differences was found. Management factors such as the limited exchange of animals between farmers (estimated gene flow Nm = 3.08) mostly due to sanitary and social constraints could be the major causes affecting Guadarrama goat population subdivision.
Genetic diversity measures revealed a good status of biodiversity in the Guadarrama goat breed. Since diseases are the first cause affecting the census in this breed, population subdivision would be an advantage for its conservation. However, to maintain private alleles present at low frequencies in such small populations minimizing the inbreeding rate, it would necessitate some mating designs of animals carrying such alleles among populations. The systematic use of molecular markers will facilitate the comprehensive management of these populations, which in combination with the actual breeding program to increase milk yield, will constitute a good strategy to preserve the breed.
Jamunapari, a dairy goat breed of India, has been gradually declining in numbers in its home tract over the years. We have analysed genetic variation and population history in Jamunapari goats based on 17 microsatellite loci, 2 milk protein loci, mitochondrial hypervariable region I (HVRI) sequencing, and three Y-chromosomal gene sequencing. We used the mitochondrial DNA (mtDNA) mismatch distribution, microsatellite data, and bottleneck tests to infer the population history and demography. The mean number of alleles per locus was 9.0 indicating that the allelic variation was high in all the loci and the mean heterozygosity was 0.769 at nuclear loci. Although the population size is smaller than 8,000 individuals, the amount of variability both in terms of allelic richness and gene diversity was high in all the microsatellite loci except ILST 005. The gene diversity and effective number of alleles at milk protein loci were higher than the 10 other Indian goat breeds that they were compared to. Mismatch analysis was carried out and the analysis revealed that the population curve was unimodal indicating the expansion of population. The genetic diversity of Y-chromosome genes was low in the present study. The observed mean M ratio in the population was above the critical significance value (Mc) and close to one indicating that it has maintained a slowly changing population size. The mode-shift test did not detect any distortion of allele frequency and the heterozygosity excess method showed that there was no significant departure from mutation-drift equilibrium detected in the population. However, the effects of genetic bottlenecks were observed in some loci due to decreased heterozygosity and lower level of M ratio. There were two observed genetic subdivisions in the population supporting the observations of farmers in different areas. This base line information on genetic diversity, bottleneck analysis, and mismatch analysis was obtained to assist the conservation decision and management of the breed.
Genetic diversity at 13 equine microsatellite loci was compared in five endangered Spanish donkey breeds: Andaluza, Catalana, Mallorquina, Encartaciones and Zamorano-Leonesa. All of the equine microsatellites used in this study were amplified and were polymorphic in the domestic donkey breeds with the exception of HMS1, which was monomorphic, and ASB2, which failed to amplify. Allele number, frequency distributions and mean heterozygosities were very similar among the Spanish donkey breeds. The unbiased expected heterozygosity (HE) over all the populations varied between 0.637 and 0.684 in this study. The low GST value showed that only 3.6% of the diversity was between breeds (P < 0.01). Significant deviations from Hardy-Weinberg equilibrium were shown for a number of locus-population combinations, except HMS5 that showed agreement in all analysed populations. The cumulative exclusion probability (PE) was 0.999 in each breed, suggesting that the loci would be suitable for donkey parentage testing. The constructed dendrogram from the DA distance matrix showed little differentiation between Spanish breeds, but great differentiation between them and the Moroccan ass and also with the horse, used as an outgroup. These results confirm the potential use of equine microsatellite loci as a tool for genetic studies in domestic donkey populations, which could also be useful for conservation plans.
donkey; endangered breed; microsatellite; diversity; genetic variability
Current research on quantitative genetics has provided efficient guidelines for the sustainable management of selected populations: genetic gain is maximized while the loss of genetic diversity is maintained at a reasonable rate. However, actual selection schemes are complex, especially for large domestic species, and they have to take into account many operational constraints. This paper deals with the actual selection of dairy goats where the challenge is to optimize diffusion of buck semen on the field. Three objectives are considered simultaneously: i) natural service buck replacement (NSR); ii) goat replacement (GR); iii) semen distribution of young bucks to be progeny-tested. An appropriate optimization method is developed, which involves five analytical steps. Solutions are obtained by simulated annealing and the corresponding algorithms are presented in detail.
The whole procedure was tested on two French goat populations (Alpine and Saanen breeds) and the results presented in the abstract were based on the average of the two breeds. The procedure induced an immediate acceleration of genetic gain in comparison with the current annual genetic gain (0.15 genetic standard deviation unit), as shown by two facts. First, the genetic level of replacement natural service (NS) bucks was predicted, 1.5 years ahead at the moment of reproduction, to be equivalent to that of the progeny-tested bucks in service, born from the current breeding scheme. Second, the genetic level of replacement goats was much higher than that of their dams (0.86 unit), which represented 6 years of selection, although dams were only 3 years older than their replacement daughters. This improved genetic gain could be achieved while decreasing inbreeding coefficients substantially. Inbreeding coefficients (%) of NS bucks was lower than that of the progeny-tested bucks (-0.17). Goats were also less inbred than their dams (-0.67).
It was possible to account for complex operational constraints while developing goat selection schemes, both efficient and sustainable. Therefore, the recommended selection and mating decisions might receive attention from goat breeders using both AI and NS.
Among the European countries, Italy counts the largest number of local goat breeds. Thanks to the recent availability of a medium-density SNP (single nucleotide polymorphism) chip for goat, the genetic diversity of Italian goat populations was characterized by genotyping samples from 14 Italian goat breeds that originate from different geographical areas with more than 50 000 SNPs evenly distributed on the genome.
Analysis of the genotyping data revealed high levels of genetic polymorphism and an underlying North–south geographic pattern of genetic diversity that was highlighted by both the first dimension of the multi-dimensional scaling plot and the Neighbour network reconstruction. We observed a moderate and weak population structure in Northern and Central-Southern breeds, respectively, with pairwise FST values between breeds ranging from 0.013 to 0.164 and 7.49 % of the total variance assigned to the between-breed level. Only 2.11 % of the variance explained the clustering of breeds into geographical groups (Northern, Central and Southern Italy and Islands).
Our results indicate that the present-day genetic diversity of Italian goat populations was shaped by the combined effects of drift, presence or lack of gene flow and, to some extent, by the consequences of traditional management systems and recent demographic history. Our findings may constitute the starting point for the development of marker-assisted approaches, to better address future breeding and management policies in a species that is particularly relevant for the medium- and long-term sustainability of marginal regions.
Electronic supplementary material
The online version of this article (doi:10.1186/s12711-015-0140-6) contains supplementary material, which is available to authorized users.
Identification of true to breed type animal for conservation purpose is imperative. Breed dilution is one of the major problems in sustainability except cases of commercial crossbreeding under controlled condition. Breed descriptor has been developed to identify breed but such descriptors cover only “pure breed” or true to the breed type animals excluding undefined or admixture population. Moreover, in case of semen, ova, embryo and breed product, the breed cannot be identified due to lack of visible phenotypic descriptors. Advent of molecular markers like microsatellite and SNP have revolutionized breed identification from even small biological tissue or germplasm. Microsatellite DNA marker based breed assignments has been reported in various domestic animals. Such methods have limitations viz. non availability of allele data in public domain, thus each time all reference breed has to be genotyped which is neither logical nor economical. Even if such data is available but computational methods needs expertise of data analysis and interpretation.
We found Bayesian Networks as best classifier with highest accuracy of 98.7% using 51850 reference allele data generated by 25 microsatellite loci on 22 goat breed population of India. The FST values in the study were seen to be low ranging from 0.051 to 0.297 and overall genetic differentiation of 13.8%, suggesting more number of loci needed for higher accuracy. We report here world’s first model webserver for breed identification using microsatellite DNA markers freely accessible at http://cabin.iasri.res.in/gomi/.
Higher number of loci is required due to less differentiable population and large number of breeds taken in this study. This server will reduce the cost with computational ease. This methodology can be a model for various other domestic animal species as a valuable tool for conservation and breed improvement programmes.
Bayesian network; Breed; Goat; Microsatellite; Prediction; Webserver
Genetic diversity, introgression and relationships were studied in 521 individuals from 9 African Bos indicus and 3 Bos taurus cattle breeds in Cameroon and Nigeria using genotype information on 28 markers (16 microsatellite, 7 milk protein and 5 blood protein markers). The genotypes of 13 of the 16 microsatellite markers studied on three European (German Angus, German Simmental and German Yellow) and two Indian (Nelore and Ongole) breeds were used to assess the relationships between them and the African breeds. Diversity levels at microsatellite loci were higher in the zebu than in the taurine breeds and were generally similar for protein loci in the breeds in each group. Microsatellite allelic distribution displayed groups of alleles specific to the Indian zebu, African taurine and European taurine. The level of the Indian zebu genetic admixture proportions in the African zebus was higher than the African taurine and European taurine admixture proportions, and ranged from 58.1% to 74.0%. The African taurine breed, Muturu was free of Indian zebu genes while its counter Namchi was highly introgressed (30.2%). Phylogenic reconstruction and principal component analysis indicate close relationships among the zebu breeds in Cameroon and Nigeria and a large genetic divergence between the main cattle groups – African taurine, European taurine and Indian zebu, and a central position for the African zebus. The study presents the first comprehensive information on the hybrid composition of the individual cattle breeds of Cameroon and Nigeria and the genetic relationships existing among them and other breeds outside of Africa. Strong evidence supporting separate domestication events for the Bos species is also provided.
genetic diversity; introgression; relationship; cattle
This study investigated the geographic and pairwise distances of nine Chinese local Cashmere goat populations through the analysis of 20 microsatellite DNA markers. Fluorescence PCR was used to identify the markers, which were selected based on their significance as identified by the Food and Agriculture Organization of the United Nations (FAO) and the International Society for Animal Genetics (ISAG). In total, 206 alleles were detected; the average allele number was 10.30; the polymorphism information content of loci ranged from 0.5213 to 0.7582; the number of effective alleles ranged from 4.0484 to 4.6178; the observed heterozygosity was from 0.5023 to 0.5602 for the practical sample; the expected heterozygosity ranged from 0.5783 to 0.6464; and Allelic richness ranged from 4.7551 to 8.0693. These results indicated that Chinese Cashmere goat populations exhibited rich genetic diversity. Further, the Wright’s F-statistics of subpopulation within total (FST) was 0.1184; the genetic differentiation coefficient (GST) was 0.0940; and the average gene flow (Nm) was 2.0415. All pairwise FST values among the populations were highly significant (p<0.01 or p<0.001), suggesting that the populations studied should all be considered to be separate breeds. Finally, the clustering analysis divided the Chinese Cashmere goat populations into at least four clusters, with the Hexi and Yashan goat populations alone in one cluster. These results have provided useful, practical, and important information for the future of Chinese Cashmere goat breeding.
Cashmere Goat; Geographic Distance; Pairwise Distance; Genetic Diversity
Basic understanding of linkage disequilibrium (LD) and population structure, as well as the consistency of gametic phase across breeds is crucial for genome-wide association studies and successful implementation of genomic selection. However, it is still limited in goats. Therefore, the objectives of this research were: (i) to estimate genome-wide levels of LD in goat breeds using data generated with the Illumina Goat SNP50 BeadChip; (ii) to study the consistency of gametic phase across breeds in order to evaluate the possible use of a multi-breed training population for genomic selection and (iii) develop insights concerning the population history of goat breeds.
Average r2 between adjacent SNP pairs ranged from 0.28 to 0.11 for Boer and Rangeland populations. At the average distance between adjacent SNPs in the current 50 k SNP panel (~0.06 Mb), the breeds LaMancha, Nubian, Toggenburg and Boer exceeded or approached the level of linkage disequilibrium that is useful (r2 > 0.2) for genomic predictions. In all breeds LD decayed rapidly with increasing inter-marker distance. The estimated correlations for all the breed pairs, except Canadian and Australian Boer populations, were lower than 0.70 for all marker distances greater than 0.02 Mb. These results are not high enough to encourage the pooling of breeds in a single training population for genomic selection. The admixture analysis shows that some breeds have distinct genotypes based on SNP50 genotypes, such as the Boer, Cashmere and Nubian populations. The other groups share higher genome proportions with each other, indicating higher admixture and a more diverse genetic composition.
This work presents results of a diverse collection of breeds, which are of great interest for the implementation of genomic selection in goats. The LD results indicate that, with a large enough training population, genomic selection could potentially be implemented within breed with the current 50 k panel, but some breeds might benefit from a denser panel. For multi-breed genomic evaluation, a denser SNP panel also seems to be required.
Electronic supplementary material
The online version of this article (doi:10.1186/s12863-015-0220-1) contains supplementary material, which is available to authorized users.
Effective population size; Genomic selection; Goat breeds; Goat 50 k panel; LD
A set of eleven pig breeds originating from six European countries, and including a small sample of wild pigs, was chosen for this study of genetic diversity. Diversity was evaluated on the basis of 18 microsatellite markers typed over a total of 483 DNA samples collected. Average breed heterozygosity varied from 0.35 to 0.60. Genotypic frequencies generally agreed with Hardy-Weinberg expectations, apart from the German Landrace and Schwäbisch-Hällisches breeds, which showed significantly reduced heterozygosity. Breed differentiation was significant as shown by the high among-breed fixation index (overall FST = 0.27), and confirmed by the clustering based on the genetic distances between individuals, which grouped essentially all individuals in 11 clusters corresponding to the 11 breeds. The genetic distances between breeds were first used to construct phylogenetic trees. The trees indicated that a genetic drift model might explain the divergence of the two German breeds, but no reliable phylogeny could be inferred among the remaining breeds. The same distances were also used to measure the global diversity of the set of breeds considered, and to evaluate the marginal loss of diversity attached to each breed. In that respect, the French Basque breed appeared to be the most "unique" in the set considered. This study, which remains to be extended to a larger set of European breeds, indicates that using genetic distances between breeds of farm animals in a classical taxonomic approach may not give clear resolution, but points to their usefulness in a prospective evaluation of diversity.
genetic diversity; molecular marker; conservation; pig; European breed
Domestic goats (Capra hircus) have been selected to play an essential role in agricultural production systems, since being domesticated from their wild progenitor, bezoar (Capra aegagrus). A detailed understanding of the genetic consequences imparted by the domestication process remains a key goal of evolutionary genomics.
We constructed the reference genome of bezoar and sequenced representative breeds of domestic goats to search for genomic changes that likely have accompanied goat domestication and breed formation. Thirteen copy number variation genes associated with coat color were identified in domestic goats, among which ASIP gene duplication contributes to the generation of light coat-color phenotype in domestic goats. Analysis of rapidly evolving genes identified genic changes underlying behavior-related traits, immune response and production-related traits.
Based on the comparison studies of copy number variation genes and rapidly evolving genes between wild and domestic goat, our findings and methodology shed light on the genetic mechanism of animal domestication and will facilitate future goat breeding.
Electronic supplementary material
The online version of this article (doi:10.1186/s12864-015-1606-1) contains supplementary material, which is available to authorized users.
Wild goat genome; Animal domestication; Artificial selection; Coat color evolution; Behavioral evolution
SNP data of goats of three Mediterranean countries were used for population studies and reconstruction of geographical patterning. 496 individuals belonging to Italian, Albanian and Greek breeds were genotyped to assess the basic population parameters.
A total of 26 SNPs were used, for a total of 12,896 genotypes assayed. Statistical analysis revealed that breeds are not so similar in terms of genetic variability, as reported in studies performed using different markers. The Mantel test showed a strongly significant correlation between genetic and geographic distance. Also, PCA analysis revealed that breeds are grouped according to geographical origin, with the exception of the Greek Skopelos breed.
Our data point out that the use of SNP markers to analyze a wider breed sample could help in understanding the recent evolutionary history of domestic goats. We found correlation between genetic diversity and geographic distance. Also PCA analysis shows that the breeds are well differentiated, with good correspondence to geographical locations, thus confirming the correlation between geographical and genetic distances. This suggests that migration history of the species played a pivotal role in the present-day structure of the breeds and a scenario in which coastal routes were easier for migrating in comparison with inland routes. A westward coastal route to Italy through Greece could have led to gene flow along the Northern Mediterranean.
West African Dwarf (WAD) goats are extremely important in the rural village economy of West Africa, but still little is known about their biology, ecology and capacity to cope with gastrointestinal nematode (GIN) infections. Here, we summarise the history of this breed and explain its economic importance in rural West Africa. We review recent work showing that Nigerian WAD goats are highly trypanotolerant and resist infections with Haemonchus contortus more effectively than other breeds of domestic goat (haemonchotolerance). We believe that haemonchotolerance is largely responsible for the generally low level GIN infections and absence of clinical haemonchosis in WADs under field conditions, and has contributed to the relatively successful and sustainable, anthelmintics-free, small-scale system of goat husbandry in Nigeria’s humid zone, and is immunologically based and genetically controlled. If haemonchotolerance can be shown to be genetically controlled, it should be possible to exploit the underlying genes to improve GIN resistance among productive fibre and milk producing breeds of goats, most of which are highly susceptible to nematode infections. Genetic resistance to GIN and trypanosome infections would obviate the need for expensive chemotherapy, mostly unaffordable to small-holder farmers in Africa, and a significant cost of goat husbandry in more developed countries. Either introgression of resistance alleles into susceptible breeds by conventional breeding, or transgenesis could be used to develop novel parasite-resistant, but highly productive breeds, or to improve the resistance of existing breeds, benefitting the local West African rural economy as well as global caprine livestock agriculture.
Nematoda; Haemonchotolerance; Trypanotolerance; West African Dwarf goats; Genetics; Immune responses
Indian agriculture is an economic symbiosis of crop and livestock production with cattle as the foundation. Sadly, the population of indigenous cattle (Bos indicus) is declining (8.94 % in last decade) and needs immediate scientific management. Genetic characterization is the first step in the development of proper management strategies for preserving genetic diversity and preventing undesirable loss of alleles. Thus, in this study we investigated genetic diversity and relationship among eleven Indian cattle breeds using 21 microsatellite markers and mitochondrial D loop sequence.
The analysis of autosomal DNA was performed on 508 cattle which exhibited sufficient genetic diversity across all the breeds. Estimates of mean allele number and observed heterozygosity across all loci and population were 8.784 ± 0.25 and 0.653 ± 0.014, respectively. Differences among breeds accounted for 13.3 % of total genetic variability. Despite high genetic diversity, significant inbreeding was also observed within eight populations. Genetic distances and cluster analysis showed a close relationship between breeds according to proximity in geographic distribution. The genetic distance, STRUCTURE and Principal Coordinate Analysis concluded that the Southern Indian Ongole cattle are the most distinct among the investigated cattle populations. Sequencing of hypervariable mitochondrial DNA region on a subset of 170 cattle revealed sixty haplotypes with haplotypic diversity of 0.90240, nucleotide diversity of 0.02688 and average number of nucleotide differences as 6.07407. Two major star clusters for haplotypes indicated population expansion for Indian cattle.
Nuclear and mitochondrial genomes show a similar pattern of genetic variability and genetic differentiation. Various analyses concluded that the Southern breed ‘Ongole’ was distinct from breeds of Northern/ Central India. Overall these results provide basic information about genetic diversity and structure of Indian cattle which should have implications for management and conservation of indicine cattle diversity.
Electronic supplementary material
The online version of this article (doi:10.1186/s12863-015-0221-0) contains supplementary material, which is available to authorized users.
Conservation; Diversity; Genetic relationship; Indian cattle; Microsatellite markers; Mitochondrial DNA; Population structure