Microsatellite analysis revealed genetic diversity and population structure among Chinese cashmere goats

Most cashmere goats are found in northern China and Mongolia. They are regarded as precious resources for their production of high quality natural fibre for the textile industry. It was the first time that the genetic diversity and population structure of nine Chinese cashmere populations has been assessed using 14 ISAG/FAO microsatellite markers. In addition, two Iranian populations and one West African goat population were genotyped for comparison. Results indicated that the genetic diversity of Chinese cashmere goats was rich, but less than those of the Iranian goat populations. All pairwise F_ST values between the Chinese cashmere goat populations reached a highly significant level (P < 0.001), suggesting that they should all be considered as separate breeds. Finally, clustering analysis divided Chinese cashmere goats into at least two clusters, with the Tibetan Hegu goats alone in one cluster. An extensive admixture was detected among the Chinese goat breeds (except the Hegu), which have important implications for breeding management.     

Variation and expression of KAP9.2 gene affecting cashmere trait in goats

Keratin-associated proteins 9.2 (KAP9.2) gene encodes one of the ultra high sulfur KAPs. Variation in KAP genes may affect the structure of KAPs and hence cashmere characteristics. In order to test the association between the polymorphism of KAP9.2 gene and cashmere trait, DNA sequencing was used to detect a novel C/T polymorphism of KAP9.2 gene from a genomic DNA pool. The mutation could be recognized by Pst I restriction enzyme. To Shanbei white cashmere goat, Inner Mongolia white cashmere goat and Guanzhong dairy goat, the genotypic frequencies of TT, TC and CC from total 1,236 animals were as follows: 0.047, 0.519 and 0.434; 0.180, 0.592 and 0.228; 0.431, 0.544 and 0.025. The allelic frequencies of T and C were 0.307 and 0.693; 0.476 and 0.524; 0.703 and 0.297, respectively, in breeds mentioned above. The frequency of C allele between cashmere and dairy goat was significant (P?<?0.01). To provide support for the hypothesis that SNP 586 was responsible for KAP9.2 expression, quantitative real-time PCR analysis revealed that the expression level of KAP9.2 was reduced in individuals bearing genotype CC compared with TT individuals, suggesting that C was the nucleotide causing decreased expression of KAP9.2 or was in linkage disequilibrium with the causative SNP. The 586C/T SNP found in this study might control translation or stability of KAP9.2 mRNA, which would be beneficial for marker assistant selection in cashmere goat breeding.     

A study on the Changthangi pashmina and the Bakerwali goat breeds in Kashmir: I. Analysis of blood protein polymorphisms and genetic variability within and between the populations

Blood samples of 483 Pashmina goats and 392 Bakerwali goats were taken from the Ladakh and Jammu provinces, respectively, for characterisation of the breeds by polymorphic enzymes and proteins. Furthermore, the sex, age, body weight and hematocrit of both breeds and the pashmina yield of Pashmina goats were recorded. Polymorphisms of 12 enzymes and proteins (albumin (Al), alkaline phosphatase (Alp), amylase (Amy), NADH-diaphorase 1 (Dia1), vitamin-D-binding protein (Gc), haemoglobin (Hb), hemopexin (Hpx), nucleoside phosphorylase (Np), malic enzyme (ME), phosphohexose-isomerase (PHI), transferrin (Tf), X-protein (X)) in blood plasma and hemolysate were determined using gel electrophoresis. Out of 12 protein systems, 10 were found to be polymorphic. In four systems (Al, Amy, Dia1, Hpx) new phenotypes were detected. To estimate the genetic variability within breeds, the degrees of heterozygosity, deviations from the Hardy–Weinberg-equilibrium (HWE) and F_IS-, F_ST-, and F_IT-values were estimated. In comparison with literature data both breeds show slightly higher variability than other goat breeds, with degrees of heterozygosity ranging from 15 to 25% and 24 to 26% in the Pashmina and Bakerwali goat populations, respectively, and the percentage of polymorphic loci ranging from 50 to 67%. On the other hand, a decrease of variability at some loci can be observed in both breeds. Deviations from HWE in combination with a deficiency of heterozygote genotypes were observed in all sub-populations apart from the Pashmina sub-population ‘Likir’. Genetic differences between the goat breeds could be quantified through calculation of Nei’s genetic distances ranging from 0.002 to 0.080. With exception of the Pashmina sub-population ‘Likir’, lower distance values were found between sub-populations within the respective breeds.     

微卫星DNA标记在绒山羊群体中的初步研究

利用7个微卫星标记对4个绒山羊品种共计18个个体的遗传多样性进行了研究。计算了有效等位基因数、遗传杂合度、遗传距离等,分析了群体相关的遗传变异。结果表明：辽宁多绒山羊的有效等位基因数最大,杂合度最高,而辽宁绒山羊的有效等位基因数最小,杂合度最低;奈氏遗传距离说明：库布旗杂种绒山羊和辽宁多绒山羊的亲缘关系最近,而和阿尔巴斯绒山羊的亲缘关系最远。     

The PCR-SSCP and DNA sequencing methods detecting a large deletion mutation at KAP6.2 locus in the cashmere goat

It is known that keratin-associated proteins 6.2 gene (KAP6.2) is an important structural gene responsible for the cashmere. So in this study, the polymorphism of KAP6.2 gene was firstly detected by PCR-SSCP and DNA sequencing methods in 1052 cashmere goat samples. The results showed that two alleles were detected at the locus, which were named as allele O and X. Frequencies of KAP6.2-O allele in Inner Mongolia White cashmere (n = 847) and Shaanbei White cashmere goat breeds (n = 205) were 1.000 and 0.856, respectively. Inner Mongolia White cashmere goat was monomorphic at this locus, while Shaanbei White cashmere goat was at low polymorphic level. Therefore, the allele KAP6.2-X was considered as the breed characterization of Shaanbei White cashmere goat at KAP6.2 locus. Further sequencing analysis showed that a 24-bp deletion mutation was described for the first time in Shaanbei, while Inner Mongolia White cashmere goat did not have the deletion. The study indicated that deletion mutation was predicted as a possible cause for the multiple pattern cashmere in Shaanbei White cashmere goat.     

利用微卫星DNA标记研究绒山羊群体遗传多样性

利用7个微卫星标记对4个绒山羊品种共计18个个体的遗传多样性进行了分析和研究。计算了有效等位基因数、遗传杂合度、遗传距离等,分析了群体相关的遗传变异。结果表明,辽宁多绒山羊的有效等位基因数最大,杂合度最高;而辽宁绒山羊的有效等位基因数最小,杂合度最低。奈氏遗传距离表明,库布齐杂种绒山羊和辽宁多绒山羊的亲缘关系最近,而和阿尔巴斯绒山羊的亲缘关系最远。     

Genome‐wide population structure and admixture analysis reveals weak differentiation among Ugandan goat breeds

Uganda has a large population of goats, predominantly from indigenous breeds reared in diverse production systems, whose existence is threatened by crossbreeding with exotic Boer goats. Knowledge about the genetic characteristics and relationships among these Ugandan goat breeds and the potential admixture with Boer goats is still limited. Using a medium‐density single nucleotide polymorphism (SNP) panel, we assessed the genetic diversity, population structure and admixture in six goat breeds in Uganda: Boer, Karamojong, Kigezi, Mubende, Small East African and Sebei. All the animals had genotypes for about 46 105 SNPs after quality control. We found high proportions of polymorphic SNPs ranging from 0.885 (Kigezi) to 0.928 (Sebei). The overall mean observed (H_O) and expected (H_E) heterozygosity across breeds was 0.355 ± 0.147 and 0.384 ± 0.143 respectively. Principal components, genetic distances and admixture analyses revealed weak population sub‐structuring among the breeds. Principal components separated Kigezi and weakly Small East African from other indigenous goats. Sebei and Karamojong were tightly entangled together, whereas Mubende occupied a more central position with high admixture from all other local breeds. The Boer breed showed a unique cluster from the Ugandan indigenous goat breeds. The results reflect common ancestry but also some level of geographical differentiation. admixture and f₄ statistics revealed gene flow from Boer and varying levels of genetic admixture among the breeds. Generally, moderate to high levels of genetic variability were observed. Our findings provide useful insights into maintaining genetic diversity and designing appropriate breeding programs to exploit within‐breed diversity and heterozygote advantage in crossbreeding schemes.     

A genome-wide perspective on the diversity and selection signatures in indigenous goats using 53 K single nucleotide polymorphism array

Tibetan goats, Taihang goats, Jining grey goats, and Meigu goats are the representative indigenous goats in China, found in Qinghai-Tibet Plateau, Western pastoral area, Northern and Southern agricultural regions. Very few studies have conducted a comprehensive analysis of the genomic diversity and selection of these breeds. We genotyped 96 unrelated individuals, using goat 53 K Illumina BeadChip array, of the following goat breeds: Tibetan (TG), Taihang (THG), Jining grey (JGG), and Meigu (MGG). A total of 45 951 single nucleotide polymorphisms were filtered to estimate the genetic diversity and selection signatures. All breeds had a high proportion (over 95%) of polymorphic loci. The observed and excepted heterozygosity ranged from 0.338 (MGG) to 0.402 (JGG) and 0.339 (MGG) to 0.395 (JGG), respectively. Clustering analysis displayed a genetically distinct lineage for each breed, and their Fst were greater than 0.25, indicating that they had a higher genetic differentiation between groups. Furthermore, effective population size reduced in all four populations, indicating a loss of genetic diversity. In addition, runs of homozygosity were mainly distributed in 5–10 Mb. Lastly, we identified signature genes, which were closely related to high-altitude adaptation (ADIRF) and prolificity (CNTROB, SMC3, and PTEN). This study provides a valuable resource for future studies on genome-wide perspectives on the diversity and selection signatures of Chinese indigenous goats.     

Genetic Variability of PRNP in Chinese Indigenous Goats

Polymorphism of the prion protein gene (PRNP) is usually associated with scrapie susceptibility or resistance. To determine the variability of PRNP in Chinese indigenous goat breeds, we isolated genomic DNA from goat blood and amplified and sequenced the coding region of the gene. We identified 10 polymorphic sites that gave rise to 28 haplotypes. Clear frequency differences were found between northern and southern breeds and confirmed by genetic distance analysis, except for the Tangshan dairy goat. Phylogeographic analysis supported the idea that northern and southern breeds might be considered separate clusters, except for the Tangshan dairy goat. The finding of significant differences in allele distribution in northern and southern goats, especially if involved in modulating resistance/susceptibility, needs to be carefully considered for the feasibility of selection plans for resistance to scrapie.     

Generation of VEGF knock-in Cashmere goat via the CRISPR/Cas9 system

Cashmere is a rare and specialised animal fibre, which grows on the outer skin of goats. Owing its low yield and soft, light, and warm properties, it has a high economic value. Here, we attempted to improve existing cashmere goat breeds by simultaneously increasing their fibre length and cashmere yield. We attempted this by knocking in the vascular endothelial growth factor (VEGF) at the fibroblast growth factor 5(FGF5) site using a gene editing technology and then studying its hair growth-promoting mechanisms. We show that a combination of RS-1 and NU7441 significantly improve the efficiency of CRISPR/Cas9-mediated, homologous-directed repair without affecting the embryo cleavage rate or the percentages of embryos at different stages. In addition, we obtained a cashmere goat, which integrated the VEGF gene at the FGF5 site, and the cashmere yield and fibre length of this gene-edited goat were improved. Through next-generation sequencing, we found that the up-regulation of VEGF and the down-regulation of FGF5 affected the cell cycle, proliferation, and vascular tone through the PI3K-AKT signalling pathway and at extracellular matrix-receptor interactions. Owing to this, the gene-edited cashmere goat showed impressive cashmere performance. Overall, in this study, we generated a gene-edited cashmere goat by integrating VEGF at the FGF5 site and provided an animal model for follow-up research on hair growth mechanisms.     

Genetic Variations of 13 Indigenous Chinese Goat Breeds Based on Cytochrome b Gene Sequences

Phylogenetic relationships among and genetic variability within 13 Chinese indigenous goat breeds and Boer goat were analyzed using cytochrome b gene sequences. There were 44 variable sites found in a 642 bp sequence, and 46 Cyt b haplotypes were subsequently defined. The phylogeny analysis of haplotypes in combination with goat Cyt b sequences from GenBank shows that Chinese goats are obviously separated from wild goats and might come from Capra aegagrus. Further analysis indicated that indigenous Chinese goats might descend from at least two lineages; most of the individuals analyzed could be classified into lineage A as defined by Luikart, but five other goats were of uncertain lineage. The Tibet plateau is a possible place of origin for Chinese goats. The neighbor-joining tree based on pairwise differences among populations shows that most Tibetan goats, except the Middle Tibet type, cluster closely with North China goats, and then with South China goats. This result confirms that differences in genetic structure exist among goats in different geographic locations. Nucleotide diversity varied among populations. Tibet and North China goats had higher genetic diversity than South China goats. The fixation index (F _st=87.72%) suggested that most of the total genetic variation was due to variation within populations. In addition, the results indicate that Cyt b gene sequence information alone might not be enough for phylogeny analysis among breeds within species, as shown by fewer polymorphic sites and lower bootstrap values on the neighbor-joining tree.     

Characterization of the <Emphasis Type="Italic">GHR</Emphasis> gene genetic variation in Chinese indigenous goat breeds

The aim of the present work was to investigate single nucleotide polymorphism (SNP) of growth hormone receptor (GHR) gene exon 10, characterize the genetic variation in three Chinese indigenous goat breeds, and search for its potential association with cashmere traits. In this study, a polymerase chain reaction-single strand conformation polymorphism (PCR–SSCP) protocol has been developed for rapid genotyping of the GHR gene in goats. One hundred seventy-eight goats from Liaoning Cashmere (96), Inner Mongolia White Cashmere (40), and Chengdu Grey (42) breeds in China were genotyped at GHR locus using the protocol developed. In all goat breeds investigated, a SNP in exon 10 of GHR gene has been identified by analyzing genomic DNA. The polymorphism consists of a single nucleotide substitution A → G, resulting in two alleles named, respectively, A and G based on the nucleotide at the position. The allele A was found to be more common in the animals investigated, and seems to be more consistent with cattle and zebu at this polymorphic site found in goats. The Hardy–Weinberg equilibrium of genotype distributions of GHR locus was verified in Liaoning Cashmere, and Inner Mongolia White Cashmere breeds. According to the classification of polymorphism information content (PIC), Chengdu Grey was less polymorphic than Liaoning Cashmere and Inner Mongolia White Cashmere breeds at this locus. The phylogenetic tree of different species based on the nucleotide sequences of GHR gene exon 10 is generally in agreement with the known species relationship. No significant association was found between the polymorphism revealed and the cashmere traits analyzed in present work.     

Genetic diversity in the Churra tensina and Churra lebrijana endangered Spanish sheep breeds and relationship with other Churra group breeds and Spanish mouflon

The aim of the present study was to estimate the genetic intra-breed variability of Churra tensina and Churra lebrijana endangered breeds and to establish genetic relationships with Churra, Latxa and Merino breeds, as well as Spanish mouflon, by using 28 microsatellite markers, to provide useful information for their conservation. Allele frequencies and heterozygosity revealed high genetic variation in the two endangered breeds despite their small population size. Estimates of inbreeding coefficient (F_IS) were significant for all breeds studied, except for Churra lebrijana breed. The highest inbreeding coefficient (F_IS = 0.143) was found in the Spanish mouflon. Genetic differentiation tests (F_ST = 0.121) and assignment of individuals to populations indicated the existence of defined breed populations, and low genetic flow between these breeds. The highest pairwise Reynolds distance (D_R) values were observed between Mouflon and the domestic sheep breeds. Considering only domestic sheep breeds, the Churra lebrijana breed showed the highest pairwise D_R values_. The lowest values were found between Latxa and the other domestic sheep, except for Churra lebrijana. Results of pairwise D_R values, as well as phylogenetic tree and bottleneck analysis showed an important genetic isolation of the Churra lebrijana breed from the other Churra types, and genetic signatures of a demographic bottleneck. Finally, structure analysis of populations detected a population subdivision in the Latxa sheep breed. In conclusion, this study presents valuable insight into the existing genetic variability of two Spanish endangered breeds, as well as the first study in Spanish mouflon based on microsatellite analysis. The high degree of variability demonstrated in Churra tensina and Churra lebrijana implies that these populations are rich reservoirs of genetic diversity.     

A novel genetic variant of the goat Six6 gene and its association with production traits in Chinese goat breeds

We looked for novel genetic variations within the Six6 gene by PCR-SSCP, DNA sequencing and forced RFLP-PCR and estimated their associations with production traits in 2132 goats of eight indigenous Chinese breeds. A novel single nucleotide polymorphism (NM_001104993.1: g.232T>C) within the goat Six6 gene was identified. The frequencies of allele "C" varied from 0.8621 to 1.000, which were in Hardy-Weinberg equilibrium. Genotype and allele frequencies were found to be significantly different in the breeds. Association analysis revealed a significant relationship between genotypes and cannon circumference in Hainan black goats (P = 0.047); adult individuals with genotype CT showed higher cannon circumference than those with genotype CC. Another significant association of genotypes with five-year-old fiber length was found in Inner Mongolia white cashmere goats (P = 0.002). In addition, individuals with genotype CT had longer fiber length than those with genotype CC. The data revealed that the Six6 gene positively affects growth traits and cashmere traits. Polymorphism of Six6-PstI could be useful as a DNA marker for goat breeding and genetics via marker-assisted selection.     

None of polymorphism of ovine fecundity major genes FecB and FecX was tested in goat

The polymorphism of mutation Q249R in BMPR-IB gene (FecB) and loci FecX^I, FecX^H, FecX^G, FecX^B in BMP15 gene was analyzed by forced PCR-RFLP method in 550 individuals from 6 flocks or breeds of goats with litter size varied from 1.4 to 2.7 including Boer (209), Haimen (128), second generation of Boer goat crossed with Huanghuai goat (82), Huanghuai (71), Nubi (37) and Matou (23) goat. None of mutations was detected in these goat breeds and their crossbreed. These results suggest that fecundity of goat is not linked to the same loci in BMPR-IB and BMP15 as sheep. Therefore, it is necessary to seek for other genes or loci in order to develop marker assistance selection techniques and study the prolific mechanism of the goat.     

Genetic diversity and population structure in multiple Chinese goat populations using a SNP panel

Information about genetic diversity and population structure among goat breeds is essential for genetic improvement, understanding of environmental adaptation as well as utilization and conservation of goat breeds. Here, we measured genetic diversity and population structure in multiple Chinese goat populations, namely, Nanjiang, Qinggeda, Arbas Cashmere, Jining Grey, Luoping Yellow and Guangfeng goats. A total of 193 individuals were genotyped for about 47 401 autosomal single nucleotide polymorphisms (SNPs). We found a high proportion of informative SNPs, ranging from 69.5% in the Luoping Yellow to 93.9% in the Jining Grey goat breeds with an average mean of 84.7%. Diversity, as measured by expected heterozygosity, ranged from 0.371 in Luoping Yellow to 0.405 in Jining Grey goat populations. The average estimated pair‐wise genetic differentiation (F_ST) among the populations was 8.6%, ranging from 0.2% to 16% and indicating low to moderate genetic differentiation. Principal component analysis, genetic structure and phylogenetic tree analysis revealed a clustering of six Chinese goat populations according to geographic distribution. The results from this study can contribute valuable genetic information and can properly assist with within‐breed diversity, which provides a good opportunity for sustainable utilization of and maintenance of genetic resource improvements in the Chinese goat populations.     

Genetic diversity and relationship among southern Indian goat breeds based on microsatellite markers

Genetic diversity and relationship among goat breeds of southern India were investigated based on microsatellite markers. All five breeds of south India, namely Attappady, Osmanabadi, Sangamneri, Malabari and Kanniadu, along with Ganjam of eastern India were considered for the investigation. In total, 190 alleles were observed from 288 DNA samples analysed with 25 microsatellite loci across six breeds. The most diverse breed was Kanniadu and the least was Osmanabadi. Gene diversity for each breed ranged from 0.73 in Kanniadu to 0.61 in Osmanabadi. The genetic distance tended to be least (0.22) between Ganjam and Malabari and the widest (0.83) between Kanniadu and Malabari. The genetic differentiation between different pairs of the breeds was significantly different from zero. Majority of the loci in almost all the breeds were heterozygote deficit. The overall F_is value (0.20) was moderate and significantly different from zero. Principal-component analysis showed the clustering of the goat breeds according to their geographical origin. Therefore, the geographical origin of the breeds should be taken into consideration while deciding conservation and improvement options for these breeds.     

The polymorphism of a novel 30 bp-deletion mutation at KAP9.2 locus in the cashmere goat

The keratins and keratin-associated proteins (KAPs) are a large heterogeneous group of proteins that make up about 90% of the cashmere fiber. Keratin-associated proteins 9.2 gene (KAP9.2) is one of the ultra high sulfur KAPs, which might play an important role in the bundling of intermediate filaments. In this study, the deletion/insertion mutation of KAP9.2 gene in 997 cashmere goat samples was firstly detected, at the same time, parts of these samples were sequenced. The results showed that two alleles were detected at this KAP9.2P1 locus, named allele W and D. The frequencies of the KAP9.2-W allele in Inner Mongolia White cashmere (n = 785) and Shaanbei White cashmere goat breeds (n = 212) were 0.878 and 0.790, respectively. The χ²-test showed that the genotype distributions in these two cashmere goat breeds were not in agreement with Hardy–Weinberg equilibrium. According to the classification of polymorphism information content (PIC), Shaanbei White cashmere goat was more polymorphic at this locus. Moreover a 30 bp-deletion mutation was described at KAP9.2P2 locus for the first time and no deletion/insertion was described at KAP9.2P1 locus. The results possibly revealed that the size polymorphism existed in the two Chinese cashmere goat and the 30 bp-deletion mutation was possibly caused by variations in the number of the decapeptide repeat structures.     

Mapping of the transcription start site (TSS) and identification of SNPs in the bovine neuropeptide Y (NPY) gene

Background

The domestic goat is one of the important livestock species of India. In the present study we assess genetic diversity of Indian goats using 17 microsatellite markers. Breeds were sampled from their natural habitat, covering different agroclimatic zones.

Results

The mean number of alleles per locus (NA) ranged from 8.1 in Barbari to 9.7 in Jakhrana goats. The mean expected heterozygosity (He) ranged from 0.739 in Barbari to 0.783 in Jakhrana goats. Deviations from Hardy-Weinberg Equilibrium (HWE) were statistically significant (P < 0.05) for 5 loci breed combinations. The D_A measure of genetic distance between pairs of breeds indicated that the lowest distance was between Marwari and Sirohi (0.135). The highest distance was between Pashmina and Black Bengal. An analysis of molecular variance indicated that 6.59% of variance exists among the Indian goat breeds. Both a phylogenetic tree and Principal Component Analysis showed the distribution of breeds in two major clusters with respect to their geographic distribution.

Conclusion

Our study concludes that Indian goat populations can be classified into distinct genetic groups or breeds based on the microsatellites as well as mtDNA information.