Advancing Dairy and Beef Genetics Through Genomic Technologies

The US beef and dairy industries have made remarkable advances in sustainability and productivity through technological advancements, including selective breeding. Yet, challenges persist due to the complex nature of quantitative traits. While the beef industry has progressed in adopting genomic technologies, the availability of phenotypic data remains an obstacle. To meet the need for sustainable production systems, novel traits are being targeted for selection. Additionally, emerging approaches such as genome editing and high-throughput phenotyping hold promise for further genetic progress. Future research should address the challenges of translating functional genomic findings into practical applications, while simultaneously harnessing analytical methods.

Machine learning to identify endometrial biomarkers predictive of pregnancy success following artificial insemination in dairy cows†

The objective was to identify a set of genes whose transcript abundance is predictive of a cow's ability to become pregnant following artificial insemination. Endometrial epithelial cells from the uterine body were collected for RNA sequencing using the cytobrush method from 193 first-service Holstein cows at estrus prior to artificial insemination (day 0). A group of 253 first-service cows not used for cytobrush collection were controls. There was no effect of cytobrush collection on pregnancy outcomes at day 30 or 70 or on pregnancy loss between days 30 and 70. There were 2 upregulated and 214 downregulated genes (false discovery rate < 0.05, absolute fold change >2-fold) for cows pregnant at day 30 versus those that were not pregnant. Functional terms overrepresented in the downregulated genes included those related to immune and inflammatory responses. Machine learning for fertility biomarkers with the R package BORUTA resulted in identification of 57 biomarkers that predicted pregnancy outcome at day 30 with an average accuracy of 77%. Thus, machine learning can identify predictive biomarkers of pregnancy in endometrium with high accuracy. Moreover, sampling of endometrial epithelium using the cytobrush can help understand functional characteristics of the endometrium at artificial insemination without compromising cow fertility. Functional characteristics of the genes comprising the set of biomarkers is indicative that a major determinant of cow fertility, at least for first insemination after calving, is immune status of the uterus, which, in turn, is likely to reflect the previous history of uterine disease.

Impact of inbreeding on production, fertility, and health traits in German Holstein dairy cattle utilizing various inbreeding estimators

In dairy cattle production, it is important to understand how inbreeding affects production, fertility, and health traits. However, there is still limited use of genomic information to estimate inbreeding, despite advancements in genotyping technologies. To address this gap, we investigated the effect of inbreeding on German Holstein dairy cattle using both pedigree-based and genomic-based inbreeding estimators. We employed one method based on pedigree information (Fped) together with 6 genomic-based methods, including 3 genome-wide complex trait analysis software estimators (Fhat1, Fhat2, Fhat3), VanRaden's first method (FVR1, with observed allele frequencies, and FVR0.5, when allele frequencies are set to 0.5), and one based on runs of homozygosity (Froh). Data from 24,489 cows with both phenotypes and genotypes were used, with a pedigree including 232,780 animals born between 1970 and 2018. We analyzed the effects of inbreeding depression on production, fertility, and health traits separately, using single-trait linear animal models as well as threshold models to account for the binary nature of the health traits. For the health traits, we transformed solutions from the liability scale to a probability scale for easier interpretation. Our results showed that the mean inbreeding coefficients from all estimators ranged from -0.003 to 0.243, with negative values observed for most genomic-based methods. We found out that a 1% increase in inbreeding caused a depression ranging from 25.94 kg (Fhat1) to 40.62 kg (Fhat3), 1.18 kg (Fhat2) to 1.70 kg (Fhat3), 0.90 kg (Fhat2) to 1.45 kg (Froh and Fhat3), 0.19 (Fped) to 0.34 d (Fhat3) for 305-d milk yield, fat, protein, and calving interval, respectively. The health traits showed very slight gradual changes when inbreeding was increased steadily from 0% to 50%, with digital dermatitis showing a rather contrasting trend to that of mastitis, which increased the more an animal was inbred. Overall, our study highlights the importance of considering both pedigree-based and genomic-based inbreeding estimators when assessing the impact on inbreeding, emphasizing that not all inbreeding is harmful.

Using expression data to fine map QTL associated with fertility in dairy cattle

BACKGROUND

Female fertility is an important trait in dairy cattle. Identifying putative causal variants associated with fertility may help to improve the accuracy of genomic prediction of fertility. Combining expression data (eQTL) of genes, exons, gene splicing and allele specific expression is a promising approach to fine map QTL to get closer to the causal mutations. Another approach is to identify genomic differences between cows selected for high and low fertility and a selection experiment in New Zealand has created exactly this resource. Our objective was to combine multiple types of expression data, fertility traits and allele frequency in high- (POS) and low-fertility (NEG) cows with a genome-wide association study (GWAS) on calving interval in Australian cows to fine-map QTL associated with fertility in both Australia and New Zealand dairy cattle populations.

RESULTS

Variants that were significantly associated with calving interval (CI) were strongly enriched for variants associated with gene, exon, gene splicing and allele-specific expression, indicating that there is substantial overlap between QTL associated with CI and eQTL. We identified 671 genes with significant differential expression between POS and NEG cows, with the largest fold change detected for the CCDC196 gene on chromosome 10. Our results provide numerous candidate genes associated with female fertility in dairy cattle, including GYS2 and TIGAR on chromosome 5 and SYT3 and HSD17B14 on chromosome 18. Multiple QTL regions were located in regions with large numbers of copy number variants (CNV). To identify the causal mutations for these variants, long read sequencing may be useful.

CONCLUSIONS

Variants that were significantly associated with CI were highly enriched for eQTL. We detected 671 genes that were differentially expressed between POS and NEG cows. Several QTL detected for CI overlapped with eQTL, providing candidate genes for fertility in dairy cattle.

Inheritance of genomic regions and genes associated with number of oocytes and embryos in Gir cattle through daughter design

Over the past decades, daughter designs, including genotyped sires and their genotyped daughters, have been used as an approach to identify QTL related to economic traits. The aim of this study was to identify genomic regions inherited by Gir sire families and genes associated with number of viable oocytes (VO), total number of oocytes (TO), and number of embryos (EMBR) based on a daughter design approach. In total, 15 Gir sire families were selected. The number of daughters per family ranged from 26 to 395, which were genotyped with different SNP panels and imputed to the Illumina BovineHD BeadChip (777K) and had phenotypes for oocyte and embryo production. Daughters had phenotypic data for VO, TO, and EMBR. The search for QTL was performed through GWAS based on GBLUP. The QTL were found for each trait among and within families based on the top 10 genomic windows with the greatest genetic variance. For EMBR, genomic windows identified among families were located on BTA4, BTA5, BTA6, BTA7, BTA8, BTA13, BTA16, and BTA17, and they were most frequent on BTA7 within families. For VO, genomic windows were located on BTA2, BTA4, BTA5, BTA7, BTA17, BTA21, BTA22, BTA23, and BTA27 among families, being most frequent on BTA8 within families. For TO, the top 10 genomic windows were identified on BTA2, BTA4, BTA5, BTA7, BTA17, BTA21, BTA22, BTA26, and BTA27, being most frequent on BTA7 and BTA8 within families. Considering all results, the greatest number of genomic windows was found on BTA7, where the VCAN, XRCC4, TRNAC-ACA, HAPLN1, and EDIL3 genes were identified in the common regions. In conclusion, 15 Gir sire families with 26 to 395 daughters per family with phenotypes for oocyte and embryo production helped to identify the inheritance of several genomic regions, especially on BTA7, where the EDIL3, HAPLN1, and VCAN candidate genes were associated with number of oocytes and embryos in Gir cattle families.

Effect of postpartum body condition score change on the pregnancy outcomes of lactating Jersey cows inseminated at first service with sexed Jersey or conventional beef semen after a synchronized estrus versus a synchronized ovulation

Our objective was to compare insemination rate and pregnancies per artificial insemination (P/AI) of lactating Jersey cows inseminated at first service with sexed Jersey or conventional beef semen after submission to a Double-Ovsynch protocol for timed artificial insemination (TAI) versus a protocol to synchronize estrus at similar days in milk (DIM). Secondary objectives were to determine the effect of protocol synchrony and postpartum body condition score (BCS) change on P/AI. Lactating Jersey cows (n = 1,272) were allocated by odd versus even ear tag number, which was randomly allocated within the herd, within parity and semen type for submission to a Double-Ovsynch protocol (DO; n = 707) or a protocol to synchronize estrus (ED; n = 565). All ED cows detected in estrus were inseminated (EDAI; n = 424), with undetected cows receiving TAI after an Ovsynch protocol (EDTAI; n = 141). There was a treatment by parity interaction on insemination rate with 100% of DO cows receiving TAI, but a tendency for fewer primiparous ED cows to be detected in estrus and AI than multiparous cows (69.5% ± 0.04% vs. 77.1% ± 0.02%, respectively). For cows inseminated with sexed Jersey or conventional beef semen, DO cows tended to have and had more P/AI than EDAI cows (sexed, 49.2% ± 0.03% vs. 43.6% ± 0.03%; beef, 64.2% ± 0.04% vs. 56.3% ± 0.05%, respectively) and had more P/AI than EDAI+EDTAI cows (sexed, 49.1% ± 0.03% vs. 40.6% ± 0.03%; beef, 65.5% ± 0.04% vs. 56.2% ± 0.04%, respectively). Overall, 29.1% of DO cows expressed estrus with 5.0% and 24.2% of cows detected in estrus ≥24 h before and at TAI, respectively, and there was no difference in P/AI 61 ± 4 d after AI based on expression of estrus at TAI. The synchronization rate was greater for DO than EDAI cows (92.1% ± 0.01% vs. 79.2% ± 0.02%, respectively); however, synchronized DO cows had more P/AI than synchronized EDAI cows (55.0% ± 0.02% vs. 49.2% ± 0.03%, respectively). There was an interaction between BCS change from 7 to 39 ± 2 DIM and treatment on P/AI 61 ± 4 d after AI with no difference between DO and EDAI cows that lost = 0.25 (49.8% ± 0.04% vs. 51.0% ± 0.05%, respectively) or maintained or gained (55.6% ± 0.04% vs. 50.8% ± 0.05%, respectively) BCS, but within cows that lost ≥0.5 BCS, DO cows had more P/AI than EDAI cows (54.1% ± 0.04% vs. 36.1% ± 0.04%, respectively). In conclusion, submission of lactating Jersey cows to a Double-Ovsynch protocol for first insemination increased insemination rate and fertility to first insemination compared with AI after a detected estrus regardless of semen type and expression of estrus, particularly for cows with excessive postpartum BCS loss.

Genetic variation in bovine LAP3 and SIRT1 genes associated with fertility traits in dairy cattle

BACKGROUND

The genetic progress of fertility and reproduction traits in dairy cattle has been constrained by the low heritability of these traits. Identifying candidate genes and variants associated with fertility and reproduction could enhance the accuracy of genetic selection and expedite breeding process of dairy cattle with low-heritability traits. While the bovine LAP3 and SIRT1 genes exhibit well-documented associations with milk production traits in dairy cattle, their effect on cow fertility have not yet been explored. Eleven single nucleotide polymorphisms (SNPs), comprising five in the promoter (rs717156555: C > G, rs720373055: T > C, rs516876447: A > G, rs461857269: C > T and rs720349928: G > A), two in 5'UTR (rs722359733: C > T and rs462932574: T > G), two in intron 12 (rs110932626: A > G and rs43702363: C > T), and one in 3'UTR of exon 13 (rs41255599: C > T) in LAP3 and one in SIRT1 (rs718329990:T > C) genes, have previously been reported to be associated with various traits of milk production and clinical mastitis in Sahiwal and Karan Fries dairy cattle. In this study, the analysis primarily aimed to assess the impact of SNPs within LAP3 and SIRT1 genes on fertility traits in Sahiwal and Karan Fries cattle. Association studies were conducted using mixed linear models, involving 125 Sahiwal and 138 Karan Fries animals in each breed. The analysis utilized a designated PCR-RFLP panel.

RESULTS

In the promoter region of the LAP3 gene, all variants demonstrated significant (P < 0.05) associations with AFC, except for rs722359733: C > T. However, specific variants with the LAP3 gene's promoter region, namely rs722359733: C > T, rs110932626: A > G, rs43702363: C > T, and rs41255599: C > T, showed significant associations with CI and DO in Sahiwal and Karan Fries cows, respectively. The SNP rs718329990: T > C in the promoter region of SIRT1 gene exhibited a significant association with CI and DO in Sahiwal cattle. Haplotype-based association analysis revealed significant associations between haplotype combinations and AFC, CI and DO in the studied dairy cattle population. Animals with H2H3 and H2H4 haplotype combination exhibited higher AFC, CI and DO than other combinations.

CONCLUSIONS

These results affirm the involvement of the LAP3 and SIRT1 genes in female fertility traits, indicating that polymorphisms within these genes are linked to the studied traits. Overall, the significant SNPs and haplotypes identified in this study could have the potential to enhance herd profitability and ensure long-term sustainability on dairy farms by enabling the selection of animals with early age first calving and enhance reproductive performance in the dairy cattle breeding program.

Effects of maternal antral follicle count in Bos taurus indicus cattle on the genetic merit of male offspring and antral follicle count of female offspring

This study evaluated the effects of antral follicle count (AFC) in female cattle on offspring characteristics. Recently calved multiparous Bos indicus cows (Nelore; n = 222) were evaluated using ultrasonography on random days of their estrous cycle to determine the AFC and were classified into "low" (≤15 follicles), "intermediate" (≥16 and ≤ 29 follicles), and "high" (≥30 follicles) AFC groups. Weight and scrotal circumference (SC) of male offspring from these cows (n = 127) were determined from 20 to 27 months, and the data were added to a genetic evaluation program (economic total genetic merit, MGTe and TOP value) that uses the kinship matrix to evaluate the genetic relationship between animals. The AFC of female offspring from these cows (n = 95) was evaluated to analyze the relationship between the AFC of mothers and daughters. The effects of maternal AFC on the genetic merit of male and female offspring were analyzed using GLIMMIX and GLM, respectively. Correlations were assessed using the Pearson's coefficient. Male offspring of cows with high AFC had superior MGTe (P = 0.005) and TOP values (P = 0.01) than those from cows with low AFC. Additionally, the AFC of mothers was positively correlated with MGTe (R = 0.33; P < 0.0001) and negatively correlated with TOP values (R = -0.32; P < 0.0001). The SC (P = 0.01), but not body weight of the offspring (P = 0.46) was affected by maternal AFC. The daughters' AFC were correlated (R = 0.29; P = 0.004) with mothers' AFC and were influenced by maternal (P = 0.05) but not paternal (P = 0.77) effect. In conclusion, cows with high AFC produced males with greater MGTe, superior TOP values and higher SC. Maternal AFC did not influence the weight of male offspring but was correlated with the AFC of daughters.

Multi-omics integration identifies regulatory factors underlying bovine subclinical mastitis

BACKGROUND

Mastitis caused by multiple factors remains one of the most common and costly disease of the dairy industry. Multi-omics approaches enable the comprehensive investigation of the complex interactions between multiple layers of information to provide a more holistic view of disease pathogenesis. Therefore, this study investigated the genomic and epigenomic signatures and the possible regulatory mechanisms underlying subclinical mastitis by integrating RNA sequencing data (mRNA and lncRNA), small RNA sequencing data (miRNA) and DNA methylation sequencing data of milk somatic cells from 10 healthy cows and 20 cows with naturally occurring subclinical mastitis caused by Staphylococcus aureus or Staphylococcus chromogenes.

RESULTS

Functional investigation of the data sets through gene set analysis uncovered 3458 biological process GO terms and 170 KEGG pathways with altered activities during subclinical mastitis, provided further insights into subclinical mastitis and revealed the involvement of multi-omics signatures in the altered immune responses and impaired mammary gland productivity during subclinical mastitis. The abundant genomic and epigenomic signatures with significant alterations related to subclinical mastitis were observed, including 30,846, 2552, 1276 and 57 differential methylation haplotype blocks (dMHBs), differentially expressed genes (DEGs), lncRNAs (DELs) and miRNAs (DEMs), respectively. Next, 5 factors presenting the principal variation of differential multi-omics signatures were identified. The important roles of Factor 1 (DEG, DEM and DEL) and Factor 2 (dMHB and DEM), in the regulation of immune defense and impaired mammary gland functions during subclinical mastitis were revealed. Each of the omics within Factors 1 and 2 explained about 20% of the source of variation in subclinical mastitis. Also, networks of important functional gene sets with the involvement of multi-omics signatures were demonstrated, which contributed to a comprehensive view of the possible regulatory mechanisms underlying subclinical mastitis. Furthermore, multi-omics integration enabled the association of the epigenomic regulatory factors (dMHBs, DELs and DEMs) of altered genes in important pathways, such as 'Staphylococcus aureus infection pathway' and 'natural killer cell mediated cytotoxicity pathway', etc., which provides further insights into mastitis regulatory mechanisms. Moreover, few multi-omics signatures (14 dMHBs, 25 DEGs, 18 DELs and 5 DEMs) were identified as candidate discriminant signatures with capacity of distinguishing subclinical mastitis cows from healthy cows.

CONCLUSION

The integration of genomic and epigenomic data by multi-omics approaches in this study provided a better understanding of the molecular mechanisms underlying subclinical mastitis and identified multi-omics candidate discriminant signatures for subclinical mastitis, which may ultimately lead to the development of more effective mastitis control and management strategies.

Bovine FRAS1: mRNA Expression Profile, Genetic Variations, and Significant Correlations with Ovarian Morphological Traits, Mature Follicle, and Corpus Luteum

The amelioration of bovine fertility caused by a multi-factorial problem has always been a hot topic, among which the detection of available target genes is the most crucial. It was hypothesized that the Fraser extracellular matrix complex subunit 1 (FRAS1) gene detected by GWAS is involved in physiological activities such as ovarian development. Herein, unilateral ovaries from 2111 cows were used to examine the mRNA expression profile and polymorphisms of bovine FRAS1 and their associations with fertility-related characteristics. Firstly, it was confirmed that FRAS1 gene transcripts are expressed in various bovine tissues. Then, among five potential insertion-deletion (indel) loci, the 20 bp (named P3-D20-bp) and 15 bp (P4-D15-bp) deletion mutations were confirmed to be polymorphic with linkage equilibrium. Secondly, the P3-D20-bp polymorphism was significantly associated with ovarian weight and corpus luteum diameter in the metaestrus phase and ovarian length in the dioestrum stage. Additionally, both ovarian length and mature follicle diameter in metaestrus are significantly correlated with different genotypes of P4-D15-bp. Thirdly, the transcriptional expression of the FRAS1 gene in groups with a minimum value of ovarian weight or volume was significantly higher than the expression in groups with a maximum value. Instead of that, the more corpus luteum and mature follicles there are, the higher the transcription expression of the FRAS1 gene is. Furthermore, FRAS1 expression in cows with a heterozygous genotype (ID) of P3-D20-bp was significantly higher than others. Eventually, P3-D20-bp deletion could disturb the binding efficiency of WT1-I and Sox2 to FRAS1 sequence according to binding prediction, indicating that mutation may affect gene expression and traits by influencing the binding of transcription factors. Overall, the polymorphisms of P3-D20-bp and P4-D15-bp of the bovine FRAS1 gene significantly correlated to follicle or ovarian traits that could be applied in optimizing female fertility in cow MAS breeding programs.

Quantitative trait loci exploration and characterization of gestation length in Holstein cattle

Gestation length (GL) is a moderately heritable trait in cattle with economic and management implications. This study aimed to characterize the gestation length of an Argentinian Holstein cattle population, understand contributing factors, and explore the GL effect on production performance. Further objectives were to estimate direct and maternal heritabilities for this trait and to identify genomic regions affecting it. Data consisted of GL records from 45,738 births corresponding to 17,004 Holstein cows and heifers. The effects of age and calving season over GL were analyzed using a Student's t-test for homoscedastic samples. The effects of the GL category (GL shorter than 1.5 SD, within ±1.5 SD, and longer than 1.5 SD from the mean) on production performance were studied by analysis of variance. A single-step genome-wide association study was performed using the BLUPF90 suite of programs with genotypes from 654 Holstein animals on 40,339 SNP. The results showed that the younger the age at calving, the shorter the GL. Moreover, gestations ending in warmer seasons were, in general, statistically shorter than those ending in colder seasons for both heifers and cows. Regarding the effect of GL on production performance, cows with gestation periods within ±1.5 SD from the population mean exhibited the highest 305-day cumulative milk, fat, and protein productions. Direct and maternal heritabilities for GL were 0.42 and 0.03, respectively. We detected a SNP suggestively associated with direct gestation length at 57.7 Mb on Bos taurus autosome 18, a locus included in a region described in the literature as associated with the trait. The information obtained on the environmental and genetic factors affecting GL in Argentinian Holstein cows contributes to characterizing the population in pursuit of improving the performance of national dairy cattle breeding systems.

Genetic parameters for reproductive traits in the Italian Mediterranean buffalo using milk yield as a correlated trait

Until now, the genetic evaluation of the Italian Mediterranean Buffalo has been mainly focused on production traits. However, female fertility affects the efficiency of the dairy industry as it is essential to maintain the profitability of dairy farms. Indeed, the estimation of its genetic component is crucial for its improvement. In this study, 3 measures of buffalo's fertility were analyzed: the age at first calving (AFC), the interval between first and second calving (CIV1), and the interval between second and successive calvings (CIV2_12). Milk yield at 270 d (MY270) was used as a correlated trait. First, genetic parameters were estimated using 7,915 buffalo cows with first calving from 1991 to 2018, then breeding values were calculated from 236,087 buffalo cows. Genetic parameters were estimated by Bayesian inference fitting a multiple-trait animal model using the GIBBS1F90 program, and BLUPF90 was used for estimation of breeding value. The heritability and repeatability estimates of fertility traits were low. The genetic correlations among fertility traits ranged from 0.10 (AFC-CIV1) to 0.92 (CIV1-CIV2_12). Genetic correlation between MY270 and fertility traits was unfavorable, ranging from 0.23 to 0.48. The results from this study can be used as a basis for the future genetic improvement of fertility traits in the Italian Mediterranean Buffaloes.

k-mer-based GWAS enhances the discovery of causal variants and candidate genes in soybean

Genome-wide association studies (GWAS) are powerful statistical methods that detect associations between genotype and phenotype at genome scale. Despite their power, GWAS frequently fail to pinpoint the causal variant or the gene controlling a given trait in crop species. Assessing genetic variants other than single-nucleotide polymorphisms (SNPs) could alleviate this problem. In this study, we tested the potential of structural variant (SV)- and k-mer-based GWAS in soybean by applying these methods as well as conventional SNP/indel-based GWAS to 13 traits. We assessed the performance of each GWAS approach based on loci for which the causal genes or variants were known from previous genetic studies. We found that k-mer-based GWAS was the most versatile approach and the best at pinpointing causal variants or candidate genes. Moreover, k-mer-based analyses identified promising candidate genes for loci related to pod color, pubescence form, and resistance to Phytophthora sojae. In our dataset, SV-based GWAS did not add value compared to k-mer-based GWAS and may not be worth the time and computational resources invested. Despite promising results, significant challenges remain regarding the downstream analysis of k-mer-based GWAS. Notably, better methods are needed to associate significant k-mers with sequence variation. Our results suggest that coupling k-mer- and SNP/indel-based GWAS is a powerful approach for discovering candidate genes in crop species.

Genome-wide association analysis of heifer livability and early first calving in Holstein cattle

BACKGROUND

The survival and fertility of heifers are critical factors for the success of dairy farms. The mortality of heifers poses a significant challenge to the management and profitability of the dairy industry. In dairy farming, achieving early first calving of heifers is also essential for optimal productivity and sustainability. Recently, Council on Dairy Cattle Breeding (CDCB) and USDA have developed new evaluations of heifer health and fertility traits. However, the genetic basis of these traits has yet to be thoroughly studied.

RESULTS

Leveraging the extensive U.S dairy genomic database maintained at CDCB, we conducted large-scale GWAS analyses of two heifer traits, livability and early first calving. Despite the large sample size, we found no major QTL for heifer livability. However, we identified a major QTL in the bovine MHC region associated with early first calving. Our GO analysis based on nearby genes detected 91 significant GO terms with a large proportion related to the immune system. This QTL in the MHC region was also confirmed in the analysis of 27 K bull with imputed sequence variants. Since these traits have few major QTL, we evaluated the genome-wide distribution of GWAS signals across different functional genomics categories. For heifer livability, we observed significant enrichment in promotor and enhancer-related regions. For early calving, we found more associations in active TSS, active Elements, and Insulator. We also identified significant enrichment of CDS and conserved variants in the GWAS results of both traits. By linking GWAS results and transcriptome data from the CattleGTEx project via TWAS, we detected four and 23 significant gene-trait association pairs for heifer livability and early calving, respectively. Interestingly, we discovered six genes for early calving in the Bovine MHC region, including two genes in lymph node tissue and one gene each in blood, adipose, hypothalamus, and leukocyte.

CONCLUSION

Our large-scale GWAS analyses of two heifer traits identified a major QTL in the bovine MHC region for early first calving. Additional functional enrichment and TWAS analyses confirmed the MHC QTL with relevant biological evidence. Our results revealed the complex genetic basis of heifer health and fertility traits and indicated a potential connection between the immune system and reproduction in cattle.

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.

A comprehensive characterization of longevity and culling reasons in Canadian Holstein cattle based on various systematic factors

The decision of premature culling cows directly impacts the profitability of dairy farms. A comprehensive characterization of the primary causes of culling reasons would greatly improve both management and selection objectives in dairy cattle breeding programs. Therefore, this study aimed to analyze the temporal frequencies of 34 culling reasons in Canadian Holstein cows. After data editing and quality control, records from 3,096,872 cows culled from 9,683 herds spread across Canada were used for the analyses covering the periods from 1996 to 2020. Reproductive issues were the main culling reason accounting for 23.02%, followed by milk production (20.82%), health (20.39%), conformation problems (13.69%), economic factors (13.10%), accidents (5.67%), age-related causes (1.67%), and workability (1.63%). Nearly fifty-eight percent of cows were culled after 47 months of age. The observed frequencies of culling due to economic factors were lower than expected from 1996 to 2014 and higher than expected between 2015 and 2020. Reproduction issues had the highest culling frequencies during fall (24.54%), winter (24.02%), and spring (22.51%), while health issues were the most frequent (22.51%) culling reason in the summer season. Health issues (25.50%) and milk production (27.71%) were the most frequent culling reasons in the provinces of Quebec and Ontario, respectively. Reproductive issues showed the highest frequency across climates based on the Köppen climate classification, except for Csb (Dry-summer subtropical or Mediterranean climate) and Bsk (Middle latitude steppe climate), which correspond to small regions in Canada, where production was the most frequent culling reason (29.42% and 21.56%, respectively). Reproductive and milk performance issues were the two main culling reasons in most ecozones, except in Boreal Shield and Atlantic Marine, where health issues had the highest frequencies (25.12 and 23.75%, respectively). These results will contribute to improving management practices and selective decisions to reduce involuntary culling of Holstein cows.

Exploring molecular dynamic indicators associated with reproductive performance of Bos indicus cattle in blood plasma samples through data-independent acquisition mass spectrometry

Improving reproductive performance of cattle is of paramount importance for sustainable dairy farming. Poor reproduction performance (RP) hinders the genetic improvement of important Bos indicus cattle breeds. It is well known that incorporation of molecular information along with conventional breeding method is far better than use of conventional method alone for the genetic improvement of reproductive performance traits in cattle. Therefore, the present study sought to investigate the plasma proteome of the Deoni cows in cyclical (n = 6) and pregnant (n = 6) reproductive phases with varying reproductive performance (high and low). High-throughput data independent acquisition (DIA) based proteomics was performed to understand corresponding proteome. We identified a total of 430 plasma proteins. Among cyclic cows, twenty proteins were differentially regulated in low RP as compared to high RP. BARD1 and AFP proteins were observed upregulated in cyclical cows whose upregulation reported to affect reproductive performance in cattle. Among the pregnant cows, thirty-five proteins were differentially regulated, including the downregulation of FGL2 and ZNFX1 that modulates the maternal immune response mechanism which is required for successful implantation of the embryo. Also, proteins such as AHSG, CLU and SERPINA6 were upregulated in the pregnant cows whose upregulation reported to reduced reproductive performance. The results of this study will be helpful in establishing a framework for future research on the aspect of improving reproductive performance in Bos indicus cattle breeds. SIGNIFICANCE: The Indian subcontinent is the center of domestication for Bos indicus cattle breeds and they are known for their disease resistance, heat tolerance, ability to survive in low input regime and harsh climatic conditions. In recent times, population of many important Bos indicus breeds including Deoni cattle is declining due to various factors, especially due to reproductive performance. Traditional breeding methods are not sufficient enough to understand and improve the reproductive performance traits in important Bos indicus cattle breeds. Proteomics approach is a promising technology to understand the complex biological factors which leads to poor reproductive performance in cattle. The present study utilized DIA based LC- MS/MS analysis to identify the plasma proteins associated with reproductive performance in cyclical and pregnant cows. This study if improved further, can be used to develop potential protein markers associated with reproductive performance which is useful for the selection and genetic improvement of important Bos indicus breeds.

QTLs and Candidate Genes Associated with Semen Traits in Merino Sheep

Ram semen traits play a significant role in conception outcomes, which in turn may influence reproductive efficiency and the overall productivity and profitability of sheep enterprises. Since hundreds of ewes may be inseminated from a single ejaculate, it is important to evaluate semen quality prior to use in sheep breeding programs. Given that semen traits have been found to be heritable, genetic variation likely contributes to the variability observed in these traits. Identifying such genetic variants could provide novel insights into the molecular mechanisms underlying variability in semen traits. Therefore, this study aimed to identify quantitative trait loci (QTLs) associated with semen traits in Merino sheep. A genome-wide association study (GWAS) was undertaken using 4506 semen collection records from 246 Merino rams collected between January 2002 and May 2021. The R package RepeatABEL was used to perform a GWAS for semen volume, gross motility, concentration, and percent post-thaw motility. A total of 35 QTLs, located on 16 Ovis aries autosomes (OARs), were significantly associated with either of the four semen traits in this study. A total of 89, 95, 33, and 73 candidate genes were identified, via modified Bonferroni, within the QTLs significantly associated with volume, gross motility, concentration, and percent post-thaw motility, respectively. Among the candidate genes identified, SORD, SH2B1, and NT5E have been previously described to significantly influence spermatogenesis, spermatozoal motility, and high percent post-thaw motility, respectively. Several candidate genes identified could potentially influence ram semen traits based on existing evidence in the literature. As such, validation of these putative candidates may offer the potential to develop future strategies to improve sheep reproductive efficiency. Furthermore, Merino ram semen traits are lowly heritable (0.071-0.139), and thus may be improved by selective breeding.

Single-step genome-wide association studies and post-GWAS analyses for the number of oocytes and embryos in Gir cattle

Genome-Wide Association Studies (GWAS) are used for identification of quantitate trait loci (QTL) and genes associated with several traits. We aimed to identify genomic regions, genes, and biological processes associated with number of total and viable oocytes, and number of embryos in Gir dairy cattle. A dataset with 17,526 follicular aspirations, including the following traits: number of viable oocytes (VO), number of total oocytes (TO), and number of embryos (EMBR) from 1641 Gir donors was provided by five different stock farms. A genotype file with 2093 animals and 395,524 SNP markers was used to perform a single-step GWAS analysis for each trait. The top 10 windows with the highest percentage of additive genetic variance explained by 100 adjacent SNPs were selected. The genomic regions identified in our work were overlapped with QTLs from QTL database on chromosomes 1, 2, 5, 6, 7, 8, 9, 13, 17, 18, 20, 21, 22, 24, and 29. These QTLs were classified as External, Health, Meat and carcass, Production or Reproduction traits, and about 38% were related to Reproduction. In total, 117 genes were identified, of which 111 were protein-coding genes. Exclusively associations were observed for 42 genes with EMBR, and 1 with TO. Also, 42 genes were in common between VO and TO, 28 between VO and EMBR, and four genes were in common among all traits. In conclusion, great part of the identified genes plays a functional role in initial embryo development or general cell functions. The protein-coding genes ARNT, EGR1, HIF1A, AHR, and PAX2 are good markers for the production of oocytes and embryos in Gir cattle.

A Million-Cow Genome-Wide Association Study of Three Fertility Traits in U.S. Holstein Cows

A genome-wide association study (GWAS) of the daughter pregnancy rate (DPR), cow conception rate (CCR), and heifer conception rate (HCR) using 1,001,374-1,194,736 first-lactation Holstein cows and 75,140-75,295 SNPs identified 7567, 3798, and 726 additive effects, as well as 22, 27, and 25 dominance effects for DPR, CCR, and HCR, respectively, with log10(1/p) > 8. Most of these effects were new effects, and some new effects were in or near genes known to affect reproduction including GNRHR, SHBG, and ESR1, and a gene cluster of pregnancy-associated glycoproteins. The confirmed effects included those in or near the SLC4A4-GC-NPFFR2 and AFF1 regions of Chr06 and the KALRN region of Chr01. Eleven SNPs in the CEBPG-PEPD-CHST8 region of Chr18, the AFF1-KLHL8 region of Chr06, and the CCDC14-KALRN region of Chr01 with sharply negative allelic effects and dominance values for the recessive homozygous genotypes were recommended for heifer culling. Two SNPs in and near the AGMO region of Chr04 that were sharply negative for HCR and age at first calving, but slightly positive for the yield traits could also be considered for heifer culling. The results from this study provided new evidence and understanding about the genetic variants and genome regions affecting the three fertility traits in U.S. Holstein cows.

Review: Diving into the cow hologenome to reduce methane emissions and increase sustainability

Interest on methane emissions from livestock has increased in later years as it is an anthropogenic greenhouse gas with an important warming potential. The rumen microbiota has a large influence on the production of enteric methane. Animals harbour a second genome consisting of microbes, collectively referred to as the "microbiome". The rumen microbial community plays an important role in feed digestion, feed efficiency, methane emission and health status. This review recaps the current knowledge on the genetic control that the cow exerts on the rumen microbiota composition. Heritability estimates for the rumen microbiota composition range between 0.05 and 0.40 in the literature, depending on the taxonomical group or microbial gene function. Variables depicting microbial diversity or aggregating microbial information are also heritable within the same range. This study includes a genome-wide association analysis on the microbiota composition, considering the relative abundance of some microbial taxa previously associated to enteric methane in dairy cattle (Archaea, Dialister, Entodinium, Eukaryota, Lentisphaerae, Methanobrevibacter, Neocallimastix, Prevotella and Stentor). Host genomic regions associated with the relative abundance of these microbial taxa were identified after Benjamini-Hoschberg correction (Padj < 0.05). An in-silico functional analysis using FUMA and DAVID online tools revealed that these gene sets were enriched in tissues like brain cortex, brain amigdala, pituitary, salivary glands and other parts of the digestive system, and are related to appetite, satiety and digestion. These results allow us to have greater knowledge about the composition and function of the rumen microbiome in cattle. The state-of-the art strategies to include methane traits in the selection indices in dairy cattle populations is reviewed. Several strategies to include methane traits in the selection indices have been studied worldwide, using bioeconomical models or economic functions under theoretical frameworks. However, their incorporation in the breeding programmes is still scarce. Some potential strategies to include methane traits in the selection indices of dairy cattle population are presented. Future selection indices will need to increase the weight of traits related to methane emissions and sustainability. This review will serve as a compendium of the current state of the art in genetic strategies to reduce methane emissions in dairy cattle.

Body and milk production traits as indicators of energy requirements and efficiency of purebred Holstein and 3-breed rotational crossbred cows from Viking Red, Montbéliarde, and Holstein sires

This study aimed to compare rotational 3-breed crossbred cows of Viking Red, Montbéliarde, and Holstein breeds with purebred Holstein cows for a range of body measurements, as well as different metrics of the cows' productivity and production efficiency. The study involved 791 cows (440 crossbreds and 351 purebreds), that were managed across 2 herds. Within each herd, crossbreds and purebreds were reared and milked together, fed the same diets, and managed as one group. The heart girth, height at withers, and body length were measured, and body condition score (BCS) was determined on all the cows on a single test day. The body weight (BW) of 225 cows were used to develop an equation to predict BW from body size traits, parity, and days in milk, which was then used to estimate the BW of all the cows. Equations from the literature were used to estimate body protein and lipid contents using the predicted BW and BCS. Evidence suggests that maintenance energy requirements may be closely related to body protein mass, and Holstein and crossbred cows may be different in body composition. Therefore, we computed the requirements of net energy for maintenance (NE_M) on the basis either of the metabolic weight (NE_M-MW: 0.418 MJ/kg of metabolic BW) or of the estimated body protein mass according to a coefficient (NE_M-PM: 0.631 MJ/kg body protein mass) computed on the subset comprising the purebred Holstein. On the same day when body measurements were collected, individual test-day milk yield and fat and protein contents were retrieved once from the official Italian milk recording system, and milk was sampled to determine fresh cheese yield. Measures of NE_M were used to scale the production traits. Statistical analyses of all variables included the fixed effects of herd, days in milk, parity, and genetic group (purebred Holstein and crossbred), and the herd × genetic group interaction. External validation of the equation predicting BW yielded a correlation coefficient of 0.94 and an average bias of -4.95 ± 36.81 kg. The crossbreds had similar predicted BW and NE_M-MW compared with the Holsteins. However, NE_M-PM of crossbreds was 3.8% lower than that of the Holsteins, due to their 11% greater BCS and different estimated body composition. The crossbred cows yielded 4.8% less milk and 3.4% less milk energy than the purebred Holsteins. However, the differences between genetic groups were no longer significant when the production traits were scaled on NE_M-PM, suggesting that the crossbreds and purebreds have the same productive ability and efficiency per unit of body protein mass. In conclusion, measures of productivity and efficiency that combine the cows' production capability with traits related to body composition and the energy cost of production seem to be more effective criteria for comparing crossbred and purebred Holstein cows than just milk, fat, and protein yields.

Infrared Thermography as a Potential Non-Invasive Tool for Estrus Detection in Cattle and Buffaloes

The productivity of dairy animals has significantly increased over the past few decades due to intense genetic selection. However, the enhanced yield performance of milk animals caused a proportional increase in stress and compromised reproductive efficiency. Optimal reproductive performance is mandatory for the sustainable production of dairy animals. Reproductive efficiency is marked by proper estrus detection and precise breeding to achieve maximum pregnancies. The existing conventional methods of estrus detection are somewhat labor intensive and less efficient. Similarly, the modern automated methods that rely on detecting physical activity are expensive, and their efficiency is affected by factors such as type of housing (tie stall), flooring, and environment. Infrared thermography has recently emerged as a technique that does not depend on monitoring physical activity. Furthermore, infrared thermography is a non-invasive, user-friendly, and stress-free option that aids in the detection of estrus in dairy animals. Infrared thermography has the potential to be considered a useful non-invasive tool for detecting temperature fluctuations to generate estrus alerts without physical contact in cattle and buffaloes. This manuscript highlights the potential use of infrared thermography to understand reproductive physiology and practical implementation of this technique through discussing its advantages, limitations, and possible precautions.

Genome-Wide Association Study of Age at First Calving in U.S. Holstein Cows

A genome-wide association study (GWAS) of age at first calving (AFC) using 813,114 first lactation Holstein cows and 75,524 SNPs identified 2063 additive effects and 29 dominance effects with p-values < 10-8. Three chromosomes had highly significant additive effects in the regions of 7.86-8.12 Mb of Chr15, 27.07-27.48 Mb and 31.25-32.11 Mb of Chr19, and 26.92-32.60 Mb of Chr23. Two of the genes in those regions were reproductive hormone genes with known biological functions that should be relevant to AFC, the sex hormone binding globulin (SHBG) gene, and the progesterone receptor (PGR) gene. The most significant dominance effects were near or in EIF4B and AAAS of Chr05 and AFF1 and KLHL8 of Chr06. All dominance effects were positive overdominance effects where the heterozygous genotype had an advantage, and the homozygous recessive genotype of each SNP had a very negative dominance value. Results from this study provided new evidence and understanding about the genetic variants and genome regions affecting AFC in U.S. Holstein cows.

Genomic Regions Associated with Milk Composition and Fertility Traits in Spring-Calved Dairy Cows in New Zealand

The objective of this study was to identify genomic regions and genes that are associated with the milk composition and fertility traits of spring-calved dairy cows in New Zealand. Phenotypic data from the 2014–2015 and 2021–2022 calving seasons in two Massey University dairy herds were used. We identified 73 SNPs that were significantly associated with 58 potential candidate genes for milk composition and fertility traits. Four SNPs on chromosome 14 were highly significant for both fat and protein percentages, and the associated genes were DGAT1, SLC52A2, CPSF1, and MROH1. For fertility traits, significant associations were detected for intervals from the start of mating to first service, the start of mating to conception, first service to conception, calving to first service, and 6-wk submission, 6-wk in-calf, conception to first service in the first 3 weeks of the breeding season, and not in calf and 6-wk calving rates. Gene Ontology revealed 10 candidate genes (KCNH5, HS6ST3, GLS, ENSBTAG00000051479, STAT1, STAT4, GPD2, SH3PXD2A, EVA1C, and ARMH3) that were significantly associated with fertility traits. The biological functions of these genes are related to reducing the metabolic stress of cows and increasing insulin secretion during the mating period, early embryonic development, foetal growth, and maternal lipid metabolism during the pregnancy period.

Runs of homozygosity and signatures of selection for number of oocytes and embryos in the Gir Indicine cattle

Runs of homozygosity (ROH) and signatures of selection are the results of selection processes in livestock species that have been shown to affect several traits in cattle. The aim of the current work was to verify the profile of ROH and inbreeding depression in the number of total (TO) and viable oocytes (VO) and the number of embryos (EMBR) in Gir Indicine cattle. In addition, we aim to identify signatures of selection, genes, and enriched regions between Gir subpopulations sorted by breeding value for these traits. The genotype file contained 2093 animals and 420,718 SNP markers. Breeding values used to sort Gir animals were previously obtained. ROH and signature of selection analyses were performed using PLINK software, followed by ROH-based (F_ROH) and pedigree-based inbreeding (F_ped) and a search for genes and their functions. An average of 50 ± 8.59 ROHs were found per animal. ROHs were separated into classes according to size, ranging from 1 to 2 Mb (ROH_1-2Mb: 58.17%), representing ancient inbreeding, ROH_2-4Mb (22.74%), ROH_4-8Mb (11.34%), ROH_8-16Mb (5.51%), and ROH_>16Mb (2.24%). Combining our results, we conclude that the increase in general F_ROH and F_ped significantly decreases TO and VO; however, in different chromosomes traits can increase or decrease with F_ROH. In the analysis for signatures of selection, we identified 15 genes from 47 significant genomic regions, indicating differences in populations with high and low breeding value for the three traits.

Quantitative trait locus for calving traits on Bos taurus autosome 18 in Holstein cattle is embedded in a complex genomic region

Although the quantitative trait locus (QTL) on chromosome 18 (BTA18) associated with paternal calving ease and stillbirth in Holstein Friesian cattle and its cross has been known for over 20 years, to our knowledge, the exact causal genetic sequence has yet escaped identification. The aim of this study was to re-examine the region of the published QTL on BTA18 and to investigate the possible reasons behind this elusiveness. For this purpose, we carried out a combined linkage disequilibrium and linkage analysis using genotyping data of 2,697 German Holstein Friesian (HF) animals and subsequent whole-genome sequencing (WGS) data analyses and genome assembly of HF samples. We confirmed the known QTL in the 95% confidence interval of 1.089 Mbp between 58.34 and 59.43 Mbp on BTA18. Additionally, these 4 SNPs in the near-perfect linkage disequilibrium with the QTL haplotype were identified: rs381577268 (on 57,816,137 bp, C/T), rs381878735 (on 59,574,329 bp, A/T), rs464221818 (on 59,329,176 bp, C/T), and rs472502785 (on 59,345,689 bp, T/C). Search for the causal mutation using short and long-read sequences, and methylation data of the BTA18 QTL region did not reveal any candidates though. The assembly showed problems in the region, as well as an abundance of segmental duplications within and around the region. Taking the QTL of BTA18 in Holstein cattle as an example, the data presented in this study comprehensively characterize the genomic features that could also be relevant for other such elusive QTL in various other cattle breeds and livestock species as well.

Longitudinal genome-wide association studies of milk production traits in Holstein cattle using whole-genome sequence data imputed from medium-density chip data

Longitudinal traits, such as milk production traits in dairy cattle, are featured by having phenotypic values at multiple time points, which change dynamically over time. In this study, we first imputed SNP chip (50-100K) data to whole-genome sequence (WGS) data in a Chinese Holstein population consisting of 6,470 cows. The imputation accuracies were 0.88 to 0.97 on average after quality control. We then performed longitudinal GWAS in this population based on a random regression test-day model using the imputed WGS data. The longitudinal GWAS revealed 16, 39, and 75 quantitative trait locus regions associated with milk yield, fat percentage, and protein percentage, respectively. We estimated the 95% confidence intervals (CI) for these quantitative trait locus regions using the logP drop method and identified 581 genes involved in these CI. Further, we focused on the CI that covered or overlapped with only 1 gene or the CI that contained an extremely significant top SNP. Twenty-eight candidate genes were identified in these CI. Most of them have been reported in the literature to be associated with milk production traits, such as DGAT1, HSF1, MGST1, GHR, ABCG2, ADCK5, and CSN1S1. Among the unreported novel genes, some also showed good potential as candidate genes, such as CCSER1, CUX2, SNTB1, RGS7, OSR2, and STK3, and are worth being further investigated. Our study provided not only new insights into the candidate genes for milk production traits, but also a general framework for longitudinal GWAS based on random regression test-day model using WGS data.

Major locus on ECA18 influences effectiveness of GonaCon vaccine in feral horses

Contraceptive vaccines are used to reduce birth rates in wild and feral animal populations. While the immunocontraceptive GonaCon-Equine has proven effective in reducing fertility among female feral horses, there is individual variation in the duration of infertility following treatment. To identify genetic factors influencing the effectiveness of GonaCon-Equine, we conducted a genome-wide association study of 88 mares from a feral population genotyped using the Illumina GGP Equine 70k SNP array. Contraceptive treatment schedules and long-term foaling rates have been recorded for each individual. We used mixed linear models to control for relatedness among mares. We found a significant association (p < 5 ×10^-8) with a locus on equine chromosome 18. The most likely candidate genes in this region are STAT1 and STAT4, which are both involved in immune system function. Variation in STAT function could affect the immune response to the vaccine, leading to variation in contraceptive efficacy. Additional SNPs reaching a less stringent threshold of significance (p < 5 ×10^-6) were located on other chromosomes near known immune system genes, supporting the hypothesis that variation in immunocontraceptive efficacy can be attributed to genetic variation in immune response rather than fertility genes.

Genetic analysis of fetal loss in Holstein cattle

Pregnancy loss is recognized as one of the major factors contributing to poor reproductive performance in dairy cattle. Here, we performed a comprehensive genetic analysis of fetal loss, defined as a pregnancy loss that occurs after detection of a viable embryo around 42 d of gestation. The objectives of this study were to reveal (1) whether fetal loss is heritable and, hence, whether it will respond to selection, and (2) to what extent current fertility traits, such as daughter pregnancy rate, are associated with fetal loss. Data consisted of 59,308 confirmed pregnancy or fetal loss records distributed across nulliparous heifers and primiparous and multiparous cows. We defined fetal loss as a binary trait (yes vs. no) or as an ordinal trait (pregnancy maintenance, early fetal loss ≤150 d of gestation, and late fetal loss >150 d of gestation), and we assessed both linear and threshold models. Heritability estimates for fetal loss ranged from 1 to 18%, depending upon parity, trait definition, and statistical model used. Heritability estimates were greater for lactating cows than for nonlactating nulliparous heifers. Threshold models were able to capture more additive genetic variance and, thus, yielded higher heritability estimates than linear models. Notably, fetal loss traits were highly genetically correlated with each other but only weakly correlated with current fertility traits included in the national genetic evaluation. Overall, our study provides evidence that fetal loss is heritable enough to make genetic selection for reducing fetal loss and improving pregnancy maintenance feasible. In addition, our results suggest that fetal loss is largely independent from current traits used to select for cow fertility, and thus current breeding efforts have unfortunately little effect on reducing the incidence of fetal loss.

Peripheral blood mononuclear cell mitochondrial enzyme activity is associated with parity and lactation performance in early lactation Holstein dairy cows

Mitochondria are central to metabolism and are the primary energy producers for all biosynthesis, including lactation. The objectives of this study were to determine if high- and low-producing dairy cows exhibit differences in peripheral blood mononuclear cell mitochondrial enzyme activities of citrate synthase, complex I, complex IV, and complex V during early lactation and, thus, to determine whether those differences were related to differences in lactation performance in the dairy cow. Fifty-six Holstein cows were assigned to 1 of 4 groups: (1) primiparous high, (2) primiparous low, (3) multiparous high, or (4) multiparous low. Primiparous and multiparous cows were analyzed separately. Then, cows were divided into high or low production groups for each production parameter [peak milk, average milk, energy-corrected milk (ECM), fat-corrected milk (FCM), milk lactose, milk fat, milk protein, total solids (TS), solids-not-fat, feed efficiency, and somatic cell count (SCC)]. For all data analysis, production parameters are expressed as yields (kg/d) and SCC (10³ cells/mL). High and low production groups were defined by their respective mean production parameters for the 56 cows, with below average cows defined as low and above average cows defined as high. Whole blood samples were collected at one time point, approximately 70 d in milk at 0800 h, and processed for crude mitochondrial extracts from peripheral blood mononuclear cells to determine the activity rates of mitochondrial enzymes. Milk samples were collected 9 times (3 d, 3 times per d) during the week of blood collection and analyzed for major components (fat, protein, lactose, TS, and SCC). Multiparous cows had lower citrate synthase activity than primiparous cows across all production parameters. High-producing cows had greater complex I activity for peak milk, milk yield, ECM, FCM, milk fat, TS, and feed efficiency, and greater complex V activity for ECM, FCM, milk lactose, milk fat, and TS across parities. These findings imply that the most influential respiratory chain enzymes on the level of milk production are those responsible for electron transport chain initialization and ATP production.

Associations between anogenital distance and measures of fertility in lactating North American Holstein cows: A validation study

Anogenital distance (AGD) has been defined in dairy cows as the distance from the center of the anus to the base of the clitoris. Initial reports on nulliparous Holstein heifers and first- and second-parity Holstein cows have found inverse relationships between AGD and measures of fertility. Our primary objective was to determine the relationship between AGD and measures of fertility in a larger population of North American Holstein cows to validate our previous finding that AGD is inversely related to fertility. Secondary objectives were to determine the associations between AGD and parity, and milk yield. Using digital calipers, we measured AGD in 4,709 Holstein cows [mean ± standard deviation (SD); parity 2.3 ± 1.4; days in milk (DIM) 154 ± 94; 305-d mature equivalent (ME) milk yield 13,759 ± 2,188 kg] from 18 herds in Western Canada and 1 herd in the USA. Anogenital distance (mm) was normally distributed with a mean (±SD) of 132 ± 12, ranging from 95 to 177, and a median of 133. Anogenital distance was linearly but inversely associated with pregnancy to first artificial insemination (P/AI1). For every 1-mm increase in AGD, the estimated probability of P/AI1 decreased by 0.8%. The optimum AGD cut-point that predicted probability of P/AI1 with sensitivity and specificity of 45 and 55%, respectively, was 129 mm. Consequently, data were categorized into either short (≤129) or long (>129) AGD groups across parities, and associations between AGD, parity (first, second, and third+), and fertility measures were determined. Rates of P/AI1 were greater (36 vs. 30%) in short- than in long-AGD cows; short-AGD cows required fewer AI per conception (2.3 vs. 2.4) and had fewer days open (137 vs. 142), and a greater proportion of short-AGD cows (67 vs. 64%) was pregnant by 150 DIM compared with long-AGD cows. The rates of pregnancy up to 150 (hazard ratio of 0.91) and 250 DIM (hazard ratio of 0.93) were smaller in long- than in short-AGD cows. Anogenital distance had a weak positive association with both parity (r = 0.22) and 305-d ME milk yield (r = 0.04). Results indicate an inverse relationship between AGD and measures of fertility in lactating cows, validating our earlier report. We infer that although selecting cows for short AGD is expected to have an adverse effect on milk yield, the anticipated gain in fertility will outweigh the small decline in milk yield, strengthening the potential of AGD as a novel reproductive phenotype for use in future breeding programs to improve fertility.

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.

Genome-Wide Genomic and Functional Association Study for Workability and Calving Traits in Holstein Cattle

The goal of our study was to identify the SNPs, metabolic pathways (KEGG), and gene ontology (GO) terms significantly associated with calving and workability traits in dairy cattle. We analysed direct (DCE) and maternal (MCE) calving ease, direct (DSB) and maternal (MSB) stillbirth, milking speed (MSP), and temperament (TEM) based on a Holstein-Friesian dairy cattle population consisting of 35,203 individuals. The number of animals, depending on the trait, ranged from 22,301 bulls for TEM to 30,603 for DCE. We estimated the SNP effects (based on 46,216 polymorphisms from Illumina BovineSNP50 BeadChip Version 2) using a multi-SNP mixed model. The SNP positions were mapped to genes and the GO terms/KEGG pathways of the corresponding genes were assigned. The estimation of the GO term/KEGG pathway effects was based on a mixed model using the SNP effects as dependent variables. The number of significant SNPs comprised 59 for DCE, 25 for DSB and MSP, 17 for MCE and MSB, and 7 for TEM. Significant KEGG pathways were found for MSB (2), TEM (2), and MSP (1) and 11 GO terms were significant for MSP, 10 for DCE, 8 for DSB and TEM, 5 for MCE, and 3 for MSB. From the perspective of a better understanding of the genomic background of the phenotypes, traits with low heritabilities suggest that the focus should be moved from single genes to the metabolic pathways or gene ontologies significant for the phenotype.

Historical Evolution of Cattle Management and Herd Health of Dairy Farms in OECD Countries

This work aimed to review the important aspects of the dairy industry evolution at herd level, interrelating production with health management systems. Since the beginning of the industrialization of the dairy cattle sector (1950s), driven by the need to feed the rapidly growing urban areas, this industry has experienced several improvements, evolving in management and technology. These changes have been felt above all in the terms of milking, rearing, nutrition, reproductive management, and design of facilities. Shortage of labor, emphasis on increasing farm efficiency, and quality of life of the farmers were the driving factors for these changes. To achieve it, in many areas of the world, pasture production has been abandoned, moving to indoor production, which allows for greater nutritional and reproductive control of the animals. To keep pace with this paradigm in milk production, animal health management has also been improved. Prevention and biosecurity have become essential to control and prevent pathologies that cause great economic losses. As such, veterinary herd health management programs were created, allowing the management of health of the herd as a whole, through the common work of veterinarians and farmers. These programs address the farms holistically, from breeding to nutrition, from prevention to consultancy. In addition, farmers are now faced with a consumer more concerned on animal production, valuing certified products that respect animal health and welfare, as well as environmental sustainability.

A review of inbreeding depression in dairy cattle: current status, emerging control strategies, and future prospects

Dairy cattle breeding has historically focused on relatively small numbers of elite bulls as sires of sons. In recent years, even if generation intervals were reduced and more diverse sires of sons could have been selected, genomic selection has not fundamentally changed the fact that a large number of individuals are being analyzed. However, a relatively small number of elite bulls are still siring those animals. Therefore inbreeding-derived negative consequences in the gene pool have brought concern. The detrimental effects of non-additive genetic changes such as inbreeding depression and dominance have been widely disseminated while seriously affecting bioeconomically important parameters because of an antagonistic relationship between dairy production and reproductive traits. Therefore, the estimation of benefits and limitations of inbreeding and variance of the selection response deserves to be evaluated and discussed to preserve genetic variability, a significant concern in the selection of individuals for reproduction and production. Short-term strategies for genetic merit improvement through modern breeding programs have severely lowered high-producing dairy cattle fertility potential. Since the current selection programs potentially increase long-term costs, genetic diversity has decreased globally as a consequence. Therefore, a greater understanding of the potential that selection programs have for supporting long-term genetic sustainability and genetic diversity among dairy cattle populations should be prioritized in managing farm profitability. The present review provides a broad approach to current inbreeding-derived problems, identifying critical points to be solved and possible alternative strategies to control selection against homozygous haplotypes while maintaining sustained selection pressure. Moreover, this manuscript explores future perspectives, emphasizing theoretical applications and critical points, and strategies to avoid the adverse effects of inbreeding in dairy cattle. Finally, this review provides an overview of challenges that will soon require multidisciplinary approaches to managing dairy cattle populations, intending to combine increases in productive trait phenotypes with improvements in reproductive, health, welfare, linear conformation, and adaptability traits into the foreseeable future.

Genome-wide detection of non-additive quantitative trait loci for semen production traits in beef and dairy bulls

Semen production traits are important aspects of bull fertility, because semen quantity leads to direct profits for artificial insemination centres, and semen quality is associated with the probability of achieving a pregnancy. Most genome-wide association studies (GWASs) for semen production traits have assumed that each quantitative trait locus (QTL) has an additive effect. However, GWASs that account for non-additive effects are also important in fitness traits, such as bull fertility. Here, we performed a GWAS using models that accounted for additive and non-additive effects to evaluate the importance of non-additive effects on five semen production traits in beef and dairy bulls. A total of 65 463 records for 615 Japanese Black bulls (JB) and 50 734 records for 873 Holstein bulls (HOL), which were previously genotyped using the Illumina BovineSNP50 BeadChip, were used to estimate genetic parameters and perform GWAS. The heritability estimates were low (ranged from 0.11 to 0.23), and the repeatability estimates were low to moderate (ranged from 0.28 to 0.45) in both breeds. The estimated repeatability was approximately twice as high as the estimated heritability for all traits. In this study, only one significant region with an additive effect was detected in each breed, but multiple significant regions with non-additive effects were detected for each breed. In particular, the region at approximately 64 Mbp on Bos taurus autosome 17 had the highest significant non-additive effect on four semen production traits in HOL. The rs41843851 single nucleotide polymorphism (SNP) in the region had a much lower P-value for the non-additive effect (P-value = 1.1 × 10^-31) than for the additive effect (P-value = 1.1 × 10^-8) in sperm motility. The AA and AB genotypes on the SNP had a higher phenotype than the BB genotype in HOL, and there was no bull with the BB genotype in JB. Our results showed that non-additive QTLs affect semen production traits, and a novel QTL accounting for non-additive effects could be detected by GWAS. This study provides new insights into non-additive QTLs that affect fitness traits, such as semen production traits in beef and dairy bulls.

Sustainable Intensification of Beef Production in the Tropics: The Role of Genetically Improving Sexual Precocity of Heifers

Simple Summary

Tropical pasture-based beef production systems play a vital role in global food security. The importance of promoting sustainable intensification of such systems has been debated worldwide. Demand for beef is growing together with concerns over the impact of its production on the environment. Implementing sustainable livestock intensification programs relies on animal genetic improvement. In tropical areas, the lack of sexual precocity is a bottleneck for cattle efficiency, directly impacting the sustainability of production systems. In the present review we present and discuss the state of the art of genetic evaluation for sexual precocity in Bos indicus beef cattle, covering the definition of measurable traits, genetic parameter estimates, genomic analyses, and a case study of selection for sexual precocity in Nellore breeding programs.

Abstract

Increasing productivity through continued animal genetic improvement is a crucial part of implementing sustainable livestock intensification programs. In Zebu cattle, the lack of sexual precocity is one of the main obstacles to improving beef production efficiency. Puberty-related traits are complex, but large-scale data sets from different “omics” have provided information on specific genes and biological processes with major effects on the expression of such traits, which can greatly increase animal genetic evaluation. In addition, genetic parameter estimates and genomic predictions involving sexual precocity indicator traits and productive, reproductive, and feed-efficiency related traits highlighted the feasibility and importance of direct selection for anticipating heifer reproductive life. Indeed, the case study of selection for sexual precocity in Nellore breeding programs presented here show that, in 12 years of selection for female early precocity and improved management practices, the phenotypic means of age at first calving showed a strong decreasing trend, changing from nearly 34 to less than 28 months, with a genetic trend of almost −2 days/year. In this period, the percentage of early pregnancy in the herds changed from around 10% to more than 60%, showing that the genetic improvement of heifer’s sexual precocity allows optimizing the productive cycle by reducing the number of unproductive animals in the herd. It has a direct impact on sustainability by better use of resources. Genomic selection breeding programs accounting for genotype by environment interaction represent promising tools for accelerating genetic progress for sexual precocity in tropical beef cattle.

Genomic Estimated Breeding Value of Milk Performance and Fertility Traits in the Russian Black-and-White Cattle Population

A breakthrough in cattle breeding was achieved with the incorporation of animal genomic data into breeding programs. The introduction of genomic selection has a major impact on traditional genetic assessment systems and animal genetic improvement programs. Since 2010, genomic selection has been officially introduced in the evaluation of the breeding and genetic potential of cattle in Europe, the U.S., Canada, and many other developed countries. The purpose of this study is to develop a system for a genomic evaluation of the breeding value of the domestic livestock of Black-and-White and Russian Holstein cattle based on 3 milk performance traits: daily milk yield (kg), daily milk fat (%), and daily milk protein content (%) and 6 fertility traits: age at first calving (AFC), calving interval (CI), calving to first insemination interval (CFI), interval between first and last insemination (IFL), days open (DO), and number of services (NS). We built a unified database of breeding animals from 523 breeding farms in the Russian Federation. The database included pedigree information on 2,551,529 cows and 69,131 bulls of the Russian Holstein and Black-and-White cattle breeds, as well as information on the milk performance of 1,597,426 cows with 4,771,366 completed lactations. The date of birth of the animals included in the database was between 1975 and 2017. Genotyping was performed in 672 animals using a BovineSNP50 v3 DNA Analysis BeadChip microarray (Illumina, USA). The genomic estimated breeding value (GEBV) was evaluated only for 644 animals (427 bulls and 217 cows) using the single-step genomic best linear unbiased prediction - animal model (ssGBLUP-AM). The mean genetic potential was +0.88 and +1.03 kg for the daily milk yield, -0.002% for the milk fat content, and -0.003 and 0.001% for the milk protein content in the cows and bulls, respectively. There was negative genetic progress in the fertility traits in the studied population between 1975 and 2017. The reliability of the estimated breeding value (EBV) for genotyped bulls ranged from 89 to 93% for the milk performance traits and 85 to 90% for the fertility traits, whereas the reliability of the genomic estimated breeding value (GEBV) varied 54 to 64% for the milk traits and 23 to 60% for the fertility traits. This result shows that it is possible to use the genomic estimated breeding value with rather high reliability to evaluate the domestic livestock of Russian Holstein and Black-and-White cattle breeds for fertility and milk performance traits. This system of genomic evaluation may help bring domestic breeding in line with modern competitive practices and estimate the breeding value of cattle at birth based on information on the animal's genome.

Genetic parameters for mid-infrared spectroscopy-predicted fertility

Female fertility is a challenging trait to improve genetically because of its low heritability, its unfavorable genetic correlation with milk yield, and its relatively small number of records. The MFERT trait is the probability of conception to first insemination predicted using mid-infrared (MIR) spectroscopy of a milk sample collected during lactation as part of routine milk recording, age at calving, days in milk, and milk production. As such, MFERT could become available on many more cows than traditional fertility traits. Our objectives were (1) to estimate the heritability of MFERT; (2) to estimate genetic correlations between MFERT, traditional fertility traits, and milk production traits; and (3) to assess the potential of MFERT to be used as an indicator trait for fertility in a selection index. The MFERT trait had a heritability of 0.16, which was higher than that (0.05) obtained for traditional fertility traits. Genetic correlations between MFERT and traditional fertility traits were low to moderate. The weakest and strongest correlations (mean ± standard error) were with pregnancy at the end of the mating season (0.13 ± 0.05) and calving to first service (-0.61 ± 0.03), respectively. Based on our estimates, including MFERT in a fertility index will only substantially increase the accuracy of the index when there are many more records available for MFERT than for the traditional fertility traits. This is likely to be the case because the number of milk samples from commercial machines belonging to milk recording companies in Australia that are capable of generating MIR spectra is growing. Hence, the number of records for MFERT is expected to increase substantially in the near future.

Investigation of β-hydroxybutyrate in early lactation of Simmental cows: Genetic parameters and genomic predictions

Genomic information allows for a more accurate calculation of relationships among animals than the pedigree information, leading to an increase in accuracy of breeding values. Here, we used pedigree-based and single-step genomic approaches to estimate variance components and breeding values for β-hydroxybutyrate milk content (BHB). Additionally, we performed a genome-wide association study (GWAS) to depict its genetic architecture. BHB concentrations within the first 90 days of lactation, estimated from milk medium infrared spectra, were available for 30,461 cows (70,984 records). Genotypes at 42,152 loci were available for 9,123 animals. Low heritabilities were found for BHB using pedigree-based (0.09 ± 0.01) and genomic (0.10 ± 0.01) approaches. Genetic correlation between BHB and milk traits ranged from -0.27 ± 0.06 (BHB and protein percentage) to 0.13 ± 0.07 (BHB and fat-to-protein ratio) using pedigree and from -0.26 ± 0.05 (BHB and protein percentage) to 0.13 ± 0.06 (BHB and fat-to-protein ratio) using genomics. Breeding values were validated for 344 genotyped cows using linear regression method. The genomic EBV (GEBV) had greater accuracy (0.51 vs. 0.45) and regression coefficient (0.98 vs. 0.95) compared to EBV. The correlation between two subsequent evaluations, without and with phenotypes for validation cows, was 0.85 for GEBV and 0.82 for EBV. Predictive ability (correlation between (G)EBV and adjusted phenotypes) was greater when genomic information was used (0.38) than in the pedigree-based approach (0.31). Validation statistics in the pairwise two-trait models (milk yield, fat and protein percentage, urea, fat/protein ratio, lactose and logarithmic transformation of somatic cells count) were very similar to the ones highlighted for the single-trait model. The GWAS allowed discovering four significant markers located on BTA20 (57.5-58.2 Mb), where the ANKH gene is mapped. This gene has been associated with lactose, alpha-lactalbumin and BHB. Results of this study confirmed the usefulness of genomic information to provide more accurate variance components and breeding values, and important insights about the genomic determination of BHB milk content.

Invited review: Use of assisted reproduction techniques to accelerate genetic gain and increase value of beef production in dairy herds

The contribution of the calf enterprise to the profit of the dairy farm is generally considered small, with beef bull selection on dairy farms often not considered a high priority. However, this is likely to change in the future as the rapid rate of expansion of the dairy herd in some countries is set to plateau and improvements in dairy herd fertility combine to reduce the proportion of dairy breed calves required on dairy farms. This presents the opportunity to increase the proportion of beef breed calves born, increasing both the value of calf sales and the marketability of the calves. Beef embryos could become a new breeding tool for dairies as producers need to reassess their breeding policy as a consequence of welfare concerns and poor calf prices. Assisted reproductive technologies can contribute to accelerated genetic gain by allowing an increased number of offspring to be produced from genetically elite dams. There are the following 3 general classes of donor females of interest to an integrated dairy-beef system: (1) elite dairy dams, from which oocytes are recovered from live females using ovum pick-up and fertilized in vitro with semen from elite dairy bulls; (2) elite beef dams, where the oocytes are recovered from live females using ovum pick-up and fertilized with semen from elite beef bulls; and (3) commercial beef dams (≥50% beef genetics), where ovaries are collected from the abattoir postslaughter, and oocytes are fertilized with semen from elite beef bulls that are suitable for use on dairy cows (resulting embryo with ≥75% beef genetics). The expected benefits of these collective developments include accelerated genetic gain for milk and beef production in addition to transformation of the dairy herd calf crop to a combination of good genetic merit dairy female calves and premium-quality beef calves. The aim of this review is to describe how these technologies can be harnessed to intensively select for genetic improvement in both dairy breed and beef breed bulls suitable for use in the dairy herd.

Moving towards sustainable breeding objectives and cow welfare in dairy production: a South African perspective

Genetic advancements have resulted in improved dairy production over many decades, due to the focus of breeding objectives on production as the driving force for genetic progress and overall farm profitability. Major advancements were made in the easy-to-measure traits with moderate to high heritability, which resulted in unintended consequences on herd fertility, health, and welfare of cows. In addition, climate change and animal welfare concerns demanded balanced breeding objectives and selection approaches for sustainable production-including health and longevity. The inclusion of genomic information into genetic evaluations has been proved to benefit traits associated with welfare and sustainable production. Cow welfare traits remain complex and suitable phenotypes are not always easy to measure or readily available for genetic evaluations. The challenge for improvement of cow welfare often lies within implementation of sensitive and measurable parameters. The aim of this review was to explore the reconsideration of breeding objectives in the dairy industry towards sustainable dairy production and cow welfare with reference to selection of dairy animals in South Africa.

Use of insemination data for joint evaluation of male and female fertility in predominantly seasonal-calving dairy herds

Conception in dairy cattle is influenced by the fertility of the cow and the bull and their interaction. Despite genetic selection for female fertility in many countries, selection for male fertility is largely not practiced. The primary objective of this study was to quantify variation in male and female fertility using insemination data from predominantly seasonal-calving herds. Nonreturn rate (NRR) was derived by coding each insemination as successful (1) or failed (0) based on a minimum of at least 25 d. The NRR was treated as a trait of the bull with semen (male fertility) and the cow that is mated (female fertility). The data (805,463 cows that mated to 5,776 bulls) were used to estimate parameters using either models that only included bulls with mating data or models that fitted the genetic and permanent environmental (PE) effects of bulls and cows simultaneously. We also evaluated whether fitting genetic and PE effects of bulls as one term is better for ranking bulls based on NRR compared with a model that ignored genetic effect. The age of cows that were mated, age of the bulls with semen data, season of mating, breed of cow that mated, inbreeding of cows and bulls, and days from calving to mating date were found to have a significant effect on NRR. Only about 3% of the total variance was explained by the random effects in the model, despite fitting the genetic and PE effects of the bull and cow. The 2 components of fertility (male and fertility) were not correlated. The heritability of male fertility was low (0.001 to 0.008), and that of female fertility was also low (~0.016). The highest heritability estimate for male fertility was obtained from the model that fitted the additive genetic relationship matrix and PE component of the bull as one term. When this model was used to calculate bull solutions, the difference between bulls with at least 100 inseminations was up to 19.2% units (-9.6 to 9.6%). Bull solutions from this model were compared with bull solutions that were predicted fitting bull effects ignoring pedigree. Bull solutions that were obtained considering pedigree had (1) the highest accuracy of prediction when early insemination was used to predict yet-to-be observed insemination data of bulls, and (2) improved model stability (i.e., a higher correlation between bull solutions from 2 randomly split herds) compared with the model which fitted bull with no pedigree. For practical purposes, the model that fitted genetic and PE effect as one term can provide more accurate semen fertility values for bulls than the model without genetic effect. To conclude, insemination data from predominantly seasonal-calving herds can be used to quantify variability between bulls for male fertility, which makes their ranking on NRR feasible. Potentially this information can be used for monitoring bulls and can supplement efforts to improve herd fertility by avoiding or minimizing the use of semen from subfertile bulls.

Genetic variations of bovine PCOS-related DENND1A gene identified in GWAS significantly affect female reproductive traits

Genome-wide association studies (GWAS) have identified DENND1A as a potential candidate gene linked to the fertility-related phenotypes in dairy cows. However, to date, no studies have examined the association of the DENND1A insertion/deletions (indels) to bovine fertility on a large scale. Herein, two indel sites, including P4-del-26-bp and P8-ins-15-bp were identified in 1064 Holstein cows. The values of the minor allelic frequency (MAF) ranged between 0.471 (deletion) and 0.230 (deletion), respectively, and combined four different haplotypes by analyzing the haplotype combination. It is noteworthy that P4-del-26-bp is associated with the ovarian width (P = 0.0004) and corpus luteum diameter (P = 0.004). Meanwhile, P8-ins-15-bp was found to have a significant association with the ovarian width (P = 0.020), ovarian weight (P = 0.004), the number of mature follicles (P = 0.020), and diameter of the mature follicles (P = 0.016). Furthermore, the combinatorial analysis showed that the two indel combined-genotypes were significantly related to several reproductive traits (ovarian width, ovarian weight, etc.). Collectively, our findings indicated that these two novel indels and their combinations are correlated with the reproductive traits, and hence, they can serve in the marker-assisted selection (MAS) in cattle breeding. Nevertheless, further functional experiments are needed for understanding the mechanisms of these indels in cattle reproduction in a better way.

Unintended consequences of selection for increased production on the health and welfare of livestock

Modern farming technologies, including quantitative selection and breeding methods in farm animal species, resulted in increased production and efficiency. Selection for increased output in both intensive and extensive production systems has trade-offs and negative outcomes, often more pronounced in intensive systems. Animal welfare and health are often adversely affected and this influences sustainable production. The relative importance of animal welfare differs among developed and developing countries due to the level of economic development, food security and education, as well as religious and cultural practices which presents challenges for sound scientific research. Due to breeding goals in the past set on growth performance, traits such as fertility, welfare and health have been neglected. Fertility is the single most important trait in all livestock species. Reduced fertility and lameness, claw health and mastitis results in unnecessary culling and reduced longevity. Selection pressure for growth accompanied with inbreeding has resulted in a number of genetic defects in beef, sheep and pigs. This review demonstrated the importance of inclusion of animal welfare concepts into breeding objectives and selection strategies. Accurate phenotyping of welfare traits is a limiting factor in the implementation of mitigating strategies, which include diagnostic testing, control of inbreeding and genomic selection.