Genome-wide association for heifer reproduction and calf performance traits in beef cattle

Genomic insights into key genes and QTLs involved in cattle reproduction

From an economic standpoint, reproductive characteristics are fundamental for sustainable production, particularly for monotocous livestock like cattle. A longer inter-calving interval is indicative of low reproductive capacity. This issue changes the dynamics of current and future lactations since it necessitates more inseminations, veterinary care, and hormone interventions. Various reproductive phenotypes, including ovulation, mating, fertility, pregnancy, embryonic growth, and calving-related traits, are observed in dairy cattle, and these traits have been associated with several QTLs. Calving ease, age at puberty, scrotal circumference, and inseminations per conception have been associated with 4437, 10623, 10498, and 2476 Quantitative Trait Loci (QTLs), respectively. This data offers valuable insights into enhancing and comprehending reproductive traits in livestock breeding. Studying QTLs associated with reproductive traits has far-reaching implications across various fields, from agriculture and animal husbandry to human health, evolutionary biology, and conservation. It provides the foundation for informed breeding practices, advances in biotechnology, and a deeper understanding of the genetic underpinnings of reproduction.

Identifying Regions of the Genome Associated with Conception Rate to the First Service in Holstein Heifers Bred by Artificial Insemination and as Embryo Transfer Recipients

Heifer conception rate to the first service (HCR1) is defined as the number of heifers that become pregnant to the first breeding service compared to the heifers bred. This study aimed to identify loci associated and gene sets enriched for HCR1 for heifers that were bred by artificial insemination (AI, n = 2829) or were embryo transfer (ET, n = 2086) recipients, by completing a genome-wide association analysis and gene set enrichment analysis using SNP data (GSEA-SNP). Three unique loci, containing four positional candidate genes, were associated (p < 1 × 10-5) with HCR1 for ET recipients, while the GSEA-SNP identified four gene sets (NES ≥ 3) and sixty-two leading edge genes (LEGs) enriched for HCR1. While no loci were associated with HCR1 bred by AI, one gene set and twelve LEGs were enriched (NES ≥ 3) for HCR1 with the GSEA-SNP. This included one gene (PKD2) shared between HCR1 AI and ET services. Identifying loci associated or enriched for HCR1 provides an opportunity to use them as genomic selection tools to facilitate the selection of cattle with higher reproductive efficiency, and to better understand embryonic loss.

Identification of quantitative trait loci and associated candidate genes for pregnancy success in Angus-Brahman crossbred heifers

BACKGROUND

In beef cattle, more than 50% of the energy input to produce a unit of beef is consumed by the female that produced the calf. Development of genomic tools to identify females with high genetic merit for reproductive function could increase the profitability and sustainability of beef production.

RESULTS

Genome-wide association studies (GWAS) were performed using a single-step genomic best linear unbiased prediction approach on pregnancy outcome traits from a population of Angus-Brahman crossbred heifers. Furthermore, a validation GWAS was performed using data from another farm. Heifers were genotyped with the Bovine GGP F250 array that contains 221,077 SNPs. In the discovery population, heifers were bred in winter breeding seasons involving a single round of timed artificial insemination (AI) followed by natural mating for 3 months. Two phenotypes were analyzed: pregnancy outcome to first-service AI (PAI; n = 1,481) and pregnancy status at the end of the breeding season (PEBS; n = 1,725). The heritability was estimated as 0.149 and 0.122 for PAI and PEBS, respectively. In the PAI model, one quantitative trait locus (QTL), located between 52.3 and 52.5 Mb on BTA7, explained about 3% of the genetic variation, in a region containing a cluster of γ-protocadherin genes and SLC25A2. Other QTLs explaining between 0.5% and 1% of the genetic variation were found on BTA12 and 25. In the PEBS model, a large QTL on BTA7 was synonymous with the QTL for PAI, with minor QTLs located on BTA5, 9, 10, 11, 19, and 20. The validation population for pregnancy status at the end of the breeding season were Angus-Brahman crossbred heifers bred by natural mating. In concordance with the discovery population, the large QTL on BTA7 and QTLs on BTA10 and 12 were identified.

CONCLUSIONS

In summary, QTLs and candidate SNPs identified were associated with pregnancy outcomes in beef heifers, including a large QTL associated with a group of protocadherin genes. Confirmation of these associations with larger populations could lead to the development of genomic predictions of reproductive function in beef cattle. Moreover, additional research is warranted to study the function of candidate genes associated with QTLs.

Harnessing Genomics and Transcriptomics Approaches to Improve Female Fertility in Beef Cattle-A Review

Female fertility is the foundation of the cow-calf industry, impacting both efficiency and profitability. Reproductive failure is the primary reason why beef cows are sold in the U.S. and the cause of an estimated annual gross loss of USD 2.8 billion. In this review, we discuss the status of the genomics, transcriptomics, and systems genomics approaches currently applied to female fertility and the tools available to cow-calf producers to maximize genetic progress. We highlight the opportunities and limitations associated with using genomic and transcriptomic approaches to discover genes and regulatory mechanisms related to beef fertility. Considering the complex nature of fertility, significant advances in precision breeding will rely on holistic, multidisciplinary approaches to further advance our ability to understand, predict, and improve reproductive performance. While these technologies have advanced our knowledge, the next step is to translate research findings from bench to on-farm applications.

Identifying pleiotropic variants and candidate genes for fertility and reproduction traits in Holstein cattle via association studies based on imputed whole-genome sequence genotypes

Background

Genetic progress for fertility and reproduction traits in dairy cattle has been limited due to the low heritability of most indicator traits. Moreover, most of the quantitative trait loci (QTL) and candidate genes associated with these traits remain unknown. In this study, we used 5.6 million imputed DNA sequence variants (single nucleotide polymorphisms, SNPs) for genome-wide association studies (GWAS) of 18 fertility and reproduction traits in Holstein cattle. Aiming to identify pleiotropic variants and increase detection power, multiple-trait analyses were performed using a method to efficiently combine the estimated SNP effects of single-trait GWAS based on a chi-square statistic.

Results

There were 87, 72, and 84 significant SNPs identified for heifer, cow, and sire traits, respectively, which showed a wide and distinct distribution across the genome, suggesting that they have relatively distinct polygenic nature. The biological functions of immune response and fatty acid metabolism were significantly enriched for the 184 and 124 positional candidate genes identified for heifer and cow traits, respectively. No known biological function was significantly enriched for the 147 positional candidate genes found for sire traits. The most important chromosomes that had three or more significant QTL identified are BTA22 and BTA23 for heifer traits, BTA8 and BTA17 for cow traits, and BTA4, BTA7, BTA17, BTA22, BTA25, and BTA28 for sire traits. Several novel and biologically important positional candidate genes were strongly suggested for heifer (SOD2, WTAP, DLEC1, PFKFB4, TRIM27, HECW1, DNAH17, and ADAM3A), cow (ANXA1, PCSK5, SPESP1, and JMJD1C), and sire (ELMO1, CFAP70, SOX30, DGCR8, SEPTIN14, PAPOLB, JMJD1C, and NELL2) traits.

Conclusions

These findings contribute to better understand the underlying biological mechanisms of fertility and reproduction traits measured in heifers, cows, and sires, which may contribute to improve genomic evaluation for these traits in dairy cattle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08555-z.

Inbreeding Calculated with Runs of Homozygosity Suggests Chromosome-Specific Inbreeding Depression Regions in Line 1 Hereford

The goal of this study was to evaluate inbreeding in a closed beef cattle population and assess phenotype prediction accuracy using inbreeding information. Effects of inbreeding on average daily gain phenotype in the Line 1 Hereford cattle population were assessed in this study. Genomic data were used to calculate inbreeding based on runs of homozygosity (ROH), and pedigree information was used to calculate the probability of an allele being identical by descent. Prediction ability of phenotypes using inbreeding coefficients calculated based on pedigree information and runs of homozygosity over the whole genome was close to 0, even in the case of significant inbreeding coefficient effects. On the other hand, inbreeding calculated per individual chromosomes' ROH yielded higher accuracies of prediction. Additionally, including only ROH from chromosomes with higher predicting ability further increased prediction accuracy. Phenotype prediction accuracy, inbreeding depression, and the effects of chromosome-specific ROHs varied widely across the genome. The results of this study suggest that inbreeding should be evaluated per individual regions of the genome. Moreover, mating schemes to avoid inbreeding depression should focus more on specific ROH with negative effects. Finally, using ROH as added information may increase prediction of the genetic merit of animals in a genomic selection program.

Fertility-Associated Polymorphism within Bovine ITGβ5 and Its Significant Correlations with Ovarian and Luteal Traits

Simple Summary

The fertility of bovines is essential for cattle husbandry. ITGβ5, which is suggested to be closely related to fertility, is known to mediate cell adhesion and affect a variety of cellular activities. To investigate the relationship between the ITGβ5 gene and the fertility of bovines, 696 ovarian samples were collected and six potential indel (insertion/deletion) within ITGβ5 were analysed, from which a deletion mutation was found to be polymorphic. The genotype frequency and allele frequency of this locus in the investigated population were calculated and the population genetic parameters analyzed. In addition, this locus was found to be significantly correlated with ovarian width and corpus luteum diameter. Considering the importance of ovary and corpus luteum in reproduction, it is tempting to speculate the crucial effects of ITGβ5 on bovine fertility, which still need further validation. The results of our study might provide a theoretical basis for future breeding to enhance bovine reproduction.

Abstract

There is an urgent need to improve bovine fertility, and molecular marker-assisted selection (MAS) can accelerate this process. Genome-wide association studies suggest that Integrin β5 (ITGβ5) might affect fertility in bovines. As a member of the integrins family, ITGβ5 can bind to the extracellular matrix and mediate various cellular processes. In our study, primers spanning six potential insertion/deletion (indel) polymorphisms within the ITGβ5 gene were designed and 696 ovary samples from different individuals, the vast majority not in oestrum were collected for genetic variation detection. A deletion locus, rs522759246, namely P1-D_13-bp, was found to be polymorphic. The allele D frequency was 0.152 and the polymorphism information content (PIC) value was 0.224, indicating a low-degree PIC. This locus did not follow the Hardy–Weinberg equilibrium (p = 1.200E-23). Importantly, associations between P1-D_13-bp and ovarian morphological traits were established. Polymorphisms of this locus had significant correlations with ovarian width (p = 0.015). The corpus luteum is also linked to fertility and P1-D_13-bp was significantly correlated with corpus luteum diameter (p = 0.005). In conclusion, an indel mutation within the bovine ITGβ5 gene was identified, which was significantly associated with several ovarian and luteal traits.

Genetic Features of Reproductive Traits in Bovine and Buffalo: Lessons From Bovine to Buffalo

Bovine and buffalo are important livestock species that have contributed to human lives for more than 1000 years. Improving fertility is very important to reduce the cost of production. In the current review, we classified reproductive traits into three categories: ovulation, breeding, and calving related traits. We systematically summarized the heritability estimates, molecular markers, and genomic selection (GS) for reproductive traits of bovine and buffalo. This review aimed to compile the heritability and genome-wide association studies (GWASs) related to reproductive traits in both bovine and buffalos and tried to highlight the possible disciplines which should benefit buffalo breeding. The estimates of heritability of reproductive traits ranged were from 0 to 0.57 and there were wide differences between the populations. For some specific traits, such as age of puberty (AOP) and calving difficulty (CD), the majority beef population presents relatively higher heritability than dairy cattle. Compared to bovine, genetic studies for buffalo reproductive traits are limited for age at first calving and calving interval traits. Several quantitative trait loci (QTLs), candidate genes, and SNPs associated with bovine reproductive traits were screened and identified by candidate gene methods and/or GWASs. The IGF1 and LEP pathways in addition to non-coding RNAs are highlighted due to their crucial relevance with reproductive traits. The distribution of QTLs related to various traits showed a great differences. Few GWAS have been performed so far on buffalo age at first calving, calving interval, and days open traits. In addition, we summarized the GS studies on bovine and buffalo reproductive traits and compared the accuracy between different reports. Taken together, GWAS and candidate gene approaches can help to understand the molecular genetic mechanisms of complex traits. Recently, GS has been used extensively and can be performed on multiple traits to improve the accuracy of prediction even for traits with low heritability, and can be combined with multi-omics for further analysis.

Genome-wide association study for age at puberty and resumption of cyclicity in a crossbred dairy cattle population

Fertility is of primary economic importance in dairy cattle and the most common reason for involuntary culling. However, standard fertility traits have very low heritability that renders genetic selection slow and difficult. In this study, we explored fertility from an endocrine standpoint. A total of 1,163 crossbred Holstein-Normande females in a 3-generation familial design were studied for progesterone level measured every 10 d to determine age at puberty (PUB) and commencement of postpartum luteal activity (CPLA). Genetic parameters were estimated using REML with WOMBAT software. The heritability estimates were 0.38 ± 0.10 and 0.16 ± 0.07 for PUB and CPLA, respectively. Moreover, the 2 traits were genetically correlated (0.45 ± 0.23), suggesting a partially common determinism. Because of the family structure, a linkage disequilibrium and linkage analysis approach was preferred over standard genome-wide association study to map genomic regions associated with these traits. Ten quantitative trait loci (QTL) were detected for PUB on chromosomes 1, 3, 11, 13, 14, 21, and 29, whereas 3 QTL were associated with CPLA on chromosomes 21 and 26. Only the QTL on chromosome 21 was common to both traits. Four functional candidate genes (NCOA2, GAS2, OVOL1, and FOSL1) were identified in the detected regions. These findings will contribute to a clearer understanding of fertility determinism and enhance the value of introducing endocrinological data in fertility studies.

Loci associated with conception rate in crossbred beef heifers

The inability of beef cattle to maintain full term pregnancies has become an economic concern for the beef industry. Herd management and nutritional improvements have alleviated environmental impacts on embryonic and fetal loss, yet additional gains can be made through genomic selection. The objectives of this study were to identify loci and gene-sets in crossbred beef heifers associated with the number of services required to become pregnant (TBRD) and heifer conception rate at first service (HCR1). Heifers (n = 709) from a commercial beef operation underwent one round of artificial insemination, before exposure to bulls for natural service for 50 days. Pregnancy and time of conception was determined by ultrasound 35 days after the breeding season. Heifers were genotyped using the GeneSeek (Lincoln, NE) Bovine GGP50K BeadChip prior to genome-wide association analyses (GWAA) conducted using an EIGENSTRAT-like model to identify loci associated (P < 1 × 10⁻⁵) with TBRD and HCR1. One locus was associated (P = 8.97 × 10⁻⁶) with TBRD on BTA19 and included the positional candidate gene ASIC2, which is differentially expressed in the endometrium of fertility classified heifers, and the positional candidate gene, SPACA3. Gene-set enrichment analyses using SNP (GSEA-SNP) data, was performed and identified one gene-set, oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen as enriched (NES = 3.15) with TBRD and contained nine leading edge genes that contributed to the enrichment of the gene set. The enriched gene-set is involved in catalyzing oxidation-reduction reactions, which have been associated with oxidative stressors impacting pregnancy success. No loci were associated nor gene-sets enriched with HCR1. Identification of loci, positional candidate genes, gene-sets and leading edge genes enriched for fertility facilitate genomic selection that allows producers to select for reproductively superior cattle, reduce costs associated with infertility, and increase percent calf crop.

Fine mapping of genomic regions associated with female fertility in Nellore beef cattle based on sequence variants from segregating sires

BACKGROUND

Impaired fertility in cattle limits the efficiency of livestock production systems. Unraveling the genetic architecture of fertility traits would facilitate their improvement by selection. In this study, we characterized SNP chip haplotypes at QTL blocks then used whole-genome sequencing to fine map genomic regions associated with reproduction in a population of Nellore (Bos indicus) heifers.

METHODS

The dataset comprised of 1337 heifers genotyped using a GeneSeek® Genomic Profiler panel (74677 SNPs), representing the daughters from 78 sires. After performing marker quality control, 64800 SNPs were retained. Haplotypes carried by each sire at six previously identified QTL on BTAs 5, 14 and 18 for heifer pregnancy and BTAs 8, 11 and 22 for antral follicle count were constructed using findhap software. The significance of the contrasts between the effects of every two paternally-inherited haplotype alleles were used to identify sires that were heterozygous at each QTL. Whole-genome sequencing data localized to the haplotypes from six sires and 20 other ancestors were used to identify sequence variants that were concordant with the haplotype contrasts. Enrichment analyses were applied to these variants using KEGG and MeSH libraries.

RESULTS

A total of six (BTA 5), six (BTA 14) and five (BTA 18) sires were heterozygous for heifer pregnancy QTL whereas six (BTA 8), fourteen (BTA 11), and five (BTA 22) sires were heterozygous for number of antral follicles' QTL. Due to inadequate representation of many haplotype alleles in the sequenced animals, fine mapping analysis could only be reliably performed for the QTL on BTA 5 and 14, which had 641 and 3733 concordant candidate sequence variants, respectively. The KEGG "Circadian rhythm" and "Neurotrophin signaling pathway" were significantly associated with the genes in the QTL on BTA 5 whereas 32 MeSH terms were associated with the QTL on BTA 14. Among the concordant sequence variants, 0.2% and 0.3% were classified as missense variants for BTAs 5 and 14, respectively, highlighting the genes MTERF2, RTMB, ENSBTAG00000037306 (miRNA), ENSBTAG00000040351, PRKDC, and RGS20. The potential causal mutations found in the present study were associated with biological processes such as oocyte maturation, embryo development, placenta development and response to reproductive hormones.

CONCLUSIONS

The identification of heterozygous sires by positionally phasing SNP chip data and contrasting haplotype effects for previously detected QTL can be used for fine mapping to identify potential causal mutations and candidate genes. Genomic variants on genes MTERF2, RTBC, miRNA ENSBTAG00000037306, ENSBTAG00000040351, PRKDC, and RGS20, which are known to have influence on reproductive biological processes, were detected.

Inbreeding depression in line 1 Hereford cattle population using pedigree and genomic information

This study aimed at assessing inbreeding and its effect on growth and fertility traits using the longtime closed line 1 Hereford cattle population. Inbreeding was estimated based on pedigree (FPED) and genomic information. For the latter, three estimates were derived based on the diagonal elements of the genomic relationship matrix using estimated (FGRM) or fixed (FGRM0.5) minor allele frequencies or runs of homozygosity (ROH) (FROH). A pedigree containing 10,186 animals was used to calculate FPED. Genomic inbreeding was evaluated using 785 animals genotyped for 30,810 SNP. Traits analyzed were birth weight (BWT), weaning weight (WWT), yearling weight (YWT), ADG, and age at first calving (AFC). The number of ROH per animal ranged between 6 and 119 segments with an average of 83. The shortest and longest segments were 1.36 and 64.86 Mb long, respectively, reflecting both ancient and recent inbreeding occurring in the last 30 to 40 generations. The average inbreeding was 29.2%, 16.1%, 30.2%, and 22.9% for FPED, FGRM, FGRM0.5, and FROH, respectively. FROH provided the highest correlations with FPED (r = 0.66). Across paternal half-sib families, with minimal variation in FPED, there were substantial variations in their genomic inbreeding. Inbreeding depression analyses showed that a 1% increase in an animal's FPED resulted in a decrease of 1.20 kg, 2.03 kg, and 0.004 kg/d in WWT, YWT, and ADG, respectively. Maternal inbreeding showed significantly negative effects on progeny growth performance. AFC increased by 1.4 and 0.8 d for each 1% increase in FPED of the cow and her dam, respectively. Using genomic inbreeding, similar impact on growth traits was observed although the magnitude of the effect varied between methods. Across all genomic measures, WWT, YWT, and ADG decreased by 0.21 to 0.53 kg, 0.46 to 1.13 kg, and 0.002 to 0.006 kg/d for each 1% increase in genomic inbreeding, respectively. Four chromosomes (9, 12, 17, and 27) were identified to have a significant association between their homozygosity (FROH-CHR) and growth traits. Variability in genomic inbreeding could be useful when deciding between full and half-sib selection candidates. Despite the high level of inbreeding in this study, its negative impact on growth performance was not as severe as expected, which may be attributed to the purging of the deleterious alleles due to natural or artificial selection over time.

Genome-wide association scan for heterotic quantitative trait loci in multi-breed and crossbred beef cattle

BACKGROUND

Heterosis has been suggested to be caused by dominance effects. We performed a joint genome-wide association analysis (GWAS) using data from multi-breed and crossbred beef cattle to identify single nucleotide polymorphisms (SNPs) with significant dominance effects associated with variation in growth and carcass traits and to understand the mode of action of these associations.

METHODS

Illumina BovineSNP50 genotypes and phenotypes for 11 growth and carcass traits were available for 6796 multi-breed and crossbred beef cattle. After performing quality control, 42,610 SNPs and 6794 animals were used for further analyses. A single-SNP GWAS for the joint association of additive and dominance effects was conducted in purebred, crossbred, and combined datasets using the ASReml software. Genomic breed composition predicted from admixture analyses was included in the mixed effect model to account for possible population stratification and breed effects. A threshold of 10% genome-wide false discovery rate was applied to declare associations as significant. The significant SNPs with dominance association were mapped to their corresponding genes at 100 kb.

RESULTS

Seven SNPs with significant dominance associations were detected for birth weight, weaning weight, pre-weaning daily gain, yearling weight and marbling score across the three datasets at a false discovery rate of 10%. These SNPs were located on bovine chromosomes 1, 3, 4, 6 and 21 and mapped to six putative candidate genes: U6atac, AGBL4, bta-mir-2888-1, REPIN1, ICA1 and NXPH1. These genes have interesting biological functions related to the regulation of gene expression, glucose and lipid metabolism and body fat mass. For most of the identified loci, we observed over-dominance association with the studied traits, such that the heterozygous individuals at any of these loci had greater genotypic values for the trait than either of the homozygous individuals.

CONCLUSIONS

Our results revealed very few regions with significant dominance genetic effects across all the traits studied in the three datasets used. Regarding the SNPs that were detected with dominance associations, further investigation is needed to determine their relevance in crossbreeding programs assuming that dominance effects are the main cause of (or contribute usefully to) heterosis.

Association between single nucleotide polymorphisms and sexual precocity in Nellore heifers

The aim of this study was to determine the extent (r²) of linkage disequilibrium (LD) in the genome of Nellore cattle, and to examine associations between single nucleotide polymorphisms (SNP) and age at first calving (AFC) and early pregnancy (EP) using a panel of high-density SNPs and data from 1182 Nellore females. A total of 13 contemporary groups (CG) were used consisting of farm, season, and year of birth. For genome-wide association analysis, SNPs with a minor allele frequency (MAF)<0.05 and animals with a call rate<0.90 were excluded, totaling 431,885 SNPs. For statistical analysis, a linear model was used for AFC and a threshold model for EP. To estimate the significance of the associations for the two traits, the model included the categorical fixed effects of CG, SNPs, and sire. In addition, the polygenic effect was included in the analysis. The additive effects and dominance deviations of Bonferroni-adjusted significant SNPs for AFC and EP were estimated using orthogonal contrasts. The average estimate of r² for all autosomes was 0.18 at a distance of 4.8kb and the mean MAF was 0.25±0.13. The LD decreased as the distance between markers increased: 0.35 (1kb) to 0.12 (100kb). Eleven significant associations were detected in seven different chromosomes. Seven SNPs were associated with AFC and four were associated with EP. Three SNPs were significant for both traits. The identification of SNPs associated with AFC and EP may contribute for selecting sexually precocious animals.