Associations of a polymorphic AP-2 binding site in the 5'-flanking region of the bovine beta-lactoglobulin gene with milk proteins

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.

Genetic variation in promoter region of the bovine LAP3 gene associated with estimated breeding values of milk production traits and clinical mastitis in dairy cattle

The purpose of this study was to identify genetic variants in the promoter and 5'UTR regions of bovine leucine amino peptidase three (LAP3) gene and analysed their associations with estimated breeding values (EBVs) of milk production traits and clinical mastitis in Sahiwal and Karan Fries cattle. Eleven SNPs were identified within the region under study of the LAP3 gene, including seven promoter variants (rs717156555: C>G, rs720373055: T>C, rs715189731: A>G, rs516876447: A>G, rs461857269: C>T, rs136548163: C>T, and rs720349928: G>A) and four 5'UTR variants (rs717884982: C>T, rs722359733: C>T, rs481631804: C>T and rs462932574: T>G). Out of them, 10 SNPs variants were found in both Sahiwal and Karan Fries cattle, with one SNP variant (rs481631804: C>T) being unique to Karan Fries cattle. Seven of these identified SNPs were chosen for association analyses. Individual SNP based association analysis revealed that two SNPs (rs720373055: T>C and rs720349928: G>A) were significantly associated with EBVs of lactation milk yield (LMY), 305-day milk yield (305dMY), and one significant association of SNP rs722359733: C>T with lactation length (LL) was observed. Haplotype based association analysis indicated that diplotypes are significantly associated with EBVs of LMY, 305dMY, and LL, individuals with H1H3 (CTACGCT/GCGTACG) being linked to higher lactation performance than other diplotypes. Further logistic regression analysis revealed that, animals with diplotype H1H3 was less susceptible to the incidence of clinical mastitis than other cows, as the odds ratio for the non-incidence of clinical mastitis was found to be low. Altogether, variations in the LAP3 gene promoter could be used as a genetic marker, most notably diplotype H1H3, may greatly benefit the simultaneous improvement of mastitis resistance and milk yield traits in dairy cattle. Moreover, bioinformatics analysis predicted that the SNPs rs720373055: T>C, rs715189731:A>G and rs720349928: G>A is located in the core promoter region and in TFBs, play key role in regulation of studied phenotypes.

Novel functional mutation of the PDIA3 gene affects milk composition traits in Chinese Holstein cattle

Protein disulfide isomerase family A member 3 (PDIA3) is a multifunctional protein, and it plays a vital role in modulating various cell biological functions under physiological and pathological conditions. Our previous study on Mediterranean buffalo demonstrated that PDIA3 is a potential candidate gene associated with milk yield based on genome-wide association study analysis. However, the genetic effects of the PDIA3 gene on milk performance in dairy cattle and the corresponding mechanism have not been documented. This study aims to explore the genetic effects of PDIA3 polymorphisms on milk production traits in 362 Chinese Holstein cattle. The results showed that 4 SNPs were identified from the 5' untranslated region of the PDIA3 gene in the studied population, of which 2 SNPs (g.-1713 C>T and g.-934 G>A) were confirmed to be significantly associated with milk protein percentage, whereas g.-434 C>T was significantly associated with milk fat percentage. Notably, linkage disequilibrium analysis indicated that 3 SNPs (g.-1713 C>T, g.-934 G>A, and g.-695 A>C) formed one haplotype block, which was found to be significantly associated with milk protein percentage. The luciferase assay demonstrated that allele C of g.-434 C>T exhibited a higher promotor activity compared with allele T, suggesting that g.-434 C>T might be a potential functional mutation affecting PDIA3 expression. Furthermore, overexpression of the PDIA3 gene was found to induce higher levels of triglyceride and BODIPY fluorescence intensity. In addition, PDIA3 overexpression was also found to positively regulate the synthesis and secretion of α-casein, β-casein, and κ-casein, whereas knockdown of this gene showed the opposite effects. In summary, our findings revealed significant genetic effects of PDIA3 on milk composition traits, and the identified SNP and the haplotype block might be used as genetic markers for dairy cow selected breeding.

Screening for phenotypic outliers identifies an unusually low concentration of a β-lactoglobulin B protein isoform in bovine milk caused by a synonymous SNP

Background

Milk samples from 10,641 dairy cattle were screened by a mass spectrometry method for extreme concentrations of the A or B isoforms of the whey protein, β-lactoglobulin (BLG), to identify causative genetic variation driving changes in BLG concentration.

Results

A cohort of cows, from a single sire family, was identified that produced milk containing a low concentration of the BLG B protein isoform. A genome-wide association study (GWAS) of BLG B protein isoform concentration in milk from AB heterozygous cows, detected a group of highly significant single nucleotide polymorphisms (SNPs) within or close to the BLG gene. Among these was a synonymous G/A variation at position + 78 bp in exon 1 of the BLG gene (chr11:103256256G > A). The effect of the A allele of this SNP (which we named B’) on BLG expression was evaluated in a luciferase reporter assay in transfected CHO-K1 and MCF-7 cells. In both cell types, the presence of the B’ allele in a plasmid containing the bovine BLG gene from -922 to + 898 bp (relative to the transcription initiation site) resulted in a 60% relative reduction in mRNA expression, compared to the plasmid containing the wild-type B sequence allele. Examination of a mammary RNAseq dataset (n = 391) identified 14 heterozygous carriers of the B’ allele which were homozygous for the BLG B protein isoform (BB’). The level of expression of the BLG B’ allele was 41.9 ± 1.0% of that of the wild-type BLG B allele. Milk samples from three cows, homozygous for the A allele at chr11:103,256,256 (B’B’), were analysed (HPLC) and showed BLG concentrations of 1.04, 1.26 and 1.83 g/L relative to a mean of 4.84 g/L in milk from 16 herd contemporaries of mixed (A and B) BLG genotypes. The mechanism by which B’ downregulates milk BLG concentration remains to be determined.

Conclusions

High-throughput screening and identification of outliers, enabled the discovery of a synonymous G > A mutation in exon 1 of the B allele of the BLG gene (B’), which reduced the milk concentration of β-lactoglobulin B protein isoform, by more than 50%. Milk from cows carrying the B’ allele is expected to have improved processing characteristics, particularly for cheese-making.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00711-z.

agReg-SNPdb: A Database of Regulatory SNPs for Agricultural Animal Species

Transcription factors (TFs) govern transcriptional gene regulation by specifically binding to short DNA motifs, known as transcription factor binding sites (TFBSs), in regulatory regions, such as promoters. Today, it is well known that single nucleotide polymorphisms (SNPs) in TFBSs can dramatically affect the level of gene expression, since they can cause a change in the binding affinity of TFs. Such SNPs, referred to as regulatory SNPs (rSNPs), have gained attention in the life sciences due to their causality for specific traits or diseases. In this study, we present agReg-SNPdb, a database comprising rSNP data of seven agricultural and domestic animal species: cattle, pig, chicken, sheep, horse, goat, and dog. To identify the rSNPs, we constructed a bioinformatics pipeline and identified a total of 10,623,512 rSNPs, which are located within TFBSs and affect the binding affinity of putative TFs. Altogether, we implemented the first systematic analysis of SNPs in promoter regions and their impact on the binding affinity of TFs for livestock and made it usable via a web interface.

Complete CSN1S2 Characterization, Novel Allele Identification and Association With Milk Fatty Acid Composition in River Buffalo

The αs2-casein is one of the phosphoproteins secreted in all ruminants' milk, and it is the most hydrophilic of all caseins. However, this important gene (CSN1S2) has not been characterized in detail in buffaloes with only two alleles detected (reported as alleles A and B), and no association studies with milk traits have been carried out unlike what has been achieved for other species of ruminants. In this study, we sequenced the whole gene of two Mediterranean river buffalo homozygotes for the presence/absence of the nucleotide C (g.7539G>C) realized at the donor splice site of exon 7 and, therefore, responsible for the skipping of the same exon at mRNA level (allele B). A high genetic variability was found all over the two sequenced CSN1S2 alleles. In particular, 74 polymorphic sites were found in introns, six in the promoter, and three SNPs in the coding region (g.11072C>T, g.12803A>T, and g.14067A>G) with two of them responsible for amino acid replacements. Considering this genetic diversity, those found in the database and the SNP at the donor splice site of exon 7, it is possible to deduce at least eight different alleles (CSN1S2 A, B, B1, B2, C, D, E, and F) responsible for seven different possible translations of the buffalo αs2-casein. Haplotype data analysis suggests an evolutionary pathway of buffalo CSN1S2 gene consistent with our proposal that the published allele CSN1S2 A is the ancestral αs2-CN form, and the B2 probably arises from interallelic recombination (single crossing) between the alleles D and B (or B1). The allele CSN1S2 C is of new identification, while CSN1S2 B, B1, and B2 are deleted alleles because all are characterized by the mutation g.7539G>C. Two SNPs (g.7539G>C and g.14067A>G) were genotyped in 747 Italian buffaloes, and major alleles had a relative frequency of 0.83 and 0.51, respectively. An association study between these SNPs and milk traits including fatty acid composition was carried out. The SNP g.14067A>G showed a significant association (P < 0.05) on the content of palmitic acid in buffalo milk, thus suggesting its use in marker-assisted selection programs aiming for the improvement of buffalo milk fatty acid composition.

Single nucleotide polymorphisms in candidate genes associated with milk yield in Argentinean Holstein and Holstein x Jersey cows

BACKGROUND

Research on loci influencing milk production traits of dairy cattle is one of the main topics of investigation in livestock. Many genomic regions and polymorphisms associated with dairy production have been reported worldwide. In this context, the purpose of this study was to identify candidate loci associated with milk yield in Argentinean dairy cattle. A database of candidate genes and single nucleotide polymorphisms (SNPs) for milk production and composition was developed. Thirty-nine SNPs belonging to 22 candidate genes were genotyped on 1643 animals (Holstein and Holstein x Jersey). The genotypes obtained were subjected to association studies considering the whole population and discriminating the population by Holstein breed percentage. Phenotypic data consisted of milk production values recorded during the first lactation of 1156 Holstein and 462 Holstein x Jersey cows from 18 dairy farms located in the central dairy area of Argentina. From these records, 305-day cumulative milk production values were predicted.

RESULTS

Eight SNPs (rs43375517, rs29004488, rs132812135, rs137651874, rs109191047, rs135164815, rs43706485, and rs41255693), located on six Bos taurus autosomes (BTA4, BTA6, BTA19, BTA20, BTA22, and BTA26), showed suggestive associations with 305-day cumulative milk production (under Benjamini-Hochberg procedure with a false discovery rate of 0.1). Two of those SNPs (rs43375517 and rs135164815) were significantly associated with milk production (Bonferroni adjusted p-values < 0.05) when considering the Holstein population.

CONCLUSIONS

The results obtained are consistent with previously reported associations in other Holstein populations. Furthermore, the SNPs found to influence bovine milk production in this study may be used as possible candidate SNPs for marker-assisted selection programs in Argentinean dairy cattle.

Mammary cell-activating factor regulates the hormone-independent transcription of the early lactation protein (ELP) gene in a marsupial

The regulation of the tammar wallaby (Macropus eugenii) early lactation protein (ELP) gene is complex. ELP is responsive to the lactogenic hormones; insulin (I), hydrocortisone (HC) and prolactin (PRL) in mammary gland explants but could not be induced with lactogenic hormones in tammar primary mammary gland cells, nor in KIM-2 conditionally immortalised murine mammary epithelial cells. Similarly, ELP promoter constructs transiently-transfected into human embryonic kidney (HEK293T) cells constitutively expressing the prolactin receptor (PRLR) and Signal Transducer and Activator of Transcription (STAT)5A were unresponsive to prolactin, unlike the rat and mouse β-casein (CSN2) promoter constructs. Identification of the minimal promoter required for the hormone-independent transcription of tammar ELP in HEK293Ts and comparative analysis of the proximal promoters of marsupial ELP and the orthologous eutherian colostrum trypsin inhibitor (CTI) gene suggests that mammary cell-activating factor (MAF), an E26 transformation-specific (ETS) factor, may bind to an AGGAAG motif and activate tammar ELP.

Variability of the caprine whey protein genes and their association with milk yield, composition and renneting properties in the Sarda breed: 2. The BLG gene

The variability of the promoter region and the 3'UTR (exon-7) of the BLG gene, encoding the β-lactoglobulin, was investigated by sequencing in 263 lactating Sarda goats in order to assess its association with milk traits. Milk traits included: milk yield, fat, total protein and lactose content, pH, daily fat and protein yield (DFPY), freezing point, milk energy, somatic cell count, total microbial mesophilic count, rennet coagulation time (RCT), curd firming rate (k20) and curd firmness (a30). A total of 7 polymorphic sites were detected and the sequence analysed was given accession number KM817769. Only three SNPs (c.-381C>T, c.-323C>T and c.*420C>A) had minor allele frequency higher than 0.05. The effects of farm, stage of lactation and the interaction farm × stage of lactation significantly influenced all the milk traits (P T and c.*420C>A (P T (P < 0.001). The c.-381TT homozygous goats showed lower pH, RCT and k20 than c.-381CT (P < 0.05). In conclusion the polymorphism of the goat BLG gene did not affect the total protein content of the Sarda goat milk, and only weakly influenced RCT and k20. On the other hand, an interesting effect on milk yields and DFPY emerged in two SNPs. This information might be useful in dairy goat breeding programs.

Genomic architecture of bovine κ-casein and β-lactoglobulin

The objective of this study was to characterize the genetic architecture underlying the absolute concentrations of 2 important milk proteins, κ-casein (κ-CN) and β-lactoglobulin (β-LG), in a backcross population of (Holstein × Jersey) × Holstein cattle. A genome-wide association analysis was performed using a selective DNA pooling strategy and the Illumina BovineHD BeadChip assay [777,000 (777K) SNP markers; Illumina Inc., San Diego, CA]. After correction for multiple testing, 25 single nucleotide polymorphisms were found to be associated with κ-CN and 36 single nucleotide polymorphisms were associated with β-LG. A pathway association analysis revealed 15 Gene Ontology (GO) terms associated with the κ-CN trait and 28 GO terms associated with β-LG. In addition, several GO terms were associated with both milk proteins. Further analysis revealed that κ-CN and β-LG production is regulated by both kinase and phosphatase activity, including mechanisms regulating the extracellular matrix. These results are in concordance with the complex multihormonal process controlling the expression of milk proteins and interactions between mammary epithelial cells and extracellular matrix components. Although κ-CN and β-LG milk proteins are expressed by single genes, the results from this study showed that many loci are involved in the regulation of the concentration of these 2 proteins.

Single nucleotide polymorphisms in the leptin-a gene and associations with growth traits in the orange-spotted grouper (Epinephelus coioides)

Leptin is a multifunctional protein involved in processes such as body weight regulation, energy expenditure, fat metabolism, food intake, and appetite regulation. Duplicate leptin genes, leptin-a and leptin-b, were previously detected in the orange-spotted grouper. In this study, we cloned the full-length open reading frame (ORF) of the leptin-a gene in the orange-spotted grouper, searched for polymorphisms, and performed association analyses between these polymorphisms and seven growth traits. Six polymorphisms, consisting of 2 SNPs in intron 1 (c.182T > G, c.183G > T) and 4 SNPs in exon 2 (c.339C > G, c.345C > T, c.447G > A, c.531C > T), were identified and genotyped in 200 individuals. The c.182T > G and c.183G > T polymorphisms showed complete linkage and were analyzed together. Association analyses revealed that the c.182 + 183TG > GT polymorphism was significantly associated with body weight (BWT) and body width (BWH), with the AB (TG/GT) genotype showing positive effects on growth traits. Additionally, the SNP c.447G > A was significantly associated with BWT, BWH, overall length (OL), trunk width (TW), and head length (HL), with the GA genotype displaying positive effects on growth traits. The c.531C > T SNP showed a close association between the TT genotype and decreased growth. Our results demonstrate that several polymorphisms in the leptin-a gene are associated with growth traits and can be used for marker-assisted selection (MAS) in orange-spotted grouper populations.

Functional characterization of genetic polymorphisms identified in the promoter region of the bovine PEPS gene

Peptidase S (PEPS) is a metallopeptidase that cleaves N-terminal residues from proteins and peptides. PEPS is used as a cell maintenance enzyme with critical roles in peptide turnover. The promoter region located upstream of the initiation site plays an important role in regulating gene expression. Polymorphism in the promoter region can alter gene expression and lead to biological changes. In the current study, polymorphisms in the promoter region of the PEPS gene were investigated. Polymerase chain reaction (PCR)-restriction fragment length polymorphism and DNA sequencing methods were used to screen sequence variations in the promoter region of DNA samples from 743 Chinese Holstein cattle. Two polymorphisms (g. -534 T>C and g. -2545 G>A) were identified and eight haplotypes were classified by haplotype analysis. The two genetic polymorphisms and haplotypes were associated with fat percentage and somatic cell score in Chinese Holstein cattle. The results of real-time PCR showed that cow kidneys exhibit the highest PEPS expression level. Moreover, bioinformatics analysis predicted that the single-nucleotide polymorphism g. -534 T>C is located in the core promoter region and in the transcription factor binding sites. The promoter activities of the polymorphism of -543 T>C were measured by luciferase assay in the human kidney epithelial cell line 293T. Transcriptional activity is significantly lower in cell lines transfected with the reporter construct containing 2.5 kb upstream fragments with -543 C than in those with wild-type -543 T. The results indicated that genetic variation at locus -543 influences PEPS promoter activity. The genetic variation in the promoter region of PEPS gene may regulate PEPS gene transcription and might have consequences at a regulatory level.

Polymorphisms of beta-lactoglobulin promoter region in three Sicilian goat breeds

Several beta-lactoglobulin (BLG) polymorphisms have been described within the proximal promoter region and coding region of the caprine gene, although no genetic variants affecting the protein amino acid composition and/or expression level have been characterized so far. Binding sites for several transcription factors (TFs) are present in the BLG promoter region. The aims of this work were to sequence the full-length promoter region of three Sicilian goat breeds in order to identify polymorphisms, analyze the identified haplotypes, search for differences between breeds for the presence of polymorphisms in this gene region, search for putative TFs binding sites, and check if polymorphisms lay within the identified TFs binding sites. The promoter region of BLG gene in Sicilian goat breeds showed high level of polymorphism due to the presence of 36 single nucleotide polymorphisms (SNPs). Association between polymorphic sites was computed within the whole sample analyzed and 18 haplotypes were inferred. Binding sites for three milk protein binding factors (MPBFs) and four nuclear factor-I (NF-I) were found within BLG promoter region based on the ovine sequence. The identification of some SNPs within TFs binding sites allowed hypothesizing the loss of TFs. Further studies are in progress to evaluate the effect of these mutations on binding affinity of TFs, the functional interaction of the TFs with the goat BLG promoter, and the relationship of the polymorphisms with BLG gene expression and milk production and composition.

Genome wide association studies for milk production traits in Chinese Holstein population

Genome-wide association studies (GWAS) based on high throughput SNP genotyping technologies open a broad avenue for exploring genes associated with milk production traits in dairy cattle. Motivated by pinpointing novel quantitative trait nucleotide (QTN) across Bos Taurus genome, the present study is to perform GWAS to identify genes affecting milk production traits using current state-of-the-art SNP genotyping technology, i.e., the Illumina BovineSNP50 BeadChip. In the analyses, the five most commonly evaluated milk production traits are involved, including milk yield (MY), milk fat yield (FY), milk protein yield (PY), milk fat percentage (FP) and milk protein percentage (PP). Estimated breeding values (EBVs) of 2,093 daughters from 14 paternal half-sib families are considered as phenotypes within the framework of a daughter design. Association tests between each trait and the 54K SNPs are achieved via two different analysis approaches, a paternal transmission disequilibrium test (TDT)-based approach (L1-TDT) and a mixed model based regression analysis (MMRA). In total, 105 SNPs were detected to be significantly associated genome-wise with one or multiple milk production traits. Of the 105 SNPs, 38 were commonly detected by both methods, while four and 63 were solely detected by L1-TDT and MMRA, respectively. The majority (86 out of 105) of the significant SNPs is located within the reported QTL regions and some are within or close to the reported candidate genes. In particular, two SNPs, ARS-BFGL-NGS-4939 and BFGL-NGS-118998, are located close to the DGAT1 gene (160bp apart) and within the GHR gene, respectively. Our findings herein not only provide confirmatory evidences for previously findings, but also explore a suite of novel SNPs associated with milk production traits, and thus form a solid basis for eventually unraveling the causal mutations for milk production traits in dairy cattle.

Mapping a quantitative trait locus for the concentration of beta-lactoglobulin in milk, and the effect of beta-lactoglobulin genetic variants on the composition of milk from Holstein-Friesian x Jersey crossbred cows

AIM

To identify quantitative trait loci (QTL) affecting the concentration of beta-lactoglobulin in milk, and to evaluate the effect of beta-lactoglobulin genetic variants on the concentration of fat, protein and casein in bovine milk.

METHODS

A herd of 850 F2 Holstein-Friesian x Jersey crossbred cows was produced through mating six Holstein-Friesian x Jersey F1 bulls of high genetic merit with F1 cows from the national herd. A total of 1,610 herd-test records from 556 second-parity crossbreds were analysed. The concentration of fat, protein and casein in milk was measured at peak, mid- and late lactation, during the production seasons of 2003-2004 and 2004-2005. Liveweight was measured daily. DNA from the F2 animals, their F1 dams and sires, and selected grandsires was genotyped across the genome, initially with 285 microsatellite markers, and subsequently with 6,634 single nucleotide polymorphisms (SNP).

RESULTS

A highly significant QTL for the concentration of beta-lactoglobulin in milk was identified, which coincided with the position of the beta-lactoglobulin gene on bovine Chromosome 11. No other consistently significant QTL for the concentration of beta-lactoglobulin in milk were detected. Cows with the BB beta-lactoglobulin genotype produced milk with a 30% lower concentration of beta-lactoglobulin than cows with the AA genotype. The beta-lactoglobulin polymorphism also explained variation in the proportion of casein in total protein. In addition, the percentage of fat was higher for BB than AA animals, whereas the percentage of total protein, mean daily milk yield and liveweight did not differ between AA and BB animals.

CONCLUSIONS

A significant QTL determining the concentration of beta-lactoglobulin in milk was identified. Selection of animals for the beta-lactoglobulin B-allele may enable the production of milk naturally enriched for casein, thus allowing a potential increase in the yield of cheese. There may be additional future value in production of bovine milk more like human milk, where decreasing the concentration of beta-lactoglobulin is desirable.

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.

Molecular cloning, promoter analysis and SNP identification of Italian Nicastrese and Saanen lactoferrin gene

Lactoferrin (Lf) is an iron-binding glycoprotein found in exocrine secretions including milk. High levels of lactoferrin may have a role in the prevention of microbial infection of the mammary gland. In this report we sequenced and characterized goat lactoferrin cDNA and its promoter region in two different breeds of goat. The complete cDNA comprised 2356 nucleotides, including 38 bp at the 5'-UTR and 194 bp at the 3'-UTR. The open reading frame is 2127 bp long and it encodes a mature protein of 689 aminoacids. A total of 19 nucleotide differences, 11 of them being responsible for 8 aminoacid changes, were identified through the comparison with French, Korean and Tibetan goat lactoferrin cDNAs. About 1700 bp of the lactoferrin gene promoter were sequenced. Sequence analysis revealed a non-canonical TATA box, multiple SP1/GC elements, and other putative binding sites for transcription factors, such as NF-kappaB, STAT3 and AP2. Two SNPs were identified, one of which would seem to create a new putative AP2 consensus sequence. The presence of an additional AP2 binding site could be associated with quantitative differences of such protein fraction, which could enhance all the activities related to such protein, and improve mammary gland defence against bacterial infections.

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.

Aberrant low expression level of bovine beta-lactoglobulin is associated with a C to A transversion in the BLG promoter region

Beta-lactoglobulin (beta-LG) is the major whey protein in cow's milk. It is well established that the predominant 2 genetic variants, beta-LG A and B, are differentially expressed. Extensive investigation of the genetic variation in the promoter region of the BLG gene revealed the existence of specific haplotypes associated with the A and B variants, respectively. However, the genetic basis for the differential expression of BLG A and B alleles is still elusive. We have previously reported a quantitative beta-LG B variant, characterized by a very low beta-LG protein expression level. Here, we report that the corresponding BLG allele (BLG B*) shows a correspondingly low mRNA expression level. Comparative DNA sequencing of 7,670 bp of the BLG B* allele and the established BLG B allele revealed a unique difference of a C to A transversion at position 215 bp upstream of the translation initiation site (g.-215C>A). This mutation segregated perfectly with the differential phenotypic expression in a paternal half-sib family and could be confirmed in 2 independent cases. The sequence of the BLG B allele in the region of the mutation is highly conserved among 4 related ruminant species. The site of the mutation corresponds to a putative consensus-binding sequence for the transcription factors c-Rel and Elk-1 as predicted by searching the TRANSFAC database. The beta-LG B* site might be relevant in the natural production of milk of low beta-LG content.

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.

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.