Relationship between genetic variants of POU1F1 , PROP1 , IGFBP3 genes and milk performance in Guanzhong dairy goats

Genomic analysis uncovers novel candidate genes related to adaptation to tropical climates and milk production traits in native goats

BACKGROUND

Since domestication, both evolutionary forces and human selection have played crucial roles in producing adaptive and economic traits, resulting in animal breeds that have been selected for specific climates and different breeding goals. Pakistani goat breeds have acquired genomic adaptations to their native climate conditions, such as tropical and hot climates. In this study, using next-generation sequencing data, we aimed to assess the signatures of positive selection in three native Pakistani goats, known as milk production breeds, that have been well adapted to their local climate.

RESULTS

To explore the genomic relationship between studied goat populations and their population structure, whole genome sequence data from native goat populations in Pakistan (n = 26) was merged with available worldwide goat genomic data (n = 184), resulting in a total dataset of 210 individuals. The results showed a high genetic correlation between Pakistani goats and samples from North-East Asia. Across all populations analyzed, a higher linkage disequilibrium (LD) level (- 0.59) was found in the Pakistani goat group at a genomic distance of 1 Kb. Our findings from admixture analysis (K = 5 and K = 6) showed no evidence of shared genomic ancestry between Pakistani goats and other goat populations from Asia. The results from genomic selection analysis revealed several candidate genes related to adaptation to tropical/hot climates (such as; KITLG, HSPB9, HSP70, HSPA12B, and HSPA12B) and milk production related-traits (such as IGFBP3, LPL, LEPR, TSHR, and ACACA) in Pakistani native goat breeds.

CONCLUSIONS

The results from this study shed light on the structural variation in the DNA of the three native Pakistani goat breeds. Several candidate genes were discovered for adaptation to tropical/hot climates, immune responses, and milk production traits. The identified genes could be exploited in goat breeding programs to select efficient breeds for tropical/hot climate regions.

Milk-y Way: the impact of single-nucleotide polymorphisms on milk production traits in Kilis dairy goats

This study aimed to investigate the impact of single-nucleotide polymorphisms (SNPs) on milk production traits in Kilis dairy goats by analyzing the genotypes of POU1F1, PRLR, β -Lg, GH1, and GH2 genes and their association with lactation milk yield (LMY), lactation length (LL) and average daily milk yield (ADMY). Blood samples were collected from 227 goats, and genotyping was performed using polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP). The results revealed that the frequencies of the genotypes varied among the genes. The polymorphisms were found to be significantly linked with milk production traits. These findings suggest that SNPs of POU1F1/AluI, PRLR/RsaI and β -Lg/SacII are significantly associated with LMY and that the POU1F1-TC genotype, PRLR-TT genotype and β -Lg-AB genotype are associated with higher LMY and ADMY. Additionally, the POU1F1-TC genotype was found to have a longer LL. However, no significant association was found between the GH1 and GH2 genotypes and LMY, LL and ADMY. Overall, this study provides valuable insights into the genetic factors influencing milk production traits in Kilis dairy goats, which can be utilized for the selection of high-yielding animals in breeding programs.

Copy number variation analysis reveals variants associated with milk production traits in dairy goats

Copy number variation (CNV) is a major type of genomic structural variation. We investigated their impacts on goat dairy traits using the CaprineSNP50 array. From 120 samples of five dairy goat breeds, we totally identified 42 CNVs ranging from 56,044 bp to 4,337,625 bp. We found significant associations between two CNVs (CNV5 and CNV25) and two milk production traits (mean of milk fat yield and mean of milk protein yield) after false discovery rate (FDR) correction (P < 0.05). CNV5 overlaps the ADAMTS20 gene, which is involved in the differentiation of mammary cell and plays a crucial role in lactogenic activity of bovine mammary epithelial cells. CNV25 overlaps with PAPPA2, which has been found to be associated with bovine reproduction and milk production traits. Our results revealed that CNVs overlapped with ADAMTS20 and PAPPA2 could be involved in goat dairy traits and function as candidate markers for further genetic selection.

Indel variants within the PRL and GHR genes associated with sheep litter size

Growth hormone and prolactin belong to the class of peptide hormones that have a wide range of regulatory functions. In this study, polymorphisms of growth hormone receptor (GHR) and prolactin (PRL) genes were analysed as candidate genes, which are responsible for the litter size in Australian White (AUW) sheep. According to the statistical analyses results, the polymorphism information content (PIC) values of the PRL-P1-ins-23 bp, GHR-P2-del-23 bp and GHR-P8-del-23 bp were 0.371, 0.366 and 0.375, respectively, which indicates the high genetic polymorphism in AUW sheep. Moreover, all indel loci are not conformed to the HWE (p < .05). Further, our findings revealed that the PRL-P1-ins-23 bp polymorphism in the ovine PRL gene was significantly related to the first parity litter size (p = .001) and the DD genotype displaying the highest genotypic mean. Meanwhile, the GHR-P2-del-23 bp and GHR-P8-23 bp indels in the ovine GHR gene were significantly correlated with first parity litter size (p < .05), and the individuals with the genotype II showed significantly higher litter size than others. Collectively, these results demonstrated that our findings could be useful for future sheep breeding strategies based on the molecular-assisted selection (MAS).

Associations between genetic variants of the POU1F1 gene and production traits in Saanen goats

This study was conducted to determine the polymorphisms of the POU1F1 gene and their relationships with milk yield and components, litter size, birth weight, and weaning weight in goats. For this purpose, a total of 108 Saanen goats from two different farms (Bornova and Manisa) were used as animal materials. Polymorphisms at the exon 6 and the 3' flanking region of the POU1F1 gene were determined by using PCR-RFLP with PstI and AluI restriction enzymes and DNA sequencing analyses. Two alleles and three genotypes were identified by AluI or PstI digestions of the POU1F1 gene. The genotypes frequencies of TT, TC, and CC were 64.8 %, 31.5 % and 3.7 % for the PstI locus; 54.6 %, 31.5 % and 13.9 % for the AluI locus, respectively. T allele frequencies (0.56 and 0.88 for the AluI locus, 0.80 and 0.81 for the PstI locus, respectively) were predominant in both loci at the Bornova and Manisa farms. In terms of POU1F1-AluI and POU1F1-PstI loci, two populations were found to be in Hardy-Weinberg equilibrium. In the POU1F1-AluI locus, significant associations were found between genotypes and lactation milk yield and litter size. Similarly, a significant relationship between genotypes and birth weight in the POU1F1-PstI locus ( p < 0.05 ) was determined. The TC and CC genotypes were observed to be higher than the TT genotype for lactation milk yield and litter size at the POU1F1-AluI locus. Birth weight was found to be higher in animals that have the CC genotype at the POU1F1-PstI locus. In conclusion, the POU1F1 gene can be used as a molecular marker for economic features like reproduction, growth, milk content and yield in Saanen goats.

Relationship between SNPs of POU1F1 Gene and Litter Size and Growth Traits in Shaanbei White Cashmere Goats

POU (Pit-Oct-Unc) class 1 homeobox 1 (POU1F1, or Pit-1) is a transcription factor that directly regulates pituitary hormone-related genes, as well as affects the reproduction and growth in mammals. Thus, POU1F1 gene was investigated as a candidate gene for litter size and growth performance in goats. In the current study, using direct DNA sequencing, c.682G > T, c.723T > G and c.837T > C loci were genotyped in Shaanbei white cashmere (SBWC) goats (n = 609), but c.876 + 110T > C was monomorphic. Besides, the c.682G > T locus was first identified by HinfI (Haemophilus influenzae Rf) restriction endonuclease. Association analysis results showed that the c.682G > T, c.837T > C loci and diplotypes were significantly associated with goat litter size (p < 0.05). The positive genotypes were GT and TT for the two SNPs, respectively, and the optimal diplotype was H3H7 (GTTT-TTTT). On the other hand, the c.682G > T, c.723T > G and c.837T > C strongly affected growth traits and body measurement indexes in SBWC goats (p < 0.05). The positive genotypes or allele of these SNPs were GT, G and TT, respectively. Additionally, the goats with H3H7 diplotype also had a greater growth status than others (p < 0.05). Here, individuals with same genotype had both a better litter size and growth traits, showing a positive correlation between these economic traits. Meanwhile, the positive genotypes of four SNPs were combined to obtain the optimal diplotype, which was also H3H7. These SNPs, especially the diplotype, could be used for the genomic selection of excellent individuals with a greater litter size and better growth status in goat breeding.

The novel 22 bp insertion mutation in a promoter region of the <i>PITX2</i> gene is associated with litter size and growth traits in goats

Abstract. The paired-like homeodomain 2 (PITX2) gene plays a critical role in regulating development, reproduction, and growth traits in ruminants. Hence, the objective of this study was to explore the polymorphisms of this gene and to evaluate their associations with quantitative traits. Herein, a novel insertion in the promoter region of the PITX2 gene was reported in Shaanbei white cashmere (SBWC) goats (n=1012). The genotype distributions between mothers of single-kid and multi-kid groups within SBWC goats were significantly different (P<0.01), implying that this indel mutation might affect the litter size. Furthermore, association analysis found that this indel mutation was significantly associated with litter size (P=0.001). Individuals with genotype DD had a significantly smaller litter size than those with other genotypes (P<0.01). Besides, this indel was significantly associated with the body length (P=0.042) and the chest width (P=0.031). Especially, the individuals with genotype DD had a significantly lower body length than those with genotype II (P<0.05), which was consistent with the trend in litter size. These findings suggested that the new 22 bp indel mutation within the PITX2 gene is significantly associated with litter size and growth traits; this can be utilized as a functional molecular marker in goat breeding.

Kompetitive Allele Specific PCR (KASP™) genotyping of 48 polymorphisms at different caprine loci in French Alpine and Saanen goat breeds and their association with milk composition

Using a novel and fast genotyping method called Kompetitive Allele Specific PCR (KASP™), we carried out a pilot study on 48 single nucleotide polymorphisms (SNPs) belonging to 40 genes in French Alpine (n = 24) and Saanen (n = 25) goats reared in Romania. Furthermore, the associations of the 13 polymorphic genetic variants with milk production and composition were investigated. Thirty-five SNPs did not show polymorphism in the studied populations. Polymorphic SNPs were detected in the following genes: CAST, CLEC4E, DES, GHRHR, HSP90AA1, IL15RA, IL1RN, IL8, MITF, PPRC1, SOCS3, TNF and TNFSF13. The studied Alpine population was in Hardy-Weinberg disequilibrium at the g.62894878A>G locus (rs671391101) (P < 0.05). The results showed that four SNPs rs671391101 (GHRHR), rs640582069 (IL1RN) rs635583012 (SOCS3) and rs635969404 (IL15RA) out of the 13 polymorphic markers were significantly associated with milk production, protein, fat and lactose content in the Alpine breed. However, no significant effect was recorded in the Saanen population regarding milk yield or milk chemical composition. The current results provide new insights for the development of SNP marker-assisted selection technology in the goat industry and confirm the potential of using SNPs for the GHRHR, IL1RN, SOCS3, and IL15RA genes as candidate genes for selection, highlighting the direct implications of such genes for farm production outputs. The results from this study are relevant for future goat genomic studies and the inclusion of the associated traits into up-to-date selection schemes.