The role of bovine causal genes underlying dairy traits in Spanish Churra sheep

In dairy cattle, quantitative trait nucleotides (QTNs) underlying quantitative trait loci (QTL) for milk production traits have been identified in bovine DGAT1, GHR and ABCG2 genes. The SPP1 gene has also been proposed to be a regulator of lactation. In sheep, QTL underlying milk production traits have been reported only recently, and no proven QTN has been identified. Taking into account the close phylogenetic relationship between sheep and cattle, this study examined the possible effects of the aforementioned genes on sheep milk production traits. We first studied the genetic variability of the DGAT1, GHR, ABCG2 and SPP1 genes in 15 rams of the Spanish Churra dairy sheep breed. Second, we performed an association analysis between SNPs identified in these genes and three milk production traits recorded in a commercial population of Churra sheep. This analysis revealed only three significant associations at the nominal level (P-value <0.05) involving allelic variants of the ABCG2 gene, whereas no significant association was found for the DGAT1, GHR and SPP1 genes. When the Bonferroni correction was applied to take into account the multiple tests performed, none of the associations identified at the nominal level remained significant. Nevertheless, taking into account the high level of false-negative findings that can arise when applying the stringent Bonferroni correction, we think that our results provide a valuable primary assessment of strong candidate genes for milk traits in sheep.     

Database of cattle candidate genes and genetic markers for milk production and mastitis

A cattle database of candidate genes and genetic markers for milk production and mastitis has been developed to provide an integrated research tool incorporating different types of information supporting a genomic approach to study lactation, udder development and health. The database contains 943 genes and genetic markers involved in mammary gland development and function, representing candidates for further functional studies. The candidate loci were drawn on a genetic map to reveal positional overlaps. For identification of candidate loci, data from seven different research approaches were exploited: (i) gene knockouts or transgenes in mice that result in specific phenotypes associated with mammary gland (143 loci); (ii) cattle QTL for milk production (344) and mastitis related traits (71); (iii) loci with sequence variations that show specific allele-phenotype interactions associated with milk production (24) or mastitis (10) in cattle; (iv) genes with expression profiles associated with milk production (207) or mastitis (107) in cattle or mouse; (v) cattle milk protein genes that exist in different genetic variants (9); (vi) miRNAs expressed in bovine mammary gland (32) and (vii) epigenetically regulated cattle genes associated with mammary gland function (1). Fourty-four genes found by multiple independent analyses were suggested as the most promising candidates and were further in silico analysed for expression levels in lactating mammary gland, genetic variability and top biological functions in functional networks. A miRNA target search for mammary gland expressed miRNAs identified 359 putative binding sites in 3'UTRs of candidate genes.     

The concordance test emerges as a powerful tool for identifying quantitative trait nucleotides: lessons from BTA6 milk yield QTL

The lack of conventions for confirming the discovery of quantitative trait nucleotides in livestock was evidenced by the proposals of mutations in two different genes ( SPP1 and ABCG2 ) as the underlying functional mutation for a major quantitative trait locus (QTL) for milk concentration on bovine chromosome 6 (BTA6). Of these conflicting candidates, SPP1 was excluded by follow-up studies and by the data described here. A simple test for concordance of the zygosity state between QTL segregation status and the candidate polymorphism was shown, in this case, to be a critical step towards establishing the proof. If a given sample effectively represents the genetic variation across the QTL region, haplotype-based concordance may further enhance the functionality and resolution power of this test, allowing identification of the causative gene.     

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.     

A novel approach identified the FOLR1 gene, a putative regulator of milk protein synthesis

This study has utilised comparative functional genomics to exploit animal models with extreme adaptation to lactation to identify candidate genes that specifically regulate protein synthesis in the cow mammary gland. Increasing milk protein production is valuable to the dairy industry. The lactation strategies of both the Cape fur seal (Artocephalus pusillus pusillus) and the tammar wallaby (Macropus eugenii) include periods of high rates of milk protein synthesis during an established lactation and therefore offer unique models to target genes that specifically regulate milk protein synthesis. Global changes in mammary gene expression in the Cape fur seal, tammar wallaby, and the cow (Bos taurus) were assessed using microarray analysis. The folate receptor α (FOLR1) showed the greatest change in gene expression in all three species [cow 12.7-fold (n = 3), fur seal 15.4-fold (n = 1), tammar 2.4-fold (n = 4)] at periods of increased milk protein production. This compliments previous reports that folate is important for milk protein synthesis and suggests FOLR1 may be a key regulatory point of folate metabolism for milk protein synthesis within mammary epithelial cells (lactocytes). These data may have important implications for the dairy industry to develop strategies to increase milk protein production in cows. This study illustrates the potential of comparative genomics to target genes of interest to the scientific community.     

Associations of the bovine major histocompatibility complex DRB3 (BoLA-DRB3) with production traits in Canadian dairy cattle

Associations of two alleles of the bovine major histocompatibility complex DRB3 gene (BoLA-DRB3) with lowered somatic cell score (SCS) and occurrence of disease (BoLA-DRB3.2* 16 and *23, respectively) have previously been documented. The objective of this study was to evaluate potential relationships between BoLA-DRB3 alleles with production traits, namely 305-day milk, milk fat and milk protein yield, in a population of Canadian dairy cattle (Holstein, n = 835 and Jersey, n = 66) over the course of two lactations. No significant associations were detected between BoLA alleles and production traits in Jerseys. In Holsteins, alleles *16 and *23 also did not show associations with production traits but allele *8 was significantly associated with increased 305-day milk, fat and protein yields in the previous lactation (the lactation prior to immunization with a gram negative core antigen vaccine), and with increased protein production in the subsequent (with reference to the time of immunization) lactation. Allele *22 was associated with decreased milk and protein yield in both previous and subsequent lactations. Therefore, it can be concluded that increasing or decreasing the frequency of BoLA alleles *16 and *23 to reduce SCS or increase resistance to mastitis in this population would not have adverse effects on production in this population, and that certain BoLA alleles (*8 and *22) are associated with altered production traits in Canadian Holsteins.     

Comparative genome analysis with the human genome reveals chicken genes associated with fatness and body weight

The selection of meat-type chickens (broilers) for rapid growth has been accompanied by excessive fat deposition. In this study, we analysed 53 candidate genes that are associated with obesity and obesity-related traits in humans, for which we found chicken orthologues by BLAST searches. We have identified single nucleotide polymorphisms (SNPs) with significant differences in allele frequencies between broilers and layers in each of the following six candidate genes: adrenergic, beta-2-, receptor, surface (ADRB2); melanocortin 5 receptor (MC5R); leptin receptor (LEPR), McKusick-Kaufman syndrome (MKKS), milk fat globule-EGF factor 8 protein (MFGE8) and adenylate kinase 1 (AK1). To examine associations with fatness and/or body weight, we used birds of extreme phenotypes in F(2) and backcross populations with varying levels of abdominal fat weight per cent (%AFW) and body weight. We then assessed the level of gene expression by real-time PCR. In two genes, ADRB2 and MFGE8, we found significant association with %AFW. The ADRB2 gene was found to have a significantly higher expression in the liver of lean chickens compared with those of the fat individuals. We believe that this approach can be applied for the identification of other quantitative genes.     

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.     

Molecular differentiation of Peroxysome proliferator activated receptor coactivator-1 among different breeds of Bubalus bubalis

Peroxysome proliferator activated receptor coactivator-1 gene (PPARGC1A) is a positional and functional candidate gene for milk fat yield. It has key role in energy, fat and glucose metabolism. Single nucleotide polymorphisms (SNPs) in Exon-8 of PPARGC1A are reported to be associated with milk fat yield in dairy cattle. In the present investigation PPARGC1A was partially amplified (around 767bp) by designing gene specific primer and confirm by sequencing the amplicon and its comparison with the PPARGC1A gene of bovine. Comparative study of PPARGC1A among different breeds of buffaloes reveals different level of mutations with respect to its gene sequence 0.013-1.69% and protein sequence 0.42% to 2.99%, Similarly the protein structures modeled from their sequences were compared by structural superposition that shows variations (RMSD) from 0.736 to 1.507. Furthermore, the sequences were used to generate a dendrogram. It reveals that Murrah and reference are very close to each other, similarly Toda, Bhadawari and Surti are closely related, whereas Pandharpuri is separated from both the cluster. Especially the variations are more at the binding site of this protein that may be the cause that different breeds have different percentage of milk fat. Further study is underway to detect polymorphism and associate them with milk fat related traits in buffalo.     

Genetic variations in the SPP1 promoter affect gene expression and the level of osteopontin secretion into bovine milk

Osteopontin (OPN) is now recognized as an important cytokine and extracellular integrin‐binding protein at the crossroads of inflammation and homeostasis. In a previous study, we found that OPN gene (SPP1) polymorphisms are associated with milk performance traits and somatic cell score (SCS), a parameter used to estimate the genetic value of udder health in dairy cattle. In this study, we assessed whether the genetic variations had an impact on SPP1 promoter activity, immune response and the level of OPN secreted into milk. The influence of DNA polymorphisms on the promoter activity of SPP1 was confirmed in vitro. To measure the impact of the genetic variations on OPN secretion into milk, we measured OPN levels in both plasma and milk throughout lactation. Cows were grouped by the OPN haplotypes associated with a high (H2 × H3) or low (H1 × H4) SCS. For both H2 × H3 and H1 × H4, the OPN level in plasma remained low throughout lactation, although the concentration in the milk of H1 × H4 cows increased more in late lactation. Moreover, the macrophages of H1 × H4 cows expressed a lower SPP1 and proinflammatory IL6 in response to infection. Regarding the immune cell response, cows with the genetic potential to secrete higher OPN levels during late lactation had macrophages expressing fewer proinflammatory cytokines, a situation that might explain the genetic association with low somatic cells. Although OPN's favorable roles during late lactation remain to be elucidated, the tissue remodeling properties associated with OPN may be beneficial for reducing the incidence of infection during the transition period in lactating cows.     

DGAT1, a new positional and functional candidate gene for intramuscular fat deposition in cattle

Intramuscular fat content, also assessed as marbling of meat, represents an important beef quality trait. Recent work has mapped a quantitative trait locus (QTL) with an effect on marbling to the centromeric region of bovine chromosome 14, with the gene encoding thyroglobulin (TG) being proposed as a positional and functional candidate gene for this QTL. Recently, the gene encoding diacylglycerol O-acyltransferase (DGAT1), which also has been mapped within the region of the marbling QTL, has been demonstrated to affect the fat content of milk. In the present study, the effects of a 5'-polymorphism of TG and of a lysine/alanine polymorphism of DGAT1 on the fat content of musculus (m.) semitendinosus and m. longissimus dorsi in 55 bovine animals (28 German Holstein and 27 Charolais) has been investigated. Significant effects were found for both candidate genes in both the breeds. These effects seem to be independent of one another because the alleles of the two polymorphisms showed no statistically significant disequilibrium. The DGAT1 effect is mainly on the m. semitendinosus. The TG polymorphism only affects m. longissimus dorsi. However, both intramuscular fat enhancing effects seem to be recessive. The possibility of two linked loci, acting recessively on intramuscular fat content, will require special strategies when selecting for higher marbling scores.     

The epithelial glucocorticoid receptor is required for the normal timing of cell proliferation during mammary lobuloalveolar development but is dispensable for milk production

Glucocorticoids have been shown to influence mammary gland function in vivo and to stimulate milk protein gene expression in vitro. Here, we describe the generation and analysis of a mouse model to study glucocorticoid receptor (GR, NR3C1) function in mammary epithelial cells. Using the Cre-loxP system, mutant mice were obtained in which the GR gene is specifically deleted in epithelial cells during lobuloalveolar development, leading to a complete loss of epithelial GR at the onset of lactation. Mice harboring the mammary-epithelial-specific GR mutation are able to nurse their litters until weaning. During pregnancy, however, GR deficiency delays lobuloalveolar development, leading to an incomplete epithelial penetration of the mammary fat pad that persists throughout lactation. We identified a reduced cell proliferation during lobuloalveolar development as reason for this delay. This reduction is compensated for by increased epithelial proliferation after parturition in the mutant glands. During lactation, GR-deficient mammary epithelium is capable of milk production and secretion. The expression of two milk proteins, namely whey acidic protein and beta-casein, during lactation was not critically affected in the absence of GR. We conclude that GR function is not essential for alveolar differentiation and milk production, but influences cell proliferation during lobuloalveolar development.     

Identification of SNPs in Interferon Gamma,Interleukin-22, and Their Receptors and Associations with Health and Production-Related Traits in Canadian Holstein Bulls

Genetic variants in a number of immunoregulatory genes have been previously associated with health and production traits in dairy cattle. Therefore, in the following study, the genes coding interferon gamma (IFNG), IFNG receptor 1 and 2 domains, interleukin-22 (IL22), and IL22 receptor alpha 1, were investigated for single nucleotide polymorphisms (SNPs) in Holstein bulls. These SNPs, along with SNPs previously identified in IL10, IL10 receptor, and transforming growth factor beta 1 (TGFB1) genes, were evaluated for statistical associations to estimated breeding values for milk somatic cell score (SCS), a trait highly correlated to mastitis incidence, and various production-related traits, including milk yield, protein yield, fat yield, and lactation persistency. While no significant associations were found between these SNPs and SCS, SNPs in IL10 receptor beta subunit showed a significant effect on protein yield and lactation persistency. While there is evidence that IL10 plays an important role during lactation, it is also likely that the effects of SNPs in IL10 receptor beta subunit on protein yield and lactation persistency are due to linkage disequilibrium with a neighboring QTL.     

转基因小鼠乳腺上皮细胞的体外培养及其对催乳素反应的研究

实现转基因生物乳腺反应器对外源蛋白的高效表达是目前生物制药亟待解决的难题。催乳素对泌乳期乳蛋白的合成与分泌具有重要的调控功能。通过转基因小鼠乳腺上皮细胞模型的建立,研究催乳素如何调控乳蛋白的表达,为提高乳腺反应器高效表达外源蛋白提供技术及理论支撑。应用机械破碎及胶原酶消化法,经差速贴壁纯化,成功培养含人转铁蛋白基因的小鼠乳腺上皮细胞,细胞上清液中检测到人转铁蛋白表达。细胞经牛催乳素诱导后人转铁蛋白的表达水平明显升高。利用转基因小鼠乳腺上皮细胞模型,可以进行催乳素和环境因素等对乳腺上皮细胞合成及分泌蛋白能力影响的研究。     

The influence of extracellular matrix and prolactin on global gene expression profiles of primary bovine mammary epithelial cells in vitro

An in vitro bovine mammosphere model was characterized for use in lactational biology studies using a functional genomics approach. Primary bovine mammary epithelial cells cultured on a basement membrane, Matrigel, formed three-dimensional alveoli-like structures or mammospheres. Gene expression profiling during mammosphere formation by high-density microarray analysis indicated that mammospheres underwent similar molecular and cellular processes to developing alveoli in the mammary gland. Gene expression profiles indicated that genes involved in milk protein and fat biosynthesis were expressed, however, lactose biosynthesis may have been compromised. Investigation of factors influencing mammosphere formation revealed that extracellular matrix (ECM) was responsible for the initiation of this process and that prolactin (Prl) was necessary for high levels of milk protein expression. CSN3 (encoding κ-casein) was the most highly expressed casein gene, followed by CSN1S1 (encoding αS1-casein) and CSN2 (encoding β-casein). Eighteen Prl-responsive genes were identified, including CSN1S1 , SOCS2 and CSN2, however, expression of CSN3 was not significantly increased by Prl and CSN1S2 was not expressed at detectable levels in mammospheres. A number of novel Prl responsive genes were identified, including ECM components and genes involved in differentiation and apoptosis. This mammosphere model is a useful model system for functional genomics studies of certain aspects of dairy cattle lactation.     

ABC- and SLC-Transporters in Murine and Bovine Mammary Epithelium - Effects of Prochloraz

Some chemicals are ligands to efflux transporters which may result in high concentrations in milk. Limited knowledge is available on the influence of maternal exposure to chemicals on the expression and function of transporters in the lactating mammary gland. We determined gene expression of ABC and SLC transporters in murine mammary tissue of different gestation and lactation stages, in murine mammary cells (HC11) featuring resting and secreting phenotypes and in bovine mammary tissue and cells (BME-UV). Effects on transporter expression and function of the imidazole fungicide prochloraz, previously reported to influence BCRP in mammary cells, was investigated on transporter expression and function in the two cell lines. Transporters studied were BCRP, MDR1, MRP1, OATP1A5/OATP1A2, OCTN1 and OCT1. Gene expressions of BCRP and OCT1 in murine mammary glands were increased during gestation and lactation, whereas MDR1, MRP1, OATP1A5 and OCTN1 were decreased, compared to expressions in virgins. All transporters measured in mammary glands of mice were detected in bovine mammary tissue and in HC11 cells, while only MDR1 and MRP1 were detected in BME-UV cells. Prochloraz treatment induced MDR1 gene and protein expression in both differentiated HC11 and BME-UV cells and increased protein function in HC11 cells, resulting in decreased accumulation of the MDR1 substrate digoxin. In conclusion, our results demonstrate that murine (HC11) and bovine (BME-UV) mammary epithelial cells can be applied to characterize expression and function of transporters as well as effects of contaminants on the mammary transporters. An altered expression, induced by a drug or toxic chemical, on any of the transporters expressed in the mammary epithelial cells during lactation may modulate the well-balanced composition of nutrients and/or secretion of contaminants in milk with potential adverse effects on breast-fed infants and dairy consumers.