Genetic Diversity of Five Fecundity Related Microsatellites in Six Goat Breeds
Genetic Diversity of Five Fecundity Related Microsatellites in Six Goat Breeds
Guang-Xin E1, Yong-Ju Zhao1, Zhong-Yuan Zhang1, Yan-Guo Han1, Wei-Wei Ni1, Ming-Xing Chu2, Yue-Hui Ma2, Jian-Ning He3, Huai-Zhi Jiang4, Jia-Hua Zhang1, Zhong-Quan Zhao1, Yan Zeng1, Ri-Su Na1, Gao-Fu Wang5, Hang-Xing Ren5 and Yong-Fu Huang1*
1College of Animal Science and Technology, Chongqing Key Laboratory of Forage and Herbivore, Chongqing Engineering Research Centre for Herbivores Resource Protection and Utilization, Southwest University, Chongqing 400716, China
2Institute of Animal Science, Chinese Academy of Agricultural Sciences (CAAS), Beijing, 100193, China
3College of Animal Science and Technology, Qingdao Agricultural University, Qingdao 266109, China
4College of Animal Science and Technology, Jilin Agricultural University, Changchun, Jilin, China
5Chongqing Academy of Animal Sciences, Chongqing 402460, China
Fecundity is an important trait for the reproduction of economically domestic animals. Some microsatellites related to the fecundity of domestic animals have been identified. In the current study, we investigated five fecundity related microsatellites, including BM1329, BM143, LSCV043, OarHH35, and OarHH55, which had been identified in sheep. We estimated their diversity and population structure in three high fecundity goats (Jianzhou big ear goat, Jining grey goat, Dazu black goat) and three low fecundity goats (Boer goat, Nubian goat, Inner Mongolia Cashmere goat). The results revealed a high diversity of those markers in the goat populations. The population structure did not diverge significantly between high fecundity goats and low fecundity populations, but it followed the agro-ecological zone and management history. This indicated that the frequency of those five microsatellites did not increase during the reproduction of those six breeds. Therefore, these five fecundity related microsatellites are not appropriate biomarkers for goat breeding.
Received 29 September 2016
Revised 27 October 2016
Accepted 30 October 2016
Available online 09 February 2017
Y-FH and G-XE participated in the experimental design and wrote the manuscript; Y-JZ, M-XC, Y-HM, J-HZ and Z-YZ performed the lab work and carried out bioinformatics data analysis; H-ZJ, G-FW, H-XR, W-WN, YZ, R-SN, Z-QZ and J-NH participated in the sampling; Y-FH and Y-JZ have read and approved the final manuscript.
* Corresponding author: H67738337@swu.edu.cn
0030-9923/2017/0002-0721 $ 9.00/0
Copyright 2017 Zoological Society of Pakistan
Fecundity is important for some economically domestic animals. Recently, research using new technology has focused on this trait in an attempt to understand its mechanism (; ; ) and to develop related biomarkers for further animal breeding (; ).
The domestic goat (Capra hircus) has been an economically and culturally important farm animal species since its domestication. A wide array of goat breeds with abundant phenotypic diversity exists as a result of domestication and selection. Currently, commercial lines and industrialized livestock production systems have spread globally. Indigenous populations have low growth speed but may harbor specific features, including the number of offspring, as a result of adaptation to their environment.
Many genes contributing to reproduction have been identified in sheep, including BMP15 (; ). Some microsatellites related to the number of offspring, including OarHH35 (), BM1329 (; ), OarHH55 (), OarAE101 (; ), BM143 (; ; ), and LSCV043 (; ) have been investigated in sheep and in some goat populations. However, some studies reported these marker were not related to litter size completely in some goat population, including LSCV043 (), OarHH35 () and OarHH55 (; ).
Therefore, investigating the distribution of these microsatellite polymorphisms in different goat populations or breeds may be helpful for improving the breeding of goats, particularly with respect to fecundity. In this study, we analyzed the diversity of five microsatellite markers in four Chinese indigenous breeds and two imported breeds (Boer and Nubian) and assessed the utility of these loci as biomarkers for goat fecundity and breeding.
Material and methods
We genotyped 191 individuals from 6 breeds of goats from different geographic locations, including Boer goats (30, East Longitude (EL):144.963, North Latitude (NL):-37.813 ), Nubian goats (29, EL: 142.573, NL:-32.486), Dazu black goats (28, EL:105.721, NL:29.707), Jianzhou big ear goats (30, EL:104.546, NL:30.410), Jining grey goats (55, EL:116.587, NL:35.415) and Inner Mongolia Cashmere goats (19, EL:105.728, NL:38.851). Of these breeds, Jianzhou big ear, Jining grey, and Dazu black show high fecundity (Litter size: 252% ~ 283%); whereas, Boer, Nubian, and Inner Mongolia Cashmere show low fecundity (Litter size: 105% ~ 193%) (). Blood samples were taken from the ulnar vein and stored in vacuum tubes containing EDTA 1~2mg/ml blood at -20 °C before examination. Genomic DNA was extracted from whole blood as described by .
The animals were genotyped by PCR amplification of the BM1329, BM143, LSCV043, OarHH35 and OarHH55 markers (). The amplification conditions consisted of an initial denaturation for 1 min at 94 °C, followed by 35 cycles at 92 °C for 45 s, 51-63 °C for 45 s and 72 °C for 45 s, with a 1 min extension at 72 °C. Approximately, 1-2 μL of each PCR product was diluted with 10 μL of autoclaved distilled water for use in DNA genotyping. A 2 μL aliquot of the diluted PCR product was mixed with 7.75 μL Hi DiTM Formamide and 0.25 μL Gene Scan-500 LIZTM. The mixtures were heated at 94 °C for 5 min and then immediately chilled on ice for 2 min. Genotyping was carried out on a 3130 xl Genetic Analyzer (AB Applied Bio systems, US).
The genetic diversity expected (HE), the observed heterozygosity (Ho), the mean number of alleles (NA) and the polymorphism information content (PIC) were estimated from the allele frequencies using the Microsatellite Toolkit (). For each locus-population combination of the global data set and population groupings, we used Fisher’s exact test with Bonferroni correction to test possible deviations from the Hardy-Weinberg equilibrium (HWE) using GENEPOP 3.4 (). The pairwise difference (FST) was estimated using Arlequin software 184.108.40.206 (). The phylogenetic neighbor-joining tree was derived from Reynold’s genetic distance using the PHYLIP () software package. The Bayesian clustering algorithm implemented in STRUCTURE v2.3.3 () with 50,000 iterations following a burn-in of 100,000 Markov Chain Monte Carlo replications in an admixture model was used to assess the genotypic composition of the genetic backgrounds of the populations analyzed and the proportions of mixed ancestry. Six independent simulations from K=1 to K=6 were performed to identify the most probable clustering solution by examining the modal distribution of Delta K (). Graphical representations of these statistics were obtained with STRUCTURE HARVESTER v0.68 ().
In this study, capillary electrophoresis was used for genotyping of these five fecundity related microsatellites (). In total, 115 alleles of the 5 fecundity related microsatellite loci were found in the six goat breeds. Across populations, an average of 23 alleles per locus was observed, ranging from 12 in BM1329 to 32 in OarHH55 (). The average polymorphism information content across loci was 0.740, ranging from 0.6199 (BM143) to 0.8357 (OarHH55) (). The mean expected and observed heterozygote frequency within a locus across populations were 0.779 (0.663 in BM143 to 0.8649 in OarHH55) and 0.590 (0.2965 in LSCV043 to 1.0000 in OarHH55), respectively (). This revealed that the microsatellite loci in this study showed large frequencies of polymorphism among the six breeds.
Across loci, the NA ranged from 8.40±4.45 in Boer goats to 14.00±5.24 in Jining grey goats (). The HO within each population ranged from 0.467±0.041 in Boer goats to 0.700±0.039 in Dazu black goats. The highest HE (0.824±0.045) was found within the Jining grey goat population, and the lowest (0.656±0.075) in the Boer population. However, for most populations, the HE and HO were inconsistent in this study ().
|Dazu black goat||0.700±0.039||0.788±0.038||8.60±2.61||3||4|
|Jining Gregy goat||0.596±0.030||0.824±0.045||14.00±5.24||4||19|
|Inner mongolia cashmere goat||0.558±0.051||0.789±0.063||8.60±3.21||3||9|
|Jianzhou big ear goat||0.587±0.040||
Note, Pa is number of private allele; dHWE is number of populations deviated from Hardy-Weinberg equilibrium.
Note, number of population; BE, Boer goat; DZ, Dazu black goat; JN, Jining Gregy goat; NM, Inner mongolia cashmere goat; NB, Nubia goat; JZ, Jianzhou big ear goat.
In total, 48 separate alleles were distributed across 6 populations and 5 markers. The pair-wise differences between populations (FST) are shown in and , and the consistency with the phylogenetic network based on Reynolds genetic distance is shown in .
Goats are important to the subsistence needs of many human populations because they can provide abundant, regular supplies of meat, milk, hides and cashmere. Therefore, it is important to improve the economic traits of goats, including reproductive traits. Unfortunately, the genetic mechanisms contributing to caprine fertility are still unclear.
Marker-assisted selection (MAS) can help to improve traits, including fecundity related phenotypes, in domestic animals (). We analyzed the polymorphisms of genes that have been confirmed as major genes of fecundity in sheep to ascertain whether they were responsible for fecundity in goats.
In previous studies, it was found that the 104 bp, 106 bp and 110 bp alleles of BMP134 contributed to the number of offspring born to Xiangdong black goats (). The 112 bp and the 175 bp alleles of OarHH35 made positive and negative contributions, respectively, to the number of offspring in Lezi black goats (). In addition, the 140 bp and 120 bp alleles of LSCV043 had positive effects on the number of offspring (). However, if those microsatellite loci were related to the fecundity of goats, they should have undergone strong selection during traditional breeding. Therefore, the diversity and polymorphism of these loci should diverge significantly in high fecundity goats and goats with lower fecundity. In addition, they should be independent of drift and natural selection due to geography.
In this study, the results obtained for HE, HO, PIC and NA for the five microsatellites are consistent with other studies of local Chinese goat breeds and commercial populations, including Boer and Nubian (; ). These results indicate that the five fecundity related microsatellite loci analyzed in this study are highly polymorphic and diverse in the different goat breeds.
The number of loci with exclusive alleles was highest in Jining grey goats (14.00±5.24) and lowest in Boer goats (8.40±4.45). The expected heterozygosity with the highest mean was found in the Jining grey goats (0.824±0.045) and the lowest was in the Boer goats (0.656±0.075). The observed heterozygosity per population ranged between 0.700±0.039 (Dazu black) and 0.467±0.041 (Boer). Interestingly, high inconsistent between HE and HO of most populations indicate that those locus were not under HWE. As, the proportion of loci not in HWE was highest in the Nubian goats and lowest in the Dazu black goats and the Inner Mongolia Cashmere goats. However, it could be caused by small sample size of each population in this study.
Private alleles (Pas) were observed in all the populations studied ranging from 3 (Nubian) to 19 (Jining grey). The frequency of some private alleles within certain populations was possible because they had different domestic histories. Overall, these results indicate that, although the within-population expected and observed heterozygosities were not widely different, the deficit found in within-population heterozygosity was different among the populations. In addition, the large divergence of private alleles among populations indicated that those six breeds should reveal an interesting population structure.
In the pair-wise difference (FST) analysis, the largest difference was between Boer and Inner Mongolia Cashmere goats (0.19681), and the smallest was between Jianzhou Big ear and Nubian goats (0.00802), which indicate that the population relationship corresponded to their management history. Jianzhou Big ear goats are neutral hybrids between the Nubian goat and the indigenous breed in Southwest China, and they date from the Second World War. Two decades ago, these ecotype individuals were combined, and a specific breed was constructed by artificial breeding for meat production (). In addition, the consensus neighbor-joining phylogenetic tree of six goat breeds revealed that the clear genetic separation between breeds followed the geographic location of the sampling, including the Jianzhou Big ear goat located between the Nubian and Chinese native goats. In addition, the Boer goat, which was sampled from Oceanica, was a distant genetic variant of Chinese native goats but resembled Nubian goats. An analysis by STRUCTURE (K=3) indicated a similar genetic population pattern of clustering.
In the current study, the clustering of breeds based on the five fecundity related microsatellites was consistent with the agro-ecological zone and management history. This phenomenon could be explained by drift and natural selection of these loci. Although the clustering could result from an artificial process, the associated beneficial genotypes of these loci have not been fixed in a specific goat breed.
In short, these five fecundity related microsatellites were polymorphic, with high numbers of alleles in the six breeds of this study. In addition, the population structure was not significantly divergent between the high fecundity goats and the lower fecundity goats. Instead, the population structure reflected the agro-ecological zone and management history of the breeds. The high fecundity phenotype of the five microsatellites did not increase in those six breeds as a result of breeding. Therefore, the five microsatellites cannot explain the mechanism of fecundity in goats or provide useful biomarkers for breed selection in goats.
Co-authors declare no conflict of interest. This work was supported by the Fundamental Research Funds for the Central Universities (No. SWU114023), the 2013 Innovation Team Building Program in Chongqing universities (KJTD201334).
Conflict of Interest Statement
The authors declare that there is no conflict of interests regarding the publication of the manuscript.
Felsenstein, J., 2005. PHYLIP (Phylogeny Inference Package) version 3.6 distributed by the author. Department of Genome Sciences, University of Washington, Seattle. Available: . Accessed 24 March 2014