Fertility in dairy cattle has declined as an unintended consequence of single trait selection for high milk yield. The unfavorable genetic correlation between milk yield and fertility is now well-documented, however, the underlying physiological mechanisms are still uncertain. To understand the relationship between these traits, we developed a method that clusters variants with similar patterns of effects and, after the integration of gene expression data, identifies the genes through which they are likely to act. Biological processes that are enriched in the genes of each cluster were then identified. We identified several clusters with unique patterns of effects. One of the clusters included variants associated with increased milk yield and decreased fertility, where the 'archetypal' variant (i.e., the one with the largest effect) was associated with the gene GC, while others were associated with TRIM32, LRRK2, and U6. These genes have been linked to transcription and alternative splicing, suggesting that these processes are likely contributors to the unfavorable relationship between the 2 traits. Another cluster, with archetypal variant near DGAT1 and including variants associated with CDH2, BTRC, SFRP2, ZFHX3, and SLITRK5, appeared to affect milk yield but have little effect on fertility. These genes have been linked to insulin, adipose tissue, and energy metabolism. A third cluster with archetypal variant near ZNF613 and including variants associated with ROBO1, EFNA5, PALLD, GPC6, and PTPRT were associated with fertility but not milk yield. These genes have been linked to GnRH neuronal migration, embryonic development, and/or ovarian function. The use of archetypal clustering to group variants with similar patterns of effects may assist in identifying the biological processes underlying correlated traits. The method is hypothesis-generating and requires experimental confirmation. However, we have uncovered several novel mechanisms potentially affecting milk production and fertility such as GnRH neuronal migration. We anticipate our method to be a starting point for experimental research into novel pathways which have been previously unexplored within the context of dairy production.
Read full abstract