中国黄牛Y-STRs遗传多样性与起源研究

[目的]研究中国黄牛Y染色体STRs的遗传多样性及父系起源。[方法]利用非变性聚丙烯酰胺凝胶电泳,选择2个牛Y-STRs位点INRA189和BM861,分析16个中国地方黄牛品种284头公牛与4头缅甸黄牛公牛的Y染色体遗传多样性。[结果]在中国16个黄牛品种中,2个Y-STR位点可以区分中国黄牛中的普通牛和瘤牛类型,表明中国黄牛有普通牛和瘤牛两种父系起源。4头缅甸黄牛均为瘤牛类型。在中国16个黄牛品种中,普通牛和瘤牛分布频率分别为57.0%和43.0%,其中普通牛频率在北方黄牛中占优势(98.3%),瘤牛频率在南方黄牛中占优势(76.1%),中原黄牛中普通牛频率较高为63.8%,瘤牛频率为36.2%。[结论]中国黄牛存在普通牛和瘤牛两种父系起源;普通牛频率自北向南逐渐减少,瘤牛频率自北向南逐渐增加,中原地区为普通牛和瘤牛的交汇处。     

Genetic structure and relationships of 16 Asian and European cattle populations using DigiTag2 assay

In this study, we genotyped 117 autosomal single nucleotide polymorphisms using a DigiTag2 assay to assess the genetic diversity, structure and relationships of 16 Eurasian cattle populations, including nine cattle breeds and seven native cattle. Phylogenetic and principal component analyses showed that Bos taurus and Bos indicus populations were clearly distinguished, whereas Japanese Shorthorn and Japanese Polled clustered with European populations. Furthermore, STRUCTURE analysis demonstrated the distinct separation between Bos taurus and Bos indicus (K=2), and between European and Asian populations (K=3). In addition, Japanese Holstein exhibited an admixture pattern with Asian and European cattle (K=3‐5). Mongolian (K=13‐16) and Japanese Black (K=14‐16) populations exhibited admixture patterns with different ancestries. Bos indicus populations exhibited a uniform genetic structure at K=2‐11, thereby suggesting that there are close genetic relationships among Bos indicus populations. However, the Bhutan and Bangladesh populations formed a cluster distinct from the other Bos indicus populations at K=12‐16. In conclusion, our study could sufficiently explain the genetic construction of Asian cattle populations, including: (i) the close genetic relationships among Bos indicus populations; (ii) the genetic influences of European breeds on Japanese breeds; (iii) the genetic admixture in Japanese Holstein, Mongolian and Japanese Black cattle; and (iv) the genetic subpopulations in Southeast Asia.     

Promoter region of the bovine growth hormone receptor gene: single nucleotide polymorphism discovery in cattle and association with performance in Brangus bulls

Expression of the GH receptor (GHR) gene and its binding with GH is essential for growth and fat metabolism. A GT microsatellite exists in the promoter of bovine GHR segregating short (11 bp) and long (16 to 20 bp) allele sequences. To detect SNP and complete an association study of genotype to phenotype, we resequenced a 1,195-bp fragment of DNA including the GT microsatellite and exon 1A. Resequencing was completed in 48 familialy unrelated Holstein, Jersey, Brown Swiss, Simmental, Angus, Brahman, and Brangus cattle. Nine SNP were identified. Phylogeny analyses revealed minor distance (i.e., <5%) in DNA sequence among the 5 Bos taurus breeds; however, sequence from Brahman cattle averaged 27.4 +/- 0.07% divergence from the Bos taurus breeds, whereas divergence of Brangus was intermediate. An association study of genotype to phenotype was completed with data from growing Brangus bulls (n = 553 from 96 sires) and data from 4 of the SNP flanking the GT microsatellite. These SNP were found to be in Hardy-Weinberg equilibrium and in phase based on linkage disequilibrium analyses (r(2) = 0.84 and D'= 0.92). An A/G tag SNP was identified (ss86273136) and was located in exon 1A, which began 88 bp downstream from the GT microsatellite. Minor allele frequency of the tag SNP was greater than 10%, and Mendelian segregation was verified in 3 generation pedigrees. The A allele was derived from Brahman, and the G allele was derived from Angus. This tag SNP genotype was a significant effect in analyses of rib fat data collected with ultrasound when bulls were ~365 d of age. Specifically, bulls of the GG genotype had 6.1% more (P = 0.0204) rib fat than bulls of the AA and AG genotypes, respectively. Tag SNP (ss86273136), located in the promoter of GHR, appears to be associated with a measure of corporal fat in Bos taurus x Bos indicus composite cattle.     

Allele frequencies of gene polymorphisms related to economic traits in Bos taurus and Bos indicus cattle breeds

Allele frequencies of 10 representative polymorphisms for beef and milk traits were investigated for a total of 240 animals from Bos taurus and Bos indicus breeds, including two Japanese groups (Japanese Black and Japanese Brown), two East Asian groups (Korean and Mongolian), three European groups (Holstein, Angus and Hereford) and a Bos indicus group in South Asia (Myanmar, Laos and Cambodia). The Japanese Black revealed unique genetic construction in GH, FASN and SREBP‐1 and the other Asian populations show intermediate frequencies between European and Japanese populations. The Bos indicus group showed low favorable allele frequencies in most of the genes. The study showed the variability and distribution of 10 genes affecting economic traits among world representative cattle breeds. The genetic information would contribute to elucidating the genetic background for worldwide cattle breeds and the possibility of improvement using the markers.     

中国黄牛Y-SNPs遗传多样性与起源研究

为了研究中国黄牛Y染色体SNPs的遗传多样性及父系起源,本研究利用PCR-SSCP与测序方法,选择4个牛Y-SNPs位点DDX3Y-7、UTY-19、ZFY-9和ZFY-10,分析了16个中国地方黄牛品种284头公牛与缅甸黄牛4头公牛Y染色体的遗传多样性.结果表明,在中国16个黄牛品种中,仅发现普通牛Y2和瘤牛Y3单倍型,表明只有Y2和Y3两种父系起源,尚未发现中国黄牛存在普通牛Y1单倍型的分子证据.4头缅甸黄牛均为Y3单倍型.在中国16个黄牛品种中,Y2和Y3单倍型频率分别为57.0％和43.0％,其中Y2单倍型频率在北方黄牛中占优势(98.3％),Y3单倍型频率在南方黄牛中占优势(76.1％),中原黄牛中普通牛Y2的单倍型频率较高,为63.8％0,瘤牛Y3的单倍型频率为36.2％.本研究证明,中国黄牛存在普通牛Y2和瘤牛Y3单倍型两种父系起源,Y2单倍型频率自北向南逐渐减少,Y3单倍型频率自北向南逐渐增加,中原地区为普通牛Y2和瘤牛Y3单倍型的交汇处.     

Frequency of bovine Nramp1 (Slc11a1) alleles in Holstein and Zebu breeds

Natural resistance against brucellosis in cattle is linked to the Nramp1 gene, which encodes a divalent cation transporter that localizes in the phagolysosome membrane in macrophages. Nramp1 gene in mouse plays a critical role in innate immunity favoring bacterial killing by macrophages in addition to its influence on adaptative immunity. Polymorphisms at the bovine Nramp1 3' untranslated region (3'UTR), detectable by Single Strand Conformational Analysis (SSCA), are associated with natural resistance against brucellosis. Such polymorphisms are associated with variation in the number of GT repeats. This study compared the frequency of Nramp1 3'UTR polymorphisms between Zebu and European bovine breeds. Eighty-one Holsteins (Bos taurus taurus) and 167 Zebu (Bos taurus indicus), including the following breeds: Nelore (n=95), Guzerá (n=37), and Gir (n=35), totaling 248 pure breed cattle studied. DNA extraction was performed using the guanidium protocol and genotyping was performed by SSCA. DNA from cattle considered genotypically resistant to brucellosis resulted in a single band (homozygous) with 175bp, corresponding to the 3'UTR with 13 GT pairs (GT13), whereas DNA from genotypically susceptible cattle generated one single band with 177bp (homozygous GT14) or double bands with both 175 and 177bp, or 175 and 179bp (heterozygous GT13/GT14 or GT13/GT15, respectively). A marked difference in the frequency of alleles was detected between the Zebu and Holstein cattle. Holsteins had an extremely homogeneous genotype, with 100% of the individuals with a GT13 genotype. In sharp contrast the Nelore breed had the most heterogeneous genotype with four allelic combinations, namely, homozygous GT13, homozygous GT14, heterozygous GT13/GT14, and heterozygous GT13/GT15. When the Zebu breeds were compared to each other, the only significant difference observed was the frequencies of the genotypes GT13 and GT14 between the Nelore and Guzerá breeds. The knowledge of allelic frequencies in different breeds of cattle may prove to be very useful in the future for planning breeding strategies for selection of resistant cattle.     

Bos taurus Y chromosome of Africander cattle and the development of improved breeds for the tropics.

The Africander has anatomical and other characteristics of an animal of approximately 3/4 Bos indicus heredity. The fact that these cattle carry the Bos taurus Y chromosome supports this view and indicates that the local cattle in South Africa would have been crossed with one or more Bos taurus bulls. Droughtmaster and Braford cattle retain the Bos indicus Y chromosome because Bos indicus instead of Bos taurus bulls were used to establish these taurindicus breeds. Contrary to some assumptions, an approximate tenfold increase in productivity of cattle was made during the 18th and 19th centuries due to improvements in disease control, nutrition and genetic improvement. What is now needed is the development of taurindicus breeds combining to the maximum possible extent the disease resistance and hardiness of Bos indicus with the early maturity and productivity of Bos taurus cattle. In addition, the ravages of disease and the seasonal variations of food supply need to be overcome in tropical areas.     

Bos indicus type of growth hormone receptor gene is retained in Japanese Black cattle

The growth hormone receptor (GHR) gene is responsible for growth and carcass traits, and polymorphisms associated with the variation of meat production are thought to occur in the liver‐specific promoter of the GHR gene in cattle. The aim of this study was to analyse the structure of the liver‐specific promoter of GHR in Japanese Black cattle, as the relationship between GHR polymorphism and meat production is poorly understood in this breed. Typically in European cattle, the LINE‐1 element, a family of retrotransposons, is inserted in the liver‐specific promoter. However, a short GHR promoter without the LINE‐1 sequence was found in the Japanese Black breed as in Bos indicus cattle. The frequency of the short allele was approximately 60%. In addition, 24 of 29 Holstein/Japanese Black crosses carried the short allele from their sire. The present result suggests that the short allele for GHR may be a candidate marker for improving meat production of Japanese Black cattle.     

Genetic diversity and structure in Bos taurus and Bos indicus populations analyzed by SNP markers

The purpose of this study was to assess genetic diversity, phylogenetic relationship and population structure among nine Eurasian cattle populations using 58 single nucleotide polymorphism (SNP) markers. The calculated distribution of minor allele frequencies and heterozygosities suggested that the genetic diversity of Bos indicus populations was lower than that of Bos taurus populations. Phylogenetic analyses revealed the main divergence between the Bos taurus and Bos indicus populations, and subsequently between Asian and European populations. By principal components analysis, the Bos taurus and Bos indicus populations were clearly distinguished with PC1 (61.1%); however, six Bos taurus populations clustered loosely and the partial separation between European and Asian groups was observed by PC2 (12.5%). The structure analysis was performed using the STRUCTURE program. Distinct separation between Bos taurus and Bos indicus was shown at K = 2, and that between European and Asian populations at K = 3. At K = 4, 5 and 6, Mongolian population showed an admixture pattern with different ancestry of Asian and European cattle. At K = 7, all Bos taurus populations showed each cluster with little proportion of admixture. In conclusion, 58 SNP markers in this study could sufficiently estimate the genetic diversity, relationship and structure for nine Eurasian cattle populations, especially by analyses of principal components and STRUCTURE.     

History and economic importance of cattle (Bos taurus L.) in Switzerland from Neolithic to Early Middle Ages

In Switzerland domestic cattle (Bos primigenius f. taurus resp. Bos taurus L.) first appear with the earliest Neolithic settlements (approximately 5000 BC). With the gradual deforestation of the landscape caused by human exploitation of the environment, cattle were used more intensive and in many ways. There is evidence that cattle were used as draught animal since ca. 3400 BC, probably even earlier milk was regularly used. The size of domestic cattle gradually decreased from Early Neolithic until Iron Age. Only with Roman influence larger animals are found. However, after the withdrawal of Romans the average size of cattle decreased again. Archaeogenetic studies will have to show, whether this is due to novel breeding strategies or the import of breeding stock. First genetic results showed that a female genetic type, which is rare in European breeds, is present in Swiss Evolène cattle and in one animal of Roman time cattle from Augusta Raurica. Is this a sign for influence of Roman cattle on today's Swiss breeds?     

A new single nucleotide polymorphism in CAPN1 extends the current tenderness marker test to include cattle of Bos indicus, Bos taurus, and crossbred descent

The three objectives of this study were to 1) test for the existence of beef tenderness markers in the CAPN1 gene segregating in Brahman cattle; 2) test existing CAPN1 tenderness markers in indicus-influenced crossbred cattle; and 3) produce a revised marker system for use in cattle of all subspecies backgrounds. Previously, two SNP in the CAPN1 gene have been described that could be used to guide selection in Bos taurus cattle (designated Markers 316 and 530), but neither marker segregates at high frequency in Brahman cattle. In this study, we examined three additional SNP in CAPN1 to determine whether variation in this gene could be associated with tenderness in a large, multisire American Brahman population. One marker (termed 4751) was associated with shear force on postmortem d 7 (P < 0.01), 14 (P = 0.015), and 21 (P < 0.001) in this population, demonstrating that genetic variation important for tenderness segregates in Bos indicus cattle at or near CAPN1. Marker 4751 also was associated with shear force (P < 0.01) in the same large, multisire population of cattle of strictly Bos taurus descent that was used to develop the previously reported SNP (referred to as the Germplasm Evaluation [GPE] Cycle 7 population), indicating the possibility that one marker could have wide applicability in cattle of all subspecies backgrounds. To test this hypothesis, Marker 4751 was tested in a third large, multisire cattle population of crossbred subspecies descent (including sire breeds of Brangus, Beefmaster, Bonsmara, Romosinuano, Hereford, and Angus referred to as the GPE Cycle 8 population). The highly significant association of Marker 4751 with shear force in this population (P < 0.001) confirms the usefulness of Marker 4751 in cattle of all subspecies backgrounds, including Bos taurus, Bos indicus, and crossbred descent. This wide applicability adds substantial value over previously released Markers 316 and 530. However, Marker 316, which had previously been shown to be associated with tenderness in the GPE Cycle 7 population, also was highly associated with shear force in the GPE Cycle 8 animals (P < 0.001). Thus, Marker 316 may continue to be useful in a variety of populations with a high percentage of Bos taurus backgrounds. An optimal marker strategy for CAPN1 in many cases will be to use both Markers 316 and 4751.     

Whole genome analyses revealed genomic difference between European taurine and East Asian taurine

European taurine and East Asian taurine are two main clades in Bos taurus, but their genomic differences are not clearly elucidated. Here, we sequenced 16 Mongolian cattle genomes and compared them to the 92 genomes of 10 representative breeds worldwide. We found the highest LD level in Mishima cattle and the fastest LD decay in European taurine. Phylogenetic analysis revealed that Mongolian, Hanwoo and Mishima cattle were clustered into East Asian taurine. From selective sweep, gene annotation, functional enrichment and differential expression analysis, we identified selective signals including genes and/or pathways related to rapid growth and large body size in European taurine, and superior meat quality in East Asian taurine. Our findings will help us understand the evolutionary history and formation process of the breeds and provide theoretical materials regarding the genetic mechanism underlying breed characteristics and molecular breeding programmes of the taurine clades in the future.     

温岭高峰牛线粒体DNA全基因组遗传多样性分析

[目的]通过测定温岭高峰牛线粒体DNA全基因组序列以分析温岭高峰牛的母系起源及遗传多样性。[方法]采用DNA提取、测序及生物信息学方法。[结果]通过对19头温岭高峰牛线粒体DNA全基因组序列分析,共发现263个变异位点,定义9种单倍型,单倍型多样度(Hd±SD)为0.778±0.096,核苷酸多样度(Pi±SD)为0.0017±0.0014,表明温岭高峰牛的遗传多样性较低。构建的NJ系统发育树和单倍型进化网络图表明温岭高峰牛有普通牛和瘤牛2种母系起源。[结论]温岭高峰牛线粒体DNA基因组的遗传多样性较低,有瘤牛和普通牛两个母系起源,但主要受瘤牛的影响。     

C- AND G- BANDING PATTERNS AND CHROMOSOMAL MORPHOLOGY OF SOME BREEDS OF AUSTRALIAN CATTLE

A cytogenetical study using metaphase chromosomes from cultured lymphocytes, was made of 2 Banteng (Bibos banteng) steers and 218 bulls representing 13 purebreeds (Bos taurus type, Bos indicus type and Sanga) and 7 cross-breeds. Studies were made of photographic karyotypes of Giemsa stained and C-banded chromosomes of bulls of each breed and of B-banded chromosomes from 3 breeds of Bos indicus and one cross-breed Australian Friesian Sahiwal) cattle. The relative lengths of chromosomes of Bos taurus and Bos indicus bulls were compared and significant difference in relative lengths of the X chromosomes were noted between these two species. There was a differences in morphology of the Y chromosomes; Sanga, Banteng and Bos taurus type breeds had a small submetacentric Y chromosome, except for the Jersey which had a metacentric Y chromosome. All Bos indicus type bulls had an acrocentric Y chromosome but the Droughtmaster breed had two forms of the Y chromosome (submetacentric and acrocentric). The C-banding patterns of the autosomes and X chromosomes were similar for all breeds while those of the Y chromosomes of Bos indicus type cattle allowed their accurate identification. G-banding patterns of Bos indicus resembled those of Bos taurus and enabled pairing of homologous chromosomes. Centromeres of the autosomes were unstained but those of the sex chromosomes were darkly stained.     

Analysis of a polymorphism in the DGAT1 gene in 14 cattle breeds through PCR-SSCP methods

The diacylglycerol O-acyltransferase (DGAT1) is a microsomal enzyme that catalyzes the final step of triglyceride synthesis. Recent work have evidenced a significant association between lysine at amino acid position 232 with elevated milk fat content, while an alanine at this position is associated with lowered milk fat content. The aim of the present work was to develop a simple and inexpensive PCR-SSCP assay in order to discriminate the CG/AA alleles in exon 8 of the DGAT1 gene. In addition, this method was used to analyze the polymorphism of the DGAT1 through PCR-SSCP methods in 14 populations of cattle from Argentine, Bolivia and Uruguay. The PCR primers were designed from GenBank reported sequences. In this study, we found three PCR-SSCP variants, which were denominated from "A" to "C". However, DNA sequencing analysis showed that "A" variant corresponded with the A allele, while both "B" and "C" observed pattern have the motif AA at positions 10,433-10,434 (K allele), being two alternative conformations of the same DNA sequence. Both variants were detected within each breed with the exception of Hereford, and the heterozygosity varied between 0.000 and 0.524. The gene frequency analysis evidenced significant differences among the studied breeds (F(ST) = 0.325, p = 0.000). European Bos taurus breeds, with the exception of Jersey breed, showed the lowest frequency of the K allele, while highest K allele frequencies were harboured by Bos indicus type cattle. In addition, unselected South American Creole cattle breeds and the synthetic Brangus breed had intermediate allele frequencies.     

Susceptibility of Bos indicus and Bos taurus to Anaplasma marginale and Babesia bigemina infections

The reaction of Bos taurus and pure-bred Bos indicus heifers to infection with the intraerythrocytic parasites Anaplasma marginale and Babesia bigemina was studied. B. bigemina infection at 18 months and A. marginale infection at 13 or 24 months resulted in slightly less severe reactions in pure-bred Bos indicus cattle than in Bos taurus. In both breeds, the reaction to A. marginale infection was more severe in older cattle. The severity of B. bigemina infection was not affected by a previous infection with A. marginale.     

The major histocompatibility complex of West African cattle. Typing of lymphocyte antigens (BoLA) of Baoulé cattle (Bos taurus) and Sudanese zebus (Bos indicus) of Burkina Faso (West Africa)

Lymphocyte antigen (BoLA) typing on 247 Baoule cattle (Bos taurus) and 106 Sudanese zebus (Bos indicus), allowed us to determine gene frequencies of 43 class 1 specificities, as international "W" and European "EU", as African local from Kenya "KN" and Burkina Faso "BF". In comparing these frequencies, it appears that some specificities could be considered as significant breed markers for Baoule cattle and zebus.     

Characterization and validation of bovine gonadotripin releasing hormone receptor (GNRHR) polymorphisms

Gonadotropin releasing hormone and its receptor (GNRHR) play a critical role in sexual differentiation and reproduction. Available evidence shows a strong genetic component in the timing of puberty. In bovines, there are significant differences within and among beef breeds in the time when bulls reach puberty. Despite its economic importance, there are not many SNPs or genetic markers associated with this characteristic. The aims of the study were to identify DNA polymorphism in the bovine GNRHR by re-sequencing analysis, determine haplotype phases, and perform a population study in a selected tag SNP in six breeds. Eight SNPs were detected, including: one in the Upstream Regulatory Region (URR), five in the coding regions, and two in non-coding regions. This polymorphism level corresponds to one variant every 249.4 bp and a global nucleotide diversity of 0.385. Two haplogroups comprising nine haplotypes and two linkage blocks were detected. Despite 5 tag SNPs were required to capture all variability, just one SNP allowed to define both haplogroups, and only two SNPs were needed to differentiate the most common haplotypes. An additional taq SNP was necessary to identify both URR variants. Allele-frequency analysis of a selected taq SNP among breeds showed a geographical cline. European Bos taurus breeds had lower frequencies of the C allele than B. indicus type cattle, while Creole cattle and Wagyu breeds had intermediate frequency. There was a significant correlation between frequency profile and timing of puberty among the studied breeds, which seems to suggest that genetic variation within bovine GNRHR gene could explain at least part of the reported variability.     

A novel phenotype for Lardon dwarfism in miniature Bos indicus cattle suggests that the expression of growth hormone receptor 1A in liver is required for normal growth

Mutations within the growth hormone receptor (GHR) gene that lead to an inactivated or truncated GHR protein cause abnormal growth and small adult size in a variety of species (Laron dwarfism). We studied a line of miniature Bos indicus cattle that have phenotypic (small mature size) and endocrine (increased blood growth hormone and decreased blood insulin-like growth factor-I concentrations) similarities to Laron dwarfs. Liver mRNA from miniature and control cattle was used to amplify a cDNA within the coding region of the GHR. The miniature cattle had GHR mRNA size (determined by Northern blot) and cDNA sequence that were similar to control cattle and, therefore, were unlike most Laron dwarf genotypes in which the GHR gene is mutated. Amounts of mRNA from liver as well as muscle (superficial neck and longissimus) were analyzed by ribonuclease protection assay for IGF-I, total GHR, GHR 1A (inducible, liver-specific GHR mRNA), and GHR 1B (constitutive GHR mRNA). Four control and five miniature bulls were tested. As expected, liver IGF-I mRNA was decreased in the miniature cattle (approximately 12% of control; P < 0.01). The amount of the total GHR as well as GHR 1A mRNA were also decreased in liver (17% and 19% of control, respectively; P < 0.01). Other GHR mRNA, including GHR 1B mRNA, were similar for miniature and control cattle. In muscle, there was a tendency (P < 0.10) for decreased IGF-I mRNA and increased GHR mRNA in miniature compared with control cattle. In summary, a novel phenotype for Laron dwarfism in Bos indicus cattle was associated with underexpression of GHR 1A mRNA, but not other GHR mRNA variants in liver. In addition to decreased GHR 1A mRNA, the miniature cattle had decreased liver IGF-I mRNA. Full expression of GHR 1A in liver, therefore, may be required for full liver IGF-I expression and normal growth.     

Characterization and validation of bovine Gonadotripin releasing hormone receptor (GNRHR) polymorphisms

Gonadotropin releasing hormone and its receptor (GNRHR) play a critical role in sexual differentiation and reproduction. Available evidence shows a strong genetic component in the timing of puberty. In bovines, there are significant differences within and among beef breeds in the time when bulls reach puberty. Despite its economic importance, there are not many SNPs or genetic markers associated with this characteristic. The aims of the study were to identify DNA polymorphism in the bovine GNRHR by re-sequencing analysis, determine haplotype phases, and perform a population study in a selected tag SNP in six breeds. Eight SNPs were detected, including: one in the Upstream Regulatory Region (URR), five in the coding regions, and two in non-coding regions. This polymorphism level corresponds to one variant every 249.4 bp and a global nucleotide diversity of 0.385. Two haplogroups comprising nine haplotypes and two linkage blocks were detected. Despite 5 tag SNPs were required to capture all variability, just one SNP allowed to define both haplogroups, and only two SNPs were needed to differentiate the most common haplotypes. An additional taq SNP was necessary to identify both URR variants. Allele-frequency analysis of a selected taq SNP among breeds showed a geographical cline. European Bos taurus breeds had lower frequencies of the C allele than B. indicus type cattle, while Creole cattle and Wagyu breeds had intermediate frequency. There was a significant correlation between frequency profile and timing of puberty among the studied breeds, which seems to suggest that genetic variation within bovine GNRHR gene could explain at least part of the reported variability.