Skip to main content

Understanding the underlying genetic mechanisms for age at first calving, inter-calving period and scrotal circumference in Bonsmara cattle



Reproduction is a key feature of the sustainability of a species and thus represents an important component in livestock genetic improvement programs. Most reproductive traits are lowly heritable. In order to gain a better understanding of the underlying genetic basis of these traits, a genome-wide association was conducted for age at first calving (AFC), first inter-calving period (ICP) and scrotal circumference (SC) within the South African Bonsmara breed. Phenotypes and genotypes (120,692 single nucleotide polymorphisms (SNPs) post editing) were available on 7,128 South African Bonsmara cattle; the association analyses were undertaken using linear mixed models.


Genomic restricted maximum likelihood analysis of the 7,128 SA Bonsmara cattle yielded genomic heritability’s of 0.183 (SE = 0.021) for AFC, 0.207 (SE = 0.022) for ICP and 0.209 (SE = 0.019) for SC. A total of 16, 23 and 51 suggestive (P ≤ 4 × 10-6) SNPs were associated with AFC, ICP and SC, while 11, 11 and 44 significant (P ≤ 4 × 10-7) SNPs were associated with AFC, ICP and SC respectively. A total of 11 quantitative trait loci (QTL) and 11 candidate genes were co-located with these associated SNPs for AFC, with 10 QTL harbouring 11 candidate genes for ICP and 41 QTL containing 40 candidate genes for SC. The QTL identified were close to genes previously associated with carcass, fertility, growth and milk-related traits. The biological pathways influenced by these genes include carbohydrate catabolic processes, cellular development, iron homeostasis, lipid metabolism and storage, immune response, ovarian follicle development and the regulation of DNA transcription and RNA translation.


This was the first attempt to study the underlying polymorphisms associated with reproduction in South African beef cattle. Genes previously reported in cattle breeds for numerous traits bar AFC, ICP or SC were detected in this study. Over 20 different genes have not been previously reported in beef cattle populations and may have been associated due to the unique genetic composite background of the SA Bonsmara breed.

Peer Review reports


The improvement of reproductive efficiency is of major economic importance in beef production systems and improvements in total lifetime productivity is a key metric for efficiency, which is a function of both reproduction and total output per cow [1, 2]. In South Africa (SA), the majority of beef is produced under extensive production systems with average calving percentages (i.e., the proportion of cows that give birth in comparison to the number of cows that could possibly give birth) of 62% [3]. Bonsmara cattle are a beef breed, developed during the nineteen-sixties, resulting in a composite comprised of approximately 5/8 Afrikaner and 3/8 exotic breeds [4]. It is the largest breed represented in the seed stock and commercial beef industry in SA. As the beef industry primarily relies on reproductive efficiency, a reduction in the unproductive period of a female’s life would positively impact production costs and profit as well as having a favourable impact on the carbon footprint.

Genetic improvement in reproductive performance remains challenging, hindered by the generally low associated heritability (h2), coupled with the complexity of recording these traits and/or the late expression of reproductive traits [5,6,7]. The age at onset of puberty can impact overall productivity through the heifer becoming productive at an earlier age [8]. Puberty is experienced earlier in composite and crossbred animals compared to their purebred counterparts with the same trend observed for early maturing versus late maturing beef breeds [9]. Age of puberty is however difficult to record, and therefore age at first calving (AFC) is often used as an indicator of heifer fertility. The heritability estimates reported for AFC differ depending on the breed, with [10] estimating a h2 = 0.10 in Brahman cattle, [8, 11, 12] reporting on Nellore cattle (0.01 to 0.31) and [13] estimating a h2 = 0.08 for AFC in Angus and Hereford cattle. Recent studies on South African populations are limited with [14] reporting a h2 = 0.08 in Bonsmara cross cattle, while [5] reported higher heritability for AFC in Afrikaner (0.27) and Drakensberger (0.30) Sanga cattle breeds. Inter-calving period (ICP) is a relatively easy trait to record and is included in most beef and dairy breed societies recording schemes in South Africa [15]. Estimates of ICP yield low to moderate heritability (0.01–0.10) in Brahman and composite beef cattle breeds [6, 10, 16]. As scrotal circumference (SC) in bulls is relatively easy to measure, with a moderate to high heritability [6, 12, 17,18,19], it has been suggested as an indicator trait for age at puberty, the latter being resource-intensive to measure.

Several studies have reported positive genetic correlations between SC and growth traits like mature body weight (rg = 0.37–0.40) in composite beef cattle [20], weaning and mature weight (rg = 0.60–0.72) in Bos indicus cattle [21] as well as weaning (rg = 0.312) and yearling (rg = 0.519) weight in Nellore cattle [22]. Although weaker genetic correlations between SC and weaning weight (rg = 0.15) have been reported in SA Bonsmara bulls [23], genomic analyses of SC and SC-related traits have identified genes which are also known to associate with growth traits [24].

Genome-wide association studies (GWAS) conducted on European and tropically adapted beef cattle breeds revealed several potential genes for fertility traits including AFC [8, 12, 25, 26], ICP [27], pregnancy status [28], gestation length [29], sexual precocity [30] and SC [12, 31, 32]. Some studies have combined fertility traits with body weight at puberty [33] while multi-trait meta-analyses have also been applied [12], which indicated that similar regions of the genome harbour genetic variation that potentially influence reproductive traits in both genders.

The SA Bonsmara, classified as a Sanga type, is a unique composite breed of 3/8 exotic (Milk Shorthorn, Hereford) and 5/8 Afrikaner [4]. The breed was established through a well-documented crossbreeding program, with the aim of founding a local composite breed that was well adapted to the challenges of a diverse SA climate. This study was the first attempt to apply GWAS for fertility traits in SA Bonsmara cattle to provide insight on gene regions in the SA Bonsmara. The objective of this study was to perform a genome-wide association study for three reproductive traits (AFC, ICP and SC) in a South African Bonsmara population in order to identify quantitative trait loci (QTL) for these traits.


Variance components estimation

A genetic correlation of 0.37 was estimated between AFC and ICP in the SA Bonsmara using the breeding values derived from the bivariate model. Pedigree heritability estimates were 0.22 for AFC, 0.13 for ICP and 0.38 for SC. The genomic REML analysis yielded genomic heritability and standard errors (SE) of 0.183 (SE = 0.021) for AFC, 0.207 (SE = 0.022) for ICP and 0.209 (SE = 0.019) for SC.

Genomic population quality control

Individual based quality control resulted in the omission of ninety-five SA Bonsmara genotypes across the five genotyping arrays available. Assessment of identical by state (IBS) genetic distances yielded a multidimensional scaled plot (MDS; Fig. 1). Identification of outliers as well as their genotyped progeny led to the further removal of 103 genotypes (Fig. 2) resulting in a final sample population of 7,128 SA Bonsmara animals (4,403 males, 2,725 females). The effect of filtering out animals with reliabilities smaller than 0.01 and an effective record count (ERC) of less than 0.50 culminated in genome-wide population sizes of 4,460 animals for AFC (median ERC = 1.758), 4,276 animals (median ERC = 1.717) for ICP1, and 5,452 animals for SC (median ERC = 1.004).

Fig. 1
figure 1

Multidimensional Scaling of 7,231 SA Bonsmara Genotypes

Fig. 2
figure 2

Multidimensional Scaling of 7,128 SA Bonsmara Genotypes

Age at first calving

Eleven single nucleotide polymorphisms (SNPs) at the genome wide threshold p-value of less than 4 × 10− 7 and a further 16 suggestive SNPs at a p-value of less than 4 × 10− 6 were associated with AFC (Fig. 3). Pair-wise linkage disequilibrium (LD) analysis of significant SNPs resulted in the identification of 11 QTL across eight autosomes (Table 1). A total of 11 different genes were co-located with the associated QTL, with the frequency of the major alleles being between 0.527 and 0.982, respectively. Bos taurus autosome (BTA) 7 harboured two QTL containing nine genes (ARAP3, CLINT1, FCHSD1, LSM11, PCDHGA, PCDHGB, PCDHGC, RELL2 and THG1L) with each of these QTL having seven SNPs in pair-wise LD. The most significant QTL, located on BTA 17, was a single intergenic SNP with a minor allele frequency (MAF) of 0.13. The gene, PLCB1, resides in a QTL that spans four SNPs in pair-wise LD (r2 > 0.50), 212.17 kilobase pairs (kbp) across BTA 13, with the significant SNP (BovineHD1300000449; MAF = 0.15), associating with AFC.

Fig. 3
figure 3

Significant SNP above the red line (≤ 4 × 10− 7) for Age at First Calving (4,460 animals)

Table 1 List of identified QTLs and genes associated with Age at First Calving

Inter-calving period

Eleven SNPs were significantly associated with ICP, with a further 23 SNPs being suggestively associated (Fig. 4). The most significant SNP (BTA-16,045-no-rs; P = 2.83 × 10− 16, MAF = 0.13) was on BTA 17 and was an intergenic SNP. A total of ten QTL across eight autosomes were significantly associated with ICP (Table 2). A QTL 337.43 kbp in length containing six genes (ARAP3, FCHSD1, PCDHGA, PCDHGB, PCDHGC and RELL2) and a second QTL of 158.19 kbp in length harbouring three genes (CLINT1, LSM11 and THG1L) were both positioned on BTA 7. A QTL containing the LDAH gene, is a 131.85 kbp long and contains two SNPs in high LD (r2 = 0.74).

Fig. 4
figure 4

Significant SNP above the red line (≤ 4 × 10− 7) for Inter-Calving Period (4,276 animals)

Table 2 List of identified QTLs and genes associated with Inter Calving Period

Scrotal circumference

Forty-four SNPs were significantly (P ≤ 4 × 10− 7) associated with SC with a further 51 SNPs being suggestive (P ≤ 4 × 10− 6, Fig. 5). Pair-wise LD analysis identified a total of 41 QTL across 14 autosomes (Table 3). The most significant SNP, (BovineHD2300009170; P = 5.02 × 10− 11, MAF = 0.04), which is lowly segregating in the SA Bonsmara population, was an intron variant of SLC17A3 located on BTA 23. A QTL consisting of eleven SNPs in LD (r2 > 0.50) on BTA 11 consisted of downstream variants for the gene NEURL1B. Gene rich QTL were located on BTA 2, 11, 19 and 22. A 93.2 kbp QTL on BTA 2 consisting of six SNPs harboured the gene ABCA12, a small nucleolar RNA (RF00156) and an insertion/deletion copy number variant (CNV). Five QTL spanning a 1.65 Mega basepair (Mbp) portion of BTA 11 house the genes AAK1, ANTXR1, CNRIP1, GMCL1, MXD1, PP3R1 and SNRNP27. A homeobox gene dense QTL on BTA 19 sweeping 171.78 kbp consisted of four HOXB gene variants and TTLL6. On BTA 22, a significant SNP (Hapmap33950-BES3_Contig483_1359; P = 4.02 × 10− 9, MAF = 0.425) in moderate LD with two flanking SNPs (downstream r2 = 0.524, upstream r2 = 0.527) includes the genes ABHD6, PXK and RPP14 associated with SC in this study. The QTL with the most genes, namely ALDH1L1, CHST13, C22H3orf22, SLC41A3, TXNRD3, UROC1 and ZXDC, was comprised of five SNPs over a 137.92 kbp DNA region on BTA 22.

Fig. 5
figure 5

Significant SNP above the red line (≤ 4 × 10− 7) for Scrotal Circumference (5,452 animals)

Overlapping genes across traits

In this study, three QTL, two on BTA 7 and one on BTA 13, were detected to be associated with more than one reproductive trait. The QTL on BTA 13 contains four SNPs in LD, with the SNP (BovineHD1300000449) being significantly associated with both AFC (P = 5.13 × 10–11) and ICP (P = 5.13 × 10− 8). Two large aforementioned QTL located on BTA 7 are both significantly associated with AFC and ICP. The first QTL was significantly associated with both AFC (P = 1.15 × 10− 7) and ICP (P = 2.98 × 10− 7), while the SNP (ARS-BFGL-NGS-11368), located in this same QTL was significantly associated with ICP (P = 2.98 × 10− 7). The second QTL, with seven SNPs in LD with each other, was identified to be significantly linked to both AFC (P = 1.15 × 10− 7) and ICP (P = 5.30 × 10− 8).


This study was the first attempt to gain insight into the underlying genetic mechanisms for reproductive traits in a South African beef cattle breed. The pedigree and genomic based heritability estimates for AFC, ICP and SC were similar to those identified in studies in beef cattle breeds [6, 8, 10,11,12,13,14, 17,18,19]. The low to moderate heritability estimates for most reproductive traits and limited genomic studies on indigenous breeds in Southern Africa [34,35,36] justifies the further investigation of this study. Genomic tools hold the potential to unlock invaluable information [37], that may help explain physiological conditions that currently remain unanswered [38].

Age at first calving

In this study, a total of eleven significant SNPs associated with AFC were observed; none of these have been previously reported in Bos indicus and crossbred beef cattle to be associated with AFC. Of the eleven genes and two CNVs located in QTL associated with AFC, ten of these genes have not been previously associated with AFC or sexual precocity in beef cattle populations but have been reported in human [39] and sheep [40] populations. Four of the genes observed in this study (ENSBTAG00000032764, FCHSD1, PLCB1 and RELL2) have been previously associated in beef cattle populations [41,42,43] with two of these genes involved in more than one trait. Physical growth and body composition have a developmental effect on the reproductive organs which would determine the onset of puberty and subsequent AFC [44, 45]. ENSBTAG00000032764 (pseudogene) was previously reported as a candidate gene associated with carcass traits [41]. The resultant translated protein facilitates the storage of iron in a soluble, non-toxic form, which is an essential component for iron homeostasis.

Inter-calving period

The inter calving period which is the period between two calvingss has the lowest genetic and genomic heritability of all three traits in this study and, in practice, is a function of the ability to ovulate post-calving [45], express oestrus, establish and maintain pregnancy and gestation length [29, 46]. A total of eleven genes and three CNVs were based in QTL associated with ICP which included ten genes not previously associated with ICP or reproduction in beef cattle. The 337.43 kbp QTL on BTA 7 contains three genes (PCDHGA, PCDHGB, PCDHGC) from the protocadherin (PCDH) family group. Studies have shown the PCDH genes play an integral role in ovarian follicular [47] and embryonic [48] development. LDAH is a lipid droplet-associated hydrolase that is essential in lipid storage and was previously associated with feet and leg disorders in Danish Holstein cattle [49].

Scrotal circumference

Pair-wise LD analysis identified a total of 41 QTL across 14 autosomes (Table 3). Genes not previously reported in beef cattle populations include AAK1, ACAD19, CHST13, CNRIP1, CRISP1, DERA, KIF2A, PARVB, PNPO, PXK, RPP14, SKAP1, SLC41A3, SNRNP27, SSBP2, TCOM5B, TJP2, TTLL6, UBE2Z, UROC1 and ZXDC. Most of these genes have not been previously reported to be associated with SC or sperm-related traits, but upon further investigation into the biological pathways that the genes are involved in, it becomes clear that the genes may play some role in the expression of fertility related traits. Biological pathways discussed below include carbohydrate catabolic processes (CHST13, DERA and UROC1), cellular development (HOXB5, HOXB7, HOXB9, HOXB13 and TTL13), lipid metabolism (ABCA12, ABHD6, ACAD9 and PARVB), immune response (SKAP1) and regulation of DNA transcription and RNA translation (CDKAL1, E2F3, MXD1, RPP14, SNRNP27, SSBP2 and ZXDC).

ABCA12 mediates lipid transporter activity, signals receptor binding, and has active trans-membrane transporter activity [50, 51]  identified that ABCA12 is the major gene that influences Harlequin Ichthyosis in humans and later was identified in livestock species. ABCA12 was associated with birth weight in Holstein cattle [52]. Monoacylglycerol lipase (ABHD6) is a lipase and the major enzyme for bone morphogenic protein catabolism, which plays a key role in the formation of intraluminal vesicles and in lipid sorting [53]. Although not previously associated with SC or fertility traits in cattle breeds, ABHD6 was associated with average daily gain in crossbred beef cattle [54] as well as identified as a positional candidate gene for milk cholesterol in dairy cattle [55]. ACAD9, which catalyses the rate limiting step during the beta-oxidation of fatty acyl-CoA, was identified as a CNV in Italian sheep [56] and was found to be significantly associated with intramuscular fat content in Large White pigs [57]. Parvin beta (PARVB) is essential in establishing and/or maintaining cell polarity as well as actine cytoskeleton reorganization. Although not previously reported in beef cattle, a study on ketosis in German Holstein cattle [58] associated PARVB with non-alcoholic fatty liver disease, indicating that PARVB could be a part of lipid metabolism. An investigation into selection signatures in Indian swamp buffaloes identified PARVB [59].

Two genes previously reported to be linked to fertility traits in cattle include FXN and GMCL1. Frataxin (FXN) promotes the biosynthesis of heme and plays a role in the protection against iron-catalysed oxidative stress. FXN was the second ranking SNP in a logistic regression analysis of pregnancy status in Santa Gertrudis cows [60]. GMCL1 appears to be located in an extended selection signature that shows high haplotypic homozygosity [61] and was previously associated with fertility traits in beef cattle breeds [62]. Single Stranded Binding Protein (SSBP2), located on BTA 7, is a candidate QTL for conferring resistance to Johne’s disease in cattle [63, 64]. SSBP2 is involved in the positive regulation of transcription by RNA polymerase II and transcription by RNA polymerase II. A study in Holstein cattle associated SSBP2 with body conformation traits [65]. PPP3R1 is a regulatory subunit of calcineurin, which plays a role in neuronal calcium signalling. A study on transcriptome profiling of muscle in Nelore cattle [66] identified PPP3R1 participates in mitogen-activated protein kinase (MAPK) signalling pathway which is responsible for cell proliferation, differentiation, and apoptosis.

A range of genes identified to be associated with SC have not been previously identified in cattle, but have been in water buffaloes and crossbred buffaloes. ALDH1L1 [67], TJP2 [67] and TCOM5B [68] were all linked to milk composition traits, more specifically fat and protein yield. The QTL containing ALDH1L1, CHST13, SLC41A3, TXNRD3, UROC1 and ZXDC on BTA 22 was previously reported to be associated with somatic cell score in Holstein cattle [69]. The TXNRD3 genes is known to affect adipocyte differentiation through the Wnt signalling pathway [70].

On BTA 19, an SC-associated QTL harbours four homeobox genes and TTLL6. TTLL6 is a gene in a strong candidate region, which included HOXB7 and HOXB9, for controlling skeletal tail length in sheep [71]. The four HOXB genes in this QTL have not been previously reported in beef cattle, but other HOX gene clusters are known to affect sperm quality in humans [72]. It is known that poor sperm DNA methylation is associated with decreased male fertility and low embryo quality. The hypomethylation of microRNA and HOX gene clusters play a significant role in embryonic development and is evidence of the sperm’s epigenetic contribution. KIF2A is known modulate mitotic events during spermatogenesis [73].

A limited number of studies have considered SC as a direct trait of interest for genomic investigations. Reviews of literature on bull fertility highlight the biological processes associated with sperm quality, motility, and scrotal volume [24, 74, 75]. More recent association studies revolve around sexual precocity, especially in tropical cattle [32, 74, 76,77,78,79] located in Central and Southern America.

Overlapping genes across traits

Of the genes observed to be significantly associated with both AFC and ICP in the present study, PLCB1 was identified to be linked to stay-ability in Nelore cattle [80]. PLCB1 is known to be a target of a micro interfering RNA (miR-301b), which has been associated with ovarian follicle development in cattle breeds [81]. Puberty in a heifer occurs upon ovulation of a potentially fertile oocyte [44], while [82, 83] stress the importance of proper nutrition postpartum in order to re-establish ovarian activity for a shortened ICP. PLCB1 is an enzyme that hydrolyses phospholipids into fatty acids as well as other lipophilic molecules and is involved in oxidative stress responses [84]. The regulation of adipose tissue affects the metabolic hormone leptin, known to regulate reproductive function in female animals [85, 86]. Twelve haplotype blocks for PLCB1 were identified through an association analysis of carcass traits in Hanwoo cattle [42] and through ontology of this gene linked it to lipid metabolism. ARAP3, a GTPase-activating protein, and the multiple genes that are members of the PCDH family group have been reported as a selection signature in cattle related to immune response [87]. Although immunological studies in livestock species are limited, [88] reviewed the effect immune cells have on ovarian follicle development and the establishment of pregnancy.

FCHSD1 and RELL2 are located in a long run of homozygosity (ROH) detected in multiple Alpine-based dual-purpose breeds. These two genes are involved in the MAPK14/p38 cascade as well as apoptosis. Clathrin interactor 1 (CLINT1) plays a major role in the formation of coated vesicles. This gene was associated with milk yield, fat yield and percentage as well as protein yield and percentage in dairy cattle [84], while was linked to milk fat content in Simmental cows [43]. LSM11 and THG1L have not been previously reported in cattle breeds but have been associated with milk protein yield and milk protein percentages in Valle del Belice dairy sheep [40]. LSM11 being a small nuclear RNA that has processes the mRNA 3’-end prior to translation while THG1L is involved in the regulation of tRNA processing during translation.

The number of overlapping genes co-located in QTL shared by AFC and ICP, alongside the moderate genetic correlation (0.37) in this study, indicates fertility is initiated, regulated and maintained by pleiotropic genetic mechanisms. Multiple genes in this study have no obvious direct link with fertility traits and this further demonstrates the complexity of genetic mechanisms for traits such as AFC and ICP.


In this study numerous genes, ARAP3, CLINT1, FCHSD1, LSM11, PLCB1, RELL2, SM11 and THG1L were co-located in QTL that had a significant or suggestive association with both AFC and ICP. Numerous QTL were identified across 14 autosomes for SC, the majority of which had never been previously reported to be linked to reproductive traits. The identification of different genes with similar molecular and biological characteristics for these sex-limited traits reaffirms our understanding that these lowly heritable traits are influenced by many genes each contributing a small amount to the variation in these traits’ expression. Some genes related to carbohydrate catabolic processes, cellular development, iron homeostasis, lipid metabolism and storage, immune response, ovarian follicle development and the regulation of DNA transcription and RNA translation were identified as candidate genes for reproductive traits in SA Bonsmara cattle.


Genotypic data

Genotypes from 7,326 SA Bonsmara animals originating from one of five possible genotype arrays were available. A total of 1,950 animals were genotyped on the GeneSeek Genomic Profiler (GGP) 150K (140,113 SNPs), while 597 animals were genotyped on the GGP 80K (76,883 SNPs), 2,625 animals were genotyped on the Versa 50K (49,855 SNPs), 1,326 animals were genotyped on the SASB 50K (54,394 SNPs) with the remaining 828 on the ICBF IDB v.2 platform (52,445 SNPs). Only autosomal SNPs with a known base pair position, a call rate ≥ 0.90, a MAF ≥ 0.10 and did not significantly deviate from Hardy-Weinberg equilibrium (p > 0.001) were retained. All SNP locations were based on the UMD 3.1 genome build (GCF_000003055.6; [89]). Animals had a call rate of > 90%, while individuals with ≥ 0.95 identical genotypes were discarded as were families with more than 10% Mendelian errors. Quality control of SNP data was carried out using PLINK v.1.9 [90].

Population stratification

Identical by state genetic distances between animals were computed through a MDS analysis with PLINK v.1.9 [90]. The analysis involved a total of 7,231 SA Bonsmara genotypes at a density of 24,216 SNPs that are truly genotyped across all five arrays. Visualisation of the data, reduced into two dimensions, allowed for the detection of possible population stratification as well as outliers. The remaining SA Bonsmara genotypes (4,403 males, 2,725 females) were imputed to 120,692 SNPs using FImpute v.3 [91].

Phenotypic data

The SA Bonsmara minimum breed standards [92] indicate that a heifer must calve before 39 months of age and the first ICP cannot exceed 790 days. SA Bonsmara animals occur throughout all nine of South Africa’s provinces and are mainly raised in extensive natural pasture systems. The recording of weaning weight (205-day weight) is compulsory and facilitates the selection of bulls for post-weaning growth tests. Scrotal circumference is measured on bulls participating in central and farm-based growth tests at around 12 to 18 months of age. Standardised phenotypes for AFC (days), first ICP (days), and SC (millimetres) were available on 347,749 records for AFC, 206,505 for records ICP and 238,454 records for SC in individual SA Bonsmara animals from the LOGIX Genetic Evaluation System [93]. This was accompanied by pedigree information on 2,135,235 animals dating back to 01 June 1949, as well as data on the contributing systematic environmental effects associated with these traits.

Deregression of breeding values

In order to predict estimated breeding values, a bivariate animal linear model for AFC and ICP and a univariate animal linear model for SC were defined as follows;

$$\varvec{y} = \varvec{Xb} + \varvec{Zu} + \varvec{e}$$


y is the vector of phenotypes for AFC, ICP and SC;

b is a vector of fixed effects which include sex, herd, birth month and year, age in days at measurement of the phenotype covariate (linear regression);

u is a vector representing the direct additive-genetic effects, with u ~ N(0,A \({\sigma }_{u}^{2}\)), where A is the pedigree-based matrix and \({\sigma }_{u}^{2}\) is the direct genetic variance;

e represents the residual, where e ~ N(0,I\({\sigma }_{e}^{2}\)), with \({\sigma }_{e}^{2}\) representing the residual variance and I the identity matrix;

X and Z are incidence matrices for b and u respectively.

Estimation of variance components for the animal model stated above was calculated using restricted estimated maximised likelihood (REML) optimised with quasi-Newton procedure using analytical gradients in Variance Component Estimation (VCE) [94] software. MiX99 [95] was used to predict breeding values for AFC, ICP and SC using the same model in the estimation of variance components. Effective record contributions (ERCs) for each animal and trait were generated as described in [96] using the reversed reliability approximation method in APaX99 [97]. The EBVs of the genotyped animals for each trait were then deregressed using the Secant method [98] in MiX99 [95] alongside the generated ERC. Deregressed EBVs (DEBVs) were weighted using the formula set out by [99];



w is the weighting factor of the ith animal with a DEBV;

h2 is the heritability estimate for the respective traits,

r2 is the reliability of the DEBV for the ith animal for a specific trait and,

c is the proportion of genetic variance not accounted by the SNPs with a value of 0.90 being used for all weighting factors between all the traits under analysis.

Only animals with an ERC ≥ 0.50 and a reliability ≥ 0.01 were retained for each trait analysis.

Association analyses

A genomic relationship matrix (GRM) was constructed for each trait using the VanRaden method 1 [100]. Additive and residual genetic variances for each trait were computed via genomic REML (GREML) using GCTA v1.94 [101]. Weighted DEBVs were regressed on each SNP individually using a linear mixed model in WOMBAT [102].

$$\varvec{y} = \varvec{\mu} + \varvec{SNP} + \varvec{a} + \varvec{e}$$


y is the vector of phenotypes, the weighted DEBV;

µ is the fixed effect of the population mean;

SNP is the fixed effect of allele dosage for each SNP (coded as 0, 1 or 2);

a is the random effect of the animal, where a ~ (0,\({\varvec{G}\sigma }_{a}^{2}\)), with \({\sigma }_{a}^{2}\) representing the additive genetic variance of the animal;

G is the genomic relationship matrix among animals,

e represents the residual, where e ~ N(0,I\({\sigma }_{e}^{2}\)),

with \({\sigma }_{e}^{2}\) representing the residual variance and I the identity matrix.

The t-test statistics for all SNPs were obtained and subsequently transformed into lower tail p-values. To minimise false positives, the Benjamini-Hochberg False Discovery Rate (B-H FDR) method was applied to each SNP. SNPs with a P ≤ 4 × 10− 7 were considered to be genome-wide significant as per Bonferroni correction, with SNPs with a P ≤ 4 × 10− 6 being deemed suggestive. Manhattan plots, Figs. 3, 4 and 5, were generated in R using the qqman [103] package.

Defining QTLs and candidate genes

The extent of LD among significant SNPs (P ≤ 4 × 10− 7) was estimated, as was the pairwise LD among all SNPs within 5 Mb up and downstream of the significant SNP [104]. The start and end of each QTL was defined by SNPs furthest up and downstream of the significant SNP and had an r2 > 0.50 with other significant SNPs. If any QTL were deemed to be overlapping, these were consolidated into one large QTL. If no SNPs were in LD with the significant SNP, that SNP was deemed a quantitative trait nucleotide. Identified QTL were then explored using ENSEMBL ( according to the UMD 3.1 genome build in order to detect candidate genes residing within and Panther [105] was used to list the biological and metabolic functions and/or processes of possible genes.

Table 3 List of identified QTL and genes significantly associated with Scrotal Circumference

Data Availability

The datasets generated and/or analysed during the current study are available at



Age at first calving


Beef Genomics Program


Bos taurus autosome


Estimated breeding value


Deregressed estimated breeding value


Effective record contribution


Geneseek genomic profiler


Genomic relationship matrix


Genome-wide association studies

h2 :



Identical by state


Inter-calving period


Linkage disequilibrium


Minor allele frequency


Multidimensional scaling


Quantitative trait loci


Restricted estimated maximum likelihood


Run of homozygosity


South Africa


South African Stud Book


Scrotal circumference


Standard error


Single nucleotide polymorphism


  1. Burns BM, Fordyce G, Holroyd RG. A review of factors that impact on the capacity of beef cattle females to conceive, maintain a pregnancy and wean a calf-implications for reproductive efficiency in northern Australia. Anim Reprod Sci. 2010;122.

  2. Summers AF, Rosasco SL, Scholljegerdes EJ. Beef species-ruminant nutrition cactus beef symposium: influence of management decisions during heifer development on enhancing reproductive success and cow longevity. J Anim Sci. 2019;97:1407–14.

    PubMed  Google Scholar 

  3. Grobler SM, Scholtz MM, Greyling JPC, Neser FWC. Reproduction performance of beef cattle mated naturally following synchronization in the Central Bushveld bioregion of South Africa. S Afr J Anim Sci. 2014;44:70–4.

    Google Scholar 

  4. Bonsma J. Livestock production - a global approach. Cape Town, South Africa: Tafelberg Publishers Ltd; 1983.

    Google Scholar 

  5. Meyer K, Hammond K, Parnell PF, MacKinnon MJ, Sivarajasingam S. Estimates of heritability and repeatability for reproductive traits in australian beef cattle. Livest Prod Sci. 1990;25:15–30.

    Google Scholar 

  6. Cammack KM, Thomas MG, Enns RM. Reproductive Traits and their heritabilities in beef cattle. Prof Anim Sci. 2009;25:517–28.

    Google Scholar 

  7. Hawken RJ, Zhang YD, Fortes MRS, Collis E, Barris WC, Corbet NJ, et al. Genome-wide association studies of female reproduction in tropically adapted beef cattle. J Anim Sci. 2012;90:1398–410.

    CAS  PubMed  Google Scholar 

  8. Costa RB, Camargo GMF, Diaz IDPS, Irano N, Dias MM, Carvalheiro R, et al. Genome-wide association study of reproductive traits in Nellore heifers using bayesian inference. Genet Selection Evol. 2015;47:1–9.

    Google Scholar 

  9. Gupta SK, Singh P, Shinde KP, Lone SA, Kumar N, Kumar A. Strategies for attaining early puberty in cattle and buffalo: a review. Agricultural Reviews. 2016;37.

  10. Cavani L, Garcia DA, Carreno LOD, Ono RK, Pires MP, Farah MM, et al. Estimates of genetic parameters for reproductive traits in Brahman cattle breed. J Anim Sci. 2015;93:3287–91.

    CAS  PubMed  Google Scholar 

  11. Kluska S, Olivieri BF, Bonamy M, Chiaia HLJ, Feitosa FLB, Berton MP, et al. Estimates of genetic parameters for growth, reproductive, and carcass traits in Nelore cattle using the single step genomic BLUP procedure. Livest Sci. 2018;216:203–9.

    Google Scholar 

  12. Melo TP, Fortes MRS, Bresolin T, Mota LFM, Albuquerque LG, Carvalheiro R. Multitrait meta-analysis identified genomic regions associated with sexual precocity in tropical beef cattle. J Anim Sci. 2018;96:4087–99.

    PubMed  PubMed Central  Google Scholar 

  13. Pardo AM, Elzo MA, Gama LT, Melucci LM. Genetic parameters for growth and cow productivity traits in Angus, Hereford and crossbred cattle. Livest Sci. 2020;233:103952.

    Google Scholar 

  14. Corbet NJ, Shepherd RK, Burrow HM, Prayaga KC, van der Westhuizen J, Bosman DJ. Evaluation of Bonsmara and Belmont Red cattle breeds in South Africa. 2. Genetic parameters for growth and fertility. Aust J Exp Agric. 2006;46:213–23.

    Google Scholar 

  15. van Marle-Köster E, Visser C, Berry DP. A review of genomic selection - implications for the south african beef and dairy cattle industries. S Afr J Anim Sci. 2013;43.

  16. Van Der Westhuizen RR, Schoeman SJ, Jordaan GF, Van Wyk JB. Heritability estimates derived from threshold analyses for reproduction and stayability traits in a beef cattle herd. South Afr J Anim Sci. 2001;31:25–32.

    Google Scholar 

  17. Meyer K, Hammond K, Mackinnon MJ, Parnell PF. Estimates of covariances between reproduction and growth in australian beef cattle. J Anim Sci. 1991;69:3533–43.

    CAS  PubMed  Google Scholar 

  18. Martínez-Velázquez G, Gregory KE, Bennett GL, Van Vleck LD. Genetic relationships between scrotal circumference and female reproductive traits. J Anim Sci. 2003;81:395–401.

    PubMed  Google Scholar 

  19. Corbet NJ, Burns BM, Johnston DJ, Wolcott ML, Corbet DH, Venus BK, et al. Male traits and herd reproductive capability in tropical beef cattle. 2. Genetic parameters of bull traits. Anim Prod Sci. 2013;53:101–13.

    CAS  Google Scholar 

  20. Burrow HM. Variances and covariances between productive and adaptive traits and temperament in a composite breed of tropical beef cattle. Livest Prod Sci. 2001;70:213–33.

    Google Scholar 

  21. Abreu LRA, Martins PGMA, Mota LFM, Ferreira TA, Ribeiro VMP, Villela SDJ, et al. Genetic correlations between body weight, scrotal circumference and visual evaluation scores in Bos indicus cattle. Anim Sci J. 2018;89:1223–9.

    PubMed  Google Scholar 

  22. Schmidt PI, Campos GS, Roso VM, Souza FRP, Boligon AA. Genetic analysis of female reproductive efficiency, scrotal circumference and growth traits in Nelore cattle. Theriogenology. 2019;128:47–53.

    CAS  PubMed  Google Scholar 

  23. Van Der Westhuizen RR, Van Der Westhuizen J, Schoeman SJ. Genetic variance components for residual feed intake and feed conversion ratio and their correlations with other production traits in beef bulls. S Afr J Anim Sci. 2004;34:257–64.

    Google Scholar 

  24. Lirón JP, Prando AJ, Fernández ME, Ripoli MV, Rogberg-Muñoz A, Goszczynski DE, et al. Association between GNRHR, LHR and IGF1 polymorphisms and timing of puberty in male Angus cattle. BMC Genet. 2012;13:26.

    PubMed  PubMed Central  Google Scholar 

  25. Regatieri IC, Boligon AA, Costa RB, de Souza FRP, Baldi F, Takada L, et al. Association between single nucleotide polymorphisms and sexual precocity in Nellore heifers. Anim Reprod Sci. 2017;177:88–96.

    CAS  PubMed  Google Scholar 

  26. Mota RR, Guimarães SEF, Fortes MRS, Hayes B, Silva FF, Verardo LL, et al. Genome-wide association study and annotating candidate gene networks affecting age at first calving in Nellore cattle. J Anim Breed Genet. 2017;134:484–92.

    CAS  PubMed  Google Scholar 

  27. Berry DP, Kearney JF, Twomey K, Evans RD. Genetics of reproductive performance in seasonal calving dairy cattle production systems. Ir J Agricultural Food Res. 2012;52.

  28. McDaneld TG, Kuehn LA, Thomas MG, Snelling WM, Smith TPL, Pollak EJ, et al. Genomewide association study of reproductive efficiency in female cattle. J Anim Sci. 2014;92:1945–57.

    CAS  PubMed  Google Scholar 

  29. Purfield DC, Evans RD, Carthy TR, Berry DP. Genomic regions Associated with Gestation length detected using whole-genome sequence data differ between dairy and beef cattle. Front Genet. 2019;10:1068.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Nascimento AV, Matos MC, Seno LO, Romero ARS, Garcia JF, Grisolia AB. Genome wide association study on early puberty in Bos indicus. Genet Mol Res. 2016;15.

  31. Santana ML, Eler JP, Bignardi AB, Ferraz JBS. Genetic associations among average annual productivity, growth traits, and stayability: a parallel between Nelore and composite beef cattle. J Anim Sci. 2013;91:2566–74.

    CAS  PubMed  Google Scholar 

  32. da Silva Romero AR, Siqueira F, Santiago GG, de Almeida Regitano LC, de Souza Júnior MD et al. Almeida Torres Júnior RA,. Prospecting genes associated with navel length, coat and scrotal circumference traits in Canchim cattle. Livest Sci. 2018;210 September 2017:33–8.

  33. Casas E, Lunstra DD, Stone RT. Quantitative trait loci for male reproductive traits in beef cattle. Anim Genet. 2004;35:451–3.

    CAS  PubMed  Google Scholar 

  34. Makina SO, Muchadeyi FC, van Marle-Köster E, MacNeil MD, Maiwashe A. Genetic diversity and population structure among six cattle breeds in South Africa using a whole genome SNP panel. Front Genet. 2014;5:333.

    PubMed  PubMed Central  Google Scholar 

  35. Makina SO, Whitacre LK, Decker JE, Taylor JF, MacNeil MD, Scholtz MM, et al. Insight into the genetic composition of south african Sanga cattle using SNP data from cattle breeds worldwide. Genet Selection Evol. 2016;48:88.

    Google Scholar 

  36. Lashmar SF, Visser C, Muchadeyi FC. Factors influencing imputation accuracy for the South African Drakensberger beef cattle breed. In: World Congress on Genetics Applied to Livestock Production. 2018. p. 11.472.

  37. Meuwissen TH, Hayes BJ, Goddard ME. Prediction of total genetic value using genome-wide dense marker maps. Genetics. 2001;157:1819–29.

    CAS  PubMed  PubMed Central  Google Scholar 

  38. Dekkers JCM, Hospital F. The use of molecular genetics in the improvement of agricultural populations. Nat Rev Genet. 2002;3:22–32.

    CAS  PubMed  Google Scholar 

  39. Moua P, Checketts M, Xu LG, Shu HB, Reyland ME, Cusick JK. RELT family members activate p38 and induce apoptosis by a mechanism distinct from TNFR1. Biochem Biophys Res Commun. 2017;491:25–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  40. Mohammadi H, Farahani AHK, Moradi MH, Mastrangelo S, Di Gerlando R, Sardina MT et al. Weighted single-step genome-wide Association Study uncovers known and novel candidate genomic regions for milk production traits and somatic cell score in Valle del belice dairy Sheep. Animals. 2022;12.

  41. Purfield DC, Evans RD, Berry DP. Reaffirmation of known major genes and the identification of novel candidate genes associated with carcass-related metrics based on whole genome sequence within a large multi-breed cattle population. BMC Genomics. 2019;20.

  42. Srivastava S, Srikanth K, Won S, Son JH, Park JE, Park W et al. Haplotype-based genome-wide association study and identification of candidate genes associated with carcass traits in Hanwoo cattle. Genes (Basel). 2020;11.

  43. Macciotta NPP, Gaspa G, Bomba L, Vicario D, Dimauro C, Cellesi M, et al. Genome-wide association analysis in italian simmental cows for lactation curve traits using a low-density (7K) SNP panel. J Dairy Sci. 2015;98:8175–85.

    CAS  PubMed  Google Scholar 

  44. Perry GA. Harnessing basic knowledge of factors controlling puberty to improve synchronization of estrus and fertility in heifers. J Anim Sci. 2012;90:1172–82.

    CAS  PubMed  Google Scholar 

  45. Wathes DC, Pollott GE, Johnson KF, Richardson H, Cooke JS. Heifer fertility and carry over consequences for life time production in dairy and beef cattle. Animal. 2014;8:91–104.

    PubMed  Google Scholar 

  46. Berry DP, Wall E, Pryce JE. Genetics and genomics of reproductive performance in dairy and beef cattle. Animal. 2014;8:105–21.

    PubMed  Google Scholar 

  47. Hou Y, Wang Y, Xu S, Qi G, Wu X. Bioinformatics identification of microRNAs involved in polycystic ovary syndrome based on microarray data. Mol Med Rep. 2019;20:281–91.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Yagi T, Takeichi M. Cadherin superfamily genes: functions, genomic organization, and neurologic diversity. Genes Dev. 2000;14:1169–80.

    CAS  PubMed  Google Scholar 

  49. Wu X, Guldbrandtsen B, Lund MS, Sahana G. Association analysis for feet and legs disorders with whole-genome sequence variants in 3 dairy cattle breeds. J Dairy Sci. 2016;99:7221–31.

    CAS  PubMed  Google Scholar 

  50. Thomas AC, Cullup T, Norgett EE, Hill T, Barton S, Dale BA, et al. ABCA12 is the major harlequin ichthyosis gene. J Invest Dermatology. 2006;126:2408–13.

    CAS  Google Scholar 

  51. Charlier C, Coppieters W, Rollin F, Desmecht D, Agerholm JS, Cambisano N, et al. Highly effective SNP-based association mapping and management of recessive defects in livestock. Nat Genet. 2008;40:449–54.

    CAS  PubMed  Google Scholar 

  52. Cole JB, Waurich B, Wensch-Dorendorf M, Bickhart DM, Swalve HH. A genome-wide association study of calf birth weight in Holstein cattle using single nucleotide polymorphisms and phenotypes predicted from auxiliary traits. J Dairy Sci. 2014;97:3156–72.

    CAS  PubMed  Google Scholar 

  53. Navia-Paldanius D, Savinainen JR, Laitinen JT. Biochemical and pharmacological characterization of human α/β-hydrolase domain containing 6 (ABHD6) and 12 (ABHD12). J Lipid Res. 2012;53:2413–24.

    CAS  PubMed  PubMed Central  Google Scholar 

  54. Schweer KR, Kachman SD, Kuehn LA, Freetly HC, Pollak JE, Spangler ML. Genome-wide association study for feed efficiency traits using SNP and haplotype models. J Anim Sci. 2018;96:2086–98.

    PubMed  PubMed Central  Google Scholar 

  55. Do DN, Schenkel FS, Miglior F, Zhao X, Ibeagha-Awemu EM. Genome wide association study identifies novel potential candidate genes for bovine milk cholesterol content. Sci Rep. 2018;8:13239.

    PubMed  PubMed Central  Google Scholar 

  56. Di Gerlando R, Mastrangelo S, Tolone M, Rizzuto I, Sutera AM, Moscarelli A et al. Identification of Copy Number Variations and genetic diversity in italian Insular Sheep Breeds. Animals. 2022;12.

  57. Puig-Oliveras A, Revilla M, Castelló A, Fernández AI, Folch JM, Ballester M. Expression-based GWAS identifies variants, gene interactions and key regulators affecting intramuscular fatty acid content and composition in porcine meat. Sci Rep. 2016;6.

  58. Klein SL, Scheper C, Brügemann K, Swalve HH, König S. Phenotypic relationships, genetic parameters, genome-wide associations, and identification of potential candidate genes for ketosis and fat-to-protein ratio in german holstein cows. J Dairy Sci. 2019;102:6276–87.

    CAS  PubMed  Google Scholar 

  59. Ravi Kumar D, Joel Devadasan M, Surya T, Vineeth MR, Choudhary A, Sivalingam J, et al. Genomic diversity and selection sweeps identified in Indian swamp buffaloes reveals it’s uniqueness with riverine buffaloes. Genomics. 2020;112:2385–92.

    CAS  PubMed  Google Scholar 

  60. Li Y, Kijas J, Henshall JM, Lehnert S, Mcculloch R, Reverter A. Using Random forests (RF) to prescreen candidate genes: a new prospective for GWAS. Proc 10th World Congress Genet Appl Livest Prod. 2014.

    Article  Google Scholar 

  61. Rothammer S, Seichter D, Förster M, Medugorac I. A genome-wide scan for signatures of differential artificial selection in ten cattle breeds. BMC Genomics. 2013;14.

  62. Khatkar MS, Randhawa IAS, Raadsma HW. Meta-assembly of genomic regions and variants associated with female reproductive efficiency in cattle. Livest Sci. 2014;166:144–57.

    Google Scholar 

  63. Pant SD, Schenkel FS, Verschoor CP, You Q, Kelton DF, Moore SS, et al. A principal component regression based genome wide analysis approach reveals the presence of a novel QTL on BTA7 for MAP resistance in Holstein cattle. Genomics. 2010;95:176–82.

    CAS  PubMed  Google Scholar 

  64. Mallikarjunappa S, Brito LF, Pant SD, Schenkel FS, Meade KG, Karrow NA. Johne’s disease in dairy cattle: an immunogenetic perspective. Front Veterinary Sci. 2021;8.

  65. Jiang J, Ma L, Prakapenka D, VanRaden PM, Cole JB, Da Y. A large-scale genome-wide association study in U.S. Holstein cattle. Front Genet. 2019;10.

  66. Souza LL, Zorzetto MF, Ricci TJT, Canesin RC, Dias e Silva NC, Negrão JA, et al. Relationship between performance, metabolic profile, and feed efficiency of lactating beef cows. Trop Anim Health Prod. 2019;51:2045–55.

    PubMed  Google Scholar 

  67. Du C, Deng T, Zhou Y, Ye T, Zhou Z, Zhang S, et al. Systematic analyses for candidate genes of milk production traits in water buffalo (Bubalus Bubalis). Anim Genet. 2019;50:207–16.

    CAS  PubMed  Google Scholar 

  68. Deng TX, Ma XY, Lu XR, Duan AQ, Shokrollahi B, Shang JH. Signatures of selection reveal candidate genes involved in production traits in chinese crossbred buffaloes. J Dairy Sci. 2022;105:1327–37.

    CAS  PubMed  Google Scholar 

  69. Durán Aguilar M, Román Ponce SI, Ruiz López FJ, González Padilla E, Vásquez Peláez CG, Bagnato A, et al. Genome-wide association study for milk somatic cell score in Holstein cattle using copy number variation as markers. J Anim Breed Genet. 2017;134:49–59.

    PubMed  Google Scholar 

  70. Kipp AP, Müller MF, Göken EM, Deubel S, Brigelius-Flohé R. The selenoproteins GPx2, TrxR2 and TrxR3 are regulated by wnt signalling in the intestinal epithelium. Biochimica et Biophysica Acta (BBA) -. Gen Subj. 2012;1820:1588–96.

    CAS  Google Scholar 

  71. Ahbara AM, Musa HH, Robert C, Abebe A, Al-Jumaili AS, Kebede A et al. Natural adaptation and human selection of northeast african sheep genomes. Genomics. 2022;114.

  72. Hammoud SS, Nix DA, Zhang H, Purwar J, Carrell DT, Cairns BR. Distinctive chromatin in human sperm packages genes for embryo development. Nature. 2009;460:473–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  73. Ma DD, Wang DH, Yang WX. Kinesins in spermatogenesis. Biol Reprod. 2017;96:267–76.

    PubMed  Google Scholar 

  74. Fortes MRS, DeAtley KL, Lehnert SA, Burns BM, Reverter A, Hawken RJ, et al. Genomic regions associated with fertility traits in male and female cattle: advances from microsatellites to high-density chips and beyond. Anim Reprod Sci. 2013;141:1–19.

    CAS  PubMed  Google Scholar 

  75. Taylor JF, Schnabel RD, Sutovsky P. Identification of genomic variants causing sperm abnormalities and reduced male fertility. Anim Reprod Sci. 2018;194:57–62.

    CAS  PubMed  PubMed Central  Google Scholar 

  76. Buzanskas ME, Grossi D, do A, Ventura RV, Schenkel FS, Chud TCS, Stafuzza NB et al. Candidate genes for male and female reproductive traits in Canchim beef cattle. J Anim Sci Biotechnol. 2017;8.

  77. Sweett H, Miglior F, Livernois A, Fonseca P, Id-Lahoucine S, Troya E, et al. Genome-wide association study to identify genomic regions and single nucleotide polymorphisms functionally associated with bull fertility. J Anim Sci. 2018;96:138–9.

    PubMed Central  Google Scholar 

  78. de Melo TP, Salinas Fortes MR, Hayes B, de Albuquerque LG, Carvalheiro R. Across-breed validation study confirms and identifies new loci associated with sexual precocity in Brahman and Nellore cattle. J Anim Breed Genet. 2020;137:139–54.

    PubMed  Google Scholar 

  79. Stafuzza NB, da Silva Costa EEV, Silva RM, de O, da Costa Filho LCC, Barbosa FB, Macedo GG, et al. Genome-wide association study for age at puberty in young Nelore bulls. J Anim Breed Genet. 2020;137:234–44.

    CAS  PubMed  Google Scholar 

  80. Sbardella AP, Watanabe RN, da Costa RM, Bernardes PA, Braga LG, Rey FSB et al. Genome-wide association study provides insights into important genes for reproductive traits in nelore cattle. Animals. 2021;11.

  81. Zielak-Steciwko AE, Browne JA, McGettigan PA, Gajewska M, Dzięcioł M, Szulc T, et al. Expression of microRNAs and their target genes and pathways associated with ovarian follicle development in cattle. Physiol Genomics. 2014;46:735–45.

    CAS  PubMed  Google Scholar 

  82. Wathes DC, Pollott GE, Johnson KF, Richardson H, Cooke JS. Heifer fertility and carry over consequences for life time production in dairy and beef cattle. 2014.

  83. Mukasa-Mugerwa E. A Review of a Reproductive Performance of Female Bos Indicus (zebu) Cattle. Addis Ababa, Ethiopia: ILCA Monograph 6; 1989.

  84. Xu L, Shi L, Liu L, Liang R, Li Q, Li J et al. Analysis of liver proteome and identification of critical proteins affecting milk Fat, protein, and Lactose Metabolism in Dariy cattle with iTRAQ. Proteomics. 2019;19.

  85. D’Occhio MJ, Baruselli PS, Campanile G. Influence of nutrition, body condition, and metabolic status on reproduction in female beef cattle: a review. Theriogenology. 2019;125:277–84.

    PubMed  Google Scholar 

  86. Zieba DA, Amstalden M, Williams GL. Regulatory roles of leptin in reproduction and metabolism: a comparative review. Domest Anim Endocrinol. 2005;29:166–85.

    CAS  PubMed  Google Scholar 

  87. Porto-Neto LR, Sonstegard TS, Liu GE, Bickhart DM, Vb M, Silva D et al. Genomic divergence of zebu and taurine cattle identified through high-density SNP genotyping. BMC Genomics. 2013;14.

  88. Fair T. The contribution of the maternal immune system to the establishment of pregnancy in cattle. Front Immunol. 2015;6.

  89. Zimin AV, Delcher AL, Florea L, Kelley DR, Schatz MC, Puiu D et al. A whole-genome assembly of the domestic cow, Bos taurus. Genome Biol. 2009;10.

  90. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D et al. PLINK: a toolset for whole-genome association and population-based linkage analysis. 2007.

  91. Sargolzaei M, Chesnais JP, Schenkel FS. A new approach for efficient genotype imputation using information from relatives. BMC Genomics. 2014;15.

  92. Standards. & Rejection Codes – Bonsmara SA. Accessed 24 May 2023.

  93. SA Stud Book / SA Stamboek. Accessed 25 May 2023.

  94. Groeneveld E. VCE User’s Guide and Reference Manual Version 6.0. 2010.

  95. MiX99 Development Team. MiX99: A software package for solving large mixed model equations. Release 17.11. 2017;:Natural Resources Institute Finland (Luke). Jokioi.

  96. Harris B, Johnson D. Approximate reliability of genetic evaluations under an animal model. J Dairy Sci. 1998;81:2723–8.

    CAS  PubMed  Google Scholar 

  97. Lidauer M, Matilainen K, Mäntysaari E, Pitkänen T, Taskinen M, Strandén I. Technical Reference Guide for MiX99 Pre-Processor. 2017;:1–87.

  98. Strandén I, Mäntysaari EA. A recipe for multiple trait deregression. Interbull Bull. 2010;:21.

  99. Garrick DJ, Taylor JF, Fernando RL. Deregressing estimated breeding values and weighting information for genomic regression analyses. Genet Selection Evol. 2009;41:55.

    Google Scholar 

  100. VanRaden PM. Efficient methods to compute genomic predictions. J Dairy Sci. 2008;91:4414–23.

    CAS  PubMed  Google Scholar 

  101. Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88:76–82.

    CAS  PubMed  PubMed Central  Google Scholar 

  102. Meyer K. WOMBAT: a tool for mixed model analyses in quantitative genetics by restricted maximum likelihood (REML). J Zhejiang Univ Sci B. 2007;8:815–21.

    PubMed  PubMed Central  Google Scholar 

  103. Turner SD. qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. J Open Source Softw. 2018;3:731.

    Google Scholar 

  104. Twomey AJ, Berry DP, Evans RD, Doherty ML, Graham DA, Purfield DC. Genome-wide association study of endo-parasite phenotypes using imputed whole-genome sequence data in dairy and beef cattle. Genet Selection Evol. 2019;51.

  105. Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, et al. PANTHER version 11: expanded annotation data from Gene Ontology and Reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2017;45:D183–9.

    CAS  PubMed  Google Scholar 

Download references


Acknowledgments must be made to the South African Beef Genomics Project for the funding provided for the creation of a reference population. Further mentions must include the Bonsmara Breeders Society and SA Stud Book and Animal Improvement Association for the addition of the lower density genotypes and availability of phenotypes and pedigrees for the traits of interest in this study.


Funding was provided by the South African Beef Genomics Project for the creation of a reference population.

Author information

Authors and Affiliations



Jason J Reding was responsible for the statistical analyses, figures, tables and writing the manuscript. Robert R van der Westhuizen provided statistical support. Donagh P Berry gave expert advice on the statistical approach and interpretation. Este van Marle-Köster assisted with conceptualization, drafting and editing the paper. All authors reviewed the final manuscript for submission.

Corresponding author

Correspondence to Jason J Reding.

Ethics declarations

Ethics approval and consent to participate

Consent was provided by the Bonsmara Breed Society, for use of the phenotypic and genotypic data and ethical approval was granted by the Research/Ethics Committee (EC-180000127), Faculty of Natural and Agricultural Sciences at the University of Pretoria for the use of external data.

Consent for publication

Not applicable.

Competing interests

I declare that the authors have no competing interests as defined by BMC, or other interests that might be perceived to influence the results and/or discussion reported in this paper.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Reding, J.J., van der Westhuizen, R.R., Berry, D.P. et al. Understanding the underlying genetic mechanisms for age at first calving, inter-calving period and scrotal circumference in Bonsmara cattle. BMC Genomics 24, 480 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: