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

Wattles in goats are associated with the FMN1/GREM1 region on chromosome 10

The presence of congenital appendages (wattles) on the throat of goats is supposed to be under genetic control with a dominant mode of inheritance. Wattles contain a cartilaginous core covered with normal skin resembling early stages of extremities. To map the dominant caprine wattles (W) locus, we collected samples of 174 goats with wattles and 167 goats without wattles from nine different Swiss goat breeds. The samples were genotyped with the 53k goat SNP chip for a subsequent genome‐wide association study. We obtained a single strong association signal on chromosome 10 in a region containing functional candidate genes for limb development and outgrowth. We sequenced the whole genomes of an informative family trio containing an offspring without wattles and its heterozygous parents with wattles. In the associated goat chromosome 10 region, a total of 1055 SNPs and short indels perfectly co‐segregate with the W allele. None of the variants were perfectly associated with the phenotype after analyzing the genome sequences of eight additional goats. We speculate that the causative mutation is located in one of the numerous gaps in the current version of the goat reference sequence and/or represents a larger structural variant which influences the expression of the FMN1 and/or GREM1 genes. Also, we cannot rule out possible genetic or allelic heterogeneity. Our genetic findings support earlier assumptions that wattles are rudimentary developed extremities.     

Genetic relationship among Mongolian native goat populations estimated by blood protein polymorphism

Several local strains and populations of goats distinguished by morphogenetic and performance characteristics are kept by goat breeders in different natural climatic regions of Mongolia, namely Bayandelger, Ulgii Red, Erchim Black, Dorgon and Zavkhan Buural. The genetic relationships among eight native goat populations in Mongolia at 33 biochemical genetic loci was assessed. A total of 440 animals in eight regional zones were studied. Twelve loci, i.e. the serum transferrin, serum amylase, serum alkaline phosphatase, serum prealbumin-3, cell esterase-D, hemoglobin (Hb) β, hemoglobin (Hb) α-II, cell peptidase-B, cell tetrazolium oxidase, cell esterase-1, cell esterase-2 and cell catalase loci, were found to be polymorphic. The data indicated that Mongolian native goats are not highly differentiated (D=0.0002–0.0038) genetically. To set Mongolian native goats in a larger context, the present data were compared with those on other goat breeds and populations in east and southeast Asia that were previously reported. The average heterozygosity in the Mongolian native goats did not significantly differ from those in other Asian goat populations and breeds. A phylogenetic tree of the gene constitution of the Mongolian native goats and other Asian goat breeds and populations was constructed and revealed that genetically the Mongolian native goats had diverged slightly from the group consisting of Chinese, Japanese, Korean and Indonesian native goats, but markedly from the Indian goat group.     

Screening of indigenous goats for prolificacy associated DNA markers of sheep

The present study was undertaken to explore the genetic basis of caprine prolificacy and to screen indigenous goats for prolificacy associated markers of sheep in BMPR1B, GDF9 and BMP15 genes. To detect the associated mutations and identify novel allelic variants in the candidate genes, representative samples were collected from the breeding tract of indigenous goat breeds varying in prolificacy and geographic distribution. DNA was extracted and PCR amplification was done using primers designed or available in literature for the coding DNA sequence of candidate genes. Direct sequencing was done to identify the genetic variations. Mutations in the candidate genes associated with fecundity in sheep were not detected in Indian goats. Three non-synonymous SNPs (C818T, A959C and G1189A) were identified in exon 2 of GDF9 gene out of which mutation A959C has been associated with prolificacy in exotic goats. Two novel SNPs (G735A and C808G) were observed in exon 2 of BMP15 gene.     

Design and Characterization of a 52K SNP Chip for Goats

The success of Genome Wide Association Studies in the discovery of sequence variation linked to complex traits in humans has increased interest in high throughput SNP genotyping assays in livestock species. Primary goals are QTL detection and genomic selection. The purpose here was design of a 50–60,000 SNP chip for goats. The success of a moderate density SNP assay depends on reliable bioinformatic SNP detection procedures, the technological success rate of the SNP design, even spacing of SNPs on the genome and selection of Minor Allele Frequencies (MAF) suitable to use in diverse breeds. Through the federation of three SNP discovery projects consolidated as the International Goat Genome Consortium, we have identified approximately twelve million high quality SNP variants in the goat genome stored in a database together with their biological and technical characteristics. These SNPs were identified within and between six breeds (meat, milk and mixed): Alpine, Boer, Creole, Katjang, Saanen and Savanna, comprising a total of 97 animals. Whole genome and Reduced Representation Library sequences were aligned on >10 kb scaffolds of the de novo goat genome assembly. The 60,000 selected SNPs, evenly spaced on the goat genome, were submitted for oligo manufacturing (Illumina, Inc) and published in dbSNP along with flanking sequences and map position on goat assemblies (i.e. scaffolds and pseudo-chromosomes), sheep genome V2 and cattle UMD3.1 assembly. Ten breeds were then used to validate the SNP content and 52,295 loci could be successfully genotyped and used to generate a final cluster file. The combined strategy of using mainly whole genome Next Generation Sequencing and mapping on a contig genome assembly, complemented with Illumina design tools proved to be efficient in producing this GoatSNP50 chip. Advances in use of molecular markers are expected to accelerate goat genomic studies in coming years.     

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.     

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.     

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.     

Genetic variation within and among three Indian breeds of goat using heterologous microsatellite markers.

The parameters of genetic variation, genetic distances and time of divergence in three Indian goat breeds were studied using 16 cattle microsatellite markers. The mean number of alleles and mean allele size (bp) per microsatellite marker in goats were 5.37 +/- 0.78 and 143.9 +/- 33.75 bp respectively. The average values of heterozygosity and polymorphism information content were 0.54 +/- 0.2 and 0.48 +/- 0.20, respectively. Five of the eight genetic distance methods were highly correlated, revealing a closer relationship between Jamnapari and Barbari goats. A phylogenetic tree constructed from inter-individual distances revealed that the individuals clustered according to the breed to which they belonged, and the Jamnapari and Barbari goats formed a cluster. The divergence times between Sirohi and Jamnapari, and Sirohi and Barbari were approximately 2000 years, while its value between Barbari and Jamnapari goats was approximately 1370 years.     

Genetic markers in the blood of Spanish goat breeds

Biochemical variation of four genetic markers (Hb, Al, Tf, 'X' protein) in the blood of four Spanish goat breeds was analysed with starch gel electrophoresis. The frequencies at all of these loci have been calculated for the Spanish goats and compared with some of the goat breeds studied so far by other authors.     

Assessment of population structure by single nucleotide polymorphisms (SNPs) in goat breeds

Single nucleotide polymorphisms (SNPs) may be used in biodiversity studies and commercial tasks like traceability, paternity testing and selection for suitable genotypes. Twenty-seven SNPs were characterized and genotyped on 250 individuals belonging to eight Italian goat breeds. Multilocus genotype data were used to infer population structure and assign individuals to populations. To estimate the number of groups (K) to test in population structure analysis we used likelihood values and variance of the bootstrap samples, deriving optimal K from a drop in the likelihood and a rise in the variance plots against K.     

Association between expression of reproductive seasonality and alleles of melatonin receptor 1A in goats

To determine whether a link exists between reproductive seasonality and the structure of the melatonin receptor 1A (MTNR1A) gene, the latter was studied in year-round estrous breeds (Jining Grey and Boer goats) and seasonal estrous breeds (Liaoning Cashmere, Inner Mongolia Cashmere, Wendeng milk and Beijing native goats). A large fragment of exon 2 of MTNR1A gene was amplified by PCR using sheep sense and antisense primers in 260 does of six breeds. The uniform 824 bp PCR product was digested with restriction endonucleases MnII and RsaI, and checked for the presence of restriction sites. No polymorphism at the MnII cleavage sites was detected in all six goat breeds and no relationship could be established between the MnII cleavage sites of MTNR1A gene and reproductive seasonality in goats. For polymorphic RsaI cleavage site at base position 53, only genotype RR (267 bp/267 bp) was detected in Jining Grey goats, both genotype RR and genotype Rr (267 bp/320 bp) were found in all other goat breeds, no genotype rr (320 bp/320 bp) was detected in all six goat breeds. Frequency of genotype RR was obviously higher, and frequency of genotype Rr was obviously lower in year-round estrous goat breeds than in seasonal estrous goat breeds. Sequencing revealed one mutation (G52A) in genotype Rr compared with genotype RR. For polymorphic RsaI cleavage site, the differences of genotype distributions were significant (P<0.05) between year-round estrous goat breeds and seasonal estrous goat breeds. These results preliminarily showed an association between genotype RR and year-round estrus in goats, and an association between genotype Rr and seasonal estrus in goats.     

An investigation of five genetic loci controlling polymorphic variants in the red cells of goats

Results of a joint study carried out in South Africa and England to search for new genetic markers in the blood of goats are presented. Haemoglobin (Hb) phenotypes were reinvestigated with the technique of isoelectric focusing; frequencies in different goat breeds are given. Anaemic Hb type A, AB and B goats all produced a Hb C with an identical electrophoretic pattern. All goats tested had identical carbonic anhydrase (CA) types, but showed polymorphism of 'X' protein. Preliminary results indicated that nucleoside phosphorylase (NP) may be polymorphic.     

An investigation of five genetic loci controlling polymorphic variants in the red cells of goats*

Results of a joint study carried out in South Africa and England to search for new genetic markers in the blood of goats are presented. Haemoglobin (Hb) phenotypes were reinvestigated with the technique of isoelectric focusing; frequencies in different goat breeds are given. Anaemic Hb type A, AB and B goats all produced a Hb C with an identical electrophoretic pattern. All goats tested had identical carbonic anhydrase (CA) types, but showed polymorphism of ‘X’ protein. Preliminary results indicated that nucleoside phosphorylase (NP) may be polymorphic.     

Novel polymorphisms of goat growth hormone and growth hormone receptor genes and their effects on growth traits

In this study, the polymorphisms of growth hormone (GH) gene 5' promoter region and intron 8, exons 4 and 10 of growth hormone receptor (GHR) gene were analyzed in Xinong Saanen goats (SG) and Boer goats (BG). Two alleles (A and B) and three genotypes (AA, AB and BB) were detected at P1 locus of GH gene, and two alleles (G and T) and two genotypes (GG and GT) were detected at P4 locus of GHR gene by PCR-SSCP analysis. In addition, two single nucleotide polymorphisms (SNPs)-A73C (P1 locus) and G114T (P4 locus), were identified by DNA sequencing. The frequencies of alleles A and B in the two goat breeds were 0.61-0.62, and 0.39-0.38, respectively, and the frequencies of alleles G and T in the two goat breeds were 0.82-0.86, and 0.18-0.14, respectively. The SNP loci were in Hardy-Weinberg disequilibrium in both goat breeds (P<0.05). Polymorphisms of GH and GHR genes were shown to be associated with growth traits in BG breed. AA and GG genotypes were associated with superior growth traits in 1-, 2- and 3-month old individuals. Hence, AA and GG genotypes are suggested to be a molecular marker for superior growth traits in BG breed.     

Polymorphisms of caprine <Emphasis Type="Italic">GDF</Emphasis>9 gene and their association with litter size in Jining Grey goats

The exons 1, 2 and flanking region of growth differentiation factor 9 (GDF9) gene in five randomly selected does of Jining Grey, Boer and Liaoning Cashmere goats were amplified and analyzed. Thirteen nucleotide differences were identified in GDF9 gene between sheep (AF078545) and goats. Four SNPs (G3288A in intron 1, G423A, A959C [Gln320Pro] and G1189A [Val397Ile] in exon 2) were detected in four goat breeds with different prolificacy, in which G3288A was a new SNP in goats. The results showed that loci 3288, 423 and 1189 in Boer goats, loci 3288 and 423 in Guizhou White goats, loci 423 and 1189 in Liaoning Cashmere goats were all in complete linkage disequilibrium (D′ = 1, r ² = 1), respectively. In moderate (Boer goat) and low prolificacy (Liaoning Cashmere goat) breeds, linkage analysis indicated that there were more fervent linkage disequilibrium among loci 3288, 423 and 1189 than high prolificacy (Jining Grey and Guizhou White goats) breeds. For the 959 locus, the genotype distribution showed obvious difference between high prolificacy breeds and moderate or low prolificacy breeds (P < 0.05 or P < 0.01). The Jining Grey goat does with genotype CC or AC had 0.81 (P < 0.01) or 0.63 (P < 0.01) kids more than those with genotype AA, respectively. The present study preliminarily showed an association between allele C at 959 locus of GDF9 gene and high litter size in Jining Grey goats. These results provide further evidence that the GDF9 gene may be significantly correlated with high prolificacy in goats.     

Performance of different SNP panels for parentage testing in two East Asian cattle breeds

The International Society for Animal Genetics (ISAG) proposed a panel of single nucleotide polymorphisms (SNPs) for parentage testing in cattle (a core panel of 100 SNPs and an additional list of 100 SNPs). However, markers specific to East Asian taurine cattle breeds were not included, and no information is available as to whether the ISAG panel performs adequately for these breeds. We tested ISAG's core (100 SNP) and full (200 SNP) panels on two East Asian taurine breeds: the Korean Hanwoo and the Japanese Wagyu, the latter from the Australian herd. Even though the power of exclusion was high at 0.99 for both ISAG panels, the core panel performed poorly with 3.01% false‐positive assignments in the Hanwoo population and 3.57% in the Wagyu. The full ISAG panel identified all sire–offspring relations correctly in both populations with 0.02% of relations wrongly excluded in the Hanwoo population. Based on these results, we created and tested two population‐specific marker panels: one for the Wagyu population, which showed no false‐positive assignments with either 100 or 200 SNPs, and a second panel for the Hanwoo, which still had some false‐positive assignments with 100 SNPs but no false positives using 200 SNPs. In conclusion, for parentage assignment in East Asian cattle breeds, only the full ISAG panel is adequate for parentage testing. If fewer markers should be used, it is advisable to use population‐specific markers rather than the ISAG panel.     

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.