A genome-wide association analysis for body, udder, and leg conformation traits recorded in Murciano-Granadina goats

The African introgression of Murciano Granadina goats has a Moroccan origin and displays remarkable levels of inter-individual variability

There is evidence that Murciano Granadina (MG), the most important caprine dairy breed in Spain, has been introgressed by African goats, but the precise geographic origin of such introgression has not been identified yet. Moreover, an accurate estimate of the magnitude of this African introgression is lacking, since current estimates are based on small numbers of sampled individuals. The aim of our work was to tackle these two issues by genotyping 500 MG goats with the Goat SNP50 BeadChip and comparing their genotypes with those of reference populations from Spain (Bermeya), France (Saanen), Morocco (Barcha, Draa, Ghazalia, Noire de Atlas, Nord, Moroccan), Egypt (Barki, Oasis, Saidi), Algeria (Arabia, Makatia, M'Zabite, Kabyle), Tunisia (Tunisian native breeds) and Sudan (Desert, Nilotic, Taggar). The population of 500 MG goats was subdivided into 10 datasets of 50 individuals to ensure that sample sizes of the target (MG) and reference populations are balanced. Performance of an unsupervised ADMIXTURE analysis demonstrated that MG goats have a North African ancestry, with an average proportion of 4.4 ± 2.3%. Next, we did a supervised ADMIXTURE analysis that revealed that the Moroccan genetic component reaches a proportion of 4.01 ± 3.9% in MG goats, while the Algerian (0.001 ± 0.001%), Egyptian (0.2 ± 0.1%), Sudanese (0.1 ± 0.1%) and Tunisian (0.3 ± 0.4%) components are present in extremely small proportions. The historical circumstances of this introgression event are currently unknown, but several plausible scenarios are outlined. Moreover, our results show considerable inter-individual heterogeneity regarding the magnitude of the Moroccan introgression of MG goats (0%- 12% depending on the MG data set under analysis). This result implies that reliable estimates about the introgression of autochthonous livestock by exotic breeds can only be obtained by extensively sampling target populations.

Elucidation of population stratifying markers and selective sweeps in crossbred Landlly pig population using genome-wide SNP data

The present study was aimed at the identification of population stratifying markers from the commercial porcine SNP 60K array and elucidate the genome-wide selective sweeps in the crossbred Landlly pig population. Original genotyping data, generated on Landlly pigs, was merged in various combinations with global suid breeds that were grouped as exotic (global pig breeds excluding Indian and Chinese), Chinese (Chinese pig breeds only), and outgroup pig populations. Post quality control, the genome-wide SNPs were ranked for their stratifying power within each dataset in TRES (using three different criteria) and FIFS programs and top-ranked SNPs (0.5K, 1K, 2K, 3K, and 4K densities) were selected. PCA plots were used to assess the stratification power of low-density panels. Selective sweeps were elucidated in the Landlly population using intra- and inter-population haplotype statistics. Additionally, Tajima's D-statistics were calculated to determine the status of balancing selection in the Landlly population. PCA plots showed 0.5K marker density to effectively stratify Landlly from other pig populations. The A-score in DAPC program revealed the Delta statistic of marker selection to outperform other methods (informativeness and FST methods) and that 3000-marker density was suitable for stratification of Landlly animals from exotic pig populations. The results from selective sweep analysis revealed the Landlly population to be under selection for mammary (NAV2), reproductive efficiency (JMY, SERGEF, and MAP3K20), body conformation (FHIT, WNT2, ASRB, DMGDH, and BHMT), feed efficiency (CSRNP1 and ADRA1A), and immunity (U6, MYO3B, RBMS3, and FAM78B) traits. More than two methods suggested sweeps for immunity and feed efficiency traits, thus giving a strong indication for selection in this direction. The study is the first of its kind in Indian pig breeds with a comparison against global breeds. In conclusion, 500 markers were able to effectively stratify the breeds. Different traits under selective sweeps (natural or artificial selection) can be exploited for further improvement.

Population structure and identification of genomic regions associated with productive traits in five Italian beef cattle breeds

Italy has a long history in beef production, with local breeds such as Marchigiana, Chianina, Romagnola, Maremmana, and Podolica which produce high-quality meat. Selection has improved meat production, precocity, growth ability and muscle development, but the genetic determinism of such traits is mostly unknown. Using 33K SNPs-data from young bulls (N = 4064) belonging to these five Italian breeds, we demonstrated that the Maremmana and Podolica rustic breeds are closely related, while the specialised Marchigiana, Chianina, and Romagnola breeds are more differentiated. A genome-wide association study for growth and muscle development traits (average daily gain during the performance test, weight at 1 year old, muscularity) was conducted in the five Italian breeds. Results indicated a region on chromosome 2, containing the myostatin gene (MSTN), which displayed significant genome-wide associations with muscularity in Marchigiana cattle, a breed in which the muscle hypertrophy phenotype is segregating. Moreover, a significant SNP on chromosome 14 was associated, in the Chianina breed, to muscularity. The identification of diverse genomic regions associated with conformation traits might increase our knowledge about the genomic basis of such traits in Italian beef cattle and, eventually, such information could be used to implement marker-assisted selection of young bulls tested in the performance test.

Genome-Wide Association Study of Body Conformation Traits in Tashi Goats (Capra hircus)

Identifying genetic markers of economically valuable traits has practical benefits for the meat goat industry. To better understand the genomic variations influencing body conformation traits, a genome-wide association study was performed on Tashi goats, an indigenous Chinese goat breed. A total of 155 Tashi goats were phenotyped for eight body conformation traits: body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Then, 100 Tashi goats were randomly selected for whole-genome sequencing and genotyped. We obtained 1676.4 Gb of raw data with an average sequencing depth of 6.2X. Clean reads were aligned to the ARS1.2 reference genome, and 11,257,923 single nucleotide polymorphisms (SNPs) were identified. The structure analysis showed that these Tashi goats were almost not genetically related. The 109, 20, 52, 14, 62, 51, 70, and 7 SNPs were significantly associated with body height, body length, chest depth, chest width, chest girth, rump width, rump height, and cannon bone circumference. Within the ±500 kb region of significant SNPs, 183 genes were annotated. The most significantly enriched KEGG pathway was "olfactory transduction", and the most significantly enriched gene ontology (GO) terms were "cellular process", "cellular anatomical entity", and "molecular transducer activity". Interestingly, we found several SNPs on chromosomes 10 and 11 that have been identified multiple times for all eight body conformation traits located in two fragments (114 kb and 1.03 Mb). In chr.10:25988403-26102739, the six SNPs were tightly linked, the TACTAG genotype was the highest at 91.8%, and the FNTB (Farnesyltransferase, CAAX Box Beta) and CHURC1 (Churchill Domain Containing 1) genes were located. In chr.11:88216493-89250659, ten SNPs were identified with several dependent linkage disequilibrium (LD) blocks, and seven related genes were annotated, but no significant SNP was located in them. Our results provide valuable biological information for improving growth performance with practical applications for genomic selection in goats.

Genome-Wide Association and Pathway Analysis of Carcass and Meat Quality Traits in Karachai Young Goats

Goats with diverse economic phenotypic traits play an important role in animal husbandry. However, the genetic mechanisms underlying complex phenotypic traits are unclear in goats. Genomic studies of variations provided a lens to identify functional genes. The work aimed to search for candidate genes related to body measurements and body weight of Karachai goats and develop an experimental PCR-RV test system for genotyping significant SNPs. Comparison of GWAS results for ages 4 and 8 months revealed 58 common SNPs for significant genotypes. 11 common SNPs were identified for body weight, 4 SNPs-for group of traits withers height, rump height, body length, 2 SNPs-for withers height and rump height, 1 SNP-for body length and chest depth. Structural annotation of genomic regions covering a window of ±0.20 Mb showed the presence of 288 genes; 52 of them had the described functions in accordance with gene ontology. The main molecular functions of proteins encoded by these genes are the regulation of transcription, cell proliferation, angiogenesis, body growth, fatty acid and lipid metabolism, nervous system development, and spermatogenesis. SNPs common to body weight and localized within a window of ±200 kb from the structural genes CRADD, HMGA2, MSRB3, FUT8, MAX, and RAB15 were selected to create a test system. The study of meat productivity after slaughter and chemical analysis of muscle tissue in Karachai goats at the age of 8 months of different genotypes according to the identified SNPs revealed that rs268269710 is the most promising for further research and use in breeding. The GG genotype is associated with a larger live weight of animals, a larger carcass yield, the content of the boneless part in it, and the ratio of protein and adipose tissue in meat preferred for dietary nutrition. These results will contribute to the genetic improvement of Karachai goats.

Genome-wide association analysis of nine reproduction and morphological traits in three goat breeds from Southern China

OBJECTIVE

This study aimed to investigate the significant single nucleotide polymorphisms (SNPs) and genes associated with nine reproduction and morphological traits in three breed populations of Chinese goats.

METHODS

The genome-wide association of nine reproduction and morphological traits (litter size, nipple number, wattle, skin color, coat color, black dorsal line, beard, beard length, and hind leg hair) were analyzed in three Chinese native goat breeds (n = 336) using an Illumina Goat SNP50 Beadchip.

RESULTS

A total of 17 genome-wide or chromosome-wide significant SNPs associated with one reproduction trait (litter size) and six morphological traits (wattle, coat color, black dorsal line, beard, beard length, and hind leg hair) were identified in three Chinese native goat breeds, and the candidate genes were annotated. The significant SNPs and corresponding putative candidate genes for each trait are as follows: two SNPs located on chromosomes 6 (CSN3) and 24 (TCF4) for litter size trait; two SNPs located on chromosome 9 (KATNA1) and 1 (UBASH3A) for wattle trait; three SNPs located on chromosome 26 (SORCS3), 24 (DYM), and 20 (PDE4D) for coat color trait; two SNPs located on chromosome 18 (TCF25) and 15 (CLMP) for black dorsal line trait; four SNPs located on chromosome 8, 2 (PAX3), 5 (PIK3C2G), and 28 (PLA2G12B and OIT3) for beard trait; one SNP located on chromosome 18 (KCNG4) for beard length trait; three SNPs located on chromosome 17 (GLRB and GRIA2), 28 (PGBD5), and 4 for hind leg hair trait. In contrast, there were no SNPs identified for nipple number and skin color.

CONCLUSION

The significant SNPs or genes identified in this study provided novel insights into the genetic mechanism underlying important reproduction and morphological traits of three local goat breeds in Southern China as well as further potential applications for breeding goats.

Whole-genome selective scans detect genes associated with important phenotypic traits in goat (Capra hircus)

Goats with diverse economic phenotypic traits play an important role in animal husbandry. However, the genetic mechanisms underlying complex phenotypic traits are unclear in goats. Genomic studies of variations provided a lens to identify functional genes. In this study, we focused on the worldwide goat breeds with outstanding traits and used whole-genome resequencing data in 361 samples from 68 breeds to detect genomic selection sweep regions. We identified 210-531 genomic regions with six phenotypic traits, respectively. Further gene annotation analysis revealed 332, 203, 164, 300, 205, and 145 candidate genes corresponding with dairy, wool, high prolificacy, poll, big ear, and white coat color traits. Some of these genes have been reported previously (e.g., KIT, KITLG, NBEA, RELL1, AHCY, and EDNRA), while we also discovered novel genes, such as STIM1, NRXN1, LEP, that may be associated with agronomic traits like poll and big ear morphology. Our study found a set of new genetic markers for genetic improvement in goats and provided novel insights into the genetic mechanisms of complex traits.

Identification of Genomic Regions and Candidate Genes Associated with Body Weight and Body Conformation Traits in Karachai Goats

The objective of this study was to identify the SNPs and candidate genes related to body weight and seven body conformation traits at the age of 8 months in the Russian aboriginal Karachai goats (n = 269) by conducting genome-wide association studies (GWAS), using genotypes generated by Goat SNP BeadChip (Illumina Inc., USA). We identified 241 SNPs, which were significantly associated with the studied traits, including 47 genome-wide SNPs (p < 10−5) and 194 suggestive SNPs (p < 10−4), distributed among all goat autosomes except for autosome 23. Fifty-six SNPs were common for two and more traits (1 SNP for six traits, 2 SNPs for five traits, 12 SNPs for four traits, 20 SNPs for three traits, and 21 SNPs for two traits), while 185 SNPs were associated with single traits. Structural annotation within a window of 0.4 Mb (±0.2 Mb from causal SNPs) revealed 238 candidate genes. The largest number of candidate genes was identified at Chr13 (33 candidate genes for the five traits). The genes identified in our study were previously reported to be associated with growth-related traits in different livestock species. The most significant genes for body weight were CRADD, HMGA2, MSRB3, MAX, HACL1 and RAB15, which regulate growth processes, body sizes, fat deposition, and average daily gains. Among them, the HMGA2 gene is a well-known candidate for prenatal and early postnatal development, and the MSRB3 gene is proposed as a candidate gene affecting the growth performance. APOB, PTPRK, BCAR1, AOAH and ASAH1 genes associated with withers height, rump height and body length, are involved in various metabolic processes, including fatty acid metabolism and lipopolysaccharide catabolism. In addition, WDR70, ZBTB24, ADIPOQ, and SORCS3 genes were linked to chest width. KCNG4 was associated with rump height, body length and chest perimeter. The identified candidate genes can be proposed as molecular markers for growth trait selection for genetic improvement in Karachai goats.

Genome-Wide Association Study of Body Conformation Traits by Whole Genome Sequencing in Dazu Black Goats

Simple Summary

Body conformation traits are economically important in the goat meat industry. Good growth performance in goats, including an accelerated growth rate, can improve carcass weight and meat yield. The identification of genetic variants associated with these traits provides a basis for the genetic improvement of growth performance. In this study, we measured six body conformation traits, including body height, body length, cannon circumference, chest depth, chest width, and heart girth. By a genome-wide association study of a Chinese meat goat breed, 53 significant single nucleotide polymorphisms and 42 candidate genes associated with these traits were detected. These findings improve our understanding of the genetic basis of body conformation traits in goats.

Abstract

Identifying associations between genetic markers and economic traits has practical benefits for the meat goat industry. To better understand the genomic regions and biological pathways contributing to body conformation traits of meat goats, a genome-wide association study was performed using Dazu black goats (DBGs), a Chinese indigenous goat breed. In particular, 150 DBGs were genotyped by whole-genome sequencing, and six body conformation traits, including body height (BH), body length (BL), cannon circumference (CC), chest depth (CD), chest width (CW), and heart girth (HG), were examined. In total, 53 potential SNPs were associated with these body conformation traits. A bioinformatics analysis was performed to evaluate the genes located close to the significant SNPs. Finally, 42 candidate genes (e.g., PSTPIP2, C7orf57, CCL19, FGF9, SGCG, FIGN, and SIPA1L) were identified as components of the genetic architecture underlying body conformation traits. Our results provide useful biological information for the improvement of growth performance and have practical applications for genomic selection in goats.

Association of milk production and udder type traits with polymorphism of phosphorylase kinase regulatory subunit alpha-2 gene in Sahiwal cattle

The study was conducted in Sahiwal cattle in tropical region in the Indian sub-continent where the production of dairy animals is yet to be optimised. The present investigation was executed to identify SNPs in PHKA2 gene and to explore its effect on udder type and milk production traits. The study was based on the hypothesis that the PHKA2 gene has highly variable exons that could be related with udder traits and eventually milk production. PHKA2 gene regulates glycogen phosphorylase a, a catalyst in breakdown of glycogen. Milk production traits were recorded; 9 udder type, 5 teat type and 8 visual traits were measured for 100 animals. Five highly variable targeted regions of PHKA2 gene were amplified using PCR and sequenced. The association analysis was carried out using general linear model (SAS) to study the fixed effect of genotype on studied traits. The synonymous type SNP g.124556852C>T was found for the first time in Sahiwal cattle and possibly associated with udder type traits. The genotypic frequencies with respect to targeted loci g.124556852C>T indicated that homozygote CC (0.58) were highest in our resource population. The chi-square 2 test showed an agreement to Hardy–Weinberg equilibrium. The association analysis revealed significant association of genotypes with udder width, udder balance and 305 days milk yield. The attempt to find significant association with the visual udder traits was also done, however no significant alliance was observed. Homozygote CC animal were desirable as they favoured the selection of animal with superior udder width, udder balance and 305 day’s Milk Yield values.