Transcriptome profiling of bovine milk oligosaccharide metabolism genes using RNA-sequencing

This study examines the genes coding for enzymes involved in bovine milk oligosaccharide metabolism by comparing the oligosaccharide profiles with the expressions of glycosylation-related genes. Fresh milk samples (n = 32) were collected from four Holstein and Jersey cows at days 1, 15, 90 and 250 of lactation and free milk oligosaccharide profiles were analyzed. RNA was extracted from milk somatic cells at days 15 and 250 of lactation (n = 12) and gene expression analysis was conducted by RNA-Sequencing. A list was created of 121 glycosylation-related genes involved in oligosaccharide metabolism pathways in bovine by analyzing the oligosaccharide profiles and performing an extensive literature search. No significant differences were observed in either oligosaccharide profiles or expressions of glycosylation-related genes between Holstein and Jersey cows. The highest concentrations of free oligosaccharides were observed in the colostrum samples and a sharp decrease was observed in the concentration of free oligosaccharides on day 15, followed by progressive decrease on days 90 and 250. Ninety-two glycosylation-related genes were expressed in milk somatic cells. Most of these genes exhibited higher expression in day 250 samples indicating increases in net glycosylation-related metabolism in spite of decreases in free milk oligosaccharides in late lactation milk. Even though fucosylated free oligosaccharides were not identified, gene expression indicated the likely presence of fucosylated oligosaccharides in bovine milk. Fucosidase genes were expressed in milk and a possible explanation for not detecting fucosylated free oligosaccharides is the degradation of large fucosylated free oligosaccharides by the fucosidases. Detailed characterization of enzymes encoded by the 92 glycosylation-related genes identified in this study will provide the basic knowledge for metabolic network analysis of oligosaccharides in mammalian milk. These candidate genes will guide the design of a targeted breeding strategy to optimize the content of beneficial oligosaccharides in bovine milk.     

Regulation of peroxisome proliferator-activated receptor gamma on milk fat synthesis in dairy cow mammary epithelial cells

Peroxisome proliferator-activated receptor gamma (PPARγ) participates in lipogenesis in rats, goats, and humans. However, the exact mechanism of PPARγ regulation on milk fat synthesis in dairy cow mammary epithelial cells (DCMECs) remains largely unexplored. The aim of this study was to investigate the role of PPARγ regarding milk fat synthesis in DCMECs and to ascertain whether milk fat precursor acetic acid and palmitic acid could interact with PPARγ signaling to regulate milk fat synthesis. For this study, we examined the effects of PPARγ overexpression and gene silencing on cell growth, triacylglycerol synthesis, and the messenger RNA (mRNA) and protein expression levels of genes involved in milk fat synthesis in DCMECs. In addition, we investigated the influences of acetic acid and palmitic acid on the mRNA and protein levels of milk lipogenic genes and triacylglycerol synthesis in DCMECs transfected with PPARγ small interfering RNA (siRNA) and PPARγ expression vector. The results showed that when PPARγ was silenced, cell viability, proliferation, and triacylglycerol secretion were obviously reduced. Gene silencing of PPARγ significantly downregulated the expression levels of milk fat synthesis-related genes in DCMECs. PPARγ overexpression improved cell viability, proliferation, and triacylglycerol secretion. The expression levels of milk lipogenic genes were significantly increased when PPARγ was overexpressed. Acetic acid and palmitic acid could markedly improve triacylglycerol synthesis and upregulate the expression levels of PPARγ and other lipogenic genes in DCMECs. These results suggest that PPARγ is a positive regulator of milk fat synthesis in DCMECs and that acetic acid and palmitic acid could partly regulate milk fat synthesis in DCMECs via PPARγ signaling.     

Development of a RNA extraction method from milk for gene expression study in the mammary gland of sheep

The aim of the study was to develop a reliable method for the RNA extraction from milk of Sarda sheep breed and to highlight if the extracted RNA can be used for expression study on mammary genes involved in milk fat synthesis using RT-qPCR. The main result is that a sample of 150 ml of milk provides an optimal amount of RNA (73.5 μg/ml). The highest RNA concentration has been found in the samples analysed within 4 h after collection. The RNA extracted was positively correlated to the number of somatic cells (P < 0.001). The efficiency of the extraction method was confirmed by the results obtained from qPCR which showed a Ct value, for SREBPF1 gene of 26.8 ± 0.15. This research demonstrated that the high-quality of the RNA obtained is suited to use for studies of mammary genes expression in sheep, avoiding any damage caused by mammary gland biopsy.     

Suitable internal control microRNA genes for measuring miRNA abundance in pig milk during different lactation periods

Determination of an optimal set/number of internal control microRNA (miRNA) genes is a critical, but often undervalued, detail of quantitative gene expression analysis. No validated internal genes for miRNA quantitative PCR (q-PCR) in pig milk were available. We compared the expression stability of six porcine internal control miRNA genes in pig milk from different lactation periods (1 h, 3 days, 7 days, 14 days, 21 days, and 28 days postpartum), using an EvaGreen q-PCR approach. We found that using the three most stable internal control genes to calculate the normalization factor is sufficient for producing reliable q-PCR expression data. We also found that miRNAs are superior to ribosomal RNA (rRNA) and snRNA, which are commonly used as internal controls for normalizing miRNA q-PCR data. In terms of economic and experimental feasibility, we recommend the use of the three most stable internal control miRNA genes (miR-17, -107 and -103) for calculating the normalization factors for pig milk samples from different lactation periods. These results can be applied to future studies aimed at measuring miRNA abundance in porcine milk.     

Gelsolin expression in sheep milk somatic cells during lactation

The identification of genes involved in phenotypes related to milk quality is important for both economic and health aspects in livestock production. The aim of this study was to assess the level of gelsolin gene expression in two breeds of dairy sheep – Sarda and Gentile – with pronounced differences in quantitative and qualitative milk traits. Gelsolin, a type of actin-modulating proteins is involved in the processes of actin remodeling during cell growth and apoptosis; therefore a role of this protein in mammary changes during lactation was here hypothesized. Individual milk samples were collected three times during lactation from 26 ewes of the two breeds. The differential gene expression of gelsolin in the two breeds and the three lactation times was estimated by quantitative PCR on RNA extracted from milk somatic cells. Correlations of gelsolin gene expression with milk yield and quality and days of lactation were also estimated. The results showed that gelsolin gene expression was significantly higher in the Sarda compared to the Gentile at each lactation stage, in agreement with the longer lactation duration and the higher daily milk yield of the first breed. Significant correlations of gelsolin gene expression were found with milk fat content in Sarda breed (−0.46, P<0.05). Gelsolin expression analysis confirmed the link between gelsolin gene function and milk fat content of sheep.     

Differential expression of immune response genes associated with subclinical mastitis in dairy buffaloes

Buffalo milk production has become of significant importance on the world scale, however, there are few studies involving biotechnological tools specifically for buffalo. To verify the effects caused by subclinical mastitis on the components of milk and to study the innate immune system in the udder of dairy buffaloes with subclinical mastitis, we evaluated the levels of expression of the lactoferrin (LTF), tumor necrosis factor alpha (TNF-α), interleukin-1 beta (IL-1β), interleukin-8 (IL-8), and toll-like receptors 2 (TLR-2) and 4 (TLR-4) genes in buffaloes with and without subclinical mastitis. Milk samples were collected for the determination of milk components: somatic cell score (SCS), fat, protein, lactose, total solids and solids-not-fat (SNF), as well as for RNA extraction of milk cells, complementary DNA synthesis, and expression profile quantification by quantitative real-time PCR. For gene expression, the ΔΔCt was estimated using contrasts of the target genes expression adjusted for the expression of the housekeeping genes between both groups. Linear regression analysis was performed to determine the relationship between the genes studied and the milk components. Subclinical mastitis induced changes in the fat, lactose and SNF in milk of buffaloes, and the messenger RNA abundance was upregulated for TLR-2, TLR-4, TNF-α, IL-1β and IL-8 genes in milk cells of buffaloes with subclinical mastitis, whereas the LTF gene was not differentially expressed. Results of linear regression analysis showed that TLR-2 gene expression most explains the variation in SCS, and the change in a unit of ΔCt of the TNF-α gene would result in a higher increase in SCS. The study of these immune function genes that are active in the mammary gland is important to characterize the action mechanism of the innate immunity that occurs in subclinical mastitis in dairy buffaloes and may aid the development of strategies to preserve the health of the udder.     

奶牛产奶性状与乳房炎相关基因CpG含量及分布特征的比较分析

DNA甲基化是表观遗传学的重要组成部分,基因启动子区及第一外显子区的CpG甲基化通常抑制该基因的表达,而去甲基化则促进基因表达。已有的研究发现荷斯坦牛的乳房炎指标SCC(Somatic cell count)与产奶量呈较强负相关。文章分析并比较了这两类性状的相关基因的启动子区、第一外显子、下游2 000 bp序列中CpG含量及分布特征。结果表明,乳房炎相关基因的启动子、第一外显子中CpG含量显著低于产奶性状相关基因,而两类性状基因下游2 000 bp序列中CpG含量无显著性差异。另外,文中提出了两个量化基因序列中CpG特征的指标,一个是CpG平均距离,用来衡量序列中的CpG分布;另一个是条件概率p(G|C),用以量化序列中二核苷酸CpG随碱基C出现的可能性,并对两类基因的启动子和第一外显子区域的这两个指标做了统计检验。研究结果对产奶性状与乳房炎相关基因的DNA甲基化调控研究奠定了基础。