中国绒山羊遗传多样性现状和系统发生关系的微卫星分析

为了调查中国绒山羊遗传资源现状, 作者应用联合国粮农组织和国际家畜研究所推荐的19对微卫星引物并结合荧光PCR技术, 对9个中国地方产绒山羊群体和1个西非山羊品种进行了遗传多样性检测。14个微卫星座位在10个山羊群体中显示为高度多态, 可作为山羊遗传多样性分析的有效标记。多态信息含量和遗传杂合度等数据表明: 目前中国地方产绒山羊群体的遗传多样性较为丰富, 并且大部分保种场较好地保存了这些地方资源。采用非加权配对算术平均法构建的聚类图和采用主成分分析法得到的散点图均显示, 中国山羊与西非山羊为不同的2类; 中国产绒山羊中河谷山羊、河西绒山羊与其他山羊的遗传距离较远; 其他山羊又大致分为2类: 一类由辽宁绒山羊、新疆山羊、柴达木山羊、陕北山羊组成, 另一类由内蒙古绒山羊组成。此研究结果为开展我国地方绒山羊种质特性研究及资源保护和利用提供了科学依据。     

Genetic Structure of Chinese Indigenous Goats and the Special Geographical Structure in the Southwest China as a Geographic Barrier Driving the Fragmentation of a Large Population

Background

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.

Methodology/Principal Findings

Forty Chinese indigenous populations containing 2078 goats were sampled from different geographic regions of China. Moderate genetic diversity at the population level (H_S of 0.644) and high population diversity at the species level (H_T value of 0.737) were estimated. Significant moderate population differentiation was detected (F_ST value of 0.129). Significant excess homozygosity (F_IS 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.

Conclusions/Significance

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.     

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 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.     

Investigation of the genetic diversity of domestic Capra hircus breeds reared within an early goat domestication area in Iran

Background

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.

Results

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 (G_ST = 5.9%). Positive and highly significant F_IS 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.

Conclusions

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.     

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

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

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

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

Mitochondrial diversity and phylogeographic structure of Chinese domestic goats

China has numerous native domestic goat breeds, but so far there has been no extensive study on genetic diversity, population demographic history, and origin of Chinese goats. Here, we examined the genetic diversity and phylogeographic structure of Chinese domestic goats by determining a 481-bp fragment of the first hypervariable region of mitochondrial DNA (mtDNA) control region from 368 individuals representing 18 indigenous breeds. Phylogenetic analyses revealed that there were four mtDNA lineages (A-D) identified in Chinese goats, in which lineage A was predominant, lineage B was moderate, and lineages C and D were at low frequency. These results further support the multiple maternal origins of domestic goats. The pattern of genetic variation in goat mtDNA sequences indicated that the two larger lineages A and B had undergone population expansion events. In a combined analysis of previously reported sequences and our sequences belonging to lineage B, we detected two subclades, in which one was unique to eastern Asia and another was shared between eastern and southern Asia. A larger genetic variation in eastern Asia than southern Asia and the pattern of phylogeographic variation in lineage B suggest that at least one subclade of lineage B originated from eastern Asia. There was no significant geographical structuring in Chinese goat populations, which suggested that there existed strong gene flow among goat populations caused by extensive transportation of goats in history.     

Targeted high‐throughput growth hormone 1 gene sequencing reveals high within‐breed genetic diversity in South African goats

This study assessed the genetic diversity in the growth hormone 1 gene (GH1) within and between South African goat breeds. Polymerase chain reaction‐targeted gene amplification together with Illumina MiSeq next‐generation sequencing (NGS) was used to generate the full length (2.54 kb) of the growth hormone 1 gene and screen for SNPs in the South African Boer (SAB) (n = 17), Tankwa (n = 15) and South African village (n = 35) goat populations. A range of 27–58 SNPs per population were observed. Mutations resulting in amino acid changes were observed at exons 2 and 5. Higher within‐breed diversity of 97.37% was observed within the population category consisting of SA village ecotypes and the Tankwa goats. Highest pairwise F_ST values ranging from 0.148 to 0.356 were observed between the SAB and both the South African village and Tankwa feral goat populations. Phylogenetic analysis indicated nine genetic clusters, which reflected close relationships between the South African populations and the other international breeds with the exception of the Italian Sarda breeds. Results imply greater potential for within‐population selection programs, particularly with SA village goats.     

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.     

Genetic diversity analyses reveal first insights into breed‐specific selection signatures within Swiss goat breeds

We used genotype data from the caprine 50k Illumina BeadChip for the assessment of genetic diversity within and between 10 local Swiss goat breeds. Three different cluster methods allowed the goat samples to be assigned to the respective breed groups, whilst the samples of Nera Verzasca and Tessin Grey goats could not be differentiated from each other. The results of the different genetic diversity measures show that Appenzell, Toggenburg, Valais and Booted goats should be prioritized in future conservation activities. Furthermore, we examined runs of homozygosity (ROH) and compared genomic inbreeding coefficients based on ROH (F_ROH) with pedigree‐based inbreeding coefficients (F_PED). The linear relationship between F_ROH and F_PED was confirmed for goats by including samples from the three main breeds (Saanen, Chamois and Toggenburg goats). F_ROH appears to be a suitable measure for describing levels of inbreeding in goat breeds with missing pedigree information. Finally, we derived selection signatures between the breeds. We report a total of 384 putative selection signals. The 25 most significant windows contained genes known for traits such as: coat color variation (MITF, KIT, ASIP), growth (IGF2, IGF2R, HRAS, FGFR3) and milk composition (PITX2). Several other putative genes involved in the formation of populations, which might have been selected for adaptation to the alpine environment, are highlighted. The results provide a contemporary background for the management of genetic diversity in local Swiss goat breeds.     

Genetic diversity and patterns of population structure in Creole goats from the Americas

Biodiversity studies are more efficient when large numbers of breeds belonging to several countries are involved, as they allow for an in‐depth analysis of the within‐ and between‐breed components of genetic diversity. A set of 21 microsatellites was used to investigate the genetic composition of 24 Creole goat breeds (910 animals) from 10 countries to estimate levels of genetic variability, infer population structure and understand genetic relationships among populations across the American continent. Three commercial transboundary breeds were included in the analyses to investigate admixture with Creole goats. Overall, the genetic diversity of Creole populations (mean number of alleles = 5.82 ± 1.14, observed heterozygosity = 0.585 ± 0.074) was moderate and slightly lower than what was detected in other studies with breeds from other regions. The Bayesian clustering analysis without prior information on source populations identified 22 breed clusters. Three groups comprised more than one population, namely from Brazil (Azul and Graúna; Moxotó and Repartida) and Argentina (Long and shorthair Chilluda, Pampeana Colorada and Angora‐type goat). Substructure was found in Criolla Paraguaya. When prior information on sample origin was considered, 92% of the individuals were assigned to the source population (threshold q ≥ 0.700). Creole breeds are well‐differentiated entities (mean coefficient of genetic differentiation = 0.111 ± 0.048, with the exception of isolated island populations). Dilution from admixture with commercial transboundary breeds appears to be negligible. Significant levels of inbreeding were detected (inbreeding coefficient > 0 in most Creole goat populations, P < 0.05). Our results provide a broad perspective on the extant genetic diversity of Creole goats, however further studies are needed to understand whether the observed geographical patterns of population structure may reflect the mode of goat colonization in the Americas.     

Genetic relationships among twelve Chinese indigenous goat populations based on microsatellite analysis

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 F_ST (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.     

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.     

Genetic diversity and relationships of 10 Chinese goat breeds in the Middle and Western China

The genetic variability and relationships among 10 Chinese goat breeds (Shannan white goat, Shanbei white cashmere goat, Funiu goat, Huanghuai goat, Taihang goat, Guanzhong dairy goat, Xinong Saanen goat, Zhongwei goat, Tibetan goat and Inner Mongolia cashmere goat) were studied using 20 microsatellite markers. The mean heterozygosity and the mean polymorphism information content (PIC) for the 10 goat breeds varied from 0.7612 to 0.8390 and 0.7293 to 0.8181, respectively. Huanghuai goat demonstrated the highest genetic variability, and Tibetan goat showed the lowest. G_ST index was 0.052. A unweighted pair group method with arithmetic means (UPGMA) diagram obtained based on Nei standard genetic distance displayed some degree of consistency with their geographical locations and production.     

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.     

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.     

Analysis of caprine pituitary specific transcription factor-1 gene polymorphism in indigenous Chinese goats

Since mutations on POU1F1 gene possibly resulted in deficiency of GH, PRL, TSH and POU1F1, this study revealed the polymorphism of goat POU1F1-AluI locus and analyzed the distribution of alleles on 13 indigenous Chinese goat breeds. The PCR-RFLP analysis showed the predominance of TT genotype and the frequencies of allele T varied from 0.757 to 0.976 in the analyzed populations (SBWC, Bo, XH and HM). Further study, distributions of genotypic and allelic frequencies at this locus were found to be significantly different among populations based on a χ²-test (P < 0.001), suggesting that the breed factor significantly affected the molecular genetic character of POU1F1 gene. The genetic diversity analysis revealed that Chinese indigenous populations had a wide spectrum of genetic diversity in goat POU1F1-AluI locus. However, the ANOVA analysis revealed no significant differences for gene homozygosty, gene heterozygosty, effective allele numbers and PIC (polymorphism information content) among meat, dairy and cashmere utility types (P > 0.05), suggesting that goat utility types had no significant effect on the spectrum of genetic diversity. X. Y. Lan and M. J. Li equally contributed to this work.     

Mitochondrial DNA diversity and origins of domestic goats in Southwest China (excluding Tibet)

Southwest China contains about one third Chinese indigenous goat breeds representing special economic and ecological characteristics. Mitochondrial DNA (mtDNA) D-loop sequences of 312 individuals (including 109 new individuals, 203 individuals retrieved from GenBank) from 18 Chinese domestic goat breeds were used to investigate breed genetic diversity, origin and phylogeography. All goat breeds in this study proved to be extremely diverse, average haplotype diversity and nucleotide diversity being 0.9829 ± 0.0027 and 0.03615 ± 0.03257, respectively. The 312 sequences gave 148 different haplotypes. Phylogenetic analyses revealed that there were two mtDNA haplogroups identified in domestic goats in Southwest China, in which haplogroup A was predominant. Mismatch analysis showed haplogroup A had experienced population expansion events, whereas haplogroup B did not. Amova analysis showed there was no significant geographical structuring. Almost 86.23% of genetic variation was included in the within-breed variance component and only 3.5% was observed among the four geographic provinces. The results of this study contribute to the knowledge of the genetic structure and origin of domestic goats in Southwest China.