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.

Single-cell transcriptomic analyses of dairy cattle ruminal epithelial cells during weaning.

Using the 10× Genomics Chromium Controller, we obtained scRNA-seq data of 5064 and 1372 individual cells from two Holstein calf ruminal epithelial tissues before and after weaning, respectively. We detected six distinct cell clusters, designated their cell types, and reported their marker genes. We then examined these clusters' underlining cell types and relationships by performing cell cycle, pseudotime trajectory, regulatory network, weighted gene co-expression network and gene ontology analyses. By integrating these cell marker genes with Holstein GWAS signals, we found they were enriched for animal production and body conformation traits. Finally, we confirmed their cell identities by comparing them with human and mouse stomach epithelial cells. This study presents an initial effort to implement single-cell transcriptomic analysis in cattle, and demonstrates ruminal tissue epithelial cell types and their developments during weaning, opening the door for new discoveries about tissue/cell type roles in complex traits at single-cell resolution.

A single nucleotide polymorphism in COQ9 affects mitochondrial and ovarian function and fertility in Holstein cows.

A single missense mutation at position 159 of coenzyme Q9 (COQ9) (G→A; rs109301586) has been associated with genetic variation in fertility in Holstein cattle, with the A allele associated with higher fertility. COQ9 is involved in the synthesis of coenzyme COQ10, a component of the electron transport system of the mitochondria. Here we tested whether reproductive phenotype is associated with the mutation and evaluated functional consequences for cellular oxygen metabolism, body weight changes, and ovarian function. The mutation in COQ9 modifies predicted tertiary protein structure and affected mitochondrial respiration of peripheral blood mononuclear cells. The A allele was associated with low resting oxygen consumption and high electron transport system capacity. Phenotypic measurements for fertility were evaluated for up to five lactations in a population of 2273 Holstein cows. There were additive effects of the mutation (P < 0.05) in favor of the A allele for pregnancy rate, interval from calving to conception, and services per conception. There was no association of genotype with milk production or body weight changes postpartum. The mutation in COQ9 affected ovarian function; the A allele was associated with increased mitochondrial DNA copy number in oocytes, and there were overdominance effects for COQ9 expression in oocytes, follicle number, and antimullerian hormone concentrations. Overall, results show how a gene involved in mitochondrial function is associated with overall fertility, possibly in part by affecting oocyte quality.

A genome-wide association study for clinical mastitis in first parity US Holstein cows using single-step approach and genomic matrix re-weighting procedure.

Clinical mastitis (CM) is one of the health disorders with large impacts on dairy farming profitability and animal welfare. The objective of this study was to perform a genome-wide association study (GWAS) for CM in first-lactation Holstein. Producer-recorded mastitis event information for 103,585 first-lactation cows were used, together with genotype information on 1,361 bulls from the Illumina BovineSNP50 BeadChip. Single-step genomic-BLUP methodology was used to incorporate genomic data into a threshold-liability model. Association analysis confirmed that CM follows a highly polygenic mode of inheritance. However, 10-adjacent-SNP windows showed that regions on chromosomes 2, 14 and 20 have impacts on genetic variation for CM. Some of the genes located on chromosome 14 (LY6K, LY6D, LYNX1, LYPD2, SLURP1, PSCA) are part of the lymphocyte-antigen-6 complex (LY6) known for its neutrophil regulation function linked to the major histocompatibility complex. Other genes on chromosome 2 were also involved in regulating immune response (IFIH1, LY75, and DPP4), or are themselves regulated in the presence of specific pathogens (ITGB6, NR4A2). Other genes annotated on chromosome 20 are involved in mammary gland metabolism (GHR, OXCT1), antibody production and phagocytosis of bacterial cells (C6, C7, C9, C1QTNF3), tumor suppression (DAB2), involution of mammary epithelium (OSMR) and cytokine regulation (PRLR). DAVID enrichment analysis revealed 5 KEGG pathways. The JAK-STAT signaling pathway (cell proliferation and apoptosis) and the 'Cytokine-cytokine receptor interaction' (cytokine and interleukines response to infectious agents) are co-regulated and linked to the 'ABC transporters' pathway also found here. Gene network analysis performed using GeneMania revealed a co-expression network where 665 interactions existed among 145 of the genes reported above. Clinical mastitis is a complex trait and the different genes regulating immune response are known to be pathogen-specific. Despite the lack of information in this study, candidate QTL for CM were identified in the US Holstein population.