Medicine » Obstetrics and Gynecology » "Prolactin", book edited by György M. Nagy and Bela E. Toth, ISBN 978-953-51-0943-3, Published: January 23, 2013 under CC BY 3.0 license. © The Author(s).

Chapter 3

Use of the Bovine Prolactin Gene (bPRL) for Estimating Genetic Variation and Milk Production in Aboriginal Russian Breeds of Bos taurus L.

By I.V. Lazebnaya, O.E. Lazebny, S.R. Khatami and G.E. Sulimova
DOI: 10.5772/54756

Article top

Overview

a) Yakut, (b) Bestuzhev, (c) Kostroma, and (d) Yaroslavl cattle breeds [
                        12, 
                        13, 
                        14, 
                        15].
Figure 1. a) Yakut, (b) Bestuzhev, (c) Kostroma, and (d) Yaroslavl cattle breeds [ 12, 13, 14, 15].
Schematic representation of the regulatory cascade involving the prolactin gene.
Figure 2. Schematic representation of the regulatory cascade involving the prolactin gene.
Distribution of 
                            bPRL(
                            RsaI) allele frequencies in the four breeds studied and in other Russian and foreign breeds. The abscissa shows the breeds; the ordinate shows the allele frequencies (in fractions of unity).
Figure 3. Distribution of bPRL( RsaI) allele frequencies in the four breeds studied and in other Russian and foreign breeds. The abscissa shows the breeds; the ordinate shows the allele frequencies (in fractions of unity).
Distribution of 
                            bPRL genotype frequencies in the four breeds studied and in other Russian and foreign breeds. The radial axes show the breeds and the corresponding genotype frequencies.
Figure 4. Distribution of bPRL genotype frequencies in the four breeds studied and in other Russian and foreign breeds. The radial axes show the breeds and the corresponding genotype frequencies.
The observed and expected heterozygosities for the 
                            bPRL gene markers in Russian and foreign cattle breeds. The abscissa shows the breeds; the ordinate shows the heterozygosity coefficients (in fractions of unity).
Figure 5. The observed and expected heterozygosities for the bPRL gene markers in Russian and foreign cattle breeds. The abscissa shows the breeds; the ordinate shows the heterozygosity coefficients (in fractions of unity).
Dependence of (a) the milk yield (in kilograms), (b) the fat content, and (c) the protein content of milk of the third lactation on the 
                            bPRL genotypes in Yaroslavl cattle. The abscissas show the 
                            bPRL genotypes; the ordinates show the mean values of the traits; vertical bars are 0.95 confidence intervals.
Figure 6. Dependence of (a) the milk yield (in kilograms), (b) the fat content, and (c) the protein content of milk of the third lactation on the bPRL genotypes in Yaroslavl cattle. The abscissas show the bPRL genotypes; the ordinates show the mean values of the traits; vertical bars are 0.95 confidence intervals.
Dependence of (a) the milk yield (in kilograms), (b) the fat content, and (c) the protein content of milk of the first lactation on combined genotypes of the 
                            bPRL and 
                            bGH genes in Yaroslavl cattle. The abscissas show the 
                            bPRL genotypes (
                            AA, AB, and 
                            BB); the ordinates show the mean values of the traits; vertical bars are 0.95 confidence intervals. Different 
                            bGH genotypes are denoted by dots of different colors; the points corresponding to the same 
                            bGH genotype but different 
                            bPRL genotypes are connected by dotted lines.
Figure 7. Dependence of (a) the milk yield (in kilograms), (b) the fat content, and (c) the protein content of milk of the first lactation on combined genotypes of the bPRL and bGH genes in Yaroslavl cattle. The abscissas show the bPRL genotypes ( AA, AB, and BB); the ordinates show the mean values of the traits; vertical bars are 0.95 confidence intervals. Different bGH genotypes are denoted by dots of different colors; the points corresponding to the same bGH genotype but different bPRL genotypes are connected by dotted lines.

Use of the Bovine Prolactin Gene (bPRL) for Estimating Genetic Variation and Milk Production in Aboriginal Russian Breeds of Bos taurus L.

I.V. Lazebnaya1, O.E. Lazebny2, S.R. Khatami1, 3 and G.E. Sulimova1

1. Introduction

Prolactin is a protein hormone mainly, but not exclusively produced by lactotroph cells of the anterior pituitary. Its role in lactogenesis and galactopoiesis (maintenance of milk secretion) is well demonstrated [ 1, 2]. Therefore, the gene encoding it ( PRL) is considered to be one of the key links in the gene network constituting the hereditary component of milk productivity. Test systems for cattle breeding have been developed based on the associations of the PRL gene polymorphism with milk yield and quality.

Inbreeding, which decreases the genetic variation and viability of animals, is a well-known negative consequence of artificial selection. Its impact is further aggravated by the recent trend towards globalization of some cattle breeds [ 3]. Therefore, conservation of aboriginal breeds adapted to local conditions (which are not infrequently extreme) is necessary in countries with wide zonal climatic variations.

This is especially important when a breed in question has pronounced adaptive characteristics and its population is small. Yakut cattle represent one of such breeds ( Figure 1a); it is unique among Russian breeds in terms of ecological plasticity. These cattle live in the northernmost part of the Bos taurus species range, a hardly accessible region of the subarctic zone of the Republic of Sakha (Yakutia), Russia, surrounded with mountain ridges. The morphological and physiological characteristics of Yakut cattle and their biochemical and behavioral adaptations allow free grazing almost round the year despite a severe continental climate, with the mean air temperatures usually varying from –43°C in winter to +25°C in summer (the lowest and highest temperatures on record are –65°C and +38°C, respectively). These animals can live on rough foods. Their body, including the udder, is covered with long, thick hair protecting them from cold and gnats. The color of aboriginal Yakut cattle varies from black and red to a leopard pattern with white spots on the head and lower trunk. This breed is exceptionally resistant to tuberculosis, leukemia, and brucellosis. Yakut cattle are small, with the shoulder height shorter than 1 m and the live weights of bulls and cows of 500–550 and 350–400 kg, respectively. The milk yield is low (2100–2350 kg in the breeding stock), but the fat content of milk is as high as 7.3% [ 4]. Yakut cattle have long been providing local residents (mostly Yakuts) with beef and dairy products. Cattle leather is widely used in Yakut ethnic handicrafts; Yakuts traditionally make comfortable, durable, beautiful leather clothes. Excavations in the Olekminsk district of Yakutia have revealed remnants of nomad camps containing fossil bones of domestic cattle, which suggests an ancient origin of this breed [ 5].

Bestuzhev and Kostroma cattle are dual-purpose breeds. Bestuzhev cattle selected for both beef and milk production were bred from local cattle in Samara province, Russia, in the late 18th century ( Figure 1b). English Shorthorn cattle were used for its improvement, the offspring being crossed with the Holland, Shorthorn, Simmental, and some other breeds. The breed was completely formed by the mid-19th century. Bestuzhev cattle are well adapted to the continental climate of the Volga basin. The animals are red; the color intensity varies from light-red to dark-red or cherry-red. Some animals have white spots, mostly on the lower trunk, udder, and head. The mean milk yield of Bestuzhev cattle is 5502–8250 kg; the mean fat content of milk is 3.82–4.0% (the maximum content is 5.5%) [ 4]. Bestuzhev cattle are especially valuable because they are almost free of hereditary diseases and abnormalities and are resistant to tuberculosis and leukemia.

Kostroma cattle are classified with the group of brown cattle ( Figure 1c). The breed was registered in 1944. These cattle are characterized by a high growing capacity, strong constitution, and steady inheritance of commercially valuable traits, including a good milk quality. The Kostroma breed is regarded as one of the most productive dual-purpose breeds. The live weights of Kostroma bulls and cows are 800–900 and 550–650 kg, respectively; the fat and protein contents of milk are as high as 3.9 and 3.6%, respectively. The milk yield varies from 6000–8000 to 10,500 kg [ 4]. The Kostroma breed is characterized by a high total frequency of the BoLA–DRB3 gene alleles determining the leukemia resistance (on average, 35.9%) [ 6].

Yaroslavl cattle, formed as a native cattle breed in the 16th century, have the highest milk yield and the best milk composition among all native Russian breeds ( Figure 1d). Their milk contains, on average, 4.37% of fat (maximum content, 5.0%) and 3.4–3.6% of protein; the dry matter content is 13.6% (compared to 12.3–12.5% in other breeds) [ 4]. Yaroslavl cattle are usually black, except the white head with a black mask around the eyes and the white lower trunk. Yaroslavl cattle were first mentioned in the literature in the mid-19th century. Bestuzhev cattle were named after the original breeder; Yaroslavl, Yakut, and Kostroma cattle, after the region of origin.

The balance between the increase in the milk yield and quality of cattle breeds and preservation of diversity both within and between breeds of B. taurus is a complicated problem. Therefore, analysis of the variation of the gene markers that are affected by artificial selection because of their associations with milk yield and composition would be useful for breed monitoring.

The bovine prolactin gene ( bPRL) is traditionally regarded as a good candidate gene for marker-assisted selection (MAS) [ 7] for milk production parameters, because it has been located to chromosome 23 at 43 cM, close to the quantitative trait loci (QTLs) (36, 41, and 42 cM) [ 8- 11]. In addition, it is known that the binding of the bPRL gene product with its receptor (PRLR) initiates a signaling cascade that activates the transcription of a number of genes, including the genes of milk proteins (caseins and lactalbumin).

media/image1.jpeg

Figure 1.

a) Yakut, (b) Bestuzhev, (c) Kostroma, and (d) Yaroslavl cattle breeds [ 12, 13, 14, 15].

Note that this regulatory cascade involves growth hormone ( Figure 2), because it is recognized not only by its own receptor (GHR), but also by PRLR [ 16]. The prolactin and growth hormone genes are very similar to each other, because they have resulted from duplication of a common ancestral gene. They have the same general structure (five exons and four introns) [ 17] and a common positive transcription factor (PITI) [ 18].

The specifics of the bPRL and bGH genes suggest their combined effect on milk production; however, this has been paid little attention until now. Most studies on the polymorphism of milk production genes deal with isolated effects of individual genes. Many data on the polymorphisms of both bPRL and bGH genes and their relationship with the milk yield and composition have been accumulated. For example, a synonymous A–G transition has been found in the codon of amino acid residue 103 of bPRL in the third exon of its gene; this mutation results in an RsaI polymorphic site [ 19]. The AA and AB genotypes have been shown to be associated with the milk yield and protein content in Polish Black & White, Holstein Friesian, and Brown Swiss cattle [ 19- 21]. Other associations have been found in Russian Red Pied cattle [ 22]. The C–G transversion in the third exon of the bGH gene (nucleotide position 2141) is known to be related to milk production traits. This SNP entails the disappearance of an AluI restriction site and the substitution of valine for leucine (Leu → Val) at position 127 of the bGH amino acid sequence. There is evidence for a stronger dependence of the milk yield on the LL genotype than the LV genotype in Black Pied and Holstein Friesian cattle [ 23, 24]. However, Mitra et al. [ 25] and Dybus et al. [ 20] note a positive dependence of the milk yield and both fat and protein contents on the V allele in Holstein and Polish Black & White cattle. Study of the combined effect of the bPRL and bGH genes might explain some contradictions about the associations of individual markers of these genes. Although associations of the bPRL( RsaI) and bGH( AluI) polymorphisms with milk yield and quality have been extensively studied, the combined effect of SNPs of these genes has been hardly considered at all.

media/image2.jpeg

Top: an anterior pituitary cell. Chromosomes carrying the transcription factor ( Pit-1), growth hormone ( GH), and prolactin ( PRL) genes are shown in the nucleus (a green oval). The product of Pit-1 (TF Pit-1) activates the transcription of the prolactin and growth hormone genes. The corresponding hormones (GH and PRL, shown in dark blue and light blue, respectively) enter the bloodstream and reach mammary cells, where they are recognized by transmembrane receptors (GHR and PRLR, respectively). This initiates regulatory cascades activating the transcription of the milk protein (β-casein and α-lactalbumin) genes. The prolactin and growth hormone receptors are shown in yellow and green, respectively. The growth hormone competes with prolactin for the prolactin receptor.

Figure 2.

Schematic representation of the regulatory cascade involving the prolactin gene.

The goal of this study was twofold. First, we compared the Kostroma, Bestuzhev, Yakut, and Yaroslavl breeds both with one another and with some other Russian and foreign breeds (literature data) with respect to genetic variation as estimated from the RsaI polymorphism of the prolactin gene. Second, we estimated the effect of the RsaI polymorphism of the bPRL gene, as well as the effect of its combination with the AluI polymorphism of the bGH gene, on milk production in Yaroslavl cattle.

2. Materials and methods

We used PCR–RFLP analysis to study the RsaI polymorphism of the bPRL gene in Yakut cattle ( n = 41) in the Republic of Sakha (Yakutia), Bestuzhev cattle ( n = 57) in Samara region, Kostroma cattle ( n =124) in Kostroma region, and Yaroslavl cattle ( n = 113) in Yaroslavl region of Russia. The possible effects of the bPRL( RsaI) polymorphism and its combination with the bGH( AluI) polymorphism on milk production in Yaroslavl cattle were estimated.

DNA was isolated from 200-μl samples of whole blood using a Diatom TM DNA Prep kit (IsoGeneLab., Russia). Fragments of bPRL (156 bp) and bGH (223 bp) [ 25] were amplified in a Tertsik thermal cycler by the standard methods using a GenePak TMPCR Core kit (IsoGene Lab., Russia). The DNA digestion with the RsaI and AluI restriction endonucleases was performed as recommended by the manufacturer (MBI Fermentas, Lithuania). The oligonucleotides F PRL (5'-CGAGTCCTTATGAGCTTGATTCTT-3') and R PRL (5'-GCCTTCCAGAAGTCGTTTGTTTTC-3') served as primers for the bPRL gene fragments; F GH (5'-GCTGCTCCTGAGGGCCCTTCG-3') and R GH (5'-GCGGCGGCACTTCATGACCCT-3'), for the bGH gene fragments. The following amplification profiles were used: for the bPRL gene fragments, one cycle of 94°C for 4 min; 35 cycles of 94°C for 30 s, 58°C for 30 s, and 72°C for 30 s; and one cycle of 72°C for 10 min; for the bGH gene fragments, one cycle of 94°C for 4 min; 35 cycles of 94°C for 45 s, 67°C for 45 s, and 72°C for 45 s; and one cycle of 72°C for 10 min. The amplification products were detected by means of 2% agarose gel electrophoresis (0.5 μg/ml ethidium bromide in 1X TBE buffer). The restriction products were visualized in 2% agarose gel (the bGH gene) or 6% polyacrylamide gel (the bPRL gene) (0.5 μg/ml ethidium bromide in 1X TBE buffer). The results of electrophoresis were recorded by means of a UVT-1 transilluminator (312 nm) and a ViTran-1 photodocumentation system (Biokom, Russia). The alleles were identified as follows: RsaI(–) and RsaI(+) corresponded to the A and B alleles of the bPRL gene, respectively; AluI(–) and AluI(+), to the V and L alleles of the bGH gene, respectively.

The PopGene [ 26] and STATISTICA 8.0 [ 25] software packages were used for statistical treatment of the results. Pairwise comparisons of the allele and heterozygote frequencies (observed heterozygosity, H obs) in the group of breeds studied were performed using the Fisher exact test. For similar comparisons using published data on other Russian and foreign breeds, we additionally calculated the expected heterozygosity (H exp) where its values were not presented. The G 2 test was used for pairwise comparisons of the breeds with respect to the genotype frequencies and H exp. The fit of the observed frequencies of heterozygous genotypes to those expected from the Hardy–Weinberg equilibrium was tested using the PopGene software. The dependence of the milk yield (in kilograms) and fat and protein contents (in percent) on the bPRL and bGH genotypes in Yaroslavl cattle was estimated by one-way and two-way ANOVA with the use of the STATISTICA 8.0 software.

3. Results and discussion

3.1. Genetic structure

The genotype and allele frequencies of the bPRL gene studied in four Russian cattle breeds are shown in Table 1.

Breed bPRL( RsaI) genotypes RsaI alleles of the bPRL gene
Yakut
( n = 41)
( n) AA (23) AB (14) BB (4) A (60) B (22)
p±s.e. 0.561±0.0780.341±0.0740.098±0.0460.732±0.0980.268±0.098
H exp± s.e. 0.393±0.120
Yaroslavl
( n = 113)
( n) AA (48) AB (50) BB (15) A (146) B (80)
p±s.e. 0.425±0.0470.442±0.0470.133±0.0320.646±0.0640.354±0.064
H exp± s.e. 0.461±0.066
Bestuzhev
( n = 57)
( n) AA (27) AB (24) BB (6) A (78) B (36)
p±s.e. 0.474±0.0660.421±0.0650.105±0.0410.684±0.0870.316±0.087
H exp± s.e. 0.432±0.085
Kostroma
( n = 124)
( n) AA (69) AB (48) BB (7) A (186) B (62)
p±s.e. 0.556±0.0450.387±0.0440.056±0.0210.750±0.0550.250±0.055
H exp± s.e. 0.375±0.133

Table 1.

The frequencies of the RsaI restriction genotypes and alleles of the bPRL gene in four Russian cattle breeds

The G 2 test did not show significant differences between the breeds studied with respect to the genotype or allele frequencies for the bPRL gene. The genotype frequency distributions in the breeds studied fit the Hardy–Weinberg equilibrium. The A allele was prevailing in all breeds studied; its frequency was two to three times higher than the B allele frequency. Thus, these breeds had similar genetic structures in terms of the given SNP despite their different origins and directions of artificial selection.

The distribution of the RsaI alleles of the bPRL gene is characterized by a higher frequency of the A allele in most breeds studied ( Figure 3). In Indian Jersey cattle [ 28], the frequencies of the two alleles were approximately equal. The A allele frequency was lower in Polish Jersey cattle [ 29], which deserves special attention. These differences in allele frequencies may have resulted from different histories of selection for milk yield in different breeds. The breeds of B. taurus exhibit a considerable genotypic variation with respect to the bPRL gene marker used in this study. This can be seen in the diagram of the genotype frequency distribution ( Figure 4). The group of four Russian cattle breeds studied here significantly differed in the genotype frequency distribution from Russian Black & White [ 31] and Russian Red Pied [ 22], but not Lithuanian Black & White [ 30]; (p < 0.006 with the Bonferroni correction).

media/image3.jpeg

Designations: A and B, alleles of the bPRL gene; P, Polish; J, Jersey; I, Indian; L, Lithuanian; B&W, Black and White; R, Russian; Holstein Friesian; Mont, Montebeliard; R.R.P., Russian Red Pied; Kostr, Kostroma; Yar, Yaroslavl; Best, Bestuzhev.The diagram is based on published data on Polish Jersey [ 29], Indian Jersey, Indian Holstein Friesian [ 28], Lithuanian Black & White, Lithuanian Red [ 30], Russian Black & White [ 31], Montebeliard [ 32], and Russian Red Pied [ 22] cattle.

Figure 3.

Distribution of bPRL( RsaI) allele frequencies in the four breeds studied and in other Russian and foreign breeds. The abscissa shows the breeds; the ordinate shows the allele frequencies (in fractions of unity).

media/image4.jpeg

Designations: AA, AB, and BB, bPRL genotypes. For designations of the breeds, see Figure 2.

Figure 4.

Distribution of bPRL genotype frequencies in the four breeds studied and in other Russian and foreign breeds. The radial axes show the breeds and the corresponding genotype frequencies.

3.2. Genetic variation

The observed (H obs) and expected (H exp) heterozygosities ( Table 1) for the bPRL gene did not differ significantly from each other in any breed studied; nor did any two breeds differ in these parameters. This was unexpected, considering the substantial differences between these breeds in the HinfI site of the bPit-1 transcription factor gene (exon 6) and the AluI site of the bGH growth hormone gene (exon 5) (our unpublished data), which are also associated with milk productivity.

In order to compare the genetic diversity in the Russian cattle breeds studied here with those in other Russian and foreign breeds ( Figure 5), we calculated the H exp for all breeds analyzed. The observed and expected heterozygosities for the bPRL gene did not differ significantly in most breeds ( Figure 5). The Indian Jersey breed was an exception (χ 2 = 13.77, ν = 1, p = 0.00021); the excess of heterozygotes in these cattle apparently resulted from artificial selection.

media/image5.jpeg

Designations: H obs and H exp, the observed and expected heterozygosities, respectively. For designations of the breeds, see Figure 2. The diagram is based on the expected heterozygosity levels calculated from published data on Montebeliard [ 32], Russian Red Pied [ 22], Polish Jersey [ 29], Lithuanian Black & White, Lithuanian Red [ 30], Indian Jersey, and Indian Holstein Friesian [ 28] cattle.

Figure 5.

The observed and expected heterozygosities for the bPRL gene markers in Russian and foreign cattle breeds. The abscissa shows the breeds; the ordinate shows the heterozygosity coefficients (in fractions of unity).

Four Russian breeds studied did not differ from one another in the genetic variation (H exp). At the same time, as can be seen in Figure 5, this group had significantly higher H exp values compared to Russian Red Pied [ 22] Indian Holstein Friesian [ 28], and Montebeliard [ 32], cattle (G 2 =36.13, ν = 4, p =3.0∙10 –7; G 2 =39.19, ν = 4, p = 1.0∙10 –7; G 2 =16.72, ν = 4, p = 0.0022; and G 2 =19.41, ν = 4, p = 0.0007, respectively; after the Bonferroni correction).

Thus, the heterozygosity levels in the Kostroma, Bestuzhev, Yakut, and Yaroslavl breeds are significantly higher than in the breeds used for comparison. This indicates a stable state of the breeds studied. However, since artificial selection alters genetic variation, regular monitoring of cattle breeds for the given RsaI genetic marker is advisable.

3.3. Search for associations of SNPs of the prolactin and somatotropin genes with milk production traits in Yaroslavl cattle

Advances in molecular genetic analysis based on PCR have given rise to novel marker systems for direct genotyping at the DNA level that are independent of the animal's sex and age. This makes the procedures less time-consuming and more accurate.

The RsaI polymorphism of the bPRL gene ( Table 1), as well as the AluI polymorphism in exon 5 of the bGH gene, was used to search for associations of these genes with milk production traits, including the milk yield and milk fat and protein contents, in Yaroslavl cattle. We found the following frequency distribution of genotypes: LL ( n = 32), 0.283 ± 0.042; LV ( n = 63), 0.558 ± 0.047; VV ( n = 18), 0.159 ± 0.034. The genotype frequency distribution in these breeds fit the Hardy–Weinberg equilibrium. Table 1 shows the frequency distribution of bPRL genotypes. AAVL ( n = 29) and ABVL ( n = 27) are the most frequent combined genotypes of these genes; the AALL and ABLL genotypes are somewhat rarer ( n = 12 and n = 16, respectively).

We used ANOVA to estimate the isolated effect of the bPRL gene and the combined effect of the bPRL and bGH genes on milk production traits (the milk yield in kilograms and the percentage fat and protein contents for the first three lactations) in Yaroslavl cattle. Table 3 shows the dependence of these parameters for the third lactation on the bPRL( RsaI) polymorphism as estimated by one-way ANOVA. The bPRL gene was found to affect the fat content of milk ( F (2;63) = 3.18, p = 0.048) but not the milk yield or protein content.

Lactation 3
Milk Yield (kg)Fat (%)
FactorD. F.SSMSFPSSMSFP
PRL( Rsa1) 20.030.0150.50.6290.0020.0013.19 0.0478
Error632.0210.0320.02320.0004
Protein (%)
SSMSFP
Rsa1( PRL) 20.00050.00032.340.104
Error630.00730.0001

Table 2.

Results of one-way ANOVA showing the dependence of milk production traits for the third lactation on the PRL( RsaI) polymorphism in Yaroslavl cattle

Figures 6a–6c show the dependences of three milk production traits on the bPRL genotypes in Yaroslavl cattle. As can be seen in Figure 6a, cows with the AA genotype exhibited a significantly higher fat content of milk compared to BB cows ( p = 0.037).

media/image6.jpeg

Figure 6.

Dependence of (a) the milk yield (in kilograms), (b) the fat content, and (c) the protein content of milk of the third lactation on the bPRL genotypes in Yaroslavl cattle. The abscissas show the bPRL genotypes; the ordinates show the mean values of the traits; vertical bars are 0.95 confidence intervals.

Considering that quantitative traits, including milk production traits, are each determined by a number of genes, we also estimated the possible combined effect of two factors, the bPRL and bGH genes, on the milk production traits using two-way ANOVA. The analysis did not reveal a significant isolated effect of any one-locus genotype for any lactation. Still it showed that the bPRL genotype tended to affect the third-lactation milk fat content ( p = 0.064), which was in contrast to the results of one-way ANOVA showing its significant effect ( Table 2). Like one-way ANOVA, the two-way analysis demonstrated that AA cows were characterized by a higher fat content of milk than BB cows. However, the main result concerning the relationship between the bPRL and bGH gene polymorphisms and milk production was that we found combined effects of these two genes on the fat and protein contents of first-lactation milk in Yaroslavl cattle ( F (4;104) = 2.59, p = 0.041 and F (4;104) = 2.93, p = 0.024, respectively) ( Table 3).

In order to identify the genotypes whose carriers significantly differed in the mean percentage fat and protein contents of milk, we performed post-hoc pairwise comparisons of these traits in cows with different combined genotypes of the bPRL and bGH genes. Table 4 shows the results of this analysis. In terms of the milk fat content, cows with the ABLL genotype significantly differed from those with the AAVL, AALL, ABVV, and BBVL genotypes; and cows with the BBVV genotype significantly differed from those with the ABVV and BBVL genotypes. In terms of the milk protein content, AALL cows significantly differed from AAVV, ABVL, ABLL, and BBLL cows; and BBVV cows, from BBLL ones.

EffectLactation 1
TraitMilk YieldFat (%)
D.f.SSMSFPSSMSFP
PRL( RsaI) 20.0060.0030.20.8584.0E-052.0E-050.0600.939
GH( AluI) 20.0840.0422.10.1230.0012.6E-040.9200.401
PRL( RsaI)
*GH( AluI)
40.0640.0160.80.5200.0030.0012.590 0.041
Error1042.0490.0200.0292.8E-04
Lactation 1
TraitProtein (%)
EffectD.f.SSMSFP
PRL( RsaI) 22.4E-051.2E-050.1000.908
GH( AluI) 22.3E-041.1E-040.8900.412
PRL( RsaI)
*GH( AluI)
40.0013.7E-042.930 0.024
Error1040.0131.3E-04

Table 3.

Results of two-way ANOVA showing the dependence of the milk production parameters on the bPRL and bGH genes and their interaction in Yaroslavl cattle

Figure 7 graphically shows the combined effects of the bPRL and bGH genes. The milk yield (in kilograms) and the fat and protein contents of milk of the first lactation in Yaroslavl cows are plotted in Figures 7a, 7b, and 7c, respectively, against the bPRL genotypes ( AA, AB, and BB). The bGH genotypes are indicated by dots of different colors. The points corresponding to the same bGH genotype but different bPRL genotypes are connected by dotted lines. As evident from Figure 7, there was no consistent dependence of the fat or protein content on the dose of any allele of any gene. Let us consider how the mean milk fat content depended on the bPRL genotypes in combination with different bGH genotypes ( Figure 7b). The combinations of the AA genotype of the bPRL gene with different bGH genotypes did not differ significantly from one another, and neither did the combinations of the BB genotype of the bPRL gene with different bGH genotypes. In contrast, the mean percentage fat content of milk of AB cows was significantly lower if they had the LL genotype of the bGH gene than if they had the VV genotype. This indicates gene interaction, which should be taken into account because otherwise the effects of individual genes on the formation of these traits would be incorrectly estimated. Note that, among the combined genotypes of the bPRL and bGH genes shown in Figure 7, the AALL, AAVL, ABVV, and BBVL genotypes were characterized by a significantly higher milk fat content of milk compared to the ABLL gene ( Figure 7b).

Combined bPRL and bGH genotypes ( bPRL( RsaI)– bGH( AluI))
AAVV AAVL AALL ABVV ABVL ABLL BBVV BBVL BBLL
AAVV 0.9330.7250.2430.9450.1820.2720.2670.968
AAVL 0.1020.7010.1630.809 0.041 0.1750.1860.984
AALL 0.038 0.3780.3360.572 0.044 0.1390.3690.806
ABVV 0.0670.4840.9860.125 0.008 0.038 0.9550.339
ABVL 0.4320.181 0.058 0.1250.0690.2180.1430.919
ABLL 0.7330.086 0.029 0.0680.5700.882 0.009 0.259
BBVV 0.0580.3430.7250.7360.1080.064 0.042 0.312
BBVL 0.5590.3680.1530.2080.9600.7270.1590.36
BBLL 0.7130.086 0.036 0.0540.2930.491 0.045 0.387

Table 4.

Results of post-hoc comparisons of the mean first-lactation fat (above the diagonal) and protein (below the diagonal) contents of milk of Yaroslavl cows with different combined bPRLbGH genotypes

Regarding the dependence of the mean protein content of milk in Yaroslavl cattle on the bPRL genotypes combined with different bGH genotypes ( Figure 7c), the combined AALL genotype differs from AAVV, and BBVV differed from BBLL. In contrast, no combination of the AB genotype of bPRL with a bGH genotype differed from its combination with any other bGH genotype in this respect. Note that, in the given sample of Yaroslavl cattle, the AAVV, ABLL, and BBLL genotypes were preferable over the AALL genotype in terms of the percentage protein content of milk, as was the BBLL genotype over the BBVV genotype.

Thus, in Yaroslavl cattle, the AA genotype of the bPRL gene was characterized by a significantly higher percentage fat content of milk than the BB genotype, in contrast to Russian Red Pied cattle, where AB cows had a higher milk fat content than AA and BB cows ( p < 0.05) [ 22]. Our study has been the first to demonstrate the combined effect of the bPRL and bGH genes (i.e., their combined genotypes) on the milk fat and protein contents. A number of other authors have also studied the relationship of milk production with combined genotypes of SNPs in the same or different genes; however, the possibility of their combined effect has almost never been considered. One exception is the study on the effects of the combined genotypes of the AluI and MspI polymorphic sites in exon 5 and intron 3, respectively, of the growth hormone gene in Polish Black & White cattle [ 20]. However, these authors studied combined bGH genotypes as a single factor; hence, they did not consider the effect of interaction between individual SNPs on the traits studied. In addition, the marker system used by them was hardly suitable for revealing the interaction of these SNPs because they were located in the same gene.

media/image7.jpeg

Figure 7.

Dependence of (a) the milk yield (in kilograms), (b) the fat content, and (c) the protein content of milk of the first lactation on combined genotypes of the bPRL and bGH genes in Yaroslavl cattle. The abscissas show the bPRL genotypes ( AA, AB, and BB); the ordinates show the mean values of the traits; vertical bars are 0.95 confidence intervals. Different bGH genotypes are denoted by dots of different colors; the points corresponding to the same bGH genotype but different bPRL genotypes are connected by dotted lines.

3. Conclusion

We have found that four Russian cattle breeds, Yakut, Bestuzhev, Kostroma, and Yaroslavl cattle, are similar in genetic structure. All of them are characterized by a low frequency of the BB genotype of the RsaI polymorphic site in the bPRL gene (from 0.056 ± 0.021 to 0.133 ± 0.031) and high frequencies of heterozygotes and homozygotes for the A allele. The breeds of B. taurus exhibit a considerable genotypic variation with respect to the bPRL gene marker used in this study. The group of these four Russian cattle breeds significantly differs in the genotype frequency distribution from other breeds, such as Russian Black & White and Russian Red Pied. At the same time, these four breeds do not differ significantly in the observed or expected heterozygosity for the bPRL gene either from each other or from other breeds used for comparison.

We have demonstrated a combined effect of the bPRL and bGH genes on the percentage protein and fat contents of milk. Each trait has been found to be significantly positively associated with some of the combined genotypes. No genotype has been found to positively affect both traits. At the same time, some genotypes are associated positively with one trait and negatively with the other one: ABVV cows are characterized by a high fat content and low protein content of milk, while this is the other way round with ABLL cows. This could be used for selecting cattle for high individual productivity traits. Note that only one combined genotype ( BBVV) is unfavorable in terms of both traits. However, being doubly homozygous, it may serve as a reserve for obtaining genotypes that are valuable in terms of either fat ( ABVV and BBVL) or protein ( AAVV and BBLL) content of milk. Thus, the study of the combined effects of the bPRL and bGH genes and the breeding practice taking these effects into account allow the cattle productive potential to be analyzed in more detail. In addition, involvement of epistatic gene interaction in the formation of selectively valuable quantitative traits has been further confirmed.

We believe that data on the bPRL( RsaI) SNP marker and the heterozygosity estimates calculated from its allele frequencies may be used in programs for conservation of aboriginal breeds while maintaining the optimal balance between the goals of artificial selection and preservation of the genetic diversity, which is necessary for sustained reproduction of cattle breeds with all their unique characters. This is especially important because of the rapid decrease in the stocks of cattle breeds and the related threat of partial loss of the genetic resources for stockbreeding on a global scale. [ 33].

Acknowledgement

This study was supported by the Subprogram of the Presidium of the Russian Academy of Sciences "Gene Pools and Gene Diversity," State Contract no. 14.740.11.0164 in the framework of the Federal Target Program "Science and Education Professionals of Innovative Russia," and the Russian Foundation for Basic Research (grant no. 12-04-92214 Mong_a).

References

1 - S Bernichtein, P Touraine, V Goffin, REVIEW New concepts in prolactin biology. Journal of Endocrinology 2010 206 1 1 11
2 - N. D Horseman, W Zhao, E Montecino-rodriguez, M Tanaka, K Nakashima, S. J Engle, F Smith, E Markoff, K Dorshkind, Defective mammopoiesis, but normal hematopoiesis, in mice with a targeted disruption of the prolactin gene.EMBO Journal 1997 16 23 6926 6935
3 - B Rischkowsky, D Pilling, editors. The State of the World’s Animal Genetic Resources for Food and Agriculture Food and Agriculture Organization of the United Nations. FAO. Rome. 2007 978-9-25105-762-9
4 - L. K Ernst, N. G Dmitriev, I. A Paronyan, editors. Geneticheskie resursy selskokhozyaistvennykh zhivotnykh v Rossii i sopredel’nykh stranakh (Genetic Resources of Farm Animals in Russia and Neighboring Countries). St. Petersburg, Russia: VNIIGRZh; 1994
5 - G. P Korotov, Krupnyi Rogatyi Skot Yakutskoi FSSZ i Metody Ego Uluchsheniya (Yakut Cattle and Breeding Methods of Its Improvement). Yakutsk, USSR: Yakutknigoizdat; 1983
6 - G. E Sulimova, I. V Lazebnaja, A. V Perchun, V. N Voronkova, M. N Ruzina, G. A Badin, Uniqueness of Kostroma breed of cattle from a position of molecular genetics. Advances in science and technology of Agro-Industrial Complex 2011 9 52
7 - J. I Weller, Quantitative Trait Loci Analysis in AnimalsLondon: CABI Publishing; 2001
8 - J Bennewitz, N Reinsch, V Guiard, S Fritz, H Thomsen, C Looft, C Kühn, M Schwerin, C Weimann, G Erhardt, F Reinhardt, R Reents, D Boichard, E Kalm, Multiple quantitative trait loci mapping with cofactors and application of alternative variants of the false discovery rate in an enlarged granddaughter designGenetics 2004 168 2 1019 1027
9 - M. S Ashwell, Rexroad Jr. C.E., Miller R.H., Van Raden P.M., Da Y. Detection of loci affecting milk production and health traits in an elite US Holstein population using microsatellite markers Animal Genetics 1997 28 3 216 222
10 - Q Zhang, D Boichard, I Hoeschele, C Ernst, A Eggen, B Murkv, M Pfistergenskow, L. A Witte, F. E Grignola, P Uimari, G Thaller, M. D Bishop, Mapping quantitative trait loci for milk production and health of dairy cattle in a large outbred pedigree.Genetics 1998 149 4 1959 1973
11 - Y Plante, J. P Gibson, J Nadesalingam, H Mehrabani-yeganeh, S Lefebvre, G Vandervoort, G. B Jansen, Detection of quantitative trait loci affecting milk production traits on 10 chromosomes in Holstein cattle.J. Dairy Sci. 2001 84 6 1516 1524
12 - http://fermerru/files/blog/ 2012yak4.jpg
13 - http://agrolibru/rastenievodstvo/item/f00/s00/e0000992/index.shtml
14 - http://agrobkru/kostromskaya
15 - http://frunze37ru/zhivotnovodstvo/
16 - L. A Svensson, K Bondensgaard, L Norskov-lauritsen, L Christensen, P Becker, M. D Andersen, M. J Maltesen, K. D Rand, J Breinholt, Crystal Structure of a Prolactin Receptor Antagonist Bound to the Extracellular Domain of the Prolactin Receptor. The Journal of Biological Chemistry 2008 283 27 19085 19094
17 - S. A Camper, D. N Luck, Y Yao, R. P Woychik, R. G Goodwin, R. H Lyons, F. M Rottman, Characterization of the bovine prolactin gene.DNA 1984 3 3 237 249
18 - K Inoue, H Goda, Mogi Ch., Tomida M., Tsurugano Sh.. Chapter 6. The Role of Glucocorticoids and Retinoic Acid in the Pituitary Endocrine Cell Differentiation. In: Handa R.J., Hayashi Sh., Terasawa E., Kawata M. (eds.) Neuroplasticity, Development, and Steroid Hormone ActionBoca Raton, FL: CRC Press; 2001 73 80
19 - E. R Chung, E. R Chung, T. J Rhim, S. K Han, Associations between PCR-RFLP markers of growth hormone and prolactin genes and production traits in dairy cattleKorean J. Anim. Sci. 1996X) 321 336
20 - A Dybus, W Grzesiak, I Szatkowska, P Błaszczyk, Association between the growth hormone combined genotypes and dairy traits in Polish Black-and-White cows. Animal Science Papers and Reports 2004 22 2 185 194
21 - P Chrenek, J Huba, M Oravcova, L Hetenyi, D Peskovieova, J Bulla, Genotypes of bGH and bPRL genes in relationships to milk production. 50th Annual Meeting: Book of Abstracts, 1999Eaap, Zurich. 1999.
22 - M Alipanah, L Kalashnikova, G Rodionov, Association of prolactin gene variants with milk production traits in Russian Red Pied cattle. Iranian Journal of Biotechnology 2007 5 3 158 161
23 - M. P Sabour, C. Y Lin, C Smith, Association of genetic variants of bovine growth hormone with milk production traits in Holstein cattle Journal of Animal Breeding and Genetics 1997X) 435 442
24 - L Zwierzchowski, J Krzyzewski, N Strzalkowska, E Dymnicki, Effect of polymorphism of growth hormon (GH), Pit-1, and leptin (LEP) genes,cow’s age, lactation stage and somatic cell count on milk yield and composition of Polish Black-and-White cows. Animal Science Papers and Reports 2002
25 - A Mitra, P Schlee, C. R Balakrishman, F Pirchner, Polymorphisms at Growth-Hormone and Prolactin Loci in Indian Cattle and Buffalo Journal of Animal Breeding and Genetics 112 71 74 1995
26 - F. C Yeh, R Yang, C Boyle, T. B. J Ye, Z-H Mao, J.X. PopGene, the User-Friendly Shareware For Population Genetic Analysis. Molecular Biology and Biotechnology Centre, University of Alberta, Canada. 1997
27 - StatSoftInc. STATISTICA (Data Analysis Software System), version 8.0. 2008 www.statsoft.com.
28 - R Kumari, K. M Singh, K. J Soni, R. K Patel, J. B Chauhan, K. R. S Sambasiva-rao, Genotyping of the Polymorphism Within Exon 3 of Prolactin Gene in Various Dairy Breeds by PCR RFLP (Brief report)Archiv Tierzucht / Archives Animal Breeding 2008 51 3 298 299
29 - A Dybus, W Grzesiak, H Kamieniecki, I Szatkowska, Z Sobek, P Blaszczyk, E Czerniawska- Piatkowska, S Zych, M Muszynska, Association of Genetic Variants of Bovine Prolactin With Milk Production Traits of Black-and-White and Jersey Cattle. Archiv Tierzucht / Archives Animal Breeding 2005 48 2 149 156
30 - R Skinkyte, L Zwierzchowski, L Riaubaite, L Baltrenaite, I Miceikiene, Distribution of Allele Frequencies Important to Milk Production Traits in Lithuanian Black and White and Lithuanian Red Cattle. Veterinary Medicine and Zootechnics (Veterinarija ir Zootechnika) 2005 31 53 93 96
31 - T. S Goryacheva, G. M Goncharenko, Polymorphism in κ-Casein and Prolactin Genes and Their Influence on Dairy Productivity of Cows of the Black-and-White Breed. Agricultural biology 2010 4 51
32 - N Ghasemi, M Zadehrahmani, G Rahimi, S. H Hafezian, Associations Between Prolactin Gene Polymorphism and Milk Production in Montebeliard Cows. International Journal of Genetics and Molecular Biology 2009
33 - FAO 2007 The State of the World’s Animal Genetic Resources for Food and Agricultureedited by Barbara Rischkowsky & Dafydd Pilling. Rome.