Analysis of protein-protein interaction network based on transcriptome profiling of ovine granulosa cells identifies candidate genes in cyclic recruitment of ovarian follicles.

After pubertal, cohort of small antral follicles enters to gonadotrophin-sensitive development, called recruited follicles. This study was aimed to identify candidate genes in follicular cyclic recruitment via analysis of protein-protein interaction (PPI) network. Differentially expressed genes (DEGs) in ovine granulosa cells of small antral follicles between follicular and luteal phases were accumulated among gene/protein symbols of the Ensembl annotation. Following directed graphs, PTPN6 and FYN have the highest indegree and outdegree, respectively. Since, these hubs being up-regulated in ovine granulosa cells of small antral follicles during the follicular phase, it represents an accumulation of blood immune cells in follicular phase in comparison with luteal phase. By contrast, the up-regulated hubs in the luteal phase including CDK1, INSRR and TOP2A which stimulated DNA replication and proliferation of granulosa cells, they known as candidate genes of the cyclic recruitment.

Transcriptome analysis of ovine granulosa cells reveals differences between small antral follicles collected during the follicular and luteal phases.

Ovarian follicular growth occurs in both the follicular and luteal phases of the estrous cycle but in very different endocrine contexts. In both phases, many small antral follicles with similar morphologic and histologic characteristics are present within the ovaries as a reserve for the terminal folliculogenesis. However, there are several gaps in our molecular knowledge of the gene expression profiles of small antral follicles in the follicular and luteal phases. The aim of the present study was to use RNA sequencing to compare and analyze the global transcriptional profile of ovine granulosa cells collected from small antral follicles (1-3 mm) either during the follicular or the luteal phase of the estrous cycle, with the hypothesis that they should be differential. We identified 663 genes whose mRNA was differentially expressed or accumulated in the granulosa cell layer of small antral follicles in the two phases. A comprehensive interpretation of these data was performed through integrative analyses (Gene Ontology, Ingenuity Pathway Analysis) and the exploitation of already available transcriptomic data on follicular growth and atresia. In particular, we observed that the contrasted endocrine context between follicular and luteal phases may have an impact on estradiol, follicle-stimulating hormone (FSH), and on the activin/inhibin signaling pathways. Furthermore, we reveal the possible initiation of early follicular atresia in small antral follicles during the follicular phase in interaction with the presence of immune cells. This study provides new insights into the gene expression profile in ovine granulosa cells, and we suggest that these molecular changes may have an implication at the time of follicle selection.

Genetic analysis of growth rate and Kleiber ratio in Zandi sheep.

Genetic parameters for average daily gain from birth to weaning (ADGa), birth to 6 months (ADGb), weaning to 6 months (ADGc), weaning to yearling age (ADGd), and corresponding Kleiber ratios (KRa, KRb, KRc, and KRd) were estimated by using records of 3,533 Zandi lambs, descendent of 163 sires and 1265 dams, born between 1991 and 2005 at the Zandi Sheep Breeding Station at Khojir National Park, Tehran, Iran. A derivative-free algorithm combined with a series of six single-trait linear animal models was used to estimate phenotypic variance and its direct, maternal, and residual components. In addition, bivariate analyses were done to estimate (co)variance components between traits. Estimates of direct heritability (h(2)) were 0.11, 0.15, 0.09, 0.10, 0.10, 0.10, 0.06, and 0.07 for ADGa, ADGb, ADGc, ADGd, KRa, KRb, KRc, and KRd, respectively, thereby indicating the presence of low additive genetic variation for growth rate and Kleiber ratio in this population of Zandi sheep. Maternal genetic component was found to be significant on ADGa and KRa and contributed 3% and 5%, respectively, in total phenotypic variance of ADGa and KRa. A widespread range of genetic correlations among traits studied was observed. Except for negative genetic correlations between ADGa and KRc, ADGa and KRd, and between KRa and KRc, in other cases, genetic correlations were positive and moderate to very high. Phenotypic correlations ranged from -0.49 (ADGa/KRd) to 0.94 (ADGc/KRc). These results indicate that selecting for improved growth rate or Kleiber ratio in Zandi sheep would generate a relatively slow genetic progress.