Endometrial transcriptional profiling of a bovine fertility model by Next-Generation Sequencing

Studying the multitude of molecular networks and pathways that are potentially involved in a complex trait such as fertility requires an equally complex and broad strategy. Here, we used Next-Generation Sequencing for the characterization of the transcriptional signature of the bovine endometrial tissue. Periovulatory endocrine environments were manipulated to generate two distinctly different fertility phenotypes. Cycling, non-lactating, multiparous Nelore cows were manipulated to ovulate larger (> 13 mm; LF group; high fertility phenotype) or smaller (< 12 mm; SF group) follicles. As a result, greater proestrus estrogen concentrations, corpora lutea and early diestrus progesterone concentrations were also observed in LF group in comparison to SF group. Endometrial cell proliferation was estimated by the protein marker MKI67 on tissues collected 4 (D4) and 7 (D7) days after induction of ovulation. Total RNA extracts from D7 were sequenced and compared according to the transcriptional profile of each experimental group (LF versus SF). Functional enrichment analysis revealed that LF and SF endometria were asynchronous in regards to their phenotype manifestation. Major findings indicated an LF endometrium that was switching phenotypes earlier than the SF one. More specifically, a proliferating SF endometrium was observed on D7, whereas the LF tissue, which expressed a proliferative phenotype earlier at D4, seemed to have already shifted towards a biosynthetically and metabolically active endometrium on D7. Data on MKI67 support the transcriptomic results. RNA-Seq-derived transcriptional profile of the endometrial tissue indicated a temporal effect of the periovulatory endocrine environment, suggesting that the moment of the endometrial exposure to the ovarian steroids, E2 and P4, regulates the timing of phenotype manifestation. Gene expression profiling revealed molecules that may be targeted to elucidate ovarian steroid-dependent mechanisms that regulate endometrial tissue receptivity. Data was deposited in the SRA database from NCBI (SRA Experiment SRP051330) and are associated with the Bio-Project (PRJNA270391). An overview of the gene expression data has been deposited in NCBI's Gene Expression Omnibus (GEO) and is accessible through GEO Series accession number GSE65450. Further assessment of the data in combination with other data sets exploring the transcriptional profile of the endometrial tissue during early diestrus may potentially identify novel molecular mechanisms and/or markers of the uterine receptivity.