DGAT1, GH, GHR, PRL and PRLR Polymorphism in Water Buffalo (Bubalus bubalis)

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.

The Role of PRLR Gene Polymorphisms in Milk Production in European Wild Rabbit (Oryctolagus cuniculus)

One of the problematic points of rabbit breeding is that the nutritional requirements of the kits are not fully satisfied by the does' milk production from the third week of lactation onwards. The prolactin receptor gene has a significant effect on reproductive processes, and its polymorphisms have been associated with milk production in several species (cattle, goats, sheep, and buffalo). The European wild rabbit (Oryctolagus cuniculus), has a more diverse genetic background compared to domesticated lines. In the course of our study, sequencing of the 1210 bp long segment of the PRLR gene promoter region was accomplished. We detected four point mutations (SNP1-407G > A, SNP2-496G > C, SNP3-926T> and SNP4-973A > C) and one microsatellite at position 574. In our population, the four SNPs were segregated into four genotypes: AACCCCCC, GGGGTTAA, AAGGTTAC, and GGGGTCAC. Our results show that the genotype in the homozygous form is associated with higher milk production (1564.7 ± 444.7 g) compared to the other three genotypes (AACCCCCC 1399.1 ± 326.8 g; GTGACCTT 1403.8 ± 517.1 g; GGGGTCAC 1220.0 ± 666.2 g), and the short microsatellite repeat (167 bp) also coincides with significantly higher milk production (1623.8 ± 525.1 g). These results make the marker-assisted selection possible also for domesticated lines.

Association of PRLR, IGF1, and LEP genes polymorphism with milk production and litter size in Egyptian Zaraibi goat

Studying variation in genes responsible for physiological characters is important to enhance goat productive and reproductive efficiency. This study aimed to detect specific nucleotide polymorphisms in prolactin receptor (PRLR), insulin-like growth factor (IGF1), and leptin (LEP) genes and their correlation with milk production (MP) and litter size (LS) traits in Zaraibi goat. PCR-SSCP products of different patterns of each gene were sequenced and aligned to reveal two mutations (T > C) and (G > A) in 3′UTR of PRLR gene and registered on NCBI with accession numbers OM418863 for TT and OM418864 for CT, while (G > A) variation was registered as OM418861 for GG and OM418862 for AG in exon 10. TT, CT, AG, and GG genotypes were distributed in the studied animals with frequencies 0.43, 0.57, 0.65, and 0.35, respectively. While alleles C, T, A, and G frequencies were 0.28, 0.72, 0.32, and 0.68, respectively. CT and AG genotypes associated significantly (P < 0.05) with higher MP and LS, respectively. By studying the haplotypes of PRLR, C-A and T-A were associated with the highest and the lowest level of MP, respectively. For LS, T-A and C-G showed significant correlation with the highest and the lowest rate, respectively. Regarding IGF1 gene, two polymorphisms were detected; T74C at exon 4 which registered on NCBI as OM418860, and combined mutations as ins. G470, A531G, and T534C (PP genotype) at 5′ flanking region that registered as OM418859. For LEP, only one polymorphism was found in intron 2 (G281A) which submitted to NCBI as OM418855. All detected polymorphisms have shown to be involved in regulating the MP or LS as reproductive traits in goat.

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.

Association analysis of prolactin and prolactin receptor genes with selected productive and reproductive traits in Egyptian buffalo

A total of 266 records of buffalo raised in two experimental herds in Egypt were assessed to detect prolactin (PRL) and prolactin receptor (PRLR) genes' polymorphism using PCR-Single Strand Conformational Polymorphism (SSCP) and PCR-Restricted Fragment Length Polymorphism (RFLP) techniques as well as to investigate their association with calf birth weight (BW), weaning weight (WW), lactation period (LP), total milk yield (TMY), stillbirth, calving ease (CE), gestation length (GL), postpartum interval to pregnancy (PPIP), calving interval (CI), and age at first calving (AFC). Predicted breeding values were estimated and used in the association with detected genotypes. A monomorphic pattern of the studied PRL 156 bp segment was recorded and absence of its polymorphism in buffalo was corroborated. We also determined polymorphism of PRLR reflected in three loci: PRLR2, PRLR4, and PRLR9. Significant differences among PRLP9 genotypes (AA, AB, and BB) were displayed for all studied traits as well as among PRLR2 genotypes, except for CE, while PRLR4 genotypes significantly differed only in BW, WW, TMY, stillbirth, GL, and AFC. In practice, strong associations among genotypes of the PRLR gene and the traits of interest candidate this gene to be selective in Egyptian buffalo breeding for improving both productive and reproductive traits.

Polymorphism detection of <i>DGAT1</i> and <i>Lep</i> genes in Anatolian water buffalo (<i>Bubalus bubalis</i>) populations in Turkey

Acyl-CoA: diacylglycerol–acyltransferase 1 (DGAT1) enzyme plays a key role in controlling the synthesis rate triglyceride from diacylglycerol. Leptin (LP, OB, obese) is an important hormone that synthesizes mostly from adipose tissue and regulates glucose metabolism and homeostasis. DGAT1 and Lep genes are closely related to reproduction, growth, milk yield and composition in water buffalo breeds. This study aimed to identify genetic variation in the DGAT1 and Lep gene regions in 150 water buffalo individuals from five different provinces of Turkey using DNA sequencing. A total of 38 nucleotide variations and indels have identified 761 bp long partial intron 2 and exon 3 and 5′ UTR regions of the Lep gene in Anatolian water buffalo populations; 422 bp long partial exon 7–9 and exon 8 regions of DGAT1 gene were amplified and two mutations were defined in the point of 155 and 275 nucleotide that is three genotypes for S allele and Y allele of DGAT1 gene in intron 7 in Anatolian buffalo populations, respectively. These SNPs may have an effect on reproduction, growth, milk yield and composition in water buffalo populations and may prove to be useful for water buffalo breeding.

Signatures of selection reveal candidate genes involved in production traits in Chinese crossbred buffaloes

Identification of selection signature is important for a better understanding of genetic mechanisms that affect phenotypic differentiation in livestock. However, the genome-wide selection responses have not been investigated for the production traits of Chinese crossbred buffaloes. In this study, an SNP data set of 133 buffaloes (Chinese crossbred buffalo, n = 45; Chinese local swamp buffalo, n = 88) was collected from the Dryad Digital Repository database (https://datadryad.org/stash/). Population genetics analysis showed that these buffaloes were divided into the following 2 groups: crossbred buffalo and swamp buffalo. The crossbred group had higher genetic diversity than the swamp group. Using 3 complementary statistical methods (integrated haplotype score, cross population extended haplotype homozygosity, and composite likelihood ratio), a total of 31 candidate selection regions were identified in the Chinese crossbred population. Here, within these candidate regions, 25 genes were under the putative selection. Among them, several candidate genes were reported to be associated with production traits. In addition, we identified 13 selection regions that overlapped with bovine QTLs that were mainly involved in milk production and composition traits. These results can provide useful insights regarding the selection response for production traits of Chinese crossbred buffalo, as identified candidate genes influence production performance.

Association of DGAT1 With Cattle, Buffalo, Goat, and Sheep Milk and Meat Production Traits

Milk fatty acids are essential for many dairy product productions, while intramuscular fat (IMF) is associated with the quality of meat. The triacylglycerols (TAGs) are the major components of IMF and milk fat. Therefore, understanding the polymorphisms and genes linked to fat synthesis is important for animal production. Identifying quantitative trait loci (QTLs) and genes associated with milk and meat production traits has been the objective of various mapping studies in the last decade. Consistently, the QTLs on chromosomes 14, 15, and 9 have been found to be associated with milk and meat production traits in cattle, goat, and buffalo and sheep, respectively. Diacylglycerol O-acyltransferase 1 (DGAT1) gene has been reported on chromosomes 14, 15, and 9 in cattle, goat, and buffalo and sheep, respectively. Being a key role in fat metabolism and TAG synthesis, the DGAT1 has obtained considerable attention especially in animal milk production. In addition to milk production, DGAT1 has also been a subject of interest in animal meat production. Several polymorphisms have been documented in DGAT1 in various animal species including cattle, buffalo, goat, and sheep for their association with milk production traits. In addition, the DGAT1 has also been studied for their role in meat production traits in cattle, sheep, and goat. However, very limited studies have been conducted in cattle for association of DGAT1 with meat production traits in cattle. Moreover, not a single study reported the association of DGAT1 with meat production traits in buffalo; thus, further studies are warranted to fulfill this huge gap. Keeping in view the important role of DGAT1 in animal production, the current review article was designed to highlight the major development and new insights on DGAT1 effect on milk and meat production traits in cattle, buffalo, sheep, and goat. Moreover, we have also highlighted the possible future contributions of DGAT1 for the studied species.

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.

Genetic Association of PPARGC1A Gene Single Nucleotide Polymorphism with Milk Production Traits in Italian Mediterranean Buffalo

PPARGC1A gene plays an important role in the activation of various important hormone receptors and transcriptional factors involved in the regulation of adaptive thermogenesis, gluconeogenesis, fiber-type switching in skeletal muscle, mitochondrial biogenesis, and adipogenesis, regulating the reproduction and proposed as a candidate gene for milk-related traits in cattle. This study identified polymorphisms in the PPARGC1A gene in Italian Mediterranean buffaloes and their associations to milk production and quality traits (lactation length, peak milk yield, fat and protein yield, and percentage). As a result, a total of seven SNPs (g.-78A>G, g.224651G>C, g.286986G>A, g.304050G>A, g.325647G>A, g.325817T>C, and g.325997G>A) were identified by DNA pooled sequencing. Analysis of productivity traits within the genotyped animals revealed that the g.286986G>A located at intron 4 was associated with milk production traits, but the g.325817T>C had no association with milk production. Polymorphisms in g.-78A>G was associated with peak milk yield and milk yield, while g.304050G>A and g.325997 G>A were associated with both milk yield and protein percentage. Our findings suggest that polymorphisms in the buffalo PPARGC1A gene are associated with milk production traits and can be used as a candidate gene for milk traits and marker-assisted selection in the buffalo breeding program.

Evaluation of non-synonym mutation in DGAT1 K232A as a marker for milk production traits in Ongole cattle and Murrah buffalo from Southern India

Various candidate genes have been reported to affect milk yield and composition in dairy cattle. A non-synonymous mutation in the DGAT1 gene, i.e., K232A was reported to have a strong association with milk yield and milk composition of Bos taurus. A study has been undertaken on 502 unrelated individuals belonging to indigenous Ongole cattle, crossbred cattle, and Murrah buffaloes from the Indian sub-continent with the objective to determine the polymorphism of the K232A locus and their association with milk yield and composition. Typing DGAT1 K232A allelic variation by PCR-RFLP using CfrI restriction enzyme revealed three genotypes in crossbred cattle. Genotype KK was more prevalent (0.60) in Jersey crossbred, whereas in Holstein Friesian crossbred it was KA genotype (0.48). In Ongole cattle and Murrah buffaloes, the locus did not exhibit polymorphism. The least-square mean of milk yields pooled over lactations across the DGAT1 variants was significantly (P < 0.05) higher among the homozygous (AA) genotypes, both in Jersey crossbred and HF crossbred cattle after adjusting for the effects of farm, parity, and season. The fat, SNF, and protein content values of AA genotypes were less than the KK genotypes in both the genetic groups (P > 0.05). The fixation of the DGAT1^K allele at the locus in Bos indicus cattle and Bubalus bubalis in the present study did not support its use as a reliable universal marker for milk production and composition traits.

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.

Association of GHR Polymorphisms with Milk Production in Buffaloes

Simple Summary

The present study reported two missense mutations in the buffalo GHR gene: A novel (c.380G>A) and (c.836T>A) which was described in previous studies. These two single nucleotide polymorphisms (SNPs) were found to be associated with milk yield, fat %, protein %, and 305 day-milk, fat and protein yield, with higher performance for AA haplotype animals. Therefore, selection of buffaloes with AA haplotype would more likely improve milk production traits. Consequently, this would allow breeders to take more precise selection decisions, leading to significantly higher productivity and profitability within the Egyptian buffalo herds.

Abstract

For its role in the mediation of growth hormone (GH) galactopoietic effect, growth hormone receptor (GHR) was considered a functional candidate gene for milk performance in cattle. However, its genetic variation and potential effect have not been investigated in Egyptian buffaloes. This study aimed to screen GHR for polymorphisms and study their associations with milk traits in Egyptian buffaloes. Polymerase chain reaction, single-strand conformation polymorphism, and sequencing were used to identify mutations in 4 exons (E4–E6 and E8) of the GHR gene in 400 Egyptian buffaloes. No polymorphisms were found in E4, while 2 SNPs (c.380G>A/p.Arg127Lys and c.387C>T/p.Gly129) in E5, one silent mutation (c.435A>G/p.Pro145) in E6, and another missense mutation (c.836T>A/p.Phe279Tyr) in E8 were detected. The c.380G>A SNP in the extracellular domain was associated with milk yield, fat %, protein %, and 305-day milk, fat and protein yield, with higher levels in animals carrying the mutant A allele. The c.836T>A SNP in the transmembrane domain was associated with milk yield, fat %, protein %, and 305-day milk, fat and protein yield, with higher milk yield and lower fat %, protein %, fat and protein yield in the mutant A allele-animals. Interestingly, animals with the two mutant AA alleles produced higher milk yield, fat %, protein %, fat and protein yield, accompanied with upregulated expressions of GHR, GH, insulin-like growth factor 1 (IGF1), prolactin (PRL), prolactin receptor (PRLR), β-casein (encoded by CSN2 gene), and diacylglycerol acyltransferase-1 (DGAT1) genes and proteins in milk somatic cells. Therefore, selection of Egyptian buffaloes with mutant AA haplotypes for the novel c.380G>A SNP and the well-known c.836T>A SNP could improve milk yield and quality in buffaloes.

Genome-wide identification of Diacylglycerol Acyltransferases (DGAT) family genes influencing Milk production in Buffalo

Background

The diacylglycerol acyltransferases (DGAT) are a vital group of enzymes in catalyzing triacylglycerol biosynthesis. DGAT genes like DGAT1 and DGAT2, have been identified as two functional candidate genes affecting milk production traits, especially for fat content in milk. Buffalo milk is famous for its excellent quality, which is rich in fat and protein content. Therefore, this study aimed to characterize DGAT family genes in buffalo and to find candidate markers or DGAT genes influencing lactation performance.

Results

We performed a genome-wide study and identified eight DGAT genes in buffalo. All the DGAT genes classified into two distinct clades (DGAT1 and DGAT2 subfamily) based on their phylogenetic relationships and structural features. Chromosome localization displayed eight buffalo DGAT genes distributed on five chromosomes. Collinearity analysis revealed that the DGAT family genes were extensive homologous between buffalo and cattle. Afterward, we discovered genetic variants loci within the genomic regions that DGAT genes located in buffalo. Seven haplotype blocks were constructed and were associated with buffalo milk production traits. Single marker association analyses revealed four most significant single nucleotide polymorphisms (SNPs) mainly affecting milk protein percentage or milk fat yield in buffalo. Genes functional analysis indicated that these DGAT family genes could influence lactation performance in the mammal through regulating lipid metabolism.

Conclusion

In the present study, we performed a comprehensive analysis for the DGAT family genes in buffalo, which including identification, structural characterization, phylogenetic classification, chromosomal distribution, collinearity analysis, association analysis, and functional analysis. These findings provide useful information for an in-depth study to determine the role of DGAT family gens play in the regulation of milk production and milk quality improvement in buffalo.

Systematic analyses for candidate genes of milk production traits in water buffalo (Bubalus Bubalis)

Water buffalo (Bubalus bubalis) is of great economic importance as a provider of milk and meat in many countries. However, the milk yield of buffalo is much lower than that of Holstein cows. Selection of candidate genes related to milk production traits can be applied to improve buffalo milk performance. A systematic review of studies of these candidate genes will be greatly beneficial for researchers to timely and efficiently understand the research development of molecular markers for buffalo milk production traits. Here, we identified and classified the candidate genes associated with buffalo milk production traits. A total of 517 candidate genes have been identified as being associated with milk performance in different buffalo breeds. Nineteen candidate genes containing 47 mutation sites have been identified using the candidate gene approach. In addition, 499 candidate genes have been identified in six genome-wide association studies (GWASes) including two studies performed with the bovine SNP chip and four studies with the buffalo SNP chip. Genes CTNND2 (catenin delta 2), APOB (apolipoprotein B), FHIT (fragile histidine triad) and ESRRG (estrogen related receptor gamma) were identified in at least two GWASes. These four genes, especially APOB, deserve further study to explore regulatory roles in buffalo milk production. With growth in the number of buffalo genomic studies, more candidate genes associated with buffalo milk production traits will be identified. Therefore, future studies, such as those investigating gene location and functional analyses, are necessary to facilitate the exploitation of genetic potential and the improvement of buffalo milk performance.

Remarkable genetic diversity detected at river buffalo prolactin receptor (PRLR) gene and association studies with milk fatty acid composition

Prolactin is an anterior pituitary peptide hormone involved in many different endocrine activities and is essential for reproductive performance. This action is mediated by its receptor, the prolactin receptor, encoded by the PRLR gene. In this study, we sequenced and characterized the Mediterranean river buffalo PRLR gene (from exon 3 to 10), and we found remarkable genetic diversity. In particular, we found 24 intronic polymorphisms and 13 exonic SNPs, seven of which were non-synonymous. Furthermore, the polymorphisms identified in the 3'-UTR were investigated to establish their possible influence on microRNA binding sites. Considering all the amino acid changes and the observed allelic combinations, it is possible to deduce at least six different translations of the buffalo prolactin receptor and, consequently, the presence at the PRLR gene of at least six alleles. Furthermore, we identified a deletion of a CACTACC heptamer between nucleotides 1102 and 1103 of exon 10 (3'-UTR), and we developed an allele-specific PCR to identify the carriers of this genetic marker. Finally, the SNP g.11188A>G, detected in exon 10 and responsible for the amino acid replacement p.His328Arg, was genotyped in 308 Italian Mediterranean river buffaloes, and an association study with milk fat traits was carried out. The statistical analysis showed a tendency that approached significance for the AA genotype with higher contents of odd branched-chain fatty acids. Thus, our results suggest that the PRLR gene is a good candidate for gene association studies with qualitative traits related to buffalo milk production.

Genome-wide association studies to identify quantitative trait loci affecting milk production traits in water buffalo

Water buffalo is the second largest resource of milk supply around the world, and it is well known for its distinctive milk quality in terms of fat, protein, lactose, vitamin, and mineral contents. Understanding the genetic architecture of milk production traits is important for future improvement by the buffalo breeding industry. The advance of genome-wide association studies (GWAS) provides an opportunity to identify potential genetic variants affecting important economical traits. In the present study, GWAS was performed for 489 buffaloes with 1,424 lactation records using the 90K Affymetrix Buffalo SNP Array (Affymetrix/Thermo Fisher Scientific, Santa Clara, CA). Collectively, 4 candidate single nucleotide polymorphisms (SNP) in 2 genomic regions were found to associate with buffalo milk production traits. One region affecting milk fat and protein percentage was located on the equivalent of Bos taurus autosome (BTA)3, spanning 43.3 to 43.8 Mb, which harbored the most likely candidate genes MFSD14A, SLC35A3, and PALMD. The other region on the equivalent of BTA14 at 66.5 to 67.0 Mb contained candidate genes RGS22 and VPS13B and influenced buffalo total milk yield, fat yield, and protein yield. Interestingly, both of the regions were reported to have quantitative trait loci affecting milk performance in dairy cattle. Furthermore, we suggest that buffaloes with the C allele at AX-85148558 and AX-85073877 loci and the G allele at AX-85106096 locus can be selected to improve milk fat yield in this buffalo-breeding program. Meanwhile, the G allele at AX-85063131 locus can be used as the favorable allele for improving milk protein percentage. Genomic prediction showed that the reliability of genomic estimated breeding values (GEBV) of 6 milk production traits ranged from 0.06 to 0.22, and the correlation between estimated breeding values and GEBV ranged from 0.23 to 0.35. These findings provide useful information to understand the genetic basis of buffalo milk properties and may play a role in accelerating buffalo breeding programs using genomic approaches.

Transcript analysis at DGAT1 reveals different mRNA profiles in river buffaloes with extreme phenotypes for milk fat

Buffalo DGAT1 (diacylglycerol O-acyltransferase 1) was mainly investigated for the characterization of the gene itself and for the identification of the K232A polymorphism, similar to what has been accomplished in cattle, although no information has been reported so far at the mRNA level. The importance of DGAT1 for lipid metabolism led us to investigate the transcript profiles of lactating buffaloes characterized as high (9.13 ± 0.23) and low (7.94 ± 0.29) for milk fat percentage, and to explore the genetic diversity at the RNA and DNA level. A total of 336 positive clones for the DGAT1 cDNA were analyzed by PCR and chosen for sequencing according to the differences in length. The clone assembling revealed a very complex mRNA pattern with a total of 21 transcripts differently represented in the 2 groups of animals. Apart from the correct transcript (17 exons long), the skipping of exon 12 is the most significant in terms of distribution of clones with 11.6% difference between the 2 groups, whereas a totally different mRNA profile was found in approximately 12% of clones. The sequencing of genomic DNA allowed the identification of 10 polymorphic sites at the intron level, which clarify, at least partially, the genetic events behind the production of complex mRNA. Genetic diversity was found also at the exon level. The single nucleotide polymorphism c.1053C>T represents the first example of polymorphism in a coding region for the DGAT1 in the Italian Mediterranean breed. To establish whether this polymorphism is present in other buffalo breeds, a quick method based on PCR-RFLP was set up for allelic discrimination in the Italian Mediterranean and the Romanian Murrah (200 animals in total). The alleles were equally represented in the overall population, whereas the analysis of the 2 breeds showed different frequencies, likely indicating diverse genetic structure of the 2 breeds. The T allele might be considered as the ancestral condition of the DGAT1 gene, being present in the great part of the sequenced species. These data add knowledge at the transcript and genetic levels for the buffalo DGAT1 and open the opportunity for further investigation of other genes involved in milk fat metabolism for the river buffalo, including the future possibility of selecting alleles with quantitative or qualitative favorable effects (or both).

Identification and analysis of the expression of microRNA from lactating and nonlactating mammary glands of the Chinese swamp buffalo

To study the role of microRNA (miR) in the lactation physiology of water buffalo, 2 multiparous dairy buffaloes (including an 8-yr-old buffalo that had been lactating for 3 mo, as well as a 10-yr-old nonlactating, nonpregnant buffalo) were used for miR library construction. The profile of differentially expressed miR in lactating and nonlactating mammary gland tissues of these water buffalo were investigated using Illumina-Solexa high-throughput sequencing technology (Illumina, San Diego, CA). The data identified 259 miR families, 359 mature miR, 363 pre-miR, 230 novel buffalo miR, and 5 buffalo-specific miR that were expressed in mammary tissues. Some highly significantly differentially expressed miR were explored, including bbu-miR-497, bbu-miR-30a-5p, bbu-miR-148a, bbu-miR-29a, bbu-miR-125a, bbu-miR-125b, and bbu-miR-103. The expression patterns of 18 miR were confirmed by quantitative real-time PCR in both tissues, and the expression of bbu miR-103 and novel miR-57 constituted the largest differences between lactating and nonlactating tissues. Further functional analysis indicated that the overexpression or suppression of miR-103 in buffalo mammary epithelial cells downregulated or upregulated the expression of pantothenate kinase 3, and also significantly increased the transcription factor steroid regulatory element binding protein, followed by the acceleration of de novo synthesis of fatty acids by upregulation of acetyl coenzyme A carboxylase α expression. The expression levels of 34 predicted target genes of novel-miR-57 in lactating and nonlactating mammary gland tissues were all analyzed by quantitative real-time PCR. Finally, only the expression of docking protein 4 could be upregulated or downregulated selectively by bbu-novel-miR-57 in buffalo mammary epithelial cells and the Bcap-37 cell line. This study provides an overview of the miR expression profile of water buffalo and the interaction between some key miR and their target genes, which may improve understanding of the important roles of miR in buffalo milk fat synthesis.

Biological underpinnings of breastfeeding challenges: the role of genetics, diet, and environment on lactation physiology

Lactation is a dynamic process that has evolved to produce a complex biological fluid that provides nutritive and nonnutritive factors to the nursing offspring. It has long been assumed that once lactation is successfully initiated, the primary factor regulating milk production is infant demand. Thus, most interventions have focused on improving breastfeeding education and early lactation support. However, in addition to infant demand, increasing evidence from studies conducted in experimental animal models, production animals, and breastfeeding women suggests that a diverse array of maternal factors may also affect milk production and composition. In this review, we provide an overview of our current understanding of the role of maternal genetics and modifiable factors, such as diet and environmental exposures, on reproductive endocrinology, lactation physiology, and the ability to successfully produce milk. To identify factors that may affect lactation in women, we highlight some information gleaned from studies in experimental animal models and production animals. Finally, we highlight the gaps in current knowledge and provide commentary on future research opportunities aimed at improving lactation outcomes in breastfeeding women to improve the health of mothers and their infants.

Comparative analysis of SNP candidates in disparate milk yielding river buffaloes using targeted sequencing

River buffalo (Bubalus bubalis) milk plays an important role in economy and nutritious diet in several developing countries. However, reliable milk-yield genomic markers and their functional insights remain unexposed. Here, we have used a target capture sequencing approach in three economically important buffalo breeds namely: Banni, Jafrabadi and Mehsani, belonging to either high or low milk-yield group. Blood samples were collected from the milk-yield/breed balanced group of 12 buffaloes, and whole exome sequencing was performed using Roche 454 GS-FLX Titanium sequencer. Using an innovative approach namely, MultiCom; we have identified high-quality SNPs specific for high and low-milk yield buffaloes. Almost 70% of the reported genes in QTL regions of milk-yield and milk-fat in cattle were present among the buffalo milk-yield gene candidates. Functional analysis highlighted transcriptional regulation category in the low milk-yield group, and several new pathways in the two groups. Further, the discovered SNP candidates may account for more than half of mammary transcriptome changes in high versus low-milk yielding cattle. Thus, starting from the design of a reliable strategy, we identified reliable genomic markers specific for high and low-milk yield buffalo breeds and addressed possible downstream effects.

Isolation and Identification of Prepubertal Buffalo (Bubalus bubalis) Spermatogonial Stem Cells

Isolation and culture of spermatogonial stem cells (SSCs) are attractive for production of genetic modified offspring. In the present study, buffalo spermatogonial stem-like cells were isolated, cultured and expression pattern of different germ cell marker genes were determined. To recover spermatogonia, testes from age 3 to 7 months of buffalo were decapsulated, and seminiferous tubules were enzymatically dissociated. Two types of cells, immature sertoli cell and type A spermatogonia were observed in buffalo testes in this stage. Germ cell marker genes, OCT3/4 (Pou5f1), THY-1, c-kit, PGP9.5 (UCHL-1) and Dolichos biflorus agglutinin, were determined to be expressed both in mRNA and protein level by reverse transcription polymerase chain reaction and immunostaining in buffalo testes and buffalo spermatogonial stem-like cells, respectively. In the following, when the isolated buffalo buffalo spermatogonial stem-like cells were cultured in the medium supplemented 2.5% fetal bovine serum and 40 ng/mL glial cell-derived neurotrophic factor medium, SSCs proliferation efficiency and colony number were significantly improved than those of other groups (p<0.05). These findings may help in isolation and establishing long term in vitro culture system for buffalo spermatogonial stem-like cells, and accelerating the generation of genetic modified buffaloes.