Investigation of the Functional Proteins Related to Fertility in Cattle’s Endometrium by Protein-Protein Interactions Networks

Document Type : Research Articles

Authors

1 Department of Animal Science, College of Agriculture and Natural Resources, Gonbad Kavous University, Gonbad Kavous, Iran

2 Department of Animal Science, Faculty of Agriculture, University of Birjand, Birjand, Iran

3 Department of Animal Science, Faculty of Animal and Food Science, Ramin Agricultural and Natural Resources University, Mollasani, Ahvaz, Iran

Abstract

Pregnancy loss is an important economic loss in cattle industry. This study was conducted to identify pre- and / or post-implantation genes and cellular algorithms. For this purpose, transcriptome data of endometrium tissue were analyzed. These data refer to three heifer categories: high fertile (HF), sub-fertile (SF) and infertile (IF). After gene detection, genes were divided into two groups: Up-expressed genes, which were up-regulated in every comparison of either favorable fertility cases or unfavorable fertility cases (HF vs. SF, HF vs. IF, and SF vs. IF), and down-expressed genes, which were down-regulated in the mentioned comparisons. String database was applied to construct protein-protein interaction (PPI) networks and clusterone plugin was used to determine the significant sub-network. Enrichment analysis, which involves the gene ontology and functional pathway, was performed to enrich the results. Our results suggested that over-expression of SHCBP1, NOP14, PGM5, and DHX58 genes may have positive effect on the outcome of pregnancy, and down-expression of IMP3, ATP5O, and RPL7 genes could help the reproductive efficiency. The results of the present study showed that the genes in up-regulated clusters could manipulate epithelial differentiation, fundamental biological role, glucose metabolism, and immune response, which led to reduced pregnancy loss. Also the genes in down-regulated clusters may participate in the improvement of pregnancy outcome by inducing anti-apoptotic processes. This study proposes the pregnancy-associated key genes and pathways to improve pregnancy success in cattle and other domestic animals.

Keywords


INTRODUCTION

It is known that in dairy cattle industry the improvement of reproductive performance is essential for increasing milk production and economic efficiency. Reproductive diseases (such as retained placenta, metritis, and ovarian cysts), calving season, milk production, and genetics are the most important elements affecting fertility. Furthermore, the problems related to cattle health and management, such as impairment of immune function, poor estrous symptoms, increased non-ovulation period after calving, reducing pregnancy rate, and high rate of mortality, affect the reproductive performance (Walsh et al. 2011). In other words, fertility is affected by various factors, which all of those factors have slight, but not negligible, effect on this trait. Genetic improvement of this trait is more complicated due to its low heritability (Sheikhlou et al. 2018). Identifying the most important genetic regulators controlling the reproductive trait, and using them in genetic selection or selection index, is a difficult, complex and challenging task. The reproductive genetic markers have been investigated in several studies. Walker et al. (2009) studied the expression profiling of endometrium tissue in order to evaluate the effect of genes on pregnancy. They also provided 10 molecular markers controlling the pregnancy establishment. The potential genes which cause embryo mortality, which describe infertility and subfertility phenotypes, were also investigated (Geary et al. 2016). These researchers also reported that biomarker genes related to uterine competency for pregnancy are involved in immune signaling pathway. In another study, endometrial transcriptome (6 days after artificial insemination) was used to screen the potential gene markers for pregnancy maintenance (Binelli et al. 2015). Genome wide association analysis was implemented to provide quantitative trait loci (QTL) markers related to early pregnancy loss in cattle (Neupane et al. 2017). The gene expression profiling of peripheral white blood cells was applied to introduce the potential biomarkers of pregnancy status, and it is possible to use these biomarkers not only for physiological condition, but also for classifying the methods of pregnancy (artificial insemination or natural breeding) (Dickinson et al. 2018). Bioinformatics is a powerful tool in identification of novel biomarkers related to the complex traits from OMIC data. One of the most favorite tools in this topic is protein-protein interactions (PPIs) networks. Biomarkers proposed by this approach, have more biological importance, because these networks consider all the interactions between significant markers of a given trait (Dehmer et al. 2013). PPIs were applied to recognize the most important regulatory pathways and genes involved in activation of immune system in weaning stress condition (Behdani and Bakhtiarizadeh, 2017). Huang et al. (2018b) constructed gene regulatory network form disease-related gene sets, which was screened by meta-analysis. They also detected biomarkers involved in human disease by identifying PPIs. PPIs network was applied to discover novel interactions and elicit similar functional gene sets (Tian et al. 2017). Although studies have been conducted on the identification of molecular biomarkers related to pregnancy, but considering the importance of this issue on domestic animals, especially dairy cattle, further studies should be conducted to clarify signaling pathways controlling pregnancy in the endometrium and identify more gene markers regulating different aspects of this physiological condition. In this study, the most important signaling pathway and gene markers related to the endometrial ability to establish and maintain pregnancy in cattle were described by screening PPIs network.

 

MATERIALS AND METHODS

Data set

In this study, endometrium RNA-Seq data from the study conducted by Geary et al. (2016) were used to achieve fertility-related molecular markers and signaling pathways. These data referred to endometrium gene expression of three predominantly classified groups of Angus heifers: high fertile (HF), sub-fertile (SF), and infertile (IF). Heifers' ability to maintain pregnancy from embryo transfer was the basis of the fertility classification. Endometrium tissue was biopsied 14 days after the pregnancy. These data are available on NCBI database under the code of GSE81449 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi? acc=GSE81449).

 

Gene detection

The quality control analysis of all samples reads was performed using FastQC (a quality control tool for high throughput sequence data). Then, low quality read filtering (quality score below 20) and adaptor sequence removing were performed by Trimmomatic software (Bolger et al. 2014). In the next analysis step, the rest of reads were mapped by Tophat2 software (v2.0.9) to bovine genome (UMD3.1) from Ensemble database (Trapnell et al. 2009). HTseq (v0.6.1) was also used to detect the genes using the intersection-strict mode and the Ensemble annotation (version 3.1.78) of the bovine genome (Anders et al. 2015). In this step, in each given gene the number of reads which strictly matched to gene's sequence were determined. An expression matrix was constructed by merging HTseq output from each sample. Furthermore, the EdgeR Bioconductor package was used to differentially detect gene expression (Robinson et al. 2010). In order to increase the accuracy of differential gene expression analysis, low expression genes (rather 1 count per million in at least 3 sample) were removed from the expression matrix by this package (Labaj et al. 2011). In this study, it was attempted to classify genes into the following categories by fold change: up-expression genes and down-expression genes, in comparisons of HF vs. SF, HF vs. IF, and SF vs. IF.

 

PPIs network and significant sub-network

It was assumed that a given gene could be an important fertility biomarker depending on the up-expressed or down-expressed in comparisons of either favorable fertility cases or unfavorable fertility cases, such as HF vs. SF, HF vs. IF and SF vs. IF. String website (https://string-db.org/) was applied to survey the interactions of up-expressed and down-expressed genes (Szklarczyk et al. 2014). Clusterone cytoscape plugin was used to find significant sub-network (Maraziotis et al. 2008). This plugin was applied to the network clustering and modules detecting based on graph topology and overlapping neighborhoods (Sahrawat and Bhalla, 2013; Shao-Jun, 2015).

 

Enrichment analysis

In order to predict biological gene pathways of significant sub-networks, which manage the animal's ability to pregnancy, comparative GO (www.comparativego.com) was used (Fruzangohar et al. 2013). Comparative GO is a web tool GO resource for gene classification, based on the cellular component, biological process, and molecular function. EnrichNet (www.enrichnet.org) (Glaab et al. 2012), which is a web-based application, was applied to peruse the functional network and KEGG (Kyoto Encyclopedia of Genes and Genomes) was the pathway of the significant sub-network.

 

RESULTS AND DISCUSSION

To recognize the cellular algorithm for pregnancy maintenance, the endometrium transcriptome data of three cattle category (HF, SF and IF) were analyzed (Geary et al. 2016). In this study, co-expressed genes were also considered, based on which all the genes were divided into two groups of up-expressed genes and down-expressed genes. Up-expressed genes were those with expression of a fold change above zero in all comparisons of desirable uterine capacity vs. undesirable uterine capacity (HF vs. SF, HF vs. IF, and SF vs. IF). Moreover, those genes with expression fold change under zero in all the same comparisons, were called down-expressed genes. The results showed that 2517 genes were up-expressed (supplemental file S1: Co-expressed gene which up-regulated in every comparison of desirable uterine capacity vs. undesirable) and 2865 genes were down-expression (supplemental file S2: Co-expressed gene which down-regulated in every comparison of desirable uterine capacity vs. undesirable). The Search Tool for the Retrieval of Interacting Genesb (STRING 10.0, http://string.embl.de/) database was used to confirm the co-expressed relation of up-expressed and down-expressed genes and construct a co-expression network. Co-expression of up-regulated gene network consisted of 97 genes and 165 interactions (Figure 1), also the PPIs enrichment P-value was < 1.0e-16 (supplemental file S3: PPIs of co-expressed of up-regulated genes). Co-expression of 354 down-regulated genes was approved by STRING web-tool and based on co-expressed relation by 1648 interactions (Figure 2) and PPIs enrichment p-value was < 1.0e-16 (supplemental file S4: PPIs of co-expressed of down-regulated genes). ClosterOne, which is a cytoscape plugin, was used to dissect co-expressed networks and achieve significant sub-network. Up-expressed co-expression network was analyzed to 4 clusters (P<0.001) (Table 1, detailed information on supplemental file S5: Clustering of PPIs of up-regulated co-expression genes), and down-expressed co-expression network was dissected to 3 sub-networks (P<0.001) (Table 3, detailed information on supplemental file S6: Clustering of PPIs of down-regulated co-expression genes). To find functional roles of significant up- and down-expressed co-expression sub-networks, enrichment analysis was performed by comparative GO analysis for biological process related to fertility capacity and EnrichNet for pathway. The biological processes, which overexpress in up-expressed co-expression sub-networks, were cell proliferation, RNA processing, immune response, and gluconeogenesis, respectively (Table 5, detailed information on supplemental file S7: Biological processes of up-regulated co-expression gene clusters). The biological processes down-regulated by down-expressed co-expression sub-networks were ribosomal biogenesis, mitochondrial respiratory chain, and apoptotic process (Table 6, detailed information on supplemental file S8: Biological processes of down-regulated co-expression gene clusters). The analysis showed that the first up-expressed co-expression sub-networks activated cell cycle pathway and other up-expressed co-expression sub-networks did not affect any cellular pathways (Table 2). Second and third down-expressed co-expression sub-networks affected the oxidative phosphorylation and ribosome pathways (Table 4). Geary et al. (2016) used differential gene expression analysis to find the genes effective in pregnancy maintenance. To achieve these findings, they obtained 26 and 12 differentially expressed genes respectively in comparison of HF to SF and SF to IF. They selected two thresholds (>2fold change and false discovery rate (FDR) P<0.10). In various studies, different levels of FDR and fold change were applied to select differentially expressed genes. In each levels of thresholds, it is possible to loss some effective genes, such as transcription factors, which express in basic levels but have unignorable roles in cellular functions and programs. Therefore, in this study, co-expressed genes were investigated to reduce the ignorance of potential effective genes related to uterine capacity to stablish and maintain pregnancy. Co-expressed genes are the up- or down- expressed genes in every applied comparison (HF vs. SF, HF vs. IF and SF vs. IF). In other words, in this study, it is assumed that a given up- or down-expressed gene, not only in comparison of HF vs. SF but also in HF vs. IF and SF vs. IF, potentially affected the animal's ability to maintain pregnancy.

 

Figure 1 Pregnancy-related up-regulated PPI networks

The important regulatory genes are shown in large blue nodes

Nodes’ regulatory role reduce, as they turn to red color and small size

 

Figure 2 Pregnancy-related down-regulated PPI networks

The important regulatory genes are shown in large blue nodes

Nodes’ regulatory role reduce, as they turn to red color and small size

 

Based on this, two gene categories were formed: up co-expressed genes and down co-expressed genes. Gene interactions of each category, searched by String database, constructed the up and down PPIs networks. ClusterOne, which is a cytoscape plugin, segregated up and down PPIs networks to achieve significant up and down sub-networks. Network analysis by cytoscape showed that SHCBP1, NOP14, PGM5 and DHX58 were the most important regulator genes in up co-expressed clusters (supplemental file S9: Network analysis of up-regulated co-expression gene clusters). SHCBP1, which is the major regulator gene in first cluster, is a conserved gene in eukaryotes and plays critical roles in T-cell proliferation (Peng et al. 2017) and carcinogenesis (Buckley et al. 2014). Various studies identified the important roles of SHCBP1 and its orthologue genes in spermatogenesis in cytokinesis step of vertebrates (Montembault et al. 2010; Liu et al. 2014). In a recent report by Huang et al. (2018b) RNA-Seq data illustrated that this gene expression was induced in endometrial tissue in hamster. SHCBP1 gene was responsive to estrogen and progesterone by microarray analysis (Ren et al. 2015). The SHCBP1 up-regulated greater than five folds in endometrium, detected as functional molecule in maternal-fetal communications (Pavličev et al. 2017). Biological processes and functional role of genes, involved in the first sub-network of up-regulated PPIs network, contributed to cell proliferation (Table 5). It is known that endometrial proliferation and vascular remodeling not only is necessary during pregnancy (Robb et al. 2007), but also is critical to embryo implantation (Tranguch et al. 2005; Egashira and Hirota, 2013). Our findings suggested that this sub-network had the major candidate genes affecting pregnancy maintenance through proliferation and differentiation of endometrial tissue and proper embryo–uterine communications. Network analysis showed that in second significant up-regulated co-expressed sub-network, NOP14 had a critical biological role (supplemental file S9). NOP14 is a conserved protein and it is essential for proper maturation of 18S rRNA and ribosome biogenesis (Liu and Thiele, 2001). Exome sequencing using DNA indicated that homozygous variant of this gene led to pregnancy loss and fetal mortality (Suzuki et al. 2018). It was previously reported that this gene contributes in different tissue development in zebrafish embryo (Burns et al. 2009). Fundamentally biological role of this gene in cell leads to embryonic lethal in mutation cases (Ching et al. 2010). It was illustrated that NOP14 was down-regulated in hypoxia condition and it is related to pregnancy disorder (Chakraborty, 2013). The results of this study suggest that genes of the second significant up-regulated sub-network, are involved in ribosomal processing and affect pregnancy by fundamental cellular roles. Based on this study, PGM5 was a considerable gene in the third cluster of up-expressed genes (Figure 3). Gene ontology also showed that carbohydrate metabolic pathways are significantly overexpressed (Table 5).

 

Table 1 Significant up-expressed co-expression sub-networks

 

Table 2 Functional pathways of up-expressed co-expression sub-networks

 

Table 3 Significant down-expressed co-expression sub-networks

 

Table 4 Functional pathways of down-expressed co-expression sub-networks

 

The single blastocyst transcriptome analysis, which was performed by RNA-Seq, showed the participation of PGM5 in glucose metabolism and cell adhesion biological pathways (Chitwood et al. 2013). In another study, it was stated that homeostasis and transport of lipid and glucose were up-regulated in endometrial tissue at day 17 of pregnancy (Cerri et al. 2012). These results indicated that the genes of the third up-expressed cluster are the major markers in pregnancy maintenance, along with important biological role in glucose homeostasis. In the fourth significant sub-network of up-regulated co-expressed genes, DHX58 played the most important role in comparison with other genes (Figure 3). DHX58 can identify pathogenic RNA and contribute in response to cytosolic double-stranded DNA (dsDNA) (Pollpeter et al. 2011; Zhu et al. 2014). Previous studies have indicated that DHX58 induces antiviral immune responses by producing interferons (Satoh et al. 2010; Zhu et al. 2014). Transcriptomic analysis of oviduct epithelium demonstrated that genes related to immune terms (such as CXCL6, DHX58, HERC5, ISG15, MX1, MX2, NKL, OAS2, and OASL) up-regulated and affected the equine embryo development (Smits et al. 2016). According to these data, it was declared that this cluster contributes in the immune responses in endothelial tissue and helps to increase uterus capacity (Table 5). In many studies, it was stated that endometrial lymphocytes and cytokines not only have significant effect on successful pregnancy (Red-Horse et al. 2004; Chakraborty, 2013), but also control the endometrial differentiation and implantation (Takabatake et al. 1997).

 

Table 5 Biological processes of up-expressed co-expression sub-networks

 

Table 6 Biological processes of down-expressed co-expression sub-networks

 

Studies also revealed that cytokines and steroid hormones construct integrated pathways to control gene expression involved in the preparation of uterine for implantation (Li et al. 2001; Spencer et al. 2008). Down co-expressed clusters IMP3, ATP5O and RPL7 were the most biologically important genes, controlling the gene expression of others (supplemental file S10: Network analysis of down-regulated co-expression gene clusters). IMP3, which is the most important regulatory gene in first down-regulated sub-network (Figure 4), acts as RBPs (RNA-binding proteins) and it has biological and pathological roles in mRNA fate by regulating miRNA activity through imposing a stabilizing effect on the structure of mRNA targets (Degrauwe et al. 2016). This gene is associated with cell migration and tumor invasion and it was reported in recurrent pregnancy loss (Luizon et al. 2017), oncofetal biomarker (Zheng et al. 2008), and endometrial cancer (Djordjevic et al. 2012). Furthermore, studies have indicated that IMP3 is a regulatory element for IGFII. Also, IMP3 has a negative effect on IGFII signaling and it inhibits the differentiation and proliferation induced by this signaling pathway in central nervous system (Mori et al. 2001). These results suggest that down expression of IMP3 is necessary for pregnancy maintenance and its success. Also the first cluster of down-regulated network is involved in RNA processing and ribosome biogenesis. It seems that considering the high regulatory effect of IMP3 in this cluster, these genes could cooperate in mRNA fate by interactions between miRNA and RBPs. Also, there are studies demonstrating the functional link between miRNAs and other important regulatory genes, such as WDR12 (Shen and Liu, 2017), POLR1B (Hu et al. 2012), and UTP3 (Gutierrez et al. 2010). In second down-expressed cluster, ATP5O is the most considerable regulatory element (Figure 4). ATP5O is involved in mitochondrial energy production processes and it is one of the genes introduced as biomarker of mitochondrial electron chain activity (Brüggemann et al. 2017). Overexpression of this gene was reported in the blastocyst and its development (Gad et al. 2011). It has been reported that ATP5O is one of the 35 genes down-regulated in normal endometrium tissue (Saxena, 2015). This gene was also induced in immature uterine of rats through estradiol and estrogenic compounds (Hong et al. 2006). Also, the ATP5O involvement in the implantation processes and endometrium modulation in mouse was indicated using LongSAGE method (Ding et al. 2013). It was demonstrated that down-regulation of ATP5O and other mitochondrial electron transfer genes inhibits the apoptosis by reducing reactive oxygen species (ROS) and inducing oxidative stress pathway (He et al. 2006). It seems that down-expression of this gene help the survival and existence of endometrium cell.

 

Figure 3 Significant up-regulated sub-networks associated with pregnancy maintenance

The important regulatory genes are shown in large blue nodes

Nodes’ regulatory role reduce, as they turn to red color and small size

 

Figure 4 Significant down-regulated sub-networks associated with pregnancy maintenance

The important regulatory genes are shown in large blue nodes

Nodes’ regulatory role reduce, as they turn to red color and small size

 

According to gene ontology results (Table 6), ATP5O is related the cluster participating in respiratory-related processes, which induce production of intracellular ROS and cell degeneration. The down-regulation of this cluster potentially provides possible regulatory manners for endometrial cell maintenance by inhibiting cellular oxidative stress condition. Based on the statistical analysis of network, RPL7 is the major principal node in interactions of the third down-expressed clusters (Figure 4). This gene is a ribosomal protein and it is highly conserved across eukaryotes. RPL7 is a multifunctional protein involved in growth, transformation, and aging of cell (Hemmerich et al. 1993). It was also proposed that RPL7 is one of significant genes in apoptosis and could regulate metastasis processes (Zucchini et al. 2008) and expression of this gene was reported to be in response to cellular oxidative stress (Grewal et al. 2004). In endometrial cell, down-expression of RPL7 potentially mediates cell aging and inhibits the apoptosis process. These results indicated that the reaction of mitochondrial oxidative phosphorylation and mitochondrial translation processes are the most significant biological overexpressed processes of the third down-expressed clusters (Table 6), which induce precocious aging. Therefore, this cluster induces uterine capacity by inhibiting the mitochondrial activity. There are several studies indicating mitochondrial oxidative reaction, mitochondrial protein expression interaction, and protein synthesis machinery construct complex pathways, which determine cell fate (Wadhwa et al. 2015; Ansoleaga et al. 2016).

 

CONCLUSION

The results of this paper highlight genes and molecular mechanisms, participating in the capability of endometrium cell and providing the necessary conditions for pre- and/or post-implantation of fetus. Our analysis proposes the potential biomarkers of the optimal endometrial tissue conditions in implantation, which may have positive effect on cell proliferation, gluconeogenesis, immune response, and anti-apoptotic processes.

 

ACKNOWLEDGEMENT

This work was supported by Gonbad Kavous University (Gonbad Kavous, Iran), Grant number 6.716.

Anders S., Pyl P.T. and Huber W. (2015). HTSeq-a python framework to work with high-throughput sequencing data. Bioinformatics. 31(2), 166-169.
Ansoleaga B., Garcia-Esparcia P., Llorens F., Hernández-Ortega K., Carmona M., Antonio del Rio J., Zerr I. and Ferrer I. (2016). Altered mitochondria, protein synthesis machinery, and purine metabolism are molecular contributors to the pathogenesis of creutzfeldt–jakob disease. J. Neuropathpl. Exp. Neurol. 75(8), 755-769.
Behdani E. and Bakhtiarizadeh M.R. (2017). Construction of an integrated gene regulatory network link to stress-related immune system in cattle. Genetica. 145(4), 441-454.
Binelli M., Scolari S.C., Pugliesi G., Van Hoeck V., Gonella-Diaza A.M., Andrade S.C., Gasparin G.R. and Coutinho L.L. (2015). The transcriptome signature of the receptive bovine uterus determined at early gestation. PLoS One. 10(4), e0122874.
Bolger A.M., Lohse M. and Usadel B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics. 30(15), 2114-2120.
Brüggemann M., Gromes A., Poss M., Schmidt D., Klümper N., Tolkach Y., Dietrich D., Kristiansen G., Müller S.C. and Ellinger J. (2017). Systematic analysis of the expression of the mitochondrial ATP synthase (Complex V) subunits in clear cell renal cell carcinoma. Transl. Oncol. 10(4), 661-668.
Buckley M.W., Arandjelovic S., Trampont P.C., Kim T.S., Braciale T.J. and Ravichandran K.S. (2014). Unexpected phenotype of mice lacking Shcbp1, a protein induced during T cell proliferation. PLoS One. 9(8), e105576.
Burns C.E., Galloway J.L., Smith A.C., Keefe M.D., Cashman T.J., Paik E.J., Mayhall E.A., Amsterdam A.H. and Zon L.I. (2009). A genetic screen in zebrafish defines a hierarchical network of pathways required for hematopoietic stem cell emergence. Blood. 113(23), 5776-5782.
Cerri R., Thompson I., Kim I., Ealy A., Hansen P., Staples C.R., Li J.L., Santos P. and Thatcher W. (2012). Effects of lactation and pregnancy on gene expression of endometrium of Holstein cows at day 17 of the estrous cycle or pregnancy. J. Dairy Sci. 95(10), 5657-5675.
Chakraborty D. (2013). Natural killer cells, hypoxia, and epigenetic regulation of hemochorial placentation. Ph D. Thesis. University of Kansas, Kansas, USA.
Ching Y.H., Munroe R.J., Moran J.L., Barker A.K., Mauceli E., Fennell T., Lindblad-Toh K., Abcunas L.M., Gilmour J.F. and Harris T.P. (2010). High resolution mapping and positional cloning of ENU-induced mutations in the Rw region of mouse chromosome 5. BMC Genet. 11(1), 106-113.
Chitwood J.L., Rincon G., Kaiser G.G., Medrano J.F. and Ross P.J. (2013). RNA-seq analysis of single bovine blastocysts. BMC Genom. 14(1), 350-361.
Degrauwe N., Suvà M.L., Janiszewska M., Riggi N. and Stamenkovic I. (2016). IMPs: An RNA-binding protein family that provides a link between stem cell maintenance in normal development and cancer. Gene Dev. 30(22), 2459-2474.
Dehmer M., Mueller L.A. and Emmert-Streib, F. (2013). Quantitative network measures as biomarkers for classifying prostate cancer disease states: a systems approach to diagnostic biomarkers. PLoS One. 8(11), e77602.
Dickinson S.E., Griffin B.A., Elmore M.F., Kriese-Anderson L., Elmore J.B., Dyce P.W., Rodnig S.P. and Biase F.H. (2018). Transcriptome profiles in peripheral white blood cells at the time of artificial insemination discriminate beef heifers with different fertility potential. BMC Genom. 19(1), 129-135.
Ding Y.B., He J.L., Chen X.M., Liu X.Q. and Wang Y.X. (2013). Novel differential transcript expression identified by LongSAGE in the mouse endometrium during the implantation window. Mol. Biol. Rep. 40(1), 651-663.
Djordjevic B., Westin S. and Broaddus R.R. (2012). Application of immune-histochemistry and molecular diagnostics to clinically relevant problems in endometrial cancer. Surg Pathol Clin. 5(4), 859-878.
Egashira M. and Hirota Y. (2013). Uterine receptivity and embryo-uterine interactions in embryo implantation: Lessons from mice. Reprod. Med. Biol. 12(4), 127-132.
Fruzangohar M., Ebrahimie E., Ogunniyi A.D., Mahdi L.K., Paton J.C. and Adelson D.L. (2013). Comparative GO: A web application for comparative gene ontology and gene ontology-based gene selection in bacteria. PLoS One. 8(3), e58759.
Gad A., Besenfelder U., Rings F., Ghanem N., Salilew-Wondim D., Hossain M., Tesfaye D., Lonergan P., Becker A. and Cinar U. (2011). Effect of reproductive tract environment following controlled ovarian hyperstimulation treatment on embryo development and global transcriptome profile of blastocysts: implications for animal breeding and human assisted reproduction. Hum. Reprod. 26(7), 1693-1707.
Geary T.W., Burns G.W., Moraes J.G., Moss J.I., Denicol A.C., Dobbs K.B., Ortega M.S., Hansen P.J., Wehrman M.E. and Neibergs H. (2016). Identification of beef heifers with superior uterine capacity for pregnancy. Biol. Reprod. 95(2), 1-12.
Glaab E., Baudot A., Krasnogor N., Schneider R. and Valencia A. (2012). EnrichNet: network-based gene set enrichment analysis. Bioinformatics. 28(18), 451-457.
Grewal R., Stepczynski J., Kelln R., Erickson T., Darrow R., Barsalou L., Patterson M., Organisciak D.T. and Wong P. (2004). Coordinated changes in classes of ribosomal protein gene expression is associated with light-induced retinal degeneration. Invest. Ophthalmol. Vis. Sci. 45(11), 3885-3895.
Gutierrez A., Tschumper R.C., Wu X., Shanafelt T.D., Eckel-Passow J., Huddleston P.M., Slager S.L., Kay N.E. and Jelinek D.F. (2010). LEF-1 is a prosurvival factor in chronic lymphocytic leukemia and is expressed in the preleukemic state of monoclonal B cell lymphocytosis. Blood. 116(16), 2975-2983.
He Y., Reichow S., Ramamoorthy S., Ding X., Lathigra R., Craig J.C., Sobral B.W., Schurig G.G., Sriranganathan N. and Boyle S.M. (2006). Brucella melitensis triggers time-dependent modulation of apoptosis and down- egulation of mitochondrion-associated gene expression in mouse macrophages. Infect Immun. 74(9), 5035-5046.
Hemmerich P., Mikecz A.V., Neumann F., Sözeri O., Wolff-Vorbeck G., Zoebelein R. and Krawinkel U. (1993). Structural and functional properties of ribosomal protein L7 from humans and rodents. Nucleic Acids Res. 21(2), 223-231.
Hong E.J., Park S.H., Choi K.C., Leung P.C. and Jeung E.B. (2006). Identification of estrogen-regulated genes by microarray analysis of the uterus of immature rats exposed to endocrine disrupting chemicals. Reprod. Biol. Endocrinnol. 4(1), 49-59.
Hu Y., Matkovich S.J., Hecker P.A., Zhang Y., Edwards J.R. and Dorn G.W. (2012). Epitranscriptional orchestration of genetic reprogramming is an emergent property of stress-regulated cardiac microRNAs. Proc. Natl. Acad. Sci. 109(48), 19864-19869.
Huang J.K., Carlin D.E., Yu M.K., Zhang W., Kreisberg J.F., Tamayo P. and Ideker T. (2018a). Systematic evaluation of molecular networks for discovery of disease genes. Cell Syst. 6(4), 484-495. e485.
Huang M.Y., Zhang W.Q., Zhao M., Zhu C., He J.P. and Liu J.L. (2018b). Assessment of embryo-induced transcriptomic changes in hamster uterus using RNA-Seq. Cell Physiol. Biochem. 46(5), 1868-1878.
Labaj P.P., Leparc G.G., Linggi B.E., Markillie L.M., Wiley H.S. and Kreil D.P. (2011). Characterization and improvement of RNA-Seq precision in quantitative transcript expression profiling. Bioinformatics. 27(13), 383-391.
Li Q., Zhang M., Kumar S., Zhu L.J., Chen D., Bagchi M.K. and Bagchi I.C. (2001). Identification and implantation stage-specific expression of an interferon-α-regulated gene in human and rat endometrium. Endocrinology. 142(6), 2390-2400.
Liu M., Shi X., Bi Y., Qi L., Guo X., Wang L., Zhou Z. and Sha J. (2014). SHCBP1L, a conserved protein in mammals, is predominantly expressed in male germ cells and maintains spindle stability during meiosis in testis. Mol. Hum. Reprod. 20(6), 463-475. 
Liu P.C. and Thiele D.J. (2001). Novel stress-responsive genes EMG1 and NOP14 encode conserved, interacting proteins required for 40S ribosome biogenesis. Mol. Biol. Cell. 12(11), 3644-3657.
Luizon M.R., Palei A.C., Cavalli R.C. and Sandrim V.C. (2017). Pharmacogenetics in the treatment of pre-eclampsia: current findings, challenges and perspectives. Pharmacogenomics. 18(6), 571-583.
Maraziotis I.A., Dimitrakopoulou K. and Bezerianos A. (2008). An in silico method for detecting overlapping functional modules from composite biological networks. BMC Syst. Biol. 2(1), 93-103.
Montembault E., Zhang W., Przewloka M.R., Archambault V., Sevin E.W., Laue E.D., Glover D.M. and D'Avino P.P. (2010). Nessun Dorma, a novel centralspindlin partner, is required for cytokinesis in Drosophila spermatocytes. J. Cell Biol. 191(7), 1351-1365.
Mori H., Sakakibara S.I., Imai T., Nakamura Y., Iijima T., Suzuki A., Yuasa Y., Takeda M. and Okano H. (2001). Expression of mouse IGF2 mRNA binding protein 3 and its implications for the developing central nervous system. J. Neurosci. Res. 64(2), 132-143.
Neupane M., Geary T.W., Kiser J.N., Burns G.W., Hansen P.J., Spencer T.E. and Neibergs H.L. (2017). Loci and pathways associated with uterine capacity for pregnancy and fertility in beef cattle. PLoS One. 12(12), e0188997.
Pavličev M., Wagner G.P., Chavan A.R., Owens K., Maziarz J., Dunn-Fletcher C., Kallapur S.G., Muglia L. and Jones H. (2017). Single-cell transcriptomics of the human placenta: inferring the cell communication network of the maternal-fetal interface. Genome Res. 27(3), 349-361.
Peng C., Zhao H., Song Y., Chen W., Wang X., Liu X., Zhang C., Zhao J., Li J. and Cheng G. (2017). SHCBP1 promotes synovial sarcoma cell metastasis via targeting TGF-β1/Smad signaling pathway and is associated with poor prognosis. J. Exp. Clin. Cancer Res. 36(1), 141-149.
Pollpeter D., Komuro A., Barber G.N. and Horvath C.M. (2011). Impaired cellular responses to cytosolic DNA or infection with Listeria monocytogenes and vaccinia virus in the absence of the murine LGP2 protein. PLoS One. 6(4), e18842.
Red-Horse K., Drake P.M. and Fisher S.J. (2004). Human pregnancy: the role of chemokine networks at the fetal–maternal interface. Exp. Rev. Mol. Med. 6(11), 1-14.
Ren C.E., Zhu X., Li J., Lyle C., Dowdy S., Podratz K.C., Byck D., Chen H.B. and Jiang S.W. (2015). Microarray analysis on gene regulation by estrogen, progesterone and tamoxifen in human endometrial stromal cells. Int. J. Mol. Sci. 16(3), 5864-5885.
Robb A.O., Mills N.L., Newby D.E. and Denison F.C. (2007). Endothelial progenitor cells in pregnancy. Reproduction. 133(1), 1-9.
Robinson M.D., McCarthy D.J. and Smyth G.K. (2010). edgeR: A bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 26(1), 139-140.
Sahrawat T.R. and Bhalla S. (2013). Identification of critical target protein for cystic fibrosis using systems biology network approach. Int. J. Bioautomation. 17, 227-240.
Satoh T., Kato, H., Kumagai Y., Yoneyama M., Sato S., Matsushita K., Tsujimura T., Fujita T., Akira S. and Takeuchi O. (2010). LGP2 is a positive regulator of RIG-I–and MDA5-mediated antiviral responses. Proc. Natl. Acad. Sci. 107(4), 1512-151.
Saxena M. (2015). Pain related genes in endometriosis: A meta-analysis. MS Thesis. University of Sydney, Australia.
Shao-Jun C. (2015). Drug-target networks for Tanshinone IIA identified by data mining. Chinese J. Nat. Med. 13(10), 751-759.
Sheikhlou M.R., Badere S., Bahri Binabaj F. and Safari R. (2018). Inbreeding Trend and Inbreeding Depression on Birth and Weaning Weights of Iranian Baluchi Sheep. Pp. 78-81 in Proc. 18th Asian Australasian Anim. Prod. Congr. Kuching, Malaysia.
Shen C. and Liu Z.P. (2017). Identifying module biomarkers of hepatocellular carcinoma from gene expression data. 5404-5407 in Proc. Chinese Automation Congr, Jinan, China.
Smits K., De Coninck D.I., Van Nieuwerburgh F., Govaere J., Van Poucke M., Peelman L., Deforce D. and Van Soom A. (2016). The equine embryo influences immune-related gene expression in the oviduct. Biol. Reprod. 94(2), 31-38.
Spencer T.E., Sandra O. and Wolf E. (2008). Genes involved in conceptus–endometrial interactions in ruminants: Insights from reductionism and thoughts on holistic approaches. Reproduction. 135(2), 165-179.
Suzuki T., Behnam M., Ronasian F., Salehi M., Shiina M., Koshimizu E., Fujita A., Sekiguchi F., Miyatake S. and Mizuguchi T. (2018). A homozygous NOP14 variant is likely to cause recurrent pregnancy loss. J. Hum. Genet. 63(4), 425-431.
Szklarczyk D., Franceschini A., Wyder S., Forslund K., Heller D., Huerta-Cepas J., Simonovic M., Roth A., Santos A. and Tsafou, K.P. (2014). STRING v10: Protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 43, 447-452.
Takabatake K., Fujiwara H., Goto Y., Nakayama T., Higuchi T., Fujita J., Maeda M. and Mori T. (1997). Splenocytes in early pregnancy promote embryo implantation by regulating endometrial differentiation in mice. Hum. Reprod. 12(10), 2102-2107.
Tian Z., Guo M., Wang C., Liu X. and Wang S. (2017). Refine gene functional similarity network based on interaction networks. BMC Bioinformatics. 18(16), 550-559.
Tranguch S., Cheung-Flynn J., Daikoku T., Prapapanich V., Cox M.B., Xie H., Wang H., Das S.K., Smith D.F., Dey S.K. and Roberts S.K. (2005). Cochaperone immunophilin FKBP52 is critical to uterine receptivity for embryo implantation. Proc. Natl. Acad. Sci. 102(40), 14326-14331.
Trapnell C., Pachter L. and Salzberg S.L. (2009). TopHat: Discovering splice junctions with RNA-Seq. Bioinformatics. 25(9), 1105-1111.
Wadhwa R., Ryu J., Ahn H.M., Saxena N., Chaudhary A., Yun C.O. and Kaul S.C. (2015). Functional significance of point mutations in stress chaperone mortalin and their relevance to Parkinson disease. J. Biol. Chem. 290(13), 8447-8456.
Walker C.G., Meier S., Mitchell M.D., Roche J.R. and Littlejohn M. (2009). Evaluation of real-time PCR endogenous control genes for analysis of gene expression in bovine endometrium. BMC Mol. Biol. 10(1), 100-107.
Walsh S., Williams E. and Evans A. (2011). A review of the causes of poor fertility in high milk producing dairy cows. Anim. Reprod. Sci. 123(3), 127-138.
Zheng W., Yi X., Fadare O., Liang S.X., Martel M., Schwartz P.E. and Jiang Z. (2008). The oncofetal protein IMP3: A novel biomarker for endometrial serous carcinoma. Am. J. Surg. Pathol. 32(2), 304-315.
Zhu Z., Zhang X., Wang G. and Zheng H. (2014). The laboratory of genetics and physiology 2: Emerging insights into the controversial functions of this RIG-I-like receptor. BioMed Res. Int. 2014, 960190-960197.
Zucchini C., Rocchi A., Manara M.C., De Sanctis P., Capanni C., Bianchini M., Carinci P., Scotlandi K. and Valvassori L. (2008). Apoptotic genes as potential markers of metastatic phenotype in human osteosarcoma cell lines. Int. J. Oncol. 32(1), 17-31.