Differential A1/A2 β-casein (CSN2) gene-derived allelic and genotypic frequencies across Ecuadorian exotic dairy cattle breeds

Recently, a growing interest in dairy cattle selection has been triggered due to specific genetic variations of the β-casein (CSN2) gene which was associated to several human disorders. The aim of the present study was to evaluate the A1/A2 β-casein (CSN2) gene-derived allelic and genotypic frequencies in five Ecuadorian dairy cattle breeds to determine the genetic diversity of the A1/A2 β-casein locus. Genomic Deoxyribonucleic Acid (gDNA) was extracted from hair samples collected from different dairy cattle breeds, namely Holstein Friesian (HF), Brown Swiss (BS), Jersey (J), Gyr (G), and crossbreds (CB) raised commonly in seven different regions in order to identify A1 and A2 variants. Genomic determination of the β-casein (CSN2) gene variants was carried out evaluating the Exon 7 region in order to identify and confirm the accuracy of the A1/A2 β-casein (CSN2) gene-derived polymorphisms related to the genotypes and the allelic profiling. A2 allele frequencies were determined as 0.628, 0.643, 0.727, and 0.644 in HF, BS, J, and CB, respectively. In contrast, the lowest A1 allele frequency (0.145) as well as the greatest A2 allele frequency (0.855) was detected in G. No significant differences for allelic frequency were observed among breeds (p > 0.05). Homozygous A1 genotype frequency in HF, BS, J, and CB was 0.140, 0.110, 0.050 and 0.150, respectively. Interestingly, the A1A1 genotype was absent in G (p < 0.05). A1A2 genotypes were 0.290, 0.490, 0.450, 0.290, and 0.420 in HF, BS, J, G, and CB, respectively. Finally, the A2A2 genotype was 0.400, 0.400, 0.500, and 0.460 in HF, BS, J, and CB, respectively. However, the greatest A2A2 genotype frequency was observed in G (0.710). No deviation from the Hardy–Weinberg equilibrium was found in any breed in terms of β-casein gene (p > 0.05). Moreover, the results revealed a geographically related distribution of A1/A2 allele and genotype frequencies of β-casein (CSN2) gene. In terms of genetic diversity a differential distribution of heterozygosity was observed among breeds being average Ho, He, uHe, and F values determined as 0.429, 0.438, 0.439, and 0.022, respectively. In conclusion, the β-casein (CSN2) gene was polymorphic in all the cattle breeds analyzed. The A1/A2 allele and genotype frequencies varied depending on the breed and the geographic region which may be conditioned by the specificities related to different breeding selection programmes. A2 allele and A2A2 genotype frequency were particularly greater in Gyr cattle which have a great potential for A2 milk production without compromising genetic diversity.

1 Introduction

The introduction of new exotic dairy cattle breeds in Ecuador during the second half of the 20th century opened up new opportunities for the rural population and the national dairy markets (1, 2). During the last decade, the genetic selection of dairy cattle has been crucial for increasing milk quality and productivity; however, undesirable effects such as gastrointestinal disorders resulting from the consumption of raw or processed milk have been observed in the Ecuadorian population (3–5). Several studies reported more than twenty genes associated with essential milk trait composition (6–8). The sequence analysis of the milk protein CSN2 gene (exon 7 regions in chromosome 6) revealed different genotypes and allele frequencies (e.g., A1, A2, A3, B, among others) of the β-casein variants which were intensively genotyped in dairy cattle for selection purposes (5, 9). The β-casein protein is composed of 209 aa being the 2nd most abundant protein in milk (10, 11). During the last years, the determination of genetic variations (polymorphisms) of bovine β-casein protein has become a common practice in dairy cattle being the A1 and A2 alleles the most common among the 13 alleles identified (12). One of the main reasons was related to the A1 allele variant which has been considered as risk factor for several heart diseases, diabetes mellitus (type 1), autoimmune diseases, and milk intolerance issues in humans (13). These issues are derived from the release of the bio-active opioid β-casomorphin-7 (BCM-7) during the histidine proteolytic digestion of A1 monomers (14). Several studies reported that the A1 allele variant appeared as a mutation of the A2 allele resulting in a polymorphism mutation at the position 67 of the CSN2 gene in Bos taurus cattle (15–17). Basically, there is a change of adenine for a cytosine in the amino acid chain, and as a consequence, the mutation causes the exchange of proline for histidine (5). The export of specialized genetics (semen, embryos and live animals) of different dairy cattle breeds, caused the A1 allele variant to widespread in most parts of the world (3, 18). On the contrary, the A2 allele variant has been considered as non-problematic regarding health-related issues (19). Thus, the A2 milk can be obtained from cows carrying the A2A2 genotype of the CSN2 gene, being therefore the cows owning the genotypes A1A1 and A1A2 undesirable for row milk consumption. However, although there are studies suggesting that milk from cows with the A2 genotype is easily digestible, it should be noted that, in cheese production, milk containing A2 alleles is not recommended, as it is associated with a deterioration in the technological characteristics of the milk and, consequently, with a less efficient cheese-making process (20, 21). Moreover, the allele and genotype frequencies vary according to the cattle breed (5, 17). Thus, there is an imperative need of developing genetic selection programs for breeding dairy cattle intended for the consumption of raw milk with the aim of increasing the A2 allele variant frequency irrespective of the country considered (22). In Ecuador, according to INEN-ESPAC (23) the distribution of dairy cattle breeds showed the Holstein, Jersey and Brown Swiss as the main dairy breeds in the highlands and the Jersey, Brown Swiss, and Gyr breeds in the coast and the Amazon region (23). On the other hand, the Ecuadorian Ministry of Agriculture and Livestock (24) has established through ministerial agreement No. 061–2021 the regulations for the certification of A2A2 genotype, which aims to promote A2A2 milk production through the identification of individuals and herds at national level (24).

In view of the above, the genetic characterization programs in dairy cattle are crucial not only for increasing the milk quality and yield but also for assuring the optimal milk composition and characteristics which are related to human health (25, 26). Thus, the molecular identification and analysis of milk-related genes as CSN2 including different polymorphisms, contribute to increase the efficiency of cattle selection programmes (27, 28). There are currently no studies in the Ecuadorian dairy cattle population and its relationship to different genotypes and allele frequencies of the milk protein CSN2 gene (29, 30). Therefore, the aim of the present research was to characterize the CSN2 gene to determine the allelic (A1 and A2) and genotypic (A1A1, A1A2, and A2A2) frequencies of exotic dairy cattle breeds and crossbreds in Ecuador. This will provide a baseline on the knowledge of the CSN2 gene that will be undoubtedly useful for the different dairy cattle producers and associations with the aim of improving the cattle breeding and selection programmes.

2 Materials and methods

2.1 Ethical statement

The study was conducted according to the guidelines of the Declaration of Helsinki and following the Code of Ethics for animal experiments as reflected in the ARRIVE guidelines available at http://www.nc3rs.org.uk/ARRIVEchecklist (Accessed on 10 August 2024). However, the present research did not require any ethical review and approval by a specific Bioethics Committee for the use of experimental animals, since the hair sample collection is considered a non-invasive method and were obtained during routine diagnostic and health control services in each farm. Therefore, none of the animals were handled or restricted at any time for this study. Then, the studies were conducted in accordance with the local legislation and institutional requirements.

2.2 Animals, sample collection and geographical location

Hair root samples (at least 20–25 hair roots per animal) were collected for DNA extraction from the hair follicles (bulbs found at the base of the hair follicle). The samples (n = 1,599) were collected individually from ears avoiding contamination and assuring sufficient quantity and quality of DNA. The collection was carried out by the technical staff (veterinarian) of each farm for routine diagnostic and health control services. The samples belonged to farms located in seven provinces from two geographical regions [West region: Santo Domingo, Manabí, and Guayas (Coast) and East region: Pichincha, Cotopaxi, Tungurahua, and Azuay (Highlands)] in Ecuador. Four pure-bred [Holstein Friesian (n: 701), Gyr (n: 258), Jersey (n: 86), and Brown Swiss (n: 230)] and a cross-bred [always Holstein Friesian x another dairy breed (n: 324)] populations were analyzed for the detection of A1/A2 β-casein-derived allelic and genotypic frequencies using molecular techniques.

2.3 β-casein characterization using genomic analysis

The characterization of β-casein (CSN2) gene was carried out by genomic analysis. Thus, the individuals were studied with genomic tests using the Illumina Bovine SNP50 BeadChip (Clarifide–Zoetis, San Diego, CA, USA and Igenity-Neogen, Lansing, MI, USA). For the genotyping sampling, the Gene Max kit (TSU, ALLFLEX, Kalamazoo, MI, USA) was used to take samples from each cow ear (pinna). A chip with 12,000 DNA markers (12 K) was used, which provides the most relevant information regarding production and milk-type-related characteristics with a ~ 70% reliability based on a panel of 50,000 (HD50K) genetic markers from the Animal Improvement Programs Laboratory (AIPL) (Zoetis–Clarifide, San Diego, CA, USA) and the Animal Laboratory (Neogen– Igenity, Lansing, MI, USA). The information predicted mutations of the CSN2 gene in bovine and have been widely used for the characterization of elite sires and dams in the dairy industry (Clarifide™ Plus, Zoetis Genetics, San Diego, CA, USA and Igenity® Plus-Neogen, Lansing, MI, USA). The initial genomic analysis (screening) allowed identifying and characterizing those individuals related to different variants for the β-casein (CSN2) gene (A1A1, A1A2 and A2A2) to later determine the genotype and allele frequencies. Other β-casein (CSN2) gene variants such as A3, A4, B, C, D, E, F, G, H1, H2, I, J, and K were not considered for the present study. The distinction between A1 and A2 genotype was carried out, as the β-casein variants were categorized into two groups based on the amino acid at position 67 of the protein: Type A1 (histidine) and Type A2 (proline).

2.4 Frequency-based population genetic parameters related to A1/A2 β-casein (CSN2) gene

Frequency-based population genetic parameters related to A1/A2 β-casein (CSN2) gene of different dairy cattle breeds raised in Ecuador are showed in Supplementary Table S1. Each population genetic indicator was based on methodologies carried out previously by different authors. Then, the following parameters were evaluated: (i) allelic frequency (31); (ii) genotypic frequency (Haploid Data) (32); (iii) Hardy–Weinberg disequilibrium (32); number of alleles (33); number of effective alleles (34); Shannon diversity index (34); observed heterozygosity (32); expected heterozygosity (31); unbiased expected heterozygosity (35); fixation index (31); Nei unbiased genetic distance (32); and finally, paired Fst (36). The corresponding variables, formulas, value ranges, parameter descriptions and the reference sources are displayed in Supplementary Table S1.

2.5 Statistical analyses

Statistical differences among breeds and regions were analyzed using the excel plugin GenAlEx 6.5 (35). To understand the distribution of CSN2 gene across populations, the allele and genotype frequencies were calculated. Based on the obtained results, the population differentiation and selection effects Hardy–Weinberg disequilibrium, number of alleles, number of effective alleles, Shannon diversity index, observed/expected and unbiased expected heterozygosity and fixation index were obtained. Finally, the comparison between populations was tested by paired unbiased Nei genetic distance and paired Fst. Significance level was set at p ≤ 0.05.

3 Results

3.1 A1/A2 β-casein-derived allelic and genotypic frequencies in different dairy cattle breeds

Differential pattern of the A1A1, A1A2, and A2A2 genotype frequencies was observed across Brown Swiss, Cross-bred, Gyr, Holstein-Friesian, and Jersey cattle breeds (Table 1). The greatest frequency was observed for A2A2 genotype in Gyr (0.710) and Jersey (0.500), followed by heterozygous A1A2 in Brown Swiss (0.490) and Holstein-Friesian (0.460) and homozygous A1A1 genotype in Cross-bred (0.150) and in Holstein-Friesian (0.140). Interestingly, the A1A1 genotype was absent in G (p = 0.006) (Table 1). Regarding frequency of A1 allele, the highest value was observed in Holstein Friesian (0.372) followed by Brown Swiss (0.357) and Cross-Bred (0.356). Finally, with regard to the frequency of A2 allele, the highest value was observed in Gyr (0.855), followed by Jersey (0.727). The lowest values of A2 allele were observed in Holstein Friesian (0.628), in Brown Swiss (0.643) and in crossbred (0.644), respectively. No significant differences for allelic frequency were observed among breeds (p > 0.05) (Table 2).

Table 1

Table 1. A1/A2 β-casein-derived genotypic frequencies in different dairy cattle breeds raised in Ecuador.

Table 2

Table 2. A1/A2 β-casein-derived allelic frequencies in different dairy cattle breeds raised in Ecuador.

3.2 A1/A2 β-casein-derived allele and genotype frequencies in dairy cattle breeds and geographical location

Allelic and genotypic frequency profiles were also generated from samples in different dairy cattle breeds raised in different geographical regions (Highlands and Coast; Tables 1, 2). Specific pattern of allelic and genotype frequency distribution was observed across groups based on geographical location. For the samples screened, the frequency of A2 allele and A2A2 genotype was the highest in the Coast region, being 0.784 and 0.610, respectively. On the contrary, the frequency of A1 allele and A1A1 genotype was the highest in the highlands region, being 0.361 and 0.130, respectively, compared to the Coast region where the frequencies for A1A1 genotype and A1 allele were 0.040 and 0.216, respectively (Tables 1, 2). Finally, the A1A2 genotype was predominant in the highlands region (0.450) compared to the coast region where A1A2 genotype frequency was lower (0.360). No significant differences for allelic frequency were observed between geographical regions (p > 0.05).

3.3 Genetic diversity indices in dairy cattle breeds related to A1/A2 β-casein (CSN2) gene

Genetic diversity indices of different dairy cattle breeds raised in Ecuador are showed in Table 3. The allele and genotype-derived results obtained in the current study revealed that a significant level of genetic diversity was observed in dairy cattle breeds. All the selected loci were found 100% polymorphic in all dairy cattle breeds except in Gyr, where no A1 allele and A1A1 genotype was found. The number of alleles observed at each locus serves as an indicator of genetic diversity, which directly influences the differentiation of dairy cattle breeds as each locus exhibited different number of alleles; however, in the present study the number of observed alleles (Na) detected all dairy cattle breeds was just two (A1 and A2) (Table 3). The average number of effective alleles (Ne) varied from 1.331 in Gyr breed to 1.877 in Holstein Friesian breed with mean average Ne across cattle breeds of 1.781. Among all the breeds studied, highest detected heterozygosity (Ho) was found in Brown Swiss breed (0.487) while lowest was observed in Gyr breed (0.291). Further, average observed heterozygosity across the whole breeds was found to be 0.429 (Table 3). The expected heterozygosity (He) was found either approximately similar than Ho in all the breeds studied, except for Jersey breed where was much less (0.397) compared to Ho (0.453). Finally, the highest value related to unbiased expected heterozygosity (uHe) was observed in Holstein-Friesian breed (0.467) while the lowest was detected again in Gyr breed (0.249). The fixation index (F) to detect differences in the dairy breed related to genetic structure indicated that the more fixed positive value was detected in crossbred individuals (0.092) while the less and negative F value was observed in Gyr breed (−0.170). The negative F value in Brown Swiss, Gyr and Jersey breeds indicated the absence of inbreeding in the analyzed breeds.

Table 3

Table 3. Genetic diversity indices in different dairy cattle breeds and geographical regions related to A1/A2 variants of β-casein (CSN2) gene.

3.4 Genetic diversity indices related to β-casein A1/A2 variants and geographical location

Genetic diversity indices of different geographical locations where different dairy cattle breeds are raised in Ecuador are showed in Table 3. In the present study the number of observed alleles (Na) detected was just two (A1 and A2). The average number of effective alleles (Ne) varied from 1,513 in the coast to 1,856 in the highlands region. The highest observed heterozygosity (Ho) was found in the highlands region (0.453) while lowest was detected in the coast (0.358). In the other hand, he expected heterozygosity (He) was pretty similar in the highlands compared to Ho; however, the coast region was lower (0.339) compared to Ho (0.358). Finally, the highest value related to unbiased expected heterozygosity (uHe) was observed in the highlands population (0.461) while the lowest was detected again in the coast population (0.339). The fixation index (F) indicated a positive value (0.018) in the highlands cattle population, while a negative F value was observed in the coast population (−0.055).

3.5 Differentiation and genetic relationship among dairy cattle breeds (Fst: fixation index and DA: unbiased Nei genetic distance)

The pair-wise Fst values of different dairy cattle breeds are displayed in Table 4. The Fst values, which indicate breed differentiation, demonstrated that 4.74% of the total genetic variability resulted from allelic differences among the different breeds. The remaining 95.26% of genetic variability corresponded to differences among individuals within the breed. The Fst values ranged from 0.000 to 0.067. This indicates the least differentiation among Holstein Friesian, Brown Swiss and Crossbred individuals, while the greatest divergence was observed between Holstein Friesian-Brown Swiss-Crossbred vs. Gyr breed. The pair-wise DA values for stabilizing the genetic relationship among the different dairy cattle breeds are showed in Table 4. The lowest DA value (shortest distance) was observed among the Holstein Friesian-Brown Swiss-crossbred (0.000) followed by Jersey-Brown Swiss (0.007) and Jersey-crossbred (0.008). The Gyr was the most distant breed, displaying the largest DA when compared to the Holstein Friesian (0.067), to the cross-bred (0.056) and to the Brown Swiss (0.055) (Table 4).

Table 4

Table 4. Pairwise unbiased Nei distance (below diagonal) and Fst (above diagonal) among different dairy cattle breeds.

4 Discussion

The present research aimed to study the distribution of the beta-casein gene (CSN2) genotypes and alleles in the main dairy cattle breeds in Ecuador. The first studies related to the CSN2 gene were carried out more than two decades ago; however, the market demand for milk production with specific characteristics related to A1/A2 variants began only a decade ago, mainly in countries such as New Zealand, Australia, Canada, United States, among others (37, 38). Due to the advances made in the cattle selection using molecular techniques (genotyping), there is currently great interest in different countries in selecting sires and cows for genetic improvement related to the fixation of the A2 allele in both sexes (12). In this way, AI and genetic companies for cattle breeding and selection started to select sires as carriers of the beta-casein CSN2 gene in favor of the A2 allele promoting the use of semen from sires carrying this allele (39–41).

It has been previously described that the A2 allele mutation to A1 has occurred in the Holstein Friesian breed (5). The allele A1 frequency in Holstein Friesian has been described as high by different studies (9, 40). However, the intense selection in this breed could influence the higher A1 allele frequency observed which could be related to individuals producing higher milk yields (40). In the present study, the A2 allele frequency was slightly lower in Holstein Friesian and Brown Swiss cattle compared to the other breeds, with the exception of the Jersey and Gyr breed where the A2 allele frequency was higher (42, 43). In the case of Holstein Friesian, the A2 allele frequency was similar to that observed in other countries such as Croatia, Italy and Sweden (40), but lower than in Chile, Slovakia, and Greece (38) (Supplementary Table S2). In this study the A2 allele frequency observed in Holstein Friesian was higher than that observed in Peru, Spain, Turkey, Italy (north), and China (22). Finally, the A2 allele frequency in this breed was 1.5 to 2 times higher in Ecuador than in India, Russia, and Pakistan (40). In addition, due to the great expansion of this breed worldwide, the CSN2 gene mutation has been spread as well, partly due to the use of semen from countries that maintain solid genetic improvement programs and the massive export of genetic material from live sires and embryos (44). The high A2 allele frequency in Ecuador could be due to the fact that genetic material has been imported mainly from A2A2 genotyped sires. For example, in Ecuador out of 311 imported sires born between 2016 and 2021, a total of 151 owned the A2A2 genotype and 130 sires the A1A2 genotype (45). This fact indicates that in recent years the application of new technologies in favor of the selection of individuals carrying the A2A2 genotype is steadily increasing. In addition to the genetic influence of European Holstein Friesian sires for the selection programs related to the A2A2 genotype in Ecuador, the genetic material also included several genes for improving the milk protein profile (45).

The prevalence of the A2 allele in the Guernsey and Jersey breeds was more frequent than in other European breeds (46, 47). In the present study, the A2 allele frequency was higher for all breeds analyzed, except for the Gyr (48). The A2 allele frequency obtained was similar to those observed in countries such as Chile, Mexico, Turkey, Japan and Sweden in the same breed (17) and higher than those observed in Russia and India, even tripling the value observed in the Pakistani Jersey (49) (Supplementary Table S2). This fact could be due to the fact that in developed countries, the use of ART’s using genetic material from Jersey sires of US, Canadian or European origin has been much more accentuated and, therefore, they benefit indirectly from the increased population of individuals carrying the A2 allele (50). Thus, the low use of imported semen together with the lack of local genotyping services for the identification of A2A2 genotype carriers led to the fact that there was no genetic improvement in favor of this genotype in less developed countries. Therefore, this fact could possibly favor that the A1 allele frequency has been maintained at above average levels in non-developed countries. In any case, according to the results obtained in the present study, the Ecuadorian Jersey breed seems to have a great potential to increase the A2 allele frequency, both within the breed itself, as well as within crossbred individuals using the Jersey breed. However, Ecuador, similarly to currently other less developed countries, lacks strict selection and genetic evaluation programmes for controlling and selecting individuals based on the CSN2 gene and A1/A2 variants. In the case of the Brown Swiss breed, it has been observed that the A2 allele frequency reached the value of 0.785 and the frequency of the A2A2 genotype reached the value of 0.62 (51, 52). These frequencies observed in the Brown Swiss population were greater than those obtained for the same breed in Ecuador. The difference could be due to the fact that there has not been a specific selection program to increase the A2 allele frequency in Ecuador, which also explains the high A1 allele frequencies observed in the Ecuadorian Brown Swiss cattle population. In recent years, due to the genetic influence of European Brown Swiss sires (3, 53), the allele frequency in favor of A2 could probably increase in the coming years. Regarding the A2 allele frequencies published in the scientific literature, this allele has been variable within the European cattle breeds (Bos taurus) (54). However, the A2 allele frequency in the Indian zebu breeds (Bos indicus) was predominant (55) (Supplementary Table S2). Regarding Gyr breed, in the present study the A2 allele frequency was lower than that observed in Brazil (country of origin of this breed) (48). Taking into account that the genetic resources entering Ecuador came from Brazil, the fact that in the present study the frequency values were lower could be due to the Gyr breed entered Ecuador for the first time in the 1960s (30). From that moment on, the reproduction was carried out without control using a reduced number of individuals. In addition, most of the breeders used Gyr for crossbreeding with Holstein Friesian individuals producing the Girolando (Gyr x Holstein Friesian) synthetic breed, which in many cases remained within the herd jointly (56). So, it was possible that these Girolando individuals introduced the Holstein Friesian A1 allele into the Gyr breed in a higher proportion in Ecuador. Thus, just since 2023 the Ecuadorian Gyr-Girolando Association registered the first foundational population of this breed, which allowed better control of both the A2 allele and the A2A2 genotype through controlled genetic selection programs. The results obtained in the present study indicate that the Ecuadorian Gyr breed has prospects of increasing the frequency of the A2 allele in a few generations by selecting genotyped sires and dams as happened in Australia between 2000–2023, which managed to increase the A2 allele frequency from 0.32 to 0.52 (57).

With regard to the genetic diversity, it should be taken into account that the expected heterozygosity (He) quantifies the equality of allele frequencies in the loci (58). The genetic population structure can be evaluated when He was compared to the observed heterozygosity (Ho) (59). Thus, when analyzing the Hardy–Weinberg equilibrium and the Ho value is greater than the He, the population presents an excess of heterozygotes compared to what’s expected. On the contrary, when Ho is lower than He value, the population present a deficit of heterozygotes. Thus, in the present study an excess was observed in the dairy breeds evaluated, specifically in Brown Swiss, Gyr and Jersey, while a deficit was observed in crossbred and Holstein Friesian breed. The highest heterozygosity was observed in the Brown Swiss, and the lowest in the Gyr breed. He value in Holstein Friesian was the highest and the lowest was observed in the Gyr breed. In the case of Ecuadorian Holstein Friesian, the He value was higher than that observed in the population of the same breed in Ukraine (He = 0.343) (60) and similar to the local Black-White breed in the same country (He = 0.45–0.49) (9). Moreover, the He value was similar when compared to the Holstein Friesian from Mexico (Ho = 0.464) (61) from Chile (He = 0.375) (62) and from Turkey (He = 0.499) (17). Both He and Ho values found in the Ecuadorian Jersey breed (Ho = 0.453; He = 0.397) were greater than those observed in the Turkish Jersey cattle population (Ho = 0.262; He = 0.262) (50) in which the two variables (He and Ho) were in Hardy–Weinberg equilibrium. Thus, the values obtained in the same breed in Ecuador were similar to the Chilean Jersey breed population (He = 0.394) (62) and similar to those observed in Turkish Jersey (He = 0.338) (17). Regarding the Gyr breed, the He value (0.248) was greater than that described by Khan et al. (40), which was 0.06 (the presence of A1A1 was null) while the A1A2 genotype was 0.07. Regarding the Ecuadorian Brown Swiss breed, the He value (0.459) was similar to that obtained in the Turkish Brown Swiss (He = 0.466) (17).

The highest F (Fixation Index) value observed in the present study was 0.092 in crossbred and 0.007 in the Holstein Friesian breed, which could indicate the absence of active formation processes (e.g., selection, population size, and inbreeding, among others) and therefore, corresponds completely to the Hardy–Weinberg equilibrium of these populations (63). The F (Fixation Index) values obtained in the present study from the Brown Swiss, Jersey and Gyr breed were negative, although quite close to 0 (could even be considered = 0), therefore, there were also in Hardy–Weinberg equilibrium (9). Moreover, this could also indicate an excess of heterozygotes in the evaluated populations as established in the Polish White-Backed, Polish Red, Polish Holstein-Friesian, and Simmental populations (60). The Ne values found in the present study (1.331 to 1.848) would indicate an average level of polymorphism of the locus relative to the classical biallelic system (31). In general, the Fst and Pairwise Unbiased Nei Distance values observed in the present study were relatively low which would indicate a low differentiation and genetic distance among the evaluated breeds similar to those observed in other studies (60).

When analyzing the allelic frequencies by geographic location, a high frequency of the A1 allele was observed in the highlands region. In the same way, a high frequency of heterozygous individuals (A1A2) was observed maybe due to the high presence of the Holstein Friesian breed in this area which is characterized by having this allelic and genotypic profile. However, in the coast region, the predominant allelic and genotypic frequency was A2A2. This may be due to the fact that in that area there was a predominance of Gyr, Jersey, and crossbred individuals, which were characterized, especially the first two, by a high A2A2 allele frequency. Therefore, this difference in allelic and genotypic frequencies between regions could also be due to the unequal distribution of dairy cattle breeds. Thus, in the highlands the presence of Holstein Friesian was higher (23) while on the coast region the presence of Jersey, Gyr, Girolando and crossbreds with Brahman was predominant. As a result, it could be considered that the human population in the coast region has been consuming greater quantities of milk from individuals owning the A2A2 genotype which would be very beneficial to avoid digestibility problems derived from the consumption of milk from individuals with higher A1 allele frequencies that is predominant in the highlands region. Finally, regarding the genetic diversity parameters by region, no differences were observed in terms of genetic distances, so it could be considered that the global population of individuals maintained a balance among them.

5 Conclusion

This study suggests that the allelic and genotypic frequencies of A1/A2 variants of β-casein (CSN2) gene were variable in exotic reared dairy cattle breeds in Ecuador. A relatively high frequency of β-casein A2 allele was observed Gyr and Jersey cattle for A2 milk production since they carry a favorable A2 genotype compared to the other breeds. The polymorphisms of β-casein (CSN2) gene revealed a geographically associated distribution of allelic and genotypic frequencies. The differential frequencies observed between the West (Coast) and the East (Highlands) suggests that the genetic variation may be a result of the introduction of dairy cattle adapted to different climatic conditions as well as geographic artificial selection. Interestingly, an introgression of taurine genes into indicine breeds could have contributed to the increase of frequency of β-casein A1A2 genotype in Gyr which was higher than that obtained in other countries. Finally, the information on A1/A2 variants of β-casein (CSN2) gene could be used in genetic selection programmes for different dairy cattle breeds raised in Ecuador increasing A2 allele which was associated with the production of healthier milk for human consumption.

Data availability statement

The datasets presented in this study can be found in online repositories. The data set that supports the findings of the present research is available online in Zenodo repository (10.5281/zenodo.15114168).

Ethics statement

The study was conducted according to the guidelines of the Declaration of Helsinki and following the Code of Ethics for animal experiments as reflected in the ARRIVE guidelines available at http://www.nc3rs.org.uk/ARRIVEchecklist (Accessed on 10 August 2024). The present research did not require any ethical review and approval by a specific Bioethics Committee for the use of experimental animals, since the samples were obtained for routine diagnostic and health control services in each farm. Therefore, none of the animals were handled or restricted at any time for this study. The studies were conducted in accordance with the local legislation and institutional requirements.

Author contributions

LC-M: Resources, Conceptualization, Methodology, Writing – review & editing, Writing – original draft, Formal analysis, Investigation, Data curation. JN-M: Conceptualization, Methodology, Data curation, Investigation, Writing – original draft, Formal analysis, Writing – review & editing. MG-R: Writing – review & editing, Funding acquisition, Writing – original draft, Resources. MG-H: Methodology, Writing – review & editing, Investigation, Supervision, Writing – original draft, Funding acquisition, Conceptualization, Project administration, Validation.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the Ecuadorian Holstein-Friesian Association (AHFE), as part of the research project entitled “Population structure and genetic diversity of Holstein-Friesian breed in Ecuador” with the agreement PRO-INV-001-2022 (Quito, Ecuador) and Ministerio de Economia y Finanzas/Universidad Técnica de Cotopaxi (UTC)/Proyectos Generativos/Government of Ecuador: DIRGI-CP2022-005.

Acknowledgments

The authors would like to thank the staff involved in the collection of information for the elaboration of the different databases, as well as the stall for the excellent technical assistance.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2025.1616426/full#supplementary-material

References

1. Alfonso, L, Urrutia, O, Mendizabal, JA, et al. Conversión de las explotaciones de vacuno de leche a la producción de leche A2 ante una posible demanda del mercado: posibilidades e implicaciones. Inf Téc Econ Agraria. (2019) 115:231–51. doi: 10.12706/itea.2019.001

Crossref Full Text | Google Scholar

2. Delgado, JV, Martínez, AM, Acosta, A, Álvarez, LA, Armstrong, E, Camacho, E, et al. Genetic characterization of Latin-American creole cattle using microsatellite markers. Anim Genet. (2012) 43:2–10. doi: 10.1111/j.1365-2052.2011.02207.x

PubMed Abstract | Crossref Full Text | Google Scholar

3. Cartuche-Macas, LF, Gutierrez-Reinoso, MA, Chacón, E, Larrea-Izurieta, CO, García-Flores, JM, and Garcia-Herreros, M. Ecuadorian Holstein-Friesian cattle paternal lineages: demographic structure, inbreeding evolution, and genetic diversity. PLoS One. (2025) 20:e0318730. doi: 10.1371/JOURNAL.PONE.0318730

PubMed Abstract | Crossref Full Text | Google Scholar

4. Miglior, F, Fleming, A, Malchiodi, F, Brito, LF, Martin, P, and Baes, CF. A 100-year review: identification and genetic selection of economically important traits in dairy cattle. J Dairy Sci. (2017) 100:10251–71. doi: 10.3168/JDS.2017-12968

PubMed Abstract | Crossref Full Text | Google Scholar

5. Sebastiani, C, Arcangeli, C, Ciullo, M, Torricelli, M, Cinti, G, Fisichella, S, et al. Frequencies evaluation of β-casein gene polymorphisms in dairy cows reared in Central Italy. Animals. (2020) 10:252. doi: 10.3390/ANI10020252

PubMed Abstract | Crossref Full Text | Google Scholar

6. Ayalew, W, Wu, X, Tarekegn, GM, Sisay Tessema, T, Naboulsi, R, Van Damme, R, et al. Whole genome scan uncovers candidate genes related to Milk production traits in Barka cattle. Int J Mol Sci. (2024) 25:6142. doi: 10.3390/IJMS25116142

PubMed Abstract | Crossref Full Text | Google Scholar

7. Gutiérrez-Reinoso, MA, Aponte, PM, Cabezas, J, Rodriguez-Alvarez, L, and Garcia-Herreros, M. Genomic evaluation of Primiparous high-producing dairy cows: inbreeding effects on genotypic and phenotypic production-reproductive traits. Animals. (2020) 10:1–18. doi: 10.3390/ANI10091704

PubMed Abstract | Crossref Full Text | Google Scholar

8. Jiang, J, Gao, Y, Hou, Y, Li, W, Zhang, S, Zhang, Q, et al. Whole-genome resequencing of Holstein bulls for Indel discovery and identification of genes associated with Milk composition traits in dairy cattle. PLoS One. (2016) 11:e0168946. doi: 10.1371/JOURNAL.PONE.0168946

PubMed Abstract | Crossref Full Text | Google Scholar

9. Kulibaba, R, Sakhatskyi, M, and Liashenko, Y. Comparative analysis of A1 and A2 allele detection efficiency for bovine CSN2 gene by AS-PCR methods. Acta Biochim Pol. (2023) 70:205–9. doi: 10.18388/ABP.2020_6530

PubMed Abstract | Crossref Full Text | Google Scholar

10. Cieślińska, A, Fiedorowicz, E, Rozmus, D, Sienkiewicz-Szłapka, E, Jarmołowska, B, and Kamiński, S. Does a little difference make a big difference? Bovine β-casein A1 and A2 variants and human health-an update. J Mol Sci. (2022) 23:15687. doi: 10.3390/ijms232415637

PubMed Abstract | Crossref Full Text | Google Scholar

11. Ladyka, V, Skliarenko, Y, Pavlenko, Y, and Malikova, A. Study of beta-casein gene polymorphism in dairy cattle populations of Ukraine In: Management of the genetic fund of animals – problems, solutions, outlooks (2023). 156–61. doi: 10.61562/mgfa2023.20

Crossref Full Text | Google Scholar

12. Sanchez, MP, Fritz, S, Patry, C, Delacroix-Buchet, A, and Boichard, D. Frequencies of milk protein variants and haplotypes estimated from genotypes of more than 1 million bulls and cows of 12 French cattle breeds. J Dairy Sci. (2020) 103:9124–41. doi: 10.3168/JDS.2020-18492

PubMed Abstract | Crossref Full Text | Google Scholar

13. Chia, JSJ, McRae, JL, Kukuljan, S, Woodford, K, Elliott, RB, Swinburn, B, et al. A1 beta-casein milk protein and other environmental pre-disposing factors for type 1 diabetes. Nutr Diabetes. (2017) 7:e274. doi: 10.1038/NUTD.2017.16

PubMed Abstract | Crossref Full Text | Google Scholar

14. de Vasconcelos, ML, Oliveira, LMFS, Hill, JP, and Vidal, AMC. Difficulties in establishing the adverse effects of β-Casomorphin-7 released from β-casein variants-a review. Foods. (2023) 12:3151. doi: 10.3390/FOODS12173151

PubMed Abstract | Crossref Full Text | Google Scholar

15. Daniloski, D, McCarthy, NA, Huppertz, T, and Vasiljevic, T. What is the impact of amino acid mutations in the primary structure of caseins on the composition and functionality of milk and dairy products? Curr Res Food Sci. (2022) 5:1701–12. doi: 10.1016/J.CRFS.2022.09.026

PubMed Abstract | Crossref Full Text | Google Scholar

16. Rangel, AHN, Zaros, LG, Lima, TC, Borba, LHF, Novaes, LP, Mota, LFM, et al. Polymorphism in the Beta casein gene and analysis of milk characteristicsin Gir and Guzerá dairy cattle. Genet Mol Res. (2017) 16. doi: 10.4238/GMR16029592

PubMed Abstract | Crossref Full Text | Google Scholar

17. Şahin, Ö, and Boztepe, S. Assessment of A1 and A2 variants in the CNS2 gene of some cattle breeds by using ACRS-PCR method. Anim Biotechnol. (2023) 34:1505–13. doi: 10.1080/10495398.2022.2036176

PubMed Abstract | Crossref Full Text | Google Scholar

18. Crowe, AD, Lonergan, P, and Butler, ST. Invited review: use of assisted reproduction techniques to accelerate genetic gain and increase value of beef production in dairy herds. J Dairy Sci. (2021) 104:12189–206. doi: 10.3168/JDS.2021-20281

PubMed Abstract | Crossref Full Text | Google Scholar

19. Kay, SIS, Delgado, S, Mittal, J, Eshraghi, RS, Mittal, R, and Eshraghi, AA. Beneficial effects of Milk having A2 β-casein protein: myth or reality? J Nutr. (2021) 151:1061–72. doi: 10.1093/JN/NXAA454

PubMed Abstract | Crossref Full Text | Google Scholar

20. Bisutti, V, Pegolo, S, Giannuzzi, D, Mota, LFM, Vanzin, A, Toscano, A, et al. The β-casein (CSN2) A2 allelic variant alters milk protein profile and slightly worsens coagulation properties in Holstein cows. J Dairy Sci. (2022) 105:3794–809. doi: 10.3168/jds.2021-21537

PubMed Abstract | Crossref Full Text | Google Scholar

21. Vigolo, V, Visentin, E, Ballancin, E, Lopez-Villalobos, N, Penasa, M, and De Marchi, M. β-Casein A1 and A2: effects of polymorphism on the cheese-making process. J Dairy Sci. (2023) 106:5276–87. doi: 10.3168/jds.2022-23072

PubMed Abstract | Crossref Full Text | Google Scholar

22. Dai, R, Fang, Y, Zhao, W, Liu, S, Ding, J, Xu, K, et al. Identification of alleles and genotypes of beta-casein with DNA sequencing analysis in Chinese Holstein cow. J Dairy Res. (2016) 83:312–6. doi: 10.1017/S0022029916000303

PubMed Abstract | Crossref Full Text | Google Scholar

23. ESPAC (2022). Encuesta de Superficie y Producción Agropecuaria Continua-2022. Instituto Nacional de Estadísticas y Censos. Available online at: https://www.ecuadorencifras.gob.ec/encuesta-de-superficie-y-produccion-agropecuaria-continua-2022/

Google Scholar

24. Ministerio de Agricultura y Ganadería (2022). Informe de Informe de Gestión. Available online at: https://empresa-de-servicios-publicos-de-aguazul.micolombiadigital.gov.co/sites/empresa-de-servicios-publicos-de-aguazul/content/files/000122/6062_informe-de-gestion-20182019-vfinal.pdf (Accessed April 14, 2024).

Google Scholar

25. Caroli, AM, Savino, S, Bulgari, O, and Monti, E. Detecting β-casein variation in bovine Milk. Molecules. (2016) 21:141. doi: 10.3390/MOLECULES21020141

PubMed Abstract | Crossref Full Text | Google Scholar

26. Gutiérrez-Reinoso, MA, Aponte, PM, and García-Herreros, M. Genomic and phenotypic udder evaluation for dairy cattle selection: a review. Animals. (2023) 13:1588. doi: 10.3390/ANI13101588

PubMed Abstract | Crossref Full Text | Google Scholar

27. Amalfitano, N, Macedo Mota, LF, Rosa, GJM, Cecchinato, A, and Bittante, G. Role of CSN2, CSN3, and BLG genes and the polygenic background in the cattle milk protein profile. J Dairy Sci. (2022) 105:6001–20. doi: 10.3168/JDS.2021-21421

PubMed Abstract | Crossref Full Text | Google Scholar

28. Antonopoulos, D, Vougiouklaki, D, Laliotis, GP, Tsironi, T, Valasi, I, Chatzilazarou, A, et al. Identification of polymorphisms of the CSN2 gene encoding β-casein in Greek local breeds of cattle. Vet Sci. (2021) 8:257. doi: 10.3390/VETSCI8110257

Crossref Full Text | Google Scholar

29. Gutierrez-Reinoso, MA, Aponte, PM, and Garcia-Herreros, M. Genomic analysis, Progress and future perspectives in dairy cattle selection: a review. Animals. (2021) 11:599. doi: 10.3390/ANI11030599

PubMed Abstract | Crossref Full Text | Google Scholar

30. Roman, J. (1970). Genetics of milk production in Ecuador. Available online at: https://ufdcimages.uflib.ufl.edu/AA/00/05/85/74/00001/AA00058574_00001.pdf (Accessed June 30, 2024).

Google Scholar

31. Hartl, D.L., and Clark, A.G. (1997). Principles of population Genetics, (3rd). Sunderland, Massachusetts: Sinauer associates, Inc., Oxford University Press. 519. Available online at: https://global.oup.com/academic/product/principles-of-population-genetics-9780878933082?cc=us&lang=en& (Accessed September 23, 2024).

Google Scholar

32. Hedrick, P. W. (2011). Genetics of populations, 4th. Jones and Bartlett Publishers. Available online at: https://www.jblearning.com/catalog/productdetails/9780763757373?srsltid=AfmBOoqjf8Ew_NHPmCGFLIu6fvT4CGJZfSKQHE4D9k4a5pMit-1jzda4

Google Scholar

33. Lynch, M, and Milligan, BG. Analysis of population genetic structure with RAPD markers. Mol Ecol. (1994) 3:91–9. doi: 10.1111/J.1365-294X.1994.TB00109.X

PubMed Abstract | Crossref Full Text | Google Scholar

34. Brown, AHD, and Weir, BS. Measuring genetic variability in plant populations In: SD Tanksley and TJ Orton, editors. Developments in Plant Genetics and Breeding (1983). 219–39. doi: 10.1016/b978-0-323-90583-1.00016-7

Crossref Full Text | Google Scholar

35. Peakall, R, and Smouse, PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research—an update. Bioinformatics. (2012) 28:2537–9. doi: 10.1093/BIOINFORMATICS/BTS460

PubMed Abstract | Crossref Full Text | Google Scholar

36. Slatkin, M. A measure of population subdivision based on microsatellite allele frequencies. Genetics. (1995) 139:457–62. doi: 10.1093/GENETICS/139.1.457

PubMed Abstract | Crossref Full Text | Google Scholar

37. Bodnár, Á, Hajzser, A, Egerszegi, I, Póti, P, Kuchtík, J, and Pajor, F. A2 milk and its importance in dairy production and global market. Anim Welf Ethol Husb Technol. (2018) 14:1–7. doi: 10.17205/SZIE.AWETH.2018.1.001

Crossref Full Text | Google Scholar

38. Żbik, K, Onopiuk, A, Górska-Horczyczak, E, and Wierzbicka, A. Trends and opportunities in the dairy industry: A2 Milk and processing methods. Appl Sci. (2024) 14:6513. doi: 10.3390/APP14156513

Crossref Full Text | Google Scholar

39. Arens, SC, Sharpe, KT, Schutz, MM, Hardie, LC, Dechow, CC, and Heins, BJ. Relationships of beta-casein genetics with production, fertility, and survival of purebred organic Holstein dairy cows. JDS Commun. (2023) 4:458–63. doi: 10.3168/jdsc.2022-0367

PubMed Abstract | Crossref Full Text | Google Scholar

40. Khan, R, De, S, Dewangan, R, Tamboli, R, and Gupta, R. Potential status of A1 and A2 variants of bovine beta-casein gene in milk samples of Indian cattle breeds. Anim Biotechnol. (2023) 34:4878–84. doi: 10.1080/10495398.2023.2200502

PubMed Abstract | Crossref Full Text | Google Scholar

41. Kovalyuk, N., Shiryaeva, E., Yakusheva, L., Kaishev, V., and Aisanov, Z. (2021). Frequency distribution of polymorphic variants of the beta-casein locus in Russian various breed cattle subpopulations. E3S Web of Conferences 262, 02020.

Google Scholar

42. Nuomin,, Nguyen, QD, Aodaohu,, and Nishino, N. Frequency of β-casein gene polymorphisms in Jersey cows in Western Japan. Animals. (2022) 12:2076. doi: 10.3390/ANI12162076

PubMed Abstract | Crossref Full Text | Google Scholar

43. Potočnik, K, Lusčtrek, B, and Kaić, A. Does the selection on SS-casein affect the traits important for dairy production of Slovenian Brown Swiss cattle? Acta Agric Slov. (2016) 5:89–93. doi: 10.14720/aas-s.2016.5.18873

Crossref Full Text | Google Scholar

44. Viana, J. H. M. (2023). Statistics of embryo production and transfer in domestic farm animals. International Embryo Transfer Society (IETS) 41, 1–25.

Google Scholar

45. Chacón Marcheco, E, Cartuche Macas, LF, Villavicencio Estrella, AN, Toro Molina, BM, Silva Déley, LM, Andrade Aulestia, PM, et al. Diversidad genética de la población de toros Holstein Friesian importados al Ecuador. Rev Cient Tecnol. (2023) 10:33–40. doi: 10.26423/RCTU.V10I1.717

Crossref Full Text | Google Scholar

46. Gustavsson, F, Buitenhuis, AJ, Johansson, M, Bertelsen, HP, Glantz, M, Poulsen, NA, et al. Effects of breed and casein genetic variants on protein profile in milk from Swedish red, Danish Holstein, and Danish Jersey cows. J Dairy Sci. (2014) 97:3866–77. doi: 10.3168/JDS.2013-7312

PubMed Abstract | Crossref Full Text | Google Scholar

47. Truswell, AS. The A2 milk case: a critical review. Eur J Clin Nutr. (2005) 59:623–31. doi: 10.1038/sj.ejcn.1602104

PubMed Abstract | Crossref Full Text | Google Scholar

48. Pessoa, DBM, Mangabeira, I, Borba, L, Madruga, R, Ribeiro, C, Pereira, M, et al. Kappa and beta-casein alleles in dairy zebu cattle. Anim Biotechnol. (2023) 34:3162–4. doi: 10.1080/10495398.2022.2137677

PubMed Abstract | Crossref Full Text | Google Scholar

49. Ayaz,, Suhail, SM, Ahmad, I, Zeb, MT, Khan, R, Ijaz, A, et al. Detection of A2A2 genotype of beta casein protein (CSN2) gene in local, exotic and cross bred cattle in Pakistan. Anim Biotechnol. (2023) 34:1462–73. doi: 10.1080/10495398.2022.2031204

PubMed Abstract | Crossref Full Text | Google Scholar

50. Aytekin, İ. (2024). Β-casein (CSN2) polymorphism in female Jersey calves in a commercial dairy cattle farm, in AFRICA 5th international conference on new horizons in sciences April 25–28, 2024 – CAIRO. Available online at: https://www.researchgate.net/publication/382452771

Google Scholar

51. Braunvieh (2024). Brown Swiss. Available online at: https://homepage.braunvieh.ch/en/brown-swiss/ (Accessed April 16, 2025).

Google Scholar

52. Zúñiga, D. (2022). Determinación de la frecuencia de los alelos A1 y A2 del gen CSN2 en bovinos de raza Holstein y bovinos criollos del departamento de Apurímac. Available online at: https://repositorio.ucsm.edu.pe/handle/20.500.12920/12076 (Accessed September 21, 2024).

Google Scholar

53. Ruiz Alvarez, K. P., and Salmerón Jiménez, C. E. (2024). Evaluación genética de los alelos A1 y A2 de β-caseína en 5 razas bovinas. Available online at: https://dspace.ucuenca.edu.ec/handle/123456789/43878 (Accessed July 24, 2024).

Google Scholar

54. Jann, OC, Ibeagha-Awemu, EM, Özbeyaz, C, Zaragoza, P, Williams, JL, Ajmone-Marsan, P, et al. Geographic distribution of haplotype diversity at the bovine casein locus. Genet Sel Evol. (2004) 36:243–57. doi: 10.1186/1297-9686-36-2-243

PubMed Abstract | Crossref Full Text | Google Scholar

55. Sodhi, M, Mukesh, M, Sharma, V, Kataria, RS, and Sobti, RC. Harnessing potential of A2 milk in India: an overview In: Advances in animal experimentation and modeling: understanding life phenomena (2022). 379–92.

Google Scholar

56. Vieira, MT, dos Santos Daltro, D, and Cobuci, JA. Breed and heterosis effects on reproduction and production traits of Girolando cows. Rev Bras Zootec. (2022) 51:e20200266. doi: 10.37496/RBZ5120200266

Crossref Full Text | Google Scholar

57. Scott, BA, Haile-Mariam, M, MacLeod, IM, Xiang, R, and Pryce, JE. Evaluating the potential impact of selection for the A2 milk allele on inbreeding and performance in Australian Holstein cattle. Front Anim Sci. (2023) 4:1142673. doi: 10.3389/FANIM.2023.1142673/BIBTEX

Crossref Full Text | Google Scholar

58. Harris, AM, and Degiorgio, M. An unbiased estimator of gene diversity with improved variance for samples containing related and inbred individuals of any ploidy. Genes Genomes Genetics. (2017) 7:671–91. doi: 10.1534/g3.116.037168

PubMed Abstract | Crossref Full Text | Google Scholar

59. Kelleher, MM, Berry, DP, Kearney, JF, McParland, S, Buckley, F, and Purfield, DC. Inference of population structure of purebred dairy and beef cattle using high-density genotype data. Animal. (2017) 11:15–23. doi: 10.1017/S1751731116001099

PubMed Abstract | Crossref Full Text | Google Scholar

60. Barłowska, J., Sawicka-Zugaj, W., Ślaska, B., Król, J., Brodziak, A., Teter, A., et al. (2022). Genetic analysis of CSN2 in local and international cattle breeds raised in Poland. Anim Sci Paper Rep 40, 411–422. Available online at: https://open.up.lublin.pl/Uczelnia/Publiczna/Wyszukiwarka/PublikacjaDetail/157053?profilId=655

Google Scholar

61. Zepeda-Batista, JL, Alarcón-Zúñiga, B, Ruíz-Flores, A, Núñez-Domínguez, R, and Ramírez-Valverde, R. Polymorphism of three milk protein genes in Mexican Jersey cattle. Electron J Biotechnol. (2015) 18:1–4. doi: 10.1016/J.EJBT.2014.10.002

Crossref Full Text | Google Scholar

62. Carvajal, AM, Levicoy, D, and Pizarro-Aránguiz, N. Genotyping of the beta-casein A1 and A2 variants in chilean dairy cattle. Anim Sci Genet. (2023) 19:95–101. doi: 10.5604/01.3001.0053.7705

Crossref Full Text | Google Scholar

63. Lachance, J. Hardy-Weinberg equilibrium and random mating In: RM Kliman, Editor. Encyclopedia of Evolutionary Biology. Academic Press. (2016). 208–11.

Google Scholar