The effect of single-nucleotide polymorphism in the promoter region of bovine alpha-lactalbumin (LALBA) gene on LALBA expression in milk cells and milk traits of cows

Genetic variability among and within domestic Old and New World camels at the α-lactalbumin gene (LALBA) reveals new alleles and polymorphisms responsible for differential expression

α-Lactalbumin (α-LA), which is encoded by the LALBA gene, is a major whey protein that binds to Ca2+ and facilitates lactose synthesis as a regulatory subunit of the synthase enzyme complex. In addition, it has been shown to play central roles in immune modulation, cell-growth regulation, and antimicrobial activity. In this study, a multitechnical approach was used to fully characterize the LALBA gene and its variants in both coding and regulatory regions for domestic camelids (dromedary, Bactrian camel, alpaca, and llama). The gene analysis revealed a conserved structure among the camelids, but a slight difference in size (2,012 bp on average) due to intronic variations. Promoters were characterized for the transcription factor binding sites (11 found in total). Intraspecies sequence comparison showed 36 SNPs in total (2 in the dromedary, none in the Bactrian camel, 22 in the alpaca, and 12 in the llama), whereas interspecies comparison showed 86 additional polymorphic sites. Eight SNPs were identified as trans-specific polymorphisms, and 2 of them (g.112A>G and g.1229A>G) were particularly interesting in the New World camels. The first creates a new binding site for transcription factor SP1. An enhancing effect of the g.112G variant on the expression was demonstrated by 3 independent pGL3 gene reporter assays. The latter is responsible for the p.78Ile>Val AA replacement and represents novel allelic variants (named LALBA A and B). A link to protein variants has been established by isoelectric focusing (IEF), and bioinformatics analysis revealed that carriers of valine (g.1229G) have a higher glycosylation rate. Genotyping methods based on restriction fragment length polymorphism (PCR-RFLP) were set up for both SNPs. Overall, adenine was more frequent (0.54 and 0.76) at both loci. Four haplotypes were found, and the AA and GA were the most common with a frequency of 0.403 and 0.365, respectively. Conversely, a putative biological gain characterizes the haplotype GG. Therefore, opportunities for rapid directional selection can be realized if this haplotype is associated with favorable milk protein properties. This study adds knowledge at the gene and protein level for α-LA (LALBA) in camelids and importantly contributes to a relatively unexplored research area in these species.

Milk proteomic analysis reveals differentially expressed proteins in high-yielding and low-yielding Guanzhong dairy goats at peak lactation

The aim of this experiment was to investigate the differential proteomic characteristics of milk from high- and low-yielding Guanzhong dairy goats during the peak lactation period under the same feeding conditions. Nine Guanzhong dairy goats with high yield (H: 3.5 ± 0.17 kg/d) and nine with low yield (L:1.2 ± 0.25 kg/d) were selected for milk proteomic analysis using tandem mass tag technology. A total of 78 differentially expressed proteins were identified. Compared with L, 50 proteins including HK3, HSPB1 and ANXA2 were significantly upregulated in H milk, while 28 proteins including LALBA and XDH were significantly downregulated. Bioinformatics analysis of the differentially expressed proteins showed that galactose metabolism, purine metabolism, glycolysis/gluconeogenesis, MAPK signaling pathway, regulation of actin cytoskeleton and other pathways were closely related to milk yield. HK3, HSPB1, ANXA2, LALBA and XDH were important candidate proteins associated with the milk production characteristics of Guanzhong dairy goats. Our data provide relevant biomarkers and a theoretical basis for improving milk production in Guanzhong dairy goats.

Supplementation of dietary semen vaccariae extracts to lactating sow diets: effects on the production performance, milk components, and gene expression related to mammogenesis

This study aimed to investigate the effects of dietary semen vaccariae extracts (SVE) on the production performance, colostrum components, and relative gene expression related to mammogenesis of lactating sows. 48 pregnant sows were selected and randomly allocated into four groups, with six replicates and two sows per replicate. The first group was the control (CON), while the other groups received the same diet further supplemented with 1.5, 3.0 and 4.5 g SVE per kg (SV1, SV2 and SV3, respectively). Compared with the control group, (1) the average daily gain was increased (p < 0.05) in SV1, SV2, and SV3 during the 11-21 days and 1-21 days of lactation; (2) the serum insulin-like growth factor-1, insulin, prolactin, and estrogen contents in SV1, SV2, and SV3 were increased (p < 0.05) on the 1st and 21st day of lactation; (3) The plasma Lysine, Threonine, and Tryptophan concentrations were also higher (p < 0.05) in SV1, SV2, and SV3 on the 1st and 21st day of lactation; (4) The milk Lysine, Methionine, Threonine, and Tryptophan concentrations were higher (p < 0.05) in SV1, SV2, and SV3 on the 1st and 21st day of lactation; (5) The milk lactose ratio and milk protein content were increased (p < 0.05) in the groups treated with semen vaccariae on the 1st day of lactation, while the milkfat ratio and milk protein content were increased (p < 0.05) in SV2 and SV3 on the 21st day of lactation; (6) the immunoglobulin M, A, and G contents were increased (p < 0.05) in the groups treated with the semen vaccariae on the first day of lactation; and (7) the relative PRLR, STAT5a, FcRn, CSN2, and LALBA expressions were higher (p < 0.05) in the groups treated with the semen vaccariae on the 1st and 21st day of lactation. In this study, the optimum dosage was 3.0 g/kg semen vaccariae, which increased the average daily gain of piglets, total lactation yield, and serum hormone levels, improved the amino acid levels in plasma, and facilitated the milk quality, up-regulated the relative gene expressions in the mammogenesis.

Read-depth based approach on whole genome resequencing data reveals important insights into the copy number variation (CNV) map of major global buffalo breeds

BACKGROUND

Elucidating genome-wide structural variants including copy number variations (CNVs) have gained increased significance in recent times owing to their contribution to genetic diversity and association with important pathophysiological states. The present study aimed to elucidate the high-resolution CNV map of six different global buffalo breeds using whole genome resequencing data at two coverages (10X and 30X). Post-quality control, the sequence reads were aligned to the latest draft release of the Bubaline genome. The genome-wide CNVs were elucidated using a read-depth approach in CNVnator with different bin sizes. Adjacent CNVs were concatenated into copy number variation regions (CNVRs) in different breeds and their genomic coverage was elucidated.

RESULTS

Overall, the average size of CNVR was lower at 30X coverage, providing finer details. Most of the CNVRs were either deletion or duplication type while the occurrence of mixed events was lesser in number on a comparative basis in all breeds. The average CNVR size was lower at 30X coverage (0.201 Mb) as compared to 10X (0.013 Mb) with the finest variants in Banni buffaloes. The maximum number of CNVs was observed in Murrah (2627) and Pandharpuri (25,688) at 10X and 30X coverages, respectively. Whereas the minimum number of CNVs were scored in Surti at both coverages (2092 and 17,373). On the other hand, the highest and lowest number of CNVRs were scored in Jaffarabadi (833 and 10,179 events) and Surti (783 and 7553 events) at both coverages. Deletion events overnumbered duplications in all breeds at both coverages. Gene profiling of common overlapped genes and longest CNVRs provided important insights into the evolutionary history of these breeds and indicate the genomic regions under selection in respective breeds.

CONCLUSION

The present study is the first of its kind to elucidate the high-resolution CNV map in major buffalo populations using a read-depth approach on whole genome resequencing data. The results revealed important insights into the divergence of major global buffalo breeds along the evolutionary timescale.

Comparative analyses of copy number variations between swamp buffaloes and river buffaloes

As an important source of genomic variation, copy number variation (CNV) contributes to environmental adaptation in worldwide buffaloes. Despite this importance, CNV divergence between swamp buffaloes and river buffaloes has not been studied previously. Here, we report 21 152 CNV regions (CNVRs) in 141 buffaloes of 20 breeds detected through multiple CNV calling strategies. Only 248 CNVRs were shared between river buffalo and swamp buffalo, reflecting great variation of CNVRs between the two subspecies. Population structure analysis based on CNVs successfully separated the two buffalo subspecies. We further assessed CNV divergence by calculating F_ST for genome-wide CNVs. Totally, we identified 110 significantly divergent CNV segments and 44 putatively selected genes between river buffaloes and swamp buffaloes. In particular, LALBA, a key gene controlling milk production in cattle, presented a highly differentiated CNV in the promoter region, which makes it a strong functional candidate gene for differences between swamp buffaloes and river buffaloes in traits related to milk production. Our study provides useful information of CNVs in buffaloes, which may help explain the genetic differences between the two subspecies.