QTL affecting conformation traits in Angora goats

Genetic characterization of Markhoz goat breed using microsatellite markers

The use of molecular markers can support the management of endangered populations and should be combined with appropriate breeding strategies to improve productive traits avoiding the decline of the breed. The genetic variability at 10 microsatellite loci were investigated in a sample of 100 unrelated Markhoz goats (77 females and 23 males). The investigated population was reared at the Sanandaj Markhoz goat Performance Testing Station in Kurdistan, Iran. Markhoz goat, a multipurpose breed, is one of the most valuable genetic resources in Iran. All the studied loci were found to be polymorphic and a total number of 52 alleles were identified with an average number of alleles of 5.2. Moreover, some population genetic indices, such as observed and expected heterozygosity, observed and expected number of alleles, Shannon's index, Nei's expected heterozygosity, and polymorphism information content were also calculated. Despite the decreasing population size, Markhoz goat genetic diversity is still conserved. The breed seems to have a good level of genetic variability and, as a consequence, a potential margin of adaptability to environment and for future genetic improvement.

Genome-wide association mapping for type and mammary health traits in French dairy goats identifies a pleiotropic region on chromosome 19 in the Saanen breed

Type traits and mammary health traits are important to dairy ruminant breeding because they influence animal health, milking ability, and longevity, as well as the economic sustainability of farms. The availability of the genomic sequence and a single nucleotide polymorphism chip in goats has opened up new fields of investigation to better understand the genes and mechanisms that underlie such complex traits and to be able to select them. Our objective was to perform a genome-wide association study in dairy goats for 11 type traits and somatic cell count (SCC) as proxies for mastitis resistance. A genome-wide association study was implemented using a daughter design composed of 1,941 Alpine and Saanen goats sired by 20 artificial insemination bucks, genotyped with the Illumina GoatSNP50 BeadChip (Illumina Inc., San Diego, CA). This association study was based on both linkage analyses and linkage disequilibrium using QTLmap software (http://dga7.jouy.inra.fr/qtlmap/) interval mapping was performed with the likelihood ratio test using linear regressions. Breeds were analyzed together and separately. The study highlighted 37 chromosome-wide significant quantitative trait loci (QTL) with linkage analyses and 222 genome-wide significant QTL for linkage disequilibrium, for type and SCC traits in dairy goats. Genomic control of those traits was mostly polygenic and breed-specific, suggesting that within-breed selection would be favored for those traits. Of note, Capra hircus autosome (CHI) 19 appeared to be highly enriched in single nucleotide polymorphisms associated with type and SCC, with 2 highly significant regions in the Saanen breed. One region (33-42 Mb) was significantly associated with SCC and includes candidate genes associated with response to intramammary infections (RARA, STAT3, STAT5A, and STAT5B). Another region of the CHI 19 (24.5-27 Mb) exhibited an adverse pleiotropic effect on milk production (milk, fat yield, and protein yield) and udder traits (udder floor position and rear udder attachment) that agreed with the negative genetic correlations that exist between those 2 groups of traits. These QTL were not found in the Alpine breed. In Alpine, the 2 most significant regions were associated with chest depth on CHI 6 (45.8-46.0 Mb) and CHI 8 (80.7-81.1 Mb). These results will be helpful for goat selection in the future and could lead to identification of causal mutations.

Genome-wide association study of conformation and milk yield in mixed-breed dairy goats

Identification of genetic markers that affect economically important traits is of high value from a biological point of view, enabling the targeting of candidate genes and providing practical benefits for the industry such as wide-scale genomic selection. This study is one of the first to investigate the genetic background of economically important traits in dairy goats using the caprine 50K single nucleotide polymorphism (SNP) chip. The aim of the project was to perform a genome-wide association study for milk yield and conformation of udder, teat, and feet and legs. A total of 137,235 milk yield records on 4,563 goats each scored for 10 conformation traits were available. Out of these, 2,381 goats were genotyped with the Illumina Caprine 50K BeadChip (Illumina Inc., San Diego, CA). A range of pseudo-phenotypes were used including deregressed breeding values and pseudo-estimated breeding values. Genome-wide association studies were performed using the multi-locus mixed model (MLMM) algorithm implemented in SNP & Variation Suite v7.7.8 (Golden Helix Inc., Bozeman, MT). A genome-wise significant [-log₁₀(P-value) > 5.95] SNP for milk yield was identified on chromosome 19, with additional chromosome-wise significant (-log₁₀(P-value) > 4.46] SNP on chromosomes 4, 8, 14, and 29. Three genome-wise significant SNP for conformation of udder attachment, udder depth, and front legs were identified on chromosome 19, and chromosome-wise SNP were found on chromosomes 4, 5, 6, 10, 11, 12, 13, 14, 15, 16, 17, 18, 21, 23, and 27. The proportion of variance explained by the significant SNP was between 0.4 and 7.0% for milk yield and between 0.1 and 13.8% for conformation traits. This study is the first attempt to identify SNP associated with milk yield and conformation in dairy goats. Two genome-wise significant SNP for milk yield and 3 SNP for conformation of udder attachment, udder depth, and front legs were found. Our results suggest that conformation traits have a polygenic background because, for most of them, we did not identify any quantitative trait loci with major effect.

The Application of Genomic Technologies to Investigate the Inheritance of Economically Important Traits in Goats

Goat genomics has evolved at a low pace because of a lack of molecular tools and sufficient investment. Whilst thousands and hundreds of quantitative trait loci (QTL) have been identified in cattle and sheep, respectively, about nine genome scans have been performed in goats dealing with traits as conformation, growth, fiber quality, resistance to nematodes, and milk yield and composition. In contrast, a great effort has been devoted to the characterization of candidate genes and their association with milk, meat, and reproduction phenotypes. In this regard, causal mutations have been identified in the αS1-casein gene that has a strong effect on milk composition and the PIS locus that is linked to intersexuality and polledness. In recent times, the development of massive parallel sequencing technologies has allowed to build a reference genome for goats as well as to monitor the expression of mRNAs and microRNAs in a broad array of tissues and experimental conditions. Besides, the recent design of a 52K SNP chip is expected to have a broad impact in the analysis of the genetic architecture of traits of economic interest as well as in the study of the population structure of goats at a worldwide scale.

Genome-scan analysis for genetic mapping of quantitative trait loci underlying birth weight and onset of puberty in doe kids (Capra hircus)

The objective of this study was to locate quantitative trait loci (QTL) causing variation in birth weight and age of puberty of doe kids in a population of Rayini cashmere goats. Four hundred and thirty kids from five half-sib families were genotyped for 116 microsatellite markers located on the caprine autosomes. The traits recorded were birth weight of the male and female kids, body weight at puberty, average daily gain from birth to age of puberty and age at puberty of the doe kids. QTL analysis was conducted using the least squares interval mapping approach. Linkage analysis indicated significant QTL for birth weight on Capra hircus chromosomes (CHI) 4, 5, 6, 18 and 21. Five QTL located on CHI 5, 14 and 29 were associated with age at puberty. Across-family analysis revealed evidence for overlapping QTL affecting birth weight (78 cM), body weight at puberty (72 cM), average daily gain from birth to age of puberty (72 cM) and age at puberty (76 cM) on CHI 5 and overlapping QTL controlling body weight at puberty and age at puberty on CHI 14 at 18-19 cM. The proportion of the phenotypic variance explained by the detected QTL ranged between 7.9% and 14.4%. Confirming some of the previously reported results for birth weight and growth QTL in goats, this study identified more QTL for these traits and is the first report of QTL for onset of puberty in doe kids.

An initial comparative map of copy number variations in the goat (Capra hircus) genome

BACKGROUND

The goat (Capra hircus) represents one of the most important farm animal species. It is reared in all continents with an estimated world population of about 800 million of animals. Despite its importance, studies on the goat genome are still in their infancy compared to those in other farm animal species. Comparative mapping between cattle and goat showed only a few rearrangements in agreement with the similarity of chromosome banding. We carried out a cross species cattle-goat array comparative genome hybridization (aCGH) experiment in order to identify copy number variations (CNVs) in the goat genome analysing animals of different breeds (Saanen, Camosciata delle Alpi, Girgentana, and Murciano-Granadina) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome.

RESULTS

We identified a total of 161 CNVs (an average of 17.9 CNVs per goat), with the largest number in the Saanen breed and the lowest in the Camosciata delle Alpi goat. By aggregating overlapping CNVs identified in different animals we determined CNV regions (CNVRs): on the whole, we identified 127 CNVRs covering about 11.47 Mb of the virtual goat genome referred to the bovine genome (0.435% of the latter genome). These 127 CNVRs included 86 loss and 41 gain and ranged from about 24 kb to about 1.07 Mb with a mean and median equal to 90,292 bp and 49,530 bp, respectively. To evaluate whether the identified goat CNVRs overlap with those reported in the cattle genome, we compared our results with those obtained in four independent cattle experiments. Overlapping between goat and cattle CNVRs was highly significant (P < 0.0001) suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Genes with environmental functions were over-represented in goat CNVRs as reported in other mammals.

CONCLUSIONS

We describe a first map of goat CNVRs. This provides information on a comparative basis with the cattle genome by identifying putative recurrent interspecies CNVs between these two ruminant species. Several goat CNVs affect genes with important biological functions. Further studies are needed to evaluate the functional relevance of these CNVs and their effects on behavior, production, and disease resistance traits in goats.

QTL detection for milk production traits in goats using a longitudinal model

Eight paternal half-sib families were used to identify chromosomal regions associated with variation in the lactation curves of dairy goats. DNA samples from 162 animals were amplified by PCR for 37 microsatellite markers, from Capra hircus autosomes CHI3, CHI6, CHI14 and CHI20. Milk samples were collected during 6 years, and there were 897 records for milk yield (MY) and 814 for fat (FP) and protein percentage (PP). The analysis was conducted in two stages. First, a random regression model with several fixed effects was fitted to describe the lactation function, using a scale (alpha) plus four shape parameters: beta and gamma, both associated with a decrease in the slope of the curve, and delta and phi that are related to the increase in slope. Predictions of alpha, beta, gamma, delta and phi were regressed using an interval mapping model, and F-tests were used to test for quantitative trait loci (QTL) effects. Significant (p < 0.05) QTLs were found for: (i) MY: CHI6 at 70-80 cM for all parameters; CHI14 at 14 cM for delta and phi; (ii) FP: CHI14, at 63 cM was associated with beta; CHI20, at 72 cM, showed association with alpha; (iii) PP: chromosomal regions associated with beta were found at 59 cM in CHI3 and at 55 cM in CHI20 with alpha and gamma. Analyses using more families and more animals will be useful to confirm or to reject these findings.