Association of GHR Polymorphisms with Milk Production in Buffaloes

Whole-Genome Resequencing Analysis of the Camelus bactrianus (Bactrian Camel) Genome Identifies Mutations and Genes Affecting Milk Production Traits

Milk production is an important trait that influences the economic value of camels. However, the genetic regulatory mechanisms underlying milk production in camels have not yet been elucidated. We aimed to identify candidate molecular markers that affect camel milk production. We classified Junggar Bactrian camels (9-10-year-old) as low-yield (<1.96 kg/d) or high-yield (>2.75 kg/d) based on milk production performance. Milk fat (5.16 ± 0.51 g/100 g) and milk protein (3.59 ± 0.22 g/100 g) concentrations were significantly lower in high-yielding camels than those in low-yielding camels (6.21 ± 0.59 g/100 g, and 3.93 ± 0.27 g/100 g, respectively) (p < 0.01). There were no apparent differences in gland tissue morphology between the low- and high-production groups. Whole-genome resequencing of 12 low- and 12 high-yield camels was performed. The results of selection mapping methods, performed using two methods (FST and θπ), showed that 264 single nucleotide polymorphism sites (SNPs) overlapped between the two methods, identifying 181 genes. These genes were mainly associated with the regulation of oxytocin, estrogen, ErbB, Wnt, mTOR, PI3K-Akt, growth hormone synthesis/secretion/action, and MAPK signaling pathways. A total of 123 SNPs were selected, based on significantly associated genomic regions and important pathways for SNP genotyping, for verification in 521 additional Bactrian camels. This analysis showed that 13 SNPs were significantly associated with camel milk production yield and 18 SNPs were significantly associated with camel milk composition percentages. Most of these SNPs were located in coding regions of the genome. However, five and two important mutation sites were found in the introns of CSN2 (β-casein) and CSN3 (κ-casein), respectively. Among the candidate genes, NR4A1, ADCY8, PPARG, CSN2, and CSN3 have previously been well studied in dairy livestock. These observations provide a basis for understanding the molecular mechanisms underlying milk production in camels as well as genetic markers for breeding programs aimed at improving milk production.

Detection of Polymorphisms in FASN, DGAT1, and PPARGC1A Genes and Their Association with Milk Yield and Composition Traits in River Buffalo of Bangladesh

This study aimed to identify SNPs in the intron, exon, and UTR regions of the FASN, DGAT1, and PPARGC1A genes and to investigate their possible association with milk yield and composition traits in the riverine buffalo of Bangladesh. A total of 150 DNA samples from riverine buffalo were used for PCR amplification with five pairs of primers, followed by association studies using a generalized linear model in R. SNP genotyping was performed by direct sequencing of the respective amplicon. Traits analyzed included DMY, fat%, protein%, and SNF%. This study identified 8 SNPs in FASN (g.7163G>A and g.7271C>T), DGAT1 (g.7809C>T and g.8525C>T) and PPARGC1A (g.387642C>T, g.387758A>G, g.409354A>G, and g.409452G>A). Genotypic and allelic frequencies differed significantly for each SNP genotype and did not follow the Hardy-Weinberg principle (p < 0.01 or p < 0.001) in most cases. The g.7163G>A and g.7271C>T SNP genotypes of the FASN gene were significantly associated with milk fat%, with the latter also significantly associated with SNF%. The g.8525C>T polymorphism of the DGAT1 gene significantly affected protein% (p < 0.01). Additionally, PPARGC1A gene polymorphisms showed significant associations: g.387642C>T with fat% (p < 0.05); g.387758A>G and g.409354A>G with protein% (p < 0.001) and SNF% (p < 0.01); and g.409452G>A with DMY (p < 0.001), fat% (p < 0.05), and protein% (p < 0.01). Reconstructed haplotypes of the PPARGC1A gene were significantly associated (p < 0.01) with all traits except SNF%. These findings suggest that polymorphisms in these three candidate genes have the potential as molecular markers for improving milk yield and composition traits in the riverine buffalo of Bangladesh.

DNA methylation of GDF-9 and GHR genes as epigenetic regulator of milk production in Egyptian Zaraibi goat

BACKGROUND

DNA methylation is an epigenetic mechanism that takes place at gene promoters and a potent epigenetic marker to regulate gene expression.

OBJECTIVE

The study aimed to improve the milk production of Zaraibi goats by addressing the methylation pattern of two milk production-related genes: the growth hormone receptor (GHR) and the growth differentiation factor-9 (GDF-9).

METHODS

54 and 46 samples of low and high milk yield groups, respectively, were collected. Detection of methylation was assessed in two CpG islands in the GDF-9 promoter via methylation-specific primer assay (MSP) and in one CpG island across the GHR promoter using combined bisulfite restriction analysis (COBRA).

RESULTS

A positive correlation between the methylation pattern of GDF-9 and GHR and their expression levels was reported. Breeding season was significantly effective on both peak milk yield (PMY) and total milk yield (TMY), where March reported a higher significant difference in PMY than November. Whereas single birth was highly significant on TMY than multiple births. The 3rd and 4th parities reported the highest significant difference in PMY, while the 4th parity was the most effective one on TMY.

CONCLUSION

These results may help improve the farm animals' milk productive efficiency and develop prospective epigenetic markers to improve milk yield by epigenetic marker-assisted selection (eMAS) in goat breeding programs.

TWAS revealed significant causal loci for milk production and its composition in Murrah buffaloes

Milk yield is the most complex trait in dairy animals, and mapping all causal variants even with smallest effect sizes has been difficult with the genome-wide association study (GWAS) sample sizes available in geographical regions with small livestock holdings such as Indian sub-continent. However, Transcriptome-wide association studies (TWAS) could serve as an alternate for fine mapping of expression quantitative trait loci (eQTLs). This is a maiden attempt to identify milk production and its composition related genes using TWAS in Murrah buffaloes (Bubalus bubalis). TWAS was conducted on a test (N = 136) set of Murrah buffaloes genotyped through ddRAD sequencing. Their gene expression level was predicted using reference (N = 8) animals having both genotype and mammary epithelial cell (MEC) transcriptome information. Gene expression prediction was performed using Elastic-Net and Dirichlet Process Regression (DPR) model with fivefold cross-validation and without any cross-validation. DPR model without cross-validation predicted 80.92% of the total genes in the test group of Murrah buffaloes which was highest compared to other methods. TWAS in test individuals based on predicted gene expression, identified a significant association of one unique gene for Fat%, and two for SNF% at Bonferroni corrected threshold. The false discovery rates (FDR) corrected P-values of the top ten SNPs identified through GWAS were comparatively higher than TWAS. Gene ontology of TWAS-identified genes was performed to understand the function of these genes, it was revealed that milk production and composition genes were mainly involved in Relaxin, AMPK, and JAK-STAT signaling pathway, along with CCRI, and several key metabolic processes. The present study indicates that TWAS offers a lower false discovery rate and higher significant hits than GWAS for milk production and its composition traits. Hence, it is concluded that TWAS can be effectively used to identify genes and cis-SNPs in a population, which can be used for fabricating a low-density genomic chip for predicting milk production in Murrah buffaloes.

Detection of genome-wide copy number variation in Murrah buffaloes

Riverine Buffaloes, especially the Murrah breed because of their adaptability to harsh climatic conditions, is farmed in many countries to convert low-quality feed into valuable dairy products and meat. Here, we investigated the copy number variations (CNVs) in 296 Murrah buffalo using the Axiom® Buffalo Genotyping Array 90K (Affymetrix, Santa Clara, CA, USA). The CNVs were detected on the autosomes, using the Copy Number Analysis Module (CNAM) using the univariate analysis. 7937 CNVs were detected in 279 Buffaloes, the average length of the CNVs was 119,048.87 bp that ranged between 7800 and 4,561,030 bp. These CNVs were accounting for 10.33% of the buffalo genome, which was comparable to cattle, sheep, and goat CNV analyses. Further, CNVs were merged and 1541 CNVRs were detected using the Bedtools-mergeBed command. 485 genes were annotated within 196 CNVRs that were identified in at least 10 animals of Murrah population. Out of these, 40 CNVRs contained 59 different genes that were associated with 69 different traits. Overall, the study identified a significant number of CNVs and CNVRs in the Murrah breed of buffalo, with a wide range of lengths and frequencies across the autosomes. The identified CNVRs contained genes associated with important traits related to production and reproduction, making them potentially important targets for future breeding and genetic improvement efforts.

Whole-genome resequencing reveals genetic diversity and selection characteristics of dairy goat

The dairy goat is one of the earliest dairy livestock species, which plays an important role in the economic development, especially for developing countries. With the development of agricultural civilization, dairy goats have been widely distributed across the world. However, few studies have been conducted on the specific characteristics of dairy goat. In this study, we collected the whole-genome data of 89 goat individuals by sequencing 48 goats and employing 41 publicly available goats, including five dairy goat breeds (Saanen, Nubian, Alpine, Toggenburg, and Guanzhong dairy goat; n = 24, 15, 11, 6, 6), and three goat breeds (Guishan goat, Longlin goat, Yunshang Black goat; n = 6, 15, 6). Through compared the genomes of dairy goat and non-dairy goat to analyze genetic diversity and selection characteristics of dairy goat. The results show that the eight goats could be divided into three subgroups of European, African, and Chinese indigenous goat populations, and we also found that Australian Nubian, Toggenburg, and Australian Alpine had the highest linkage disequilibrium, the lowest level of nucleotide diversity, and a higher inbreeding coefficient, indicating that they were strongly artificially selected. In addition, we identified several candidate genes related to the specificity of dairy goat, particularly genes associated with milk production traits (GHR, DGAT2, ELF5, GLYCAM1, ACSBG2, ACSS2), reproduction traits (TSHR, TSHB, PTGS2, ESR2), immunity traits (JAK1, POU2F2, LRRC66). Our results provide not only insights into the evolutionary history and breed characteristics of dairy goat, but also valuable information for the implementation and improvement of dairy goat cross breeding program.

Comparative analyses of copy number variations between swamp and river buffalo

Domestic buffalo is an important livestock in the tropical and sub-tropical region, including two types: swamp and river buffalo. The swamp buffalo is mainly used as draft animal, while the river buffalo is raised for milk production. In this study, based on the new high-quality buffalo reference genome UOA_WB_1, we firstly investigated the copy number variants in buffalo using whole-genome Illumina sequencing. A total of 3,734 CNV regions (CNVRs) were detected in 106 buffalo population with a total length of 23,429,066 bp, corresponding to ∼ 0.88% of the water buffalo genome (UOA_WB_1). Our results revealed a clear population differentiation in CNV between swamp and river buffalo. In addition, a total of 667 highly differentiated CNVRs (covering 886 genes) were detected between river and swamp buffalo population. We detected a set of CNVR-overlapping genes associated with exercise, immunity, nerve, and milk trait which exhibited different copy numbers between swamp and river buffalo population. This study provides valuable genome variation resources for buffalo and would contribute to understanding the genetic differences between swamp and river buffalo.

Meprin A1 subunit beta gene polymorphism is associated with the length of post-partum anestrus interval in Murrah buffaloes

Postpartum anestrus interval (PPAI) is the interval between parturition and the first postpartum estrus exhibition in animals. Appearance of both normal and PPA buffaloes under the same farm conditions indicates the role of possible genetic predisposition to PPA. To identify the genetic and non-genetic factors associated with PPA in buffaloes, we collected data on PPAI and other 38 non-genetic variables from 575 Murrah buffaloes in the field conditions and identified the PPA associated non-genetic factors in our previous study. To explore the genetic factors associated with the unexplained variation in PPAI residuals, the present study identified 41 single nucleotide polymorphisms (SNPs) in 13 candidate genes using Sanger sequencing. Exploration of their association with the PPAI residuals of 50 extreme PPA and 50 normal buffaloes identified the significant (P < 0.01) association of the SNP (g.37219977A>G) in the 3'-UTR region of the Meprin A 1 subunit beta (Mep1b) gene with PPAI, which was further validated (P = 0.058) in a large population sample (n = 417). Bioinformatics analysis of the 3'-UTR region has identified three miRNA, bta-miR-2420, bta-miR-2325b and bta-miR-453 that could regulate Igf-1 in the plasma of animals with different genotypes (GG, AG and AA). The higher Igf-1 levels in the GG genotypes than that of AA and AG genotypes of this SNP (g.37219977A>G) further suggest the association of Mep1b gene with PPA condition in Murrah buffaloes. As a result of this study, we propose that buffaloes with protective alleles at this SNP be selected to improve the herd's reproductive efficiency.

Exploring the relationship between polymorphisms of leptin and IGF-1 genes with milk yield in indicine and taurine crossbred cows

Leptin an adipokine plays a significant role in several physiological processes and has been indicated as a candidate gene for marker-assisted selection for high-yielding cows. Insulin-like growth factor-1 (IGF-1) hormone plays an important physiological role in growth, development, metabolism, and lactation in bovines. It is believed to be one of the main mediators of energy balance effects on the reproductive performance of dairy cows after calving. The present investigation was carried out to identify the novel polymorphisms in exon 3 region of leptin and exon 3 partial intron 3 of IGF-1 genes and their association with the milk production performance in indicine and taurine crossbred (Karan Fries) cows. Blood samples were collected from 160 apparently healthy Karan Fries (KF) cows. Four SNPs (single-nucleotide polymorphisms) at positions rs29004508 (C > T), rs29004509 (C > T), rs29004510 (T > C), and rs29004511 (T > C) in leptin gene and two SNPs at positions rs133251968 (C > A) and rs137289661 (C > T) in IGF-1 gene were found in KF cows; however, rs29004509 (C > T) had a positive correlation (r = 0.376; P < 0.05) with milk yield. The genetic variants observed in exon 3 region of leptin gene and their association with milk yield traits revealed the importance of CT genotype, which had been useful for genetic improvement of KF cow for milk production traits and can also be utilized as a potential genetic marker to select appropriate animals.

Polymorphisms of the PRLR Gene and Their Association with Milk Production Traits in Egyptian Buffaloes

Simple Summary

The two non-synonymous g.11685G>A and g.11773T>C SNPs of PRLR(L2) were significantly associated with milk yield, fat%, and protein%, and mRNA and protein levels of PRL and PRLR in milk somatic cells. GT-animals had the best milk performance; however, AC-animals had inferior milk production. Thus, the selection of buffaloes with GT haplotypes may enhance milk performance in Egyptian buffaloes.

Abstract

Prolactin (PRL) and its receptor (PRLR) were considered as potential genetic markers for milk production and quality traits in cattle. However, little information is available regarding PRLR genetic diversity and association studies with milk traits in Egyptian water buffaloes. Therefore, the present study was conducted to search for mutations in PRLR and determine their associations with milk performance in these animals. Exon3 (E3) and E10 of PRLR were screened for polymorphisms using single strand conformation polymorphism (SSCP) and sequencing in 400 buffaloes. The associations between haplotypes and milk production (fat%, protein%, lactose%, and solid%) traits as well as mRNA and protein levels of PRL and PRLR were studied. Two single nucleotide polymorphisms (SNPs) in E10 were detected: g.11685G>A (p.Ala494Thr) and g.11773T>C (p.Val523Aal). The G and T alleles were wild (ancestral) alleles, while the A and C alleles were mutant alleles. These SNPs resulted in four haplotypes; AC, AT, GC, and GT. Buffaloes with wild GT haplotypes showed significantly higher milk yield, fat% and protein%, mRNA and protein levels of PRL and PRLR in milk somatic cells than other animals. Animals carrying mutant AC haplotype had inferior milk traits and lowest levels of associated mRNAs and proteins. With these results, we could conclude that the selection of buffaloes with wild GT haplotypes for g.11685G>A and g.11773T>C SNPs of the PRLR gene might improve the milk production traits of Egyptian water buffaloes.

Whole-Genome Sequencing and Characterization of Buffalo Genetic Resources: Recent Advances and Future Challenges

The buffalo was domesticated around 3000-6000 years ago and has substantial economic significance as a meat, dairy, and draught animal. The buffalo has remained underutilized in terms of the development of a well-annotated and assembled reference genome de novo. It is mandatory to explore the genetic architecture of a species to understand the biology that helps to manage its genetic variability, which is ultimately used for selective breeding and genomic selection. Morphological and molecular data have revealed that the swamp buffalo population has strong geographical genomic diversity with low gene flow but strong phenotypic consistency, while the river buffalo population has higher phenotypic diversity with a weak phylogeographic structure. The availability of recent high-quality reference genome and genotyping marker panels has invigorated many genome-based studies on evolutionary history, genetic diversity, functional elements, and performance traits. The increasing molecular knowledge syndicate with selective breeding should pave the way for genetic improvement in the climatic resilience, disease resistance, and production performance of water buffalo populations globally.

Opportunities and Challenges for Improving the Productivity of Swamp Buffaloes in Southeastern Asia

The swamp buffalo is a domesticated animal commonly found in Southeast Asia. It is a highly valued agricultural animal for smallholders, but the production of this species has unfortunately declined in recent decades due to rising farm mechanization. While swamp buffalo still plays a role in farmland cultivation, this species' purposes has shifted from draft power to meat, milk, and hide production. The current status of swamp buffaloes in Southeast Asia is still understudied compared to its counterparts such as the riverine buffaloes and cattle. This review discusses the background of swamp buffalo, with an emphasis on recent work on this species in Southeast Asia, and associated genetics and genomics work such as cytogenetic studies, phylogeny, domestication and migration, genetic sequences and resources. Recent challenges to realize the potential of this species in the agriculture industry are also discussed. Limited genetic resource for swamp buffalo has called for more genomics work to be done on this species including decoding its genome. As the economy progresses and farm mechanization increases, research and development for swamp buffaloes are focused on enhancing its productivity through understanding the genetics of agriculturally important traits. The use of genomic markers is a powerful tool to efficiently utilize the potential of this animal for food security and animal conservation. Understanding its genetics and retaining and maximizing its adaptability to harsher environments are a strategic move for food security in poorer nations in Southeast Asia in the face of climate change.

Association of GH polymorphisms with growth traits in buffaloes

Members of the somatotrophic axis, especially GH and IGF1, are essential for growth. The association between GH polymorphisms and growth traits was numerously studied in cattle; however, no data are available for such association studies in buffalo. Therefore, this study was conducted to screen for polymorphisms in the GH gene and to study their putative association with growth traits in 200 Egyptian buffaloes. Polymerase chain reaction single-strand conformation polymorphism and sequencing were applied to look for polymorphisms in 3 loci spanning all exons and introns of buffalo GH. The C (MspI+) >T (MspI-) SNP in intron3, which is well known in cattle, was not detected in the examined buffaloes. However, 2 missense mutations were detected in exon5: one previously detected p.Leu153Val SNP, with very low frequencies for the mutant (Val) allele and one novel p.Asn174His SNP. At weaning age, the p.Leu153Val SNP was significantly associated with weaning body weight and gain with the positive effect of the wild allele (Leu) and higher GH serum, mRNA, and protein levels in animals with Leu/Asn and Leu/His haplotypes. At yearling age, the 2 SNPs associated with yearling weight and gain with positive effect for the mutant (Val and His) alleles with increased GH, and IGF1 serum, mRNA, and protein and GHR mRNA and protein levels in animals with Val/Asn haplotype. Therefore, the selection of Egyptian buffaloes with the Val/Asn haplotype could improve the growth traits of Egyptian buffaloes at yearling age which is the target age for perfect growing.