A multi-omics analysis identifies molecular features associated with fertility in heifers (Bos taurus).

Infertility or subfertility is a critical barrier to sustainable cattle production, including in heifers. The development of heifers that do not produce a calf within an optimum window of time is a critical factor for the profitability and sustainability of the cattle industry. In parallel, heifers are an excellent biomedical model for understanding the underlying etiology of infertility because well-nourished heifers can still be infertile, mostly because of inherent physiological and genetic causes. Using a high-density single nucleotide polymorphism (SNP) chip, we collected genotypic data, which were analyzed using an association analysis in PLINK with Fisher's exact test. We also produced quantitative transcriptome data and proteome data. Transcriptome data were analyzed using the quasi-likelihood test followed by the Wald's test, and the likelihood test and proteome data were analyzed using a generalized mixed model and Student's t-test. We identified two SNPs significantly associated with heifer fertility (rs110918927, chr12: 85648422, P = 6.7 × 10-7; and rs109366560, chr11:37666527, P = 2.6 × 10-5). We identified two genes with differential transcript abundance (eFDR ≤ 0.002) between the two groups (Fertile and Sub-Fertile): Adipocyte Plasma Membrane Associated Protein (APMAP, 1.16 greater abundance in the Fertile group) and Dynein Axonemal Intermediate Chain 7 (DNAI7, 1.23 greater abundance in the Sub-Fertile group). Our analysis revealed that the protein Alpha-ketoglutarate-dependent dioxygenase FTO was more abundant in the plasma collected from Fertile heifers relative to their Sub-Fertile counterparts (FDR < 0.05). Lastly, an integrative analysis of the three datasets identified a series of molecular features (SNPs, gene transcripts, and proteins) that discriminated 21 out of 22 heifers correctly based on their fertility category. Our multi-omics analyses confirm the complex nature of female fertility. Very importantly, our results also highlight differences in the molecular profile of heifers associated with fertility that transcend the constraints of breed-specific genetic background.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing.

BACKGROUND: A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation. RESULTS: In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework. CONCLUSIONS: Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

Cleavage kinetics is a better indicator of embryonic developmental competency than brilliant cresyl blue staining of oocytes.

In vitro production of embryos (IVP) is a valuable technology to produce embryos of high genetic value. Despite advances in IVP, the efficiency of culture systems remains low. One method to increase IVP success is the early selection of oocytes or embryos that may have greater developmental potential. Here, we investigated two methods of selection, namely BCB staining and cleavage kinetics, both individually and in conjunction, for improved developmental outcomes in vitro. We hypothesized that a synergistic use of both BCB staining and cleavage kinetics would result in identification of embryos of greater developmental potential. The selection of oocytes by BCB staining does select for those oocytes with higher developmental potential, as noted by a greater blastocyst development between BCB positive (32.6%) and BCB negative (22.0%) on day 8 post-fertilization. However, the utilization of BCB staining and cleavage kinetics in tandem resulted in a complete masking of the effect observed when using BCB alone. We obtained the highest proportion of blastocyst development per selection group using cleavage kinetics alone, in which 53.1% of embryos grouped as Fast produced a blastocyst, which was significantly different from the three other groups (Fast+, Slow, not cleaved). We observed, however, that the separation of embryos by cleavage kinetics did not predict their survival to cryopreservation. In conclusion, in standard culture systems, cleavage kinetics is an effective method for the selection of embryos with increased developmental potential to develop blastocysts, however, it may not be effective to select healthy embryos for transfer following cryopreservation.

Tight gene co-expression in BCB positive cattle oocytes and their surrounding cumulus cells.

BACKGROUND: Cytoplasmic and nuclear maturation of oocytes, as well as interaction with the surrounding cumulus cells, are important features relevant to the acquisition of developmental competence. METHODS: Here, we utilized Brilliant cresyl blue (BCB) to distinguish cattle oocytes with low activity of the enzyme Glucose-6-Phosphate Dehydrogenase, and thus separated fully grown (BCB positive) oocytes from those in the growing phase (BCB negative). We then analyzed the developmental potential of these oocytes, mitochondrial DNA (mtDNA) copy number in single oocytes, and investigated the transcriptome of single oocytes and their surrounding cumulus cells of BCB positive versus BCB negative oocytes. RESULTS: The BCB positive oocytes were twice as likely to produce a blastocyst in vitro compared to BCB- oocytes (P < 0.01). We determined that BCB negative oocytes have 1.3-fold more mtDNA copies than BCB positive oocytes (P = 0.004). There was no differential transcript abundance of genes expressed in oocytes, however, 172 genes were identified in cumulus cells with differential transcript abundance (FDR < 0.05) based on the BCB staining of their oocyte. Co-expression analysis between oocytes and their surrounding cumulus cells revealed a subset of genes whose co-expression in BCB positive oocytes (n = 75) and their surrounding cumulus cells (n = 108) compose a unique profile of the cumulus-oocyte complex. CONCLUSIONS: If oocytes transition from BCB negative to BCB positive, there is a greater likelihood of producing a blastocyst, and a reduction of mtDNA copies, but there is no systematic variation of transcript abundance. Cumulus cells present changes in transcript abundance, which reflects in a dynamic co-expression between the oocyte and cumulus cells.

Comparison of reproductive performance of AI- and natural service-sired beef females under commercial management.

The objective of this study was to assess differences in reproductive performance of natural service and artificial insemination (AI) sired beef females based on pregnancy outcomes, age at first calving, and calving interval. Data were sourced from 8,938 cows sired by AI bulls and 3,320 cows sired by natural service bulls between 2010 and 2017. All cows were in a commercial Angus herd with 17 management units located throughout Virginia and represented spring and fall calving seasons. All calves were born to dams managed with estrus synchronization. Pregnancy was analyzed with generalized linear mixed models and other reproductive measures with linear mixed models in R. Six models were evaluated with the dependent variables of pregnancy status at the first diagnosis, pregnancy status at the second diagnosis, pregnancy type (AI or natural service) at the first diagnosis, pregnancy type at the second diagnosis, calving interval, and age at first calving. Independent variables differed by model but included sire type of the female (AI or natural service), prebreeding measures of age, weight, and body condition score, postpartum interval, sex of the calf nursing the cow, and management group. No differences were observed between AI- and natural service-sired females based on pregnancy status at first and second pregnancy diagnosis (P > 0.05). Sire type was only found to be significant for age at first calving (P < 0.05) with AI-sired females being 26.6 ± 1.6 d older at their first calving, which was expected because AI-sired females were born early in the calving season making them older at breeding. Surprisingly, age and body condition score were not significant predictors of pregnancy (P > 0.05). Body weight at breeding was not significant for pregnancy (P > 0.05) but was significant for age at first calving (P < 0.05). These data suggested that lighter heifers calved earlier which contradicts our original hypothesis. Overall, commercial Angus females sired by AI or natural service bulls had similar reproductive performance. Factors that were commonly associated with reproductive success were not significant in this commercial Angus herd managed with estrus synchronization. Given the size of these data, the importance of body condition, age, and weight should be reassessed in modern genetics and management practices.