Polymorphism of the bovine CSN1S1 promoter: linkage mapping, intragenic haplotypes, and effects on milk production traits

Genome-wide identification and annotation of SNPs and their mapping in candidate genes related to milk production and fertility traits in Badri cattle

This study was conducted in Badri cattle using a double digest restriction-site associated DNA sequencing approach. The study aimed to identify and annotate high confidence single nucleotide polymorphisms (SNPs) and their mapping in candidate genes related to production and fertility in dairy cattle. A total of 7,168,552 genome-wide SNPs were initially identified in Badri cattle by alignment with the Bos indicus reference genome. After filtration of SNPs, 65,483 high confidence SNPs were retained and further used for downstream analysis. Annotation of high confidence SNPs revealed 99.197% SNPs had modifier impact, 0.326% SNPs were low impact, 0.036% were high impact, and 0.441% were moderate impact SNPs. Most SNPs in Badri cattle were found in intergenic, transcript and intronic regions. The candidate genes for milk production PRKCE, ABCG2, GHR, EPS8, CAST and NRXN1 were found to harbour maximum high confidence variants. Among candidate genes for fertility in cattle, ATP2B1, SOX5, WDR27, ARHGAP12, CACNA1D, ANKRD6, GRIA3, ZNF521 and CAST822 have maximum high confidence variants mapped in them. The SNPs found mapped in the candidate genes will be important genetic tools in the search for phenotype-modifying nucleotide changes and will aid in formulating relevant genetic improvement programmes for dairy cattle.

A New AS-PCR Method to Detect CSN201 Allele, Genotyping at Ca-Sensitive Caseins Loci and Milk Traits Association Studies in Autochthonous Lazio Goats

Calcium-sensitive caseins are the main protein component of milk. In the goat, they are encoded by three genes (CSN1S1, CSN2, and CSN1S2) located on chromosome 6. A high number of alleles has been discovered for these genes in the goat species, responsible for changes in the milk’s qualitative and quantitative characteristics. This study aimed to develop an Allele-Specific PCR (AS-PCR), which allowed us to unequivocally detect goat carriers of the CSN201 allele. Subsequently, the calcium-sensitive casein loci genotype was investigated in three native goat breeds of the Lazio Region (Bianca Monticellana, Capestrina, and Ciociara Grigia). No individuals were carriers of the CSN1S101, CSN1S1E, CSN201, CSN1S2D, and CSN1S20 alleles, while a high frequency of the alleles CSN1S1F and CSN1S1A*,B* was observed. Association analyses between the different genotypes at the CSN1S1 locus and some milk traits, namely the fat and protein yielded and the fat, protein, solids-not-fat, and casein percentages without an effect on the milk yield, were observed.

GWAS and genomic prediction of milk urea nitrogen in Australian and New Zealand dairy cattle

Background

Urinary nitrogen leakage is an environmental concern in dairy cattle. Selection for reduced urinary nitrogen leakage may be done using indicator traits such as milk urea nitrogen (MUN). The result of a previous study indicated that the genetic correlation between MUN in Australia (AUS) and MUN in New Zealand (NZL) was only low to moderate (between 0.14 and 0.58). In this context, an alternative is to select sequence variants based on genome-wide association studies (GWAS) with a view to improve genomic prediction accuracies. A GWAS can also be used to detect quantitative trait loci (QTL) associated with MUN. Therefore, our objectives were to perform within-country GWAS and a meta-GWAS for MUN using records from up to 33,873 dairy cows and imputed whole-genome sequence data, to compare QTL detected in the GWAS for MUN in AUS and NZL, and to use sequence variants selected from the meta-GWAS to improve the prediction accuracy for MUN based on a joint AUS-NZL reference set.

Results

Using the meta-GWAS, we detected 14 QTL for MUN, located on chromosomes 1, 6, 11, 14, 19, 22, 26 and the X chromosome. The three most significant QTL encompassed the casein genes on chromosome 6, PAEP on chromosome 11 and DGAT1 on chromosome 14. We selected 50,000 sequence variants that had the same direction of effect for MUN in AUS and MUN in NZL and that were most significant in the meta-analysis for the GWAS. The selected sequence variants yielded a genetic correlation between MUN in AUS and MUN in NZL of 0.95 and substantially increased prediction accuracy in both countries.

Conclusions

Our results demonstrate how the sharing of data between two countries can increase the power of a GWAS and increase the accuracy of genomic prediction using a multi-country reference population and sequence variants selected based on a meta-GWAS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00707-9.

First DNA Sequencing in Beninese Indigenous Cattle Breeds Captures New Milk Protein Variants

This study investigated polymorphisms in the milk protein genes CSN1S1, CSN2, CSN1S2, CSN3, LALBA, and LGB, and casein haplotypes in Beninese indigenous cattle. Considering 67 animals, DNA sequencing of the genes' exons, flanking regions and parts of the 5'-upstream regions identified 1058 genetic variants including 731 previously unknown. In addition, four novel milk protein variants were detected, including CSN3K (p.Ala66Val), LALBAF (p.Arg58Trp), LGBB1 (p.Ala134Val) and LGBK (p.Thr92Asnfs*13). CSN3K is caused by a novel SNP (BTA6:85656526C>T, exon 4) whereas LALBAF and LGBB1 are due to rs714688595C>T (exon 1) and rs109625649C>T (exon 4), respectively. Regarding LGBK, a frameshift insertion of one adenine residue at BTA11:103257980 (exon 3) induces a premature translation termination resulting in a 46% reduction of the reference protein sequence. The casein polymorphisms formed five main CSN1S1-CSN2-CSN1S2-CSN3 haplotypes including B-A1-A-B, B-A1-A-A and C-A2-A-B which are predominant in the investigated cattle breeds. Moreover, in silico analyses of polymorphisms within the 5'- and 3'- untranslated regions of all six milk proteins revealed effects on microRNA and transcription factor binding sites. This study suggests a large genetic variation of milk protein genes in Beninese cattle, which should be investigated in further studies for their effects on milk production, including quality and yield traits.

Casein (CSN) gene variants and parity affect the milk protein traits in crossbred (Bos taurus x Bos indicus) cows in sub-tropical climate

Genotypes at four casein (CSN) loci-A26181G of CSN1S1, C6227T of CSN1S2, A8101C of CSN2, and A13104C of CSN3-along with non-genetic factors were studied for their effects on various milk protein traits in 100 crossbred cows with major inheritance of Holstein Frisian (Bos taurus) and Tharparkar (Bos indicus). Results demonstrated the presence of all CSN genotypes with a predominance of heterozygotes. At CSN2 (A8101C; His67Pro) locus, the A2 allele, desirable for human health, was present in 62% as heterozygous and 29% in homozygous condition. Among non-genetic factors, parity of the cows had a significant influence on the milk protein traits in these crossbreds. The genotypes at CSN1S1, CSN2, and CSN3 loci were found to influence (p<0.05 to 0.01) the casein and whey protein yields and composition traits. The casein index and total milk yield were most influenced by the CSN1S2 locus. The AA (A1 milk) genotype of CSN2 had significantly higher yields and percentages of casein and whey proteins. Positive influence of CC genotype of CSNS3 on milk proteins of was observed similar to Bos taurus cows; however, such influence of AA genotype of CSN2 locus may be distinctive to the crossbred cows maintained in subtropical condition. Overall, the results revealed the diverse effects of CSN genotypes on milk proteins in crossbred cattle.

Differential Expression of the Alpha S1 Casein and Beta-Lactoglobulin Genes in Different Physiological Stages of the Adani Goats Mammary Glands

Background

Milk proteins genes have been the focus of the researches as the candidate target genes that play a decisive role when animal breeding is desired.

Objectives

In the present study, the transcriptional levels of Beta-lactoglobulin (BLG) and Alpha S1 casein (CSN1S1) genes were investigated during prenatal, milking and drying times in mammary glands of the Adani goats which showed high and low breeding values.

Materials and Methods

The breeding values of the animals were estimated first by applying multi-trait random regression model. Using the biopsy gun, the mammary gland samples were taken and real-time PCR was applied to search the expression of the genes. Fixed factors of the model were the breeding value groups, sampling times and their interactions.

Results

The interactions were significant for both genes. At milking time, the high breeding value group exhibited more transcriptional levels for BLG and less transcriptional levels for CSN1S1 gene compared with the low breeding value group. The expression patterns of these genes were also different between the two breeding value groups. The maximum level of BLG and CSN1S1 transcriptions were found to occur at drying time.

Conclusions

A difference in the gene expression was observed between the two groups which indicate the change in the nucleotide sequence for transcription factor binding sites, or miRNA binding sites, otherwise in the coding regions. Therefore, the variations in the coding and promoter regions of this gene should be investigated in the further studies.

Genome-wide association study for longevity with whole-genome sequencing in 3 cattle breeds

Longevity is an important economic trait in dairy production. Improvements in longevity could increase the average number of lactations per cow, thereby affecting the profitability of the dairy cattle industry. Improved longevity for cows reduces the replacement cost of stock and enables animals to achieve the highest production period. Moreover, longevity is an indirect indicator of animal welfare. Using whole-genome sequencing variants in 3 dairy cattle breeds, we carried out an association study and identified 7 genomic regions in Holstein and 5 regions in Red Dairy Cattle that were associated with longevity. Meta-analyses of 3 breeds revealed 2 significant genomic regions, located on chromosomes 6 (META-CHR6-88MB) and 18 (META-CHR18-58MB). META-CHR6-88MB overlaps with 2 known genes: neuropeptide G-protein coupled receptor (NPFFR2; 89,052,210-89,059,348 bp) and vitamin D-binding protein precursor (GC; 88,695,940-88,739,180 bp). The NPFFR2 gene was previously identified as a candidate gene for mastitis resistance. META-CHR18-58MB overlaps with zinc finger protein 717 (ZNF717; 58,130,465-58,141,877 bp) and zinc finger protein 613 (ZNF613; 58,115,782-58,117,110 bp), which have been associated with calving difficulties. Information on longevity-associated genomic regions could be used to find causal genes/variants influencing longevity and exploited to improve the reliability of genomic prediction.

Estimates of epistatic and pleiotropic effects of () and () genetic markers on beef heifer performance traits enhanced by selection

Genetic marker effects and type of inheritance are estimated with poor precision when minor marker allele frequencies are low. A stable composite population (MARC II) was subjected to marker assisted selection for 2 yr to equalize and genetic marker frequencies to evaluate the epistatic and pleiotropic effects of these markers on BW, reproduction, and first calf performance traits in replacement beef females ( = 171) managed under 2 postweaning development protocols. Traits evaluated on the heifers were birth BW, weaning BW, 11-mo BW, 12-mo BW, 13-mo BW, first breeding season pregnancy evaluation BW, first calving season BW, 11-mo puberty, 12-mo puberty, 13-mo puberty, first breeding season pregnancy, and first calf weaning rate. Additionally, heifer's first calf performance traits of ordinal calving date, first calf birth BW, and first calf weaning BW (with and without age adjustment) were analyzed. Selection to increase minor allele frequencies and balanced sampling across genotype classes enhanced the ability to detect all genetic effects except dominance × dominance epistasis. The × genotype effect was significant ( < 0.05) for 11-mo BW and 12-mo BW and tended to be significant ( = 0.08) for 13-mo BW. Consistently, for all 3 traits, the most significant effect among epistatic × genotype effects was the additive effect, with the G allele decreasing BW. There were no associations between × genotype and fertility related traits ( ≥ 0.46) in this study. Additionally, there were no × genotype associations with first progeny performance traits ( ≥ 0.14). The large effect of the additive × additive interaction on first calf weaning BW was imprecisely estimated, which may warrant further investigation.

Mediterranean river buffalo CSN1S1 gene: search for polymorphisms and association studies

The aim of the present study was to investigate the variability at CSN1S1 locus of the Italian Mediterranean river buffalo and to study possible allele effects on milk yield and its composition. Effects of parity, calving season and month of production were also evaluated. Three single-nucleotide polymorphisms were detected. The first mutation, located at position 89 of the 17th exon (c.628C>T), is responsible for the amino acid change p.Ser178 (B allele)/Leu178 (A allele). The other two polymorphisms, detected at the positions 144 (c.882G>A) and 239 (c.977A>G) of 19th exon, respectively, are silent (3ʹ UTR, untranslated region). Associations between the CSN1S1 genotypes and milk production traits were investigated using 4122 test day records of 503 lactations from 175 buffalo cows. Milk yield, fat and protein percentages were analysed using a mixed linear model. A significant association between the c.628C>T SNP and the protein percentage was found. In particular, the CC genotype showed an average value ~0.04% higher than the CT and TT genotypes. The allele substitution effect of cytosine into thymine was –0.014, with a quite low (0.3%) protein percentage contribution to total phenotypic variance. A large dominance effect was detected. Characterisation of the CSN1S1 transcripts and a method based on MboI amplification created restriction site PCR for a rapid genotyping of c.628C>T are provided.

The β-casein in camels: molecular characterization of the CSN2 gene, promoter analysis and genetic variability

The β-casein is the most abundant protein in camel milk and its encoding gene (CSN2) is considered in other species a 'major' gene for the presence of alleles associated to different level of expression. In the present paper, we report for the first time the characterization of the nucleotide sequence of the whole β-casein-encoding gene (CSN2) plus 2,141 bp at the 5'-flanking region in Camelus dromedarius. The promoter region and the complete cDNA are also provided for the first time in Camelus bactrianus. The gene is spread over 7.8 kb and consists of 9 exons varying in length from 24 bp (exon 5) to 519 bp (exon 7), and 8 introns from 95 bp (intron 5) to 1,950 bp (intron 1). The composite response element (CoRE) region was identified in the promoter, whereas the presence of mature microRNA sequences improves the knowledge on the factors putatively involved in the gene regulation. A total of 46 polymorphic sites have been detected. The transition g.2126A>G falls within the TATA-box of dromedary CSN2 promoter with a putative influence on the transcription factor binding activity. The frequency of the G allele is 0.35 in a population of 180 she-camels belonging to 4 different ecotypes. In the same population, a conservative SNP (g.4175C>A) was found at the codon 7 of the signal peptide, whereas a comparative analysis with a cDNA sequence available in the database evidenced a missense SNP (g.4180T(Leu)>G(Arg)) at exon 2. Four SNPs were found in the bactrian camel. The SNP c.666G>A is responsible for the amino acid change Met(201)→Ile and it represents the first missense allele at the β-casein in camels. Finally, five interspersed repeated elements were identified at intronic level, whereas the presence of putative bio-functional peptides belonging to ACE-inhibitor and anti-oxidative families confirms the potential protective role of the camel milk for the human nutrition.

Using genome-wide association analysis to characterize environmental sensitivity of milk traits in dairy cattle

Genotype-by-environment interaction (GxE) has been widely reported in dairy cattle. One way to analyze GxE is to apply reaction norm models. The first derivative of a reaction norm is the environmental sensitivity (ES). In the present study we conducted a large-scale, genome-wide association analysis to identify single-nucleotide polymorphisms (SNPs) that affect general production (GP) and ES of milk traits in the German Holstein population. Sire estimates for GP and for ES were calculated from approximately 13 million daughter records by the use of linear reaction norm models. The daughters were offspring from 2297 sires. Sires were genotyped for 54k SNPs. The environment was defined as the average milk energy yield performance of the herds at the time during which the daughter observations were recorded. The sire estimates were used as observations in a genome-wide association analysis, using 1797 sires. Significant SNPs were confirmed in an independent validation set (500 sires of the same population). To separate GxE scaling and other GxE effects, the observations were log-transformed in some analyses. Results from the reaction norm model revealed GxE effects. Numerous significant SNPs were validated for both GP and ES. Many SNPs that affect GP also affect ES. We showed that ES of milk traits is a typical quantitative trait, genetically controlled by many genes with small effects and few genes with larger effect. A log-transformation of the observation resulted in a reduced number of validated SNPs for ES, pointing to genes that not only caused scaling GxE effects. The results will have implications for breeding for robustness in dairy cattle.

Sequence analysis and identification of new variations in the 5'-flanking region of αS2-casein gene in Indian zebu cattle

Regulatory region of milk protein alpha S2-casein (αS2-CN) gene sequence was characterized and analyzed for nucleotide variations in animals representing 13 Indian zebu cattle (Bos indicus) breeds. A total of 15 variations; 11 in promoter region (1.56 Kb): -1481 (C>T), -1412 (C>T), -1342 (C>T), -1084 (G>A), -979 (A>G), -657 (A>T), -508 (A>G), -186 (T>C), -184 (T>C), -151 (T>C) and -135 (C>T); 1 in 5'-UTR (44 bp): 7 (C>T) while, 3 in intron-I region (73 bp): 186 (C>T), 194 (A>C) and 301 (A>T) were identified. Additionally, single deletion was observed at -975 (A>-) but not involve any known potential transcription factor binding sites (TFBS). Comparison with Bos taurus sequence revealed two additional variations -1085 (T>C) and -739 (A>G). Out of the total 18 variations observed between indicine and taurine αS2-CN regulatory region sequence, 15 were novel to B. indicus and are reported for the first time. Among these, four variations were located within the potential TFBSs; -1342 (C>T) within HNF-3beta, -739 (A>G) within C/EBP-alpha while -657 (A>T) and -508 (A>G) were found within glucocorticoid receptor TFBSs. Variations located within or in proximity to putative TFBSs could possibly influence the binding affinity of nuclear factors towards DNA binding domains, thus affecting transcriptional rate of αS2-CN gene. Phylogenetically, as expected, Indian zebu cattle were grouped close to B. taurus and were most distantly placed in comparison to human. The study indicated possible genetic variations in the regulatory regions of αS2-CN gene within Indian native cattle (B. indicus) and also its comparison with evolutionary different B. taurus breeds.

Variations in the Regulatory Region of Alpha S1-Casein Milk Protein Gene among Tropically Adapted Indian Native (Bos Indicus) Cattle

Regulatory region of milk protein alpha S1-casein (αS1-CN) gene was sequenced, characterized, and analyzed to detect variations among 13 Indian cattle (Bos indicus) breeds. Comparative analysis of 1,587 bp region comprising promoter (1,418 bp), exon-I (53 bp), and partial intron-I (116 bp) revealed 35 nucleotide substitutions (32 within promoter region, 1 in exon-I, and 2 in partial intron-I region) and 4 Indels. Within promoter, 15 variations at positions −1399 (A > G), −1288 (G > A), −1259 (T > C), −1158 (T > C), −1016 (A > T), −941 (T > G), −778 (C > T), −610 (G > A), −536 (A > G), −521 (A > G), −330 (A > C), −214 (A > G), −205 (A > T), −206 (C > A), and −175 (A > G) were located within the potential transcription factor binding sites (TFBSs), namely, NF-κE1/c-Myc, GATA-1, GATA-1/NF-E, Oct-1/POU3F2, MEF-2/YY1, GATA-1, AP-1, POU1F1a/GR, TMF, GAL4, YY1/Oct-1, HNF-1, GRalpha/AR, GRalpha/AR, and AP-1, respectively. Seventy-four percent (26/35) of the observed SNPs were novel to Indian cattle and 11 of these novel SNPs were located within one or more TFBSs. Collectively, these might influence the binding affinity towards their respective nuclear TFs thus modulating the level of transcripts in milk and affecting overall protein composition. The study provides information on several distinct variations across indicine and taurine αS1-CN regulatory domains.

Selection for genetic markers in beef cattle reveals complex associations of thyroglobulin and casein1-s1 with carcass and meat traits

Genetic markers in casein (CSN1S1) and thyroglobulin (TG) genes have previously been associated with fat distribution in cattle. Determining the nature of these genetic associations (additive, recessive, or dominant) has been difficult, because both markers have small minor allele frequencies in most beef cattle populations. This results in few animals homozygous for the minor alleles. selection to increase the frequencies of the minor alleles for 2 SNP markers in these genes was undertaken in a composite population. The objective was to obtain better estimates of genetic effects associated with these markers and determine if there were epistatic interactions. Selection increased the frequencies of minor alleles for both SNP from <0.30 to 0.45. Bulls (n = 24) heterozygous for both SNP were used in 3 yr to produce 204 steer progeny harvested at an average age of 474 d. The combined effect of the 9 CSN1S1 × TG genotypes was associated with carcass-adjusted fat thickness (P < 0.06) and meat tenderness predicted at the abattoir by visible and near-infrared reflectance spectroscopy (P < 0.04). Genotype did not affect BW from birth through harvest, ribeye area, marbling score, slice shear force, or image-based yield grade (P > 0.10). Additive, dominance, and epistatic SNP association effects were estimated from genotypic effects for adjusted fat thickness and predicted meat tenderness. Adjusted fat thickness showed a dominance association with TG SNP (P < 0.06) and an epistatic additive CSN1S1 × additive TG association (P < 0.03). For predicted meat tenderness, heterozygous TG meat was more tender than meat from either homozygote (P < 0.002). Dominance and epistatic associations can result in different SNP allele substitution effects in populations where SNP have the same linkage disequilibrium with causal mutations but have different frequencies. Although the complex associations estimated in this study would contribute little to within-population selection response, they could be important for marker-assisted management or reciprocal selection schemes.

The milk protein α-casein functions as a tumor suppressor via activation of STAT1 signaling, effectively preventing breast cancer tumor growth and metastasis

Here, we identified the milk protein α-casein as a novel suppressor of tumor growth and metastasis. Briefly, Met-1 mammary tumor cells expressing α-casein showed a ~5-fold reduction in tumor growth and a near 10-fold decrease in experimental metastasis. To identify the molecular mechanism(s), we performed genome-wide transcriptional profiling. Interestingly, our results show that α-casein upregulates gene transcripts associated with interferon/STAT1 signaling and downregulates genes associated with "stemness." These findings were validated by immunoblot and FACS analysis, which showed the upregulation and hyperactivation of STAT1 and a decrease in the number of CD44(+) "cancer stem cells." These gene signatures were also able to predict clinical outcome in human breast cancer patients. Thus, we conclude that a lactation-based therapeutic strategy using recombinant α-casein would provide a more natural and non-toxic approach to the development of novel anticancer therapies.

Invited review: milk protein polymorphisms in cattle: effect on animal breeding and human nutrition

The 6 main milk proteins in cattle are encoded by highly polymorphic genes characterized by several nonsynonymous and synonymous mutations, with up to 47 protein variants identified. Such an extensive variation was used for linkage analysis with the description of the casein cluster more than 30 yr ago and has been applied to animal breeding for several years. Casein haplotype effects on productive traits have been investigated considering information on the whole casein complex. Moreover, mutations within the noncoding sequences have been shown to affect the specific protein expression and, as a consequence, milk composition and cheesemaking. Milk protein variants are also a useful tool for breed characterization, diversity, and phylogenetic studies. In addition, they are involved in various aspects of human nutrition. First, the occurrence of alleles associated with a reduced content of different caseins might be exploited for the production of milk with particular nutritional qualities; that is, hypoallergenic milk. On the other hand, the frequency of these alleles can be decreased by selection of sires using simple DNA tests, thereby increasing the casein content in milk used for cheesemaking. Furthermore, the biological activity of peptides released from milk protein digestion can be affected by amino acid exchanges or deletions resulting from gene mutations. Finally, the gene-culture coevolution between cattle milk protein genes and human lactase genes, which has been recently highlighted, is impressive proof of the nonrandom occurrence of milk protein genetic variation over the centuries.

A single nucleotide polymorphism and sequence analysis of CSN1S1 gene promoter region in Chinese Bos grunniens (yak)

The aim of this study was to investigate the polymorphism of the CSN1S1 gene promoter region in 4 Chinese yak breeds, and compare the yak CSN1S1 gene promoter region sequences with other ruminants. A Polymerase Chain Reaction-Single Strand Conformation Polymorphism protocol was developed for rapid genotyping of the yak CSN1S1 gene. One hundred fifty-eight animals from 4 Chinese yak breeds were genotyped at the CSN1S1 locus using the protocol developed. A single nucleotide polymorphism of the CSN1S1 gene promoter region has been identified in all yak breeds investigated. The polymorphism consists of a single nucleotide substitution G-->A at position 386 of the CSN1S1 gene promoter region, resulting in two alleles named, respectively, G(386) and A(386), based on the nucleotide at position 386. The allele G(386) was found to be more common in the animals investigated. The corresponding nucleotide sequences in GenBank of yak (having the same nucleotides as allele G(386) in this study), bovine, water buffalo, sheep, and goat had similarity of 99.68%, 99.35%, 97.42%, 95.14%, and 94.19%, respectively, with the yak allele A(386.).

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.

Loss of caveolin-3 induces a lactogenic microenvironment that is protective against mammary tumor formation

Here, we show that functional loss of a single gene is sufficient to confer constitutive milk protein production and protection against mammary tumor formation. Caveolin-3 (Cav-3), a muscle-specific caveolin-related gene, is highly expressed in muscle cells. We demonstrate that Cav-3 is also expressed in myoepithelial cells within the mammary gland. To determine whether genetic ablation of Cav-3 expression affects adult mammary gland development, we studied the phenotype(s) of Cav-3(-/-)-null mice. Interestingly, Cav-3(-/-) virgin mammary glands developed lobulo-alveolar hyperplasia, akin to the changes normally observed during pregnancy and lactation. Genome-wide expression profiling revealed up-regulation of gene transcripts associated with pregnancy/lactation, mammary stem cells, and human breast cancers, consistent with a constitutive lactogenic phenotype. Expression levels of three key transcriptional regulators of lactation, namely Elf5, Stat5a, and c-Myc, were also significantly elevated. Experiments with pregnant mice directly showed that Cav-3(-/-) mice underwent precocious lactation. Finally, using orthotopic tumor cell implantation, we demonstrated that virgin Cav-3(-/-) mice were dramatically protected against mammary tumor formation. Thus, Cav-3(-/-) mice are a novel preclinical model to study the protective effects of a lactogenic microenvironment on mammary tumor onset and progression. Our current studies have broad implications for using the lactogenic microenvironment as a paradigm to discover new therapies for the prevention and/or treatment of human breast cancers.

A critical analysis of production-associated DNA polymorphisms in the genes of cattle, goat, sheep, and pig

Increasing productivity is one of the main objectives in animal production. Traditional breeding methods have led to increased gains in some traits but gains are not easily attainable in traits with low heritabilities. Exploiting the genetic variations underlying desired phenotypes is the goal of today's animal producers. Such positive genetic variants must, however, be known before possible application. Consequently, candidate genes of traits of interest have been searched for possible relationships with such traits or to explain reported quantitative trait loci (QTL) for such traits. DNA variants or polymorphisms have been identified in many such genes and their relationships with production traits determined. However, only a few genes have been evaluated, given the wealth of information on reported QTL for production traits, and in most cases genes are only partially investigated. This review presents available information on DNA variants for production traits and discusses steps that are required for effective utilization of this information for successful marker-assisted selection programs.

Molecular characterization of bovine CSN1S2*B and extensive distribution of zebu-specific milk protein alleles in European cattle

The B allele of the bovine alpha (S2)-casein gene (CSN1S2) was characterized at the molecular level and the distribution of zebu-specific milk protein alleles was determined in 26 cattle breeds originating from 3 continents. The CSN1S2*B allele is characterized by a C --> T transition affecting nucleotide 17 of exon 3, which leads to a change in the eighth amino acid of the mature protein, from Ser to Phe (i.e., TCC --> TTC). DNA-based methods were developed to identify carriers of CSN1S2*B and the other alleles (CSN1S2*A, C, and D) at the same locus. CSN1S2*B and other zebu-specific milk protein alleles and casein haplotypes are widely distributed in European cattle breeds, particularly those of southeastern origin. Alleles CSN1S2*B and CSN3*H are important in searching for zebu imprints in European cattle breeds. Diversity estimates at the milk protein loci were highest in the zebus followed by southeastern European taurines. Anatolian Black had the highest number of zebu alleles among European taurines. Common, group, and intergroup haplotypes occurred in the breeds and demonstrated relationships that concurred with developmental histories, genetic makeup, and, in particular, exposed the extent of zebu influence on southeastern European cattle.

Short Communication: Composition of Milk Protein and Milk Fatty Acids Is Stable for Cows Differing in Genetic Merit for Milk Production

Changing the composition of milk protein and of milk fatty acids alters nutritional and physical properties of dairy products and their consumer appeal. Genetic selection for milk yield decreases concentrations of milk protein and of milk fat. Little is known, however, about how the decrease affects composition of milk protein and milk fatty acids. The objective of this study was to quantify changes in composition of milk protein and of milk fatty acids in cows differing in genetic merit for milk production. Three measures of genetic merit for milk production were used for each cow: genetic line, parent average predicted transmitting ability (PTA) for milk, and cow milk PTA. Composition of milk protein and milk fatty acids were compared in 448 milk samples from 178 cows representing 2 divergent lines of Holsteins that were bred for high or average PTA for milk and combined milk protein and fat yield. High-line cows (n = 97) produced more milk that contained less fat and had higher proportions of alphaS1-casein in milk protein than did average-line cows (n = 81). We additionally obtained from 233 cows (178 cows representing the 2 genetic lines and 55 cows with ancestors from both genetic lines) the parent average milk PTA and cow milk PTA and compared composition of milk protein and of milk fatty acids in 592 milk samples. Cows whose parent average milk PTA was above or equal to the median of the 233 cows produced more milk that contained less protein and less fat and that tended to have greater proportions of alphaS1-casein in milk protein than cows whose average milk PTA was below the median. Similarly, cows with above or equal median milk PTA of the 233 cows produced more milk that contained less protein and less fat and had greater proportions of alphaS1-casein in milk protein than did cows with below-median milk PTA. Milk fatty acid composition was not consistently different between groups. Therefore, selection for milk yield decreased concentrations of milk protein and milk fat but had little effect on composition of milk protein and milk fatty acids.

Evaluation in beef cattle of six deoxyribonucleic acid markers developed for dairy traits reveals an osteopontin polymorphism associated with postweaning growth

The objectives of this study were to 1) estimate the allelic frequencies in US beef cattle of 6 DNA markers reported to be associated with variation in dairy production traits; and 2) evaluate the association of these markers with beef production traits. Several genetic markers have been associated with milk yield or composition, including polymorphisms in secreted phosphoprotein 1 (SPP1; also called osteopontin), growth hormone receptor (GHR), casein S1 (CSN1S1), diacylglycerol O-acyltransferase 1 (DGAT1), peroxisome proliferator-activated receptor gamma co-activator-1alpha (PPARGC1A), and ATP-binding cassette subfamily G (white) member 2 (ABCG2). Allelic frequencies for these 6 markers, and their association with 21 phenotypes, were evaluated in 2 crossbred beef cattle populations that sample influential industry sires. Five of 6 markers were segregating in beef cattle populations; the exception was ABCG2. The SPP1 marker was associated with yearling weight (P = 0.025), live weight at slaughter (P = 0.016), postweaning ADG (P = 0.007), and HCW (P = 0.007) in a large, multisire population representing the 7 most populous beef breeds in the United States. Postweaning growth trait associations were confirmed in an independent population of similar construction, including sires from tropically adapted breeds. The SPP1 marker was associated with yearling weight (P = 0.034), live weight at slaughter (P = 0.011), and postweaning ADG (P = 0.015) and showed a trend toward association with HCW (P = 0.083) in this population. Whereas DGAT1, GHR, and CSN1S1 polymorphisms showed association with some traits in individual populations, the lack of consistent predictive merit between populations indicates they may not be suited for beef cattle selection. No significant associations were observed for the PPARGC1A marker and any of 21 recorded traits, indicating this marker had no apparent value in selection for the beef cattle traits tested in these populations. The SPP1 marker had consistent associations and effect sizes (10.5 to 11.5 kg of live weight at slaughter) in both populations, providing strong evidence for utility of the SPP1 marker for postweaning growth in beef cattle.

An evaluation of genetic diversity indices of the Red Bororo and White Fulani cattle breeds with different molecular markers and their implications for current and future improvement options

The genetic diversity of the Red Bororo and White Fulani cattle breeds of Cameroon and Nigeria was assessed with a panel of 32 markers. Estimates for the various indices of genetic diversity, total number of alleles (TNA), mean observed number of alleles (MNA), mean effective number of alleles (MNE), observed heterozygosity (Hob) and expected heterozygosity (Hex), were higher at microsatellite loci than at protein loci. Mean Hex values were above 71% at microsatellite loci in all the breeds and ranged from 37% to 41.6% at milk protein loci and from 40.9% to 45.6% at blood protein loci. The highest TNA and MNA of microsatellites were recorded for the Nigerian White Fulani. MNE of milk protein loci was highest in the Cameroonian Red Bororo, while TNA of blood protein loci was highest in the Cameroonian White Fulani. The high genetic diversity levels indicate the presence of the necessary ingredients for improvement breeding and conservation. Multi-locus estimates of within-population inbreeding (f), total inbreeding (F) and population differentiation (theta) of the breeds were significantly different from zero, except for theta of blood proteins. A high level of gene flow was found between the breeds (5.829). The phylogenetic relationship existing among the four breeds is greatly influenced by location. The high gene flow between the breeds may lead to a loss of genetic diversity through genetic uniformity and a reduction in opportunities for future breed development. We propose an improvement scheme with aims to prevent loss of genetic diversity, improve productivity and reduce uncontrolled genetic exchanges between breeds.

Development of a Single Nucleotide Polymorphism Genotyping Microarray Platform for the Identification of Bovine Milk Protein Genetic Polymorphisms

The objective of this study was to develop and validate a fast method for typing the main mutations of bovine milk protein genes by using microarray technology. An approach based on the ligation detection reaction (LDR) and a universal array (UA) was used. Polymorphisms in both the coding and noncoding sequences of alpha(S1)-casein, beta-casein, kappa-casein, and beta-lactoglobulin genes were considered because of their well-known effects on milk composition and cheese production. A total of 22 polymorphic sites, corresponding to 21 different variants, were included in the diagnostic microarray. First, a multiplex PCR was developed to amplify all the DNA target sequences simultaneously. Second, the LDR-UA assay was implemented. The method was validated by analyzing 100 Italian Friesian DNA samples, which were also genotyped by conventional methods both at the protein level by means of milk isoelectrofocusing and at the molecular level using PCR-RFLP and PCR-single strand conformation polymorphism techniques. The genotypes obtained using the LDR-UA approach were in full agreement with those obtained by the conventional analyses. An important result of the LDR-UA assay was a more accurate genotyping of the different milk protein alleles than was found with conventional typing methods. At the kappa-casein gene, in fact, 4 samples were heterozygous (3 reference samples and 1 validation sample) for an allele coding for Thr(136) and Ala(148). This variant, which can be considered as the wild type of the genus Bos, is not usually identifiable by the conventional typing methods used. The multiplex PCR-LDR-UA approach developed provides for an accurate, inexpensive, and high-throughput assay that does not exhibit false positive or false negative signals, thus making it highly suitable for animal genotyping.

Characterization of 37 Breed-Specific Single-Nucleotide Polymorphisms in Sheep

We identified 37 single-nucleotide polymorphisms (SNPs) in sheep and screened 16 individuals from 8 different sheep breeds selected throughout Europe. Population genetic measures based on the genotyping of about 30 sheep from the same 8 breeds are reported. To date, there are no sheep SNPs documented in the National Center for Biotechnology Information dbSNP database. Therefore, the markers presented here contribute significantly to those currently available.

Characterization of the DGAT1 Mutations and the CSN1S1 Promoter in the German Angeln Dairy Cattle Population

The identification of quantitative trait loci (QTL) and genes with influence on milk production traits has been the objective of various mapping studies in the last decade. In the centromeric region of Bos taurus autosome (BTA) 14, the acyl-CoA:diacylglycerol acyltransferase1 gene (DGAT1) has been identified as the most likely causative gene underlying a QTL for milk fat yield and content. Recently, a second polymorphism in the promoter of DGAT1 emerged as an additional source of variation. In this study, the frequencies and the effects of alleles at the DGAT1 K232A and at the DGAT1 promoter variable number of tandem repeat (VNTR) locus on BTA14, and of alleles at the CSN1S1 (alpha(S1)-casein-encoding gene) promoter on BTA6 in the German Angeln dairy cattle population were investigated. Analyzed traits were milk, fat, protein, lactose, and milk energy yield, fat, protein, lactose, and milk energy content and somatic cell score. The lysine variant of the DGAT1 K232A mutation showed significant effects for most of the milk production traits. A specific allele of the DGAT1 promoter VNTR showed significant effects on the traits lactose yield and content, milk energy content, and SCS compared with the other alleles. Additionally, a regulation mechanism between the DGAT1 K232A mutation and the DGAT1 promoter VNTR was found for fat yield and content, which could be caused by an upper physiological bound for the effects of the DGAT1 gene. At the CSN1S1 promoter, 2 of 4 alleles showed significant allele substitution effects on the milk yield traits.

Detection of Quantitative Trait Loci Affecting Milk Production Traits on Bovine Chromosome 6 in a Chinese Holstein Population by the Daughter Design

Fourteen microsatellite markers with a coverage of 63.5 cM on bovine chromosome 6 were selected, and 26 sire families with 2,260 daughters were analyzed for mapping quantitative trait loci (QTL) affecting 5 milk production traits in a Chinese Holstein population. In the analyses across 26 families and within the largest significant families with a one-QTL model fitted, a QTL near BMS470 was detected that affected fat yield at the 5% experiment-wide significance level. When a 2-QTL model was fitted in the across-family analysis, it was found that there might exist 2 QTL affecting the 3 yield traits, although the exact or empirical thresholds for the significance testing were unknown. In all analyses, the results for milk yield and protein yield were generally consistent, which might have resulted from the same genetic background for milk and protein yield.

Polymorphic AP-1 Binding Site in Bovine CSN1S1 Shows Quantitative Differences in Protein Binding Associated with Milk Protein Expression

Polymorphisms in 5'-flanking regions of milk protein encoding genes can influence the binding activity of the affected response elements and thus have an impact on the expression of the gene products. However, precise quantitative data concerning the binding properties of such variable response elements have so far not been described. In this study we present the results of a quantitative fluorescent electromobility shift assay comparing the allelic variants of a polymorphic activator protein-1 binding site in the promoter region of the bovine alphas1-casein encoding gene (CSN1S1), which is affected by an A-->G exchange at -175 bp (CSN1S1(-175bp)). A supershift assay using a commercial c-jun antibody was carried out to verify the specificity of protein binding. The gel shift analysis revealed specific and significantly reduced protein binding of oligonucleotides containing the G variant of the CSN1S1(-175bp) binding site. Further investigations comprised genotyping of the variable CSN1S1(-175bp) activator protein-1 element by an NmuCl restriction fragment length polymorphism in 62 cows of the breed Simmental and 80 cows of the breed German Holstein. Single milk proteins from at least 4 milk samples per cow were quantified by alkaline urea polyacrylamide gel electrophoresis. Homozygotes for CSN1S1(-175bp)*G were not observed, and the allele frequencies were 0.19 in Simmental and 0.05 in German Holstein. Carriers of CSN1S1(-175bp)*G showed higher content (%) as well as quantity (g/d) of alphas1-casein than CSN1S1(-175bp)*A homozygotes, independent of breed. We assume that the positive association of the CSN1S1(-175bp)*G variant with CSN1S1 expression is likely to be caused by a reduced affinity of the affected response element to a c-jun-containing CSN1S1 dimer with repressor properties.

Effects of casein haplotypes on milk production traits in Italian Holstein and Brown Swiss cattle

The objective of this study was to estimate the effects of different haplotypes of the casein genes on milk production traits in Italian dairy cattle. Traits of interest were yields of milk, fat, and protein, and percentages of fat and protein in milk. The data included 728 multiparous records from 347 Holsteins and 773 records from 298 Brown Swiss cows. Records were preadjusted for effects of age and parity, season of calving, and region, and expressed as deviations from herdmate averages. Twenty half-sib families were represented in each breed. Haplotype probabilities were estimated for each animal and phenotypes were regressed on these probabilities. Nine haplotypes were observed in Holsteins and 17 were identified among the Brown Swiss. For Holsteins, significant effects were observed for protein percentage, with some indication of an effect for fat percentage. For the Brown Swiss, effects of haplotypes were significant for milk yield and fat and protein percentages. Effects were strongest for protein percentage. Correlation coefficients of solutions across breeds tended to be strong and positive, indicating that the same haplotypes had similar estimated effects in the 2 breeds. Although the data were limited (<350 cows in each study), this latter result may suggest that genes in the casein complex itself are responsible for the effects observed, rather than loci that are physically linked on either side of the casein cluster.