Next Article in Journal
Deciphering Conformational Changes of the GDP-Bound NRAS Induced by Mutations G13D, Q61R, and C118S through Gaussian Accelerated Molecular Dynamic Simulations
Next Article in Special Issue
The Recent Development of Acoustic Sensors as Effective Chemical Detecting Tools for Biological Cells and Their Bioactivities
Previous Article in Journal
Molecular Mechanism of Tanshinone against Prostate Cancer
Previous Article in Special Issue
The Effect of Combining Post-Harvest Calcium Nanoparticles with a Salicylic Acid Treatment on Cucumber Tissue Breakdown via Enzyme Activity during Shelf Life
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Genomic Study of Inhibin Molecule Production through Granulosa Cell Gene Expression in Inhibin-Deficient Mice

1
National Center for International Research on Animal Genetics, Breeding and Reproduction (NCIRAGBR), Huazhong Agricultural University, Ministry of Science and Technology of the People’s Republic of China, Wuhan 430070, China
2
Department of Animal Breeding and Genetics, Sindh Agriculture University, Tandojam 70060, Sindh, Pakistan
3
College of Veterinary Sciences, Faculty of Animal Husbandry and Veterinary Sciences, University of Agriculture, Peshawar 25120, Khyber Pakhtunkhwa, Pakistan
4
College of Biosystems Engineering and Food Science, Zhejiang University, Hangzhou 310058, China
5
Department of Nutrition & Food Science, National Research Centre, Dokki, Giza 12622, Egypt
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(17), 5595; https://doi.org/10.3390/molecules27175595
Submission received: 11 August 2022 / Revised: 23 August 2022 / Accepted: 25 August 2022 / Published: 30 August 2022

Abstract

:
Inhibin is a molecule that belongs to peptide hormones and is excreted through pituitary gonadotropins stimulation action on the granulosa cells of the ovaries. However, the differential regulation of inhibin and follicle-stimulating hormone (FSH) on granulosa cell tumor growth in mice inhibin-deficient females is not yet well understood. The objective of this study was to evaluate the role of inhibin and FSH on the granulosa cells of ovarian follicles at the premature antral stage. This study stimulated immature wild-type (WT) and Inhibin-α knockout (Inha−/−) female mice with human chorionic gonadotropin (hCG) and examined hCG-induced gene expression changes in granulosa cells. Also, screening of differentially expressed genes (DEGs) was performed in the two groups under study. In addition, related modules to external traits and key gene drivers were determined through Weighted Gene Co-Expression Network Analysis (WGCNA) algorithm. The results identified a number of 1074 and 931 DEGs and 343 overlapping DEGs (ODEGs) were shared in the two groups. Some 341 ODEGs had high relevance and consistent expression direction, with a significant correlation coefficient (r2 = 0.9145). Additionally, the gene co-expression network of selected 153 genes showed 122 nodes enriched to 21 GO biological processes (BP) and reproduction and 3 genes related to genomic pathways. By using principal component analysis (PCA), the 14 genes in the regulatory network were fixed and the cumulative proportion of fitted top three principal components was 94.64%. In conclusion, this study revealed the novelty of using ODEGs for investigating the inhibin and FSH hormone pathways that might open the way toward gene therapy for granulosa cell tumors. Also, these genes could be used as biomarkers for tracking the changes in inhibin and FSH hormone from the changes in the nutrition pattern.

Graphical Abstract

1. Introduction

Gonadotropins are glycoprotein hormones produced in the pituitary by gonadotropic cells the lightest copious ovarian cells and regulate ovarian follicle development [1,2]. These peptide hormones are also known for regulating ovarian and testicular function and are essential for normal growth, sexual development, and reproduction. Also, gonadotropins include some essential hormones like inhibin, FSH, and luteinizing hormone (LH) [2,3]. For instance, inhibin is released from the ovarian granulosa cell. It is a heterodimeric glycoprotein that consists of α-subunit linked with a β-hormone [4]. It is a member of the transforming growth factor-β superfamily that stimulate the release of FSH from the pituitary cells. Moreover, inhibin plays an important role in the FSH feedback secretion regulation during puberty in females [5]. In which, its secretion changes during female puberty have disclose correlation with adulthood chronic diseases like diabetes type 2 and heart disease which are generally considered among the common malnutrition diseases [6]. In addition, the inhibin hormone is essential as a diagnostic marker for ovarian cancer [7,8]. For instance, Rathore, et al. [9] mentioned that Inhibin deficiency in mice leads to the growth of gonadal sex cord-stromal tumors. In such cases, the mice mostly died at the age of 28 days due to testicular and ovarian tumors with cachexia-related signs [10,11,12]. A genetic method assured that gonadotropins play a key modifier role for gonadal sex cord-stromal tumor growth in inhibin-deficient mice. In that study, Nagaraja, et al. [13] reported that inhibin genetically interacts with multiple factors that influence testicular and ovarian growth and differentiation, including the pituitary gonadotropins [14]. Furthermore, they found a complex interplay among inhibins, gonadotropins, and ovarian cancer. Thus, the tracking of inhibin and FSH hormones could enhance the diagnostic performance of the health-related disease, not just ovarian female disease.
Furthermore, it is known that granulosa cells (GCs) are somatic cells of the sex rope that is mainly related to an embryonic female gamete identified to be an oocyte that exists in the animals’ ovary [15]. GCs show various phenotypes in the follicle, reliant on their position. Additionally, GCs gene expression analysis is very important to understand its functional mechanism that is related to animal growth. Hence gene association networks are very essential for expressing the relationship patterns between genes transversely microarray data, where the weighted gene co-expression network analysis (WGCNA) algorithm is a necessary tool to determine the relationship patterns among the genes microarray data [16], for which, the WGCNA algorithm is frequently used to understand the genes’ molecular processes and identification of interrelated genes and modules [17]. Therefore, it can show the co-expression structure and cluster the expression data into modules of conserved function that allow one to detect patterns of gene connectivity that can be aligned with behavioral and physiological phenotypes [18]. Meanwhile, PCA is a multivariate statistical procedure that uses an orthogonal transformation to convert a set of observations of possibly correlated variables into a set of values of linearly uncorrelated variables called principal components [19,20,21].
Therefore, this study aimed to investigate the variations in hCG-induced gene expression in WT and Inha−/− granulosa cells. It provided a novel reference for the pathway mechanisms. Initially, differentially expressed genes (DEGs) in WT/Inha−/−with and without human chorionic gonadotropin (hCG) stimulation samples were identified. Then the overlapping genes as the characteristic genes in hCG stimulation and Inha−/− female mice granulosa cells were selected and further investigated. For this process, the WGCNA algorithm was used to identify the interrelated gene patterns. Besides, PCA was used to examine the important regulated genes. The examined outcome recommended that Inhibin α knockout and hCG stimulation can down-regulate JUP expression, and up-regulate Psmc3ip expression. These genes could be used as a marker for the early prediction of granulosa cell carcinomas.

2. Results

2.1. Hierarchical Clustering and Comparison Analysis of Selected DEGs in Different Groups

We extracted DEG expressions from two groups and drew hierarchical clustering heatmaps, as shown in Figure 1. The experimental samples in each group were divided into two separate parts, indicating that DEGs had obviously different expression patterns in each group (gene expression matrix can be found in Table S1, Supplementary Materials).

2.2. GO and KEGG Pathway Enrichment Analysis for the ODEGs

Figure 2a shows the Venn diagram of DEGs in groups Inha−/− vs. WT and Inha−/− (hCG) vs. WT(hCG). In which, 343 overlapped genes were found with a total of 25 (10 BP, 10 CC, 5 MF) significant related GO annotations for 343 overlapped genes that listed in Table 1. Also, the hierarchical clustering heatmap of ODEGs showed a significant negative correlation of the down-regulated Inha DEGs compared to the WT group (Figure 2b,c). What’s more, among the 343 shared DEGs, 341 had a consistent expression direction and high relevance, with a significant correlation coefficient of 0.9145 (p < 0.000001). The shared 343 DEGs were differentially expressed in both Inha−/− vs. WT group and Inha−/− (hCG) vs. WT(hCG) group (Table S1). Meanwhile, the GO BPs were significantly related to the cell cycle process (GO: 0022402), which had the most enrichment significance (p-value < 0.0001). Also, 17 genes were involved, such as CDC6, KIFC1, MKI67, DSN1, NUF2. Besides, according to KEGG enrichment analysis, ODEGs were significantly enriched in eight pathways (RFC5, PRIM1, RPA2, RFC4, LIG1, POLD2, POLE, MCM2) participating in the most significant related pathway: DNA replication (mmu03030). Using the ggplot2 package in R, significant related GO and KEGG pathways annotations were displayed in Figure 3.

2.3. Physiological Phenotypes R Modules and Genes Identification Based on WGCNA

For clustering the ODEGs based on the physiological phenotypes, the expression data were processed by square root transformation and used to infer co-expression gene network modules with the WGCNA network construction and the module detection method. Firstly, the distances among all the samples were studied to eliminate discrete samples with no discrete samples to be removed (Figure 4a). In which, a high affinity between Inha and Inha(hCG) was found by the clustering tree. Then a proper power-law coefficient was selected using the soft-threshold method (Figure 4b). Through this model the selected soft-threshold (X-axis) was 18 when the scale-free topology model fit was signed at correlation coefficient (R2) = 0.8 (Y-axis). Then, a dynamic hierarchical tree cut algorithm was used to detect the co-expression modules, and a total of six related modules were found (Figure 4c). Moreover, the R2 between the physiological phenotypes and each module had a very high correlation coefficient (over 0.8) with the physiological phenotypes (Figure 4d and Table 2). Thus, 153 genes in the top three modules (blue, green, and brown) in total were selected as representative ODEGs for further analysis based on their significant physiological phenotypes.

2.4. Co-Expression Network Construction

The Co-expression Network had 153 DEGs in blue, green, and brown modules, as well as an expression correlation coefficient from the WGCNA algorithm (expression correlation matrix was shown in Table S2). Also, the selected gene pair was based on the expression correlation coefficient >0.8 as shown in Figure 5. The co-expression network included 122 nodes in total. These nodes had 35 down-regulated genes (14 blue, 17 brown, and 4 green) and 87 up-regulated genes (30 blue, 36 brown, and 21 green genes) with 410 edges (129 negative coefficient connections and 281 positive coefficient connections) (Table 2 and Table S3).
Moreover, GO and KEGG pathway enrichment analysis showed that a total of 21 significantly related GO BPs and 3 KEGG pathways were found for the DEGs in the co-expression network. DEGs in the co-expression network were significantly (p-value < 0.01) related to the cell cycle BPs and participated in ECM-receptor interaction (mmu04512), Focal adhesion (mmu04510), and DNA replication (mmu03030) pathways (Figure 6 and Table 3).

2.5. miRNA-DEGs-TF Regulatory Network Construction

MicroRNAs (miRNAs) assume a pivotal role in controlling inborn and versatile immunity in humans and animals [22,23]. A total number of eight miRNAs and seven TFs that have the potential to regulate ODEGs in the co-expression network are listed in Table 4 and Table 5, respectively. Integration regulatory relationships were identified and constructed as a miRNA-DEGs-TF regulatory network between miRNA and DEGs, TFs and DEGs (Figure 7). In that network, 29 nodes (8 miRNAs, 7 TFs, and 14 DEGs), 7 down (1 blue, 6 brown) and 7 up-regulated (3 blue, 2 brown, and 2 green genes) and 56 edges (20 miRNA-DEGs regulation and 36 TFs-DEGs) have been confirmed (Table S4).

2.6. PCA for Genes in Regulatory Network

To further refine the important genes, the PCA algorithm defined 14 genes in the regulatory network. The cumulative proportion of fitted top 3 principal components accounts for 94.64% of the total variance, which means that they can effectively describe the vast majority of input gene variables [20]. In general, the cumulative contribution rate of more than 80% is considered to have caught most of the input variable information. Figure 8 showed the differences in the top 3PCs between the two groups where significant differences (p < 0.01) were observed among the four different groups. Also, the gene contributions to the PCs were listed in Table 6, where there were 10 genes whose contribution rate (absolute value) was over 0.9.

3. Discussion

It is shown that Inhibin-α plays an important role in follicular development, oocyte development, cell differentiation, and finally reproduction. The Inhibin-α knockout/down female mouse can develop ovarian cancer and the LH and FSH may play a crucial role in GCs tumor development [24,25]. This study aimed to examine hCG-induced gene expression changes in different types of granulosa cells (WT and Inha−/− types). Also, it provided an important reference for the pathway mechanisms by showing that DEGs were different in WT/Inha−/− before and after hCG stimulation. The ODEGs were used as characteristic genes in hCG stimulation and Inhibin α knockout (Inha−/−) female mice granulosa cells. This observation is in agreement with FarmanUllah, Liang, Khan, Salim, Rehman, Khan, Talpur, Schreurs, Gouda, Khan and Shujun [22] who mentioned that ODEGs can effectively work as biomarkers for immune-related tumors. In which, 341 DEGs had high relevance, with a significant correlation coefficient (p < 0.000001) in both Inha−/− vs. WT group and Inha−/− (hCG) vs. WT(hCG) group which means that they could be used as characteristic genes in hCG stimulation and Inhibin α knockout (Inha−/−) female mice granulosa cells. Therefore, this study demonstrated for the first time that hCG induces the granulosa cells to excrete Inha through stimulating Fndc5, Sertad4, JUP, and Psmc3ip genes. In which, quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) on selected gene expression changes were observed in the gene array analysis verified the most important ODEGs of the knockout mice. Vasilache, et al. [26] mentioned that qRT-PCR microarray combined with modeling is an effective technique to detect the knockout mice’s important DEGs.
According to the analysis of KEGG pathway annotations, ODEGs were significantly enriched in eight biological pathways: DNA replication, Focal adhesion, and purine metabolism pathways which significantly enriched the GO term. In the biological process category, the genes were mainly enriched in GO terms associated with extracellular matrix and axon [27], for which, the R2 between the physiological phenotypes and each KEGG module for DEGs had a very high correlation coefficient (>0.8). Chen, et al. [28] mentioned that a positive regulator of the steroidogenesis pathway of FSH is essential for the granulosa cell proliferation, death, and differentiation in almost all cell types. The WGCNA algorithm was used to detect related modules and genes, significantly related to eight miRNAs; whereas, the seven TFs and regulatory networks were utilized to get regulated DEGs [29]. Finally, PCA differentiated the four groups under study to determine the importance of 14 regulated genes. Also, the top 3PCs between the two groups were significantly different (p < 0.01) among the four groups under study. Among them, the JUP gene was significantly related to cell adhesion (GO: 0007155), DNA metabolic process (GO:0006259), and biological adhesion (GO: 0022610), while Psmc3ip participated in the cell cycle (GO: 0022402) and M phase (GO: 0000279), and they all belonged to the brown module of the physiological phenotypes of WGCNA results, so they had close expression relationship. Chen, et al. [30] reported that the JUP form is a member of the catenin family that can affect various processes such as proliferation, migration, and differentiation by mediating cellular adhesion. Thus, the mutation in its gene is associated with several gene-related diseases. [30]. In addition, Psmc3ip (also known as GT198) is used as a unique tumor marker suppressor gene for the mutant cells in ovarian cancer. Psmc3ip protein has been shown as a steroid hormone receptor regulator and also as a crucial factor in DNA repair [31]. Thus, studying such genes could facilitate the complex mission of dealing with ovarian cancer.
Additionally, the result suggested that Inhibin α knockout and hCG stimulation can down-regulate the expression of JUP and up-regulate Psmc3ip. In which, the co-expression Network had 153 DEGs expression correlation coefficient from WGCNA algorithm with R2 > 0.8 with 35 down-regulated genes and 87 up-regulated genes. Moreover, JUP forms distinct complexes with cadherins and desmosomal cadherins through an amino acid motif called the armadillo repeat, which can affect the diverse processes and modulate the function of extracellular ligands [32,33]. It also showed that JUP and Psmc3ip genes had close relationships in both expression patterns and functions in the Inha−/− hCG stimulation female mice granulosa cells. Similarly, the activity of the Psmc3ip gene is revealed to have a crucial role in ovarian dysgenesis and male fertility in mammalians [34,35]. Thus, the data analysis detected DEGs and relevant biological functions after the knockdown of the Inha and associated gene expression for further research guidance in mammalian reproduction.

4. Materials and Methods

4.1. Experimental Animals

A number of 100 specific-pathogen-free (SPF) mice were grouped (25 mice WT, 25 mice Inha−/−, 25 mice WT (hCG), and 25 mice Inha−/−(hCG)) according to Hofland, et al. [36]. In which, 21 to 23-day-old WT and Inha−/− female mice were injected with 5 IU hCG for 6 h to stimulate hCG groups and granulosa cells with and without hCG stimulation which were collected from 2 genotypes (WT and Inha−/−) according to National Institutes of Health (NIH) Guidelines for the Care and use of Laboratory Animals, USA (Approval ID: SCXK Hubei 20080005).

4.2. Data and Experimental Design

A schematic diagram of the overall research procedure for data analysis is shown in Figure 9. Dataset and description.
The target gene expression profiles were downloaded from NCBI Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/ accessed on 1 February 2011) through accession number GSE20466 (Platform: GPL1261 [Mouse430_2] Affymetrix Mouse Genome 430 2.0 Array), which contained 12 samples in total [37,38].

4.3. Data Reprocessing and Differentially Expressed Genes (DEGs) Screening

The main objective of this part was to initially normalize the datasets’ differences and functions. The data before and normalization were shown in Figure 10, and the detailed normalized gene expression data can be found in Table S5. Only those genes meeting FDR < 0.05 and |log2 FC (fold change)| > 1 were chosen as DEGs from each group. In Inha−/− vs. WT and Inha−/− with hCG vs. WT groups only 1074 and 931 DEGs were identified based on the cut-off criteria and showed in volcano plots for WT (Figure 11A), and Inha−/− (Figure 11B) respectively. The list of DEGs could be found in Table S6 (Supplementary Materials).
Raw CEL files and annotation files were downloaded, and the gene expression data of all samples were preprocessed via background correction, quantile normalization, and probe summarization using the Robust Multi-array Average (RMA) algorithm (http://www.bioconductor.org/packages/release/bioc/html/affy.html, accessed on 1 February 2022) in R 3.4.1 (R Studio, USA). Linear Models of Microarray Data package (LIMMA, version 3.32.5) from the link http://www.bioconductor.org/packages/release/bioc/html/limma.html (accessed on 1 February 2022) was used to identify DEGs [39].

4.4. Hierarchical Clustering and Comparison Analysis of Selected DEGs in Different Groups

The expression of selected DEGs in Inha−/− vs. WT and Inha−/− (hCG) vs. WT (hCG), Inha−/− vs. Inha−/− (hCG), and WT vs. WT (hCG) were used to generate a hierarchical clustering image by heatmap (version 1.0.8) package in R 3.4.1 (RStudio; http://www.cran.r-project.org/web/packages/pheatmap/, accessed on 1 February 2022) [40,41]. Then, the identified DEGs were compared in the two groups and the ODEGs were selected by using VennDiagram package R 3.4.1 (http://www.cran.r-project.org/web/packages/VennDiagram/, accessed on 1 February 2022). After that, Pearson Correlation Coefficient (PCC) was used for further studying the ODEGs correlations following Huang da, et al. [42].

4.5. Enrichment Analysis for the Overlapping DEGs

To explore the functions of ODEGs and their pathways, the DAVID version 6.8 (Database for Annotation, Visualization and Integrated Discovery; http://www.david.ncifcrf.gov/, accessed on 1 February 2022) database was used to perform GO (Go Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analyses for ODEGs. The p-value < 0.05 and gene count ≥ 2 were set as the cut-off criteria. Furthermore, the category of enriched GO, KEGG terms, and the gene number were displayed as scatterplots by the ggplot2 package in R3.4.1 (http://www.cran.r-project.org/web/packages/ggplot2/, accessed on 1 February 2022).

4.6. Physiological Phenotypes-Related Modules and Genes Identification Based on WGCNA

Weighted Gene Co-expression Network Analysis (WGCNA) algorithm was used to investigate the co-expression modules and genes which were related to phenotypes through the WGCNA package (version 1.61) (http://www.cran.r-project.org/web/packages/WGCNA/index.html, accessed on 1 February 2022).

4.7. Co-Expression Network Construction

Based on the results of the WGCNA algorithm, only gene pairs with expression correlation coefficient > 0.8 were used to construct a gene co-expression network which was then built by Cytoscape3.2.0 (http://www.cytoscape.org/, accessed on 1 February 2022). Also, GO and KEGG pathway enrichment analysis for the genes in the co-expression network was made.

4.8. miRNA-DEGs-TF Target Regulatory Network Analysis

WEB-based Gene Set Analysis Toolkit (WebGestalt; http://www.webgestalt.org/option.php, accessed on 1 February 2022) was used to search Transcription Factors (TFs) and miRNAs that regulated the DEGs in co-expression networks. Besides, p-value < 0.05 was set as the significance cut-off criteria. As a result of TFs and miRNAs were integrated and then identified in the miRNA-DEGs-TF regulatory network. The regulatory network consisting of DEGs, miRNAs, and TFs was then constructed and visualized by Cytoscape3.2.0 (http://www.cytoscape.org/, accessed on 1 February 2022).

4.9. Principal Component Analysis (PCA) for Genes in the Regulatory Network

In order to refine genes and get the most specific ones, we further narrowed the gene range by using the PCA algorithm in the psych package (version 1.7.5) in R3.1.4 (http://www.cran.r-project.org/web/packages/psych/, accessed on 1 February 2022). Then scatterplot3d package (version 0.3-40) (http://www.cran.r-project.org/web/packages/scatterplot3d/, accessed on 1 February 2022) was used to display the effect of PCA based on the top 3 components: PC1, PC2, and PC3.

5. Conclusions

In this study, 1074 and 931 DEGs aggregates were identified in inha and wild-type. Through bioinformatics investigation 8 miRNAs, 7 TFs and 14 DEGs and 7 up-regulated genes with 20 miRNA-DEGs regulation and 36 TFs-DEGs were confirmed. This study provides potential key information for using ODEGs as biomarkers for granulosa cell cancer regulation. Also, further integration of the DEGs and the TF related to the immune response can facilitate the development of the target drugs for controlling the transcription pathways of the inhibin-deficient females.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules27175595/s1.

Author Contributions

Conceptualization, H.S.T., L.Y., M.G. and Z.u.R.; methodology, F.F. and M.S.H.; software and M.S.H.; formal analysis, A.L., M.G. and M.S.H.; validation, H.S.T., Z.u.R., M.G., F.F. and M.S.H.; investigation, F.F., H.S.T., M.G. and M.S.H.; resources, F.F., I.B. and M.S.H.; data curation, F.F., H.S.T., M.G., I.B. and M.S.H.; writing—original draft preparation, F.F., M.G., Z.u.R. and writing—review and editing, M.G., L.Y., F.F. and M.S.H.; visualization, M.G., F.F., I.B. and M.S.H.; supervision, L.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research and its APC were funded by Earmarked Fund for Modern Agro-industry Technology Research System (No. CARS-37-04B) and the National Natural Science Foundation of China (31272446).

Informed Consent Statement

Not applicable.

Data Availability Statement

Available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dufour, S.; Quérat, B.; Tostivint, H.; Pasqualini, C.; Vaudry, H.; Rousseau, K. Origin and Evolution of the Neuroendocrine Control of Reproduction in Vertebrates, With Special Focus on Genome and Gene Duplications. Physiol. Rev. 2020, 100, 869–943. [Google Scholar] [CrossRef] [PubMed]
  2. Filatov, M.; Khramova, Y.; Parshina, E.; Bagaeva, T.; Semenova, M. Influence of gonadotropins on ovarian follicle growth and development in vivo and in vitro. Zygote 2017, 25, 235–243. [Google Scholar] [CrossRef] [PubMed]
  3. Kandaraki, E.A.; Chatzigeorgiou, A.; Papageorgiou, E.; Piperi, C.; Adamopoulos, C.; Papavassiliou, A.G.; Koutsilieris, M.; Diamanti-Kandarakis, E. Advanced glycation end products interfere in luteinizing hormone and follicle stimulating hormone signaling in human granulosa KGN cells. Exp. Biol. Med. 2018, 243, 29–33. [Google Scholar] [CrossRef] [PubMed]
  4. Haas, J.; Bassil, R.; Meriano, J.; Samara, N.; Barzilay, E.; Gonen, N.; Casper, R.F. Does daily co-administration of letrozole and gonadotropins during ovarian stimulation improve IVF outcome? Reprod. Biol. Endocrinol. 2017, 15, 70. [Google Scholar] [CrossRef]
  5. Han, L.; Wu, C.; Riaz, H.; Bai, L.; Chen, J.; Zhen, Y.; Guo, A.; Yang, L. Characterization of the mechanism of inhibin alpha-subunit gene in mouse anterior pituitary cells by RNA interference. PLoS ONE 2013, 8, e74596. [Google Scholar]
  6. Calcaterra, V.; Cena, H.; Regalbuto, C.; Vinci, F.; Porri, D.; Verduci, E.; Chiara, M.; Zuccotti, G.V. The Role of Fetal, Infant, and Childhood Nutrition in the Timing of Sexual Maturation. Nutrients 2021, 13, 419. [Google Scholar] [CrossRef]
  7. Robertson, D.M.; Pruysers, E.; Jobling, T. Inhibin as a diagnostic marker for ovarian cancer. Cancer Lett. 2007, 249, 14–17. [Google Scholar] [CrossRef]
  8. Matzuk, M.M.; Finegold, M.J.; Su, J.G.; Hsueh, A.J.; Bradley, A. Alpha-inhibin is a tumour-suppressor gene with gonadal specificity in mice. Nature 1992, 360, 313–319. [Google Scholar] [CrossRef]
  9. Rathore, R.; Arora, D.; Agarwal, S.; Sharma, S. Correlation of foxl2 with inhibin and calretinin in the diagnosis of ovarian sex cord stromal tumors. Turk. J. Pathol. 2017, 33, 121–128. [Google Scholar] [CrossRef]
  10. Doroszko, M.; Chrusciel, M.; Belling, K.; Vuorenoja, S.; Dalgaard, M.; Leffers, H.; Nielsen, H.B.; Huhtaniemi, I.; Toppari, J.; Rahman, N.A. Novel genes involved in pathophysiology of gonadotropin-dependent adrenal tumors in mice. Mol. Cell. Endocrinol. 2017, 444, 9–18. [Google Scholar] [CrossRef]
  11. Haverfield, J.T.; Stanton, P.G.; Loveland, K.L.; Zahid, H.; Nicholls, P.K.; Olcorn, J.S.; Makanji, Y.; Itman, C.M.; Simpson, E.R.; Meachem, S.J. Suppression of Sertoli cell tumour development during the first wave of spermatogenesis in inhibin α-deficient mice. Reprod. Fertil. Dev. 2017, 29, 609–620. [Google Scholar] [CrossRef]
  12. Hetzler, K.L.; Hardee, J.P.; LaVoie, H.A.; Murphy, E.A.; Carson, J.A. Ovarian function’s role during cancer cachexia progression in the female mouse. Am. J. Physiol. Endocrinol. Metab. 2017, 312, E447–E459. [Google Scholar] [CrossRef]
  13. Nagaraja, A.K.; Agno, J.E.; Kumar, T.R.; Matzuk, M.M. Luteinizing hormone promotes gonadal tumorigenesis in inhibin-deficient mice. Mol. Cell. Endocrinol. 2008, 294, 19–28. [Google Scholar] [CrossRef]
  14. Kumar, T.R.; Palapattu, G.; Wang, P.; Woodruff, T.K.; Boime, I.; Byrne, M.C.; Matzuk, M.M. Transgenic models to study gonadotropin function: The role of follicle-stimulating hormone in gonadal growth and tumorigenesis. Mol. Endocrinol. 1999, 13, 851–865. [Google Scholar] [CrossRef]
  15. Chermuła, B.; Brązert, M.; Iżycki, D.; Ciesiółka, S.; Kranc, W.; Celichowski, P.; Ożegowska, K.; Nawrocki, M.J. New Gene Markers of Angiogenesis and Blood Vessels Development in Porcine Ovarian Granulosa Cells during Short-Term Primary Culture In Vitro. BioMed Res. Intl. 2019, 2019, 6545210. [Google Scholar] [CrossRef]
  16. Li, Z.; Wang, J.; Zhao, Y. scRNA-seq of ovarian follicle granulosa cells from different fertility goats reveals distinct expression patterns. Reprod. Domest. Anim. 2021, 56, 801–811. [Google Scholar] [CrossRef]
  17. Nia, A.M.; Chen, T.; Barnette, B.L.; Khanipov, K.; Ullrich, R.L.; Bhavnani, S.K.; Emmett, M.R. Efficient identification of multiple pathways: RNA-Seq analysis of livers from (56)Fe ion irradiated mice. BMC Bioinform. 2020, 21, 118. [Google Scholar] [CrossRef]
  18. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef]
  19. Gouda, M.; Huang, Z.; Liu, Y.; He, Y.; Li, X. Physicochemical impact of bioactive terpenes on the microalgae biomass structural characteristics. Bioresour. Technol. 2021, 334, 125232. [Google Scholar] [CrossRef]
  20. Gouda, M.; Chen, K.; Li, X.; Liu, Y.; He, Y. Detection of microalgae single-cell antioxidant and electrochemical potentials by gold microelectrode and Raman micro-spectroscopy combined with chemometrics. Sens. Actuators B Chem. 2021, 329, 129229. [Google Scholar] [CrossRef]
  21. Shankar, V.; Gouda, M.; Moncivaiz, J.; Gordon, A.; Reo, N.V.; Hussein, L.; Paliy, O. Differences in Gut Metabolites and Microbial Composition and Functions between Egyptian and U.S. Children Are Consistent with Their Diets. mSystems 2017, 2, e00169-16. [Google Scholar] [CrossRef] [Green Version]
  22. FarmanUllah; Liang, X.; Khan, F.; Salim, M.; Rehman, Z.; Khan, M.; Talpur, H.; Schreurs, N.; Gouda, M.; Khan, S.; et al. Transcriptomic in silico analysis of bovine Escherichia coli mastitis highlights its immune-related expressed genes as an effective biomarker. J. Genet. Eng. Biotechnol. 2021, 19, 00290–00291. [Google Scholar] [CrossRef]
  23. Ahmed, F.E.; Gouda, M.M.; Hussein, L.A.; Ahmed, N.C.; Vos, P.W.; Mohammad, M.A. Role of Melt Curve Analysis in Interpretation of Nutrigenomics’ MicroRNA Expression Data. Cancer Genom. Proteom. 2017, 14, 469–481. [Google Scholar]
  24. Fong, M.Y.; Kakar, S.S. Ovarian cancer mouse models: A summary of current models and their limitations. J. Ovarian Res. 2009, 2, 12. [Google Scholar] [CrossRef] [PubMed]
  25. Kadariya, I.; Wang, J.; ur Rehman, Z.; Ali, H.; Riaz, H.; He, J.; Bhattarai, D.; Liu, J.J.; Zhang, S.J. RNAi-mediated knockdown of inhibin α subunit increased apoptosis in granulosa cells and decreased fertility in mice. J. Steroid Biochem. Mol. Biol. 2015, 152, 161–170. [Google Scholar] [CrossRef]
  26. Vasilache, A.M.; Kugelberg, U.; Blomqvist, A.; Nilsberth, C. Minor changes in gene expression in the mouse preoptic hypothalamic region by inflammation-induced prostaglandin E2. J. Neuroendocrinol. 2013, 25, 635–643. [Google Scholar] [CrossRef]
  27. Huang, W.; Zhang, X.; Li, A.; Xie, L.; Miao, X. Genome-Wide Analysis of mRNAs and lncRNAs of Intramuscular Fat Related to Lipid Metabolism in Two Pig Breeds. Cell. Physiol. Biochem. 2018, 50, 2406–2422. [Google Scholar] [CrossRef]
  28. Chen, C.; Ahmad, M.J.; Ye, T.; Du, C.; Zhang, X.; Liang, A.; Yang, L. Cathepsin B Regulates Mice Granulosa Cells’ Apoptosis and Proliferation In Vitro. Int. J. Mol. Sci. 2021, 22, 11827. [Google Scholar] [CrossRef]
  29. Jing, R.; Gu, L.; Li, J.; Gong, Y. A transcriptomic comparison of theca and granulosa cells in chicken and cattle follicles reveals ESR2 as a potential regulator of CYP19A1 expression in the theca cells of chicken follicles. Comp. Biochem. Physiol. Part D Genom. Proteom. 2018, 27, 40–53. [Google Scholar] [CrossRef]
  30. Chen, Y.; Yang, L.; Qin, Y.; Liu, S.; Qiao, Y.; Wan, X.; Zeng, H.; Tang, X.; Liu, M.; Hou, Y. Effects of differential distributed-JUP on the malignancy of gastric cancer. J. Adv. Res. 2021, 28, 195–208. [Google Scholar] [CrossRef]
  31. Peng, M.; Zhang, H.; Jaafar, L.; Risinger, J.I.; Huang, S.; Mivechi, N.F.; Ko, L. Human ovarian cancer stroma contains luteinized theca cells harboring tumor suppressor gene GT198 mutations. J. Biol. Chem. 2013, 288, 33387–33397. [Google Scholar] [CrossRef] [PubMed]
  32. Leask, A.; Abraham, D.J. All in the CCN family: Essential matricellular signaling modulators emerge from the bunker. J. Cell Sci. 2006, 119, 4803–4810. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Fang, J.; Xiao, L. Junction plakoglobin, a potential prognostic marker of oral squamous cell carcinoma, promotes proliferation, migration and invasion. J. Oral Pathol. Med. 2020, 49, 30–38. [Google Scholar] [CrossRef] [PubMed]
  34. Zangen, D.; Kaufman, Y.; Zeligson, S.; Perlberg, S.; Fridman, H.; Kanaan, M.; Abdulhadi-Atwan, M.; Abu Libdeh, A.; Gussow, A.; Kisslov, I.; et al. XX ovarian dysgenesis is caused by a PSMC3IP/HOP2 mutation that abolishes coactivation of estrogen-driven transcription. Am. J. Hum. Genet. 2011, 89, 572–579. [Google Scholar] [CrossRef]
  35. Yin, T.; Getsios, S.; Caldelari, R.; Kowalczyk, A.P.; Müller, E.J.; Jones, J.C.R.; Green, K.J. Plakoglobin suppresses keratinocyte motility through both cell–cell adhesion-dependent and -independent mechanisms. Proc. Natl. Acad. Sci. USA 2005, 102, 5420–5425. [Google Scholar] [CrossRef]
  36. Hofland, J.; Steenbergen, J.; Voorsluijs, J.M.; Verbiest, M.M.; de Krijger, R.R.; Hofland, L.J.; de Herder, W.W.; Uitterlinden, A.G.; Feelders, R.A.; de Jong, F.H. Inhibin alpha-subunit (INHA) expression in adrenocortical cancer is linked to genetic and epigenetic INHA promoter variation. PLoS ONE 2014, 9, e104944. [Google Scholar] [CrossRef]
  37. Nagaraja, A.K.; Middlebrook, B.S.; Rajanahally, S.; Myers, M.; Li, Q.; Matzuk, M.M.; Pangas, S.A. Defective gonadotropin-dependent ovarian folliculogenesis and granulosa cell gene expression in inhibin-deficient mice. Endocrinology 2010, 151, 4994–5006. [Google Scholar] [CrossRef]
  38. Chen, L.; Zhang, W.; Huang, R.; Miao, X.; Li, J.; Yu, D.; Li, Y.; Hsu, W.; Qiu, M.; Zhang, Z.; et al. The function of Wls in ovarian development. Mol. Cell. Endocrinol. 2021, 522, 111142. [Google Scholar] [CrossRef]
  39. Wang, L.; Cao, C.; Ma, Q.; Zeng, Q.; Wang, H.; Cheng, Z.; Zhu, G.; Qi, J.; Ma, H.; Nian, H.; et al. RNA-seq analyses of multiple meristems of soybean: Novel and alternative transcripts, evolutionary and functional implications. BMC Plant Biol. 2014, 14, 169. [Google Scholar] [CrossRef]
  40. Chen, H.; Boutros, P.C. VennDiagram: A package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinform. 2011, 12, 35. [Google Scholar] [CrossRef]
  41. Diao, C.; Xi, Y.; Xiao, T. Identification and analysis of key genes in osteosarcoma using bioinformatics. Oncol. Lett. 2018, 15, 2789–2794. [Google Scholar] [CrossRef]
  42. Huang da, W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
Figure 1. Hierarchical clustering heatmap of DEGs in group Inha−/− vs. WT (a) and Inha−/− (hCG) vs. WT(hCG) (b).
Figure 1. Hierarchical clustering heatmap of DEGs in group Inha−/− vs. WT (a) and Inha−/− (hCG) vs. WT(hCG) (b).
Molecules 27 05595 g001
Figure 2. (a) Venn diagram of DEGs in groups Inha−/− vs. WT and Inha−/− (hCG) vs. WT(hCG). (b) Hierarchical clustering heatmap of overlapping DEGs. (c) Scatter-plot of correlation between logFC of Inha−/− vs. WT and Inha−/− (hCG) vs. WT(hCG). The red triangle and green inverted triangle refer to up and down-regulated DEGs in both Inha−/− vs. WT and Inha−/− (hCG) vs. WT(hCG).
Figure 2. (a) Venn diagram of DEGs in groups Inha−/− vs. WT and Inha−/− (hCG) vs. WT(hCG). (b) Hierarchical clustering heatmap of overlapping DEGs. (c) Scatter-plot of correlation between logFC of Inha−/− vs. WT and Inha−/− (hCG) vs. WT(hCG). The red triangle and green inverted triangle refer to up and down-regulated DEGs in both Inha−/− vs. WT and Inha−/− (hCG) vs. WT(hCG).
Molecules 27 05595 g002
Figure 3. The histogram of the category of enriched GO terms and KEGG pathways for the overlapped DEGs. The horizontal axis represents the number of blue, purple, green, and orange mean Biology Process, Cellular Component, Molecular Function, and pathways, respectively; the red dot curve means –log10 (p-value).
Figure 3. The histogram of the category of enriched GO terms and KEGG pathways for the overlapped DEGs. The horizontal axis represents the number of blue, purple, green, and orange mean Biology Process, Cellular Component, Molecular Function, and pathways, respectively; the red dot curve means –log10 (p-value).
Molecules 27 05595 g003
Figure 4. (a) Sample clustering tree. The Black and red bars mean different types of samples. (b) power-law coefficient parameter plot X-axis means soft-threshold, Y-axis means scale-free topology model fit signed R-square. (c) Modules clustering tree, different colors in the bottom mean different modules. (d) Module bar plot, X-axis means different modules, Y-axis means the significance of genes in each color module based on their different physiological phenotypes.
Figure 4. (a) Sample clustering tree. The Black and red bars mean different types of samples. (b) power-law coefficient parameter plot X-axis means soft-threshold, Y-axis means scale-free topology model fit signed R-square. (c) Modules clustering tree, different colors in the bottom mean different modules. (d) Module bar plot, X-axis means different modules, Y-axis means the significance of genes in each color module based on their different physiological phenotypes.
Molecules 27 05595 g004
Figure 5. Gene co-expression network based on 153 DEGs in blue, green, and brown modules. Triangle and inverted triangle refer to up and down-regulated DEGs; blue, green, and brown nodes mean genes from the corresponding colored module. Redline connections mean a positive correlation coefficient, and green line connections mean a negative correlation coefficient.
Figure 5. Gene co-expression network based on 153 DEGs in blue, green, and brown modules. Triangle and inverted triangle refer to up and down-regulated DEGs; blue, green, and brown nodes mean genes from the corresponding colored module. Redline connections mean a positive correlation coefficient, and green line connections mean a negative correlation coefficient.
Molecules 27 05595 g005
Figure 6. The scatterplot of the category of enriched GO terms and KEGG pathways for the DEGs in the co-expression network. Node size means gene count, color means p-value.
Figure 6. The scatterplot of the category of enriched GO terms and KEGG pathways for the DEGs in the co-expression network. Node size means gene count, color means p-value.
Molecules 27 05595 g006
Figure 7. miRNA-DEG-TF regulatory network of DEGs in the co-expression network. Triangle and inverted triangles refer to up and down-regulated DEGs; blue, green, and brown nodes mean genes from the corresponding colored module. Yellow square and diamond mean miRNAs and TFs.
Figure 7. miRNA-DEG-TF regulatory network of DEGs in the co-expression network. Triangle and inverted triangles refer to up and down-regulated DEGs; blue, green, and brown nodes mean genes from the corresponding colored module. Yellow square and diamond mean miRNAs and TFs.
Molecules 27 05595 g007
Figure 8. Three-dimensional distribution of samples based on PC1, PC2, and PC3. Black solid points and hollow points mean Inha−/− with and without hCG samples, and solid and hollow triangles mean WT with and without hCG samples.
Figure 8. Three-dimensional distribution of samples based on PC1, PC2, and PC3. Black solid points and hollow points mean Inha−/− with and without hCG samples, and solid and hollow triangles mean WT with and without hCG samples.
Molecules 27 05595 g008
Figure 9. Schematic illustration of the analysis strategy.
Figure 9. Schematic illustration of the analysis strategy.
Molecules 27 05595 g009
Figure 10. Boxplot of GSE20466 data preprocessing before (A) and after normalization (B). Pink, hot pink, gold, and orange boxes refer to Inha−/−(hCG), Inha−/−, WT(hCG), and WT samples.
Figure 10. Boxplot of GSE20466 data preprocessing before (A) and after normalization (B). Pink, hot pink, gold, and orange boxes refer to Inha−/−(hCG), Inha−/−, WT(hCG), and WT samples.
Molecules 27 05595 g010
Figure 11. Volcano plot of DEGs in group Inha−/− vs. WT (A) and Inha−/− (hCG) vs. WT (hCG) (B). The red triangle and green inverted triangle refer to up and down-regulated DEGs; the Red horizontal dot line means FDR = 0.05 cutoff line, and two red vertical dot lines mean logFC = 1 and logFC = −1 cutoff line.
Figure 11. Volcano plot of DEGs in group Inha−/− vs. WT (A) and Inha−/− (hCG) vs. WT (hCG) (B). The red triangle and green inverted triangle refer to up and down-regulated DEGs; the Red horizontal dot line means FDR = 0.05 cutoff line, and two red vertical dot lines mean logFC = 1 and logFC = −1 cutoff line.
Molecules 27 05595 g011
Table 1. Enriched GOs and KEGG pathways for overlapped genes.
Table 1. Enriched GOs and KEGG pathways for overlapped genes.
CategoryTermCountp-Value
Biology ProcessGO:0022402~cell cycle process174.05 × 10−6
GO:0006259~DNA metabolic process178.56 × 10−4
GO:0022403~cell cycle phase161.79 × 10−4
GO:0000279~M phase151.31 × 10−4
GO:0010033~response to organic substance152.53 × 10−2
GO:0006260~DNA replication134.14 × 10−6
GO:0051301~cell division111.21 × 10−2
GO:0001568~blood vessel development101.37 × 10−2
GO:0001944~vasculature development101.58 × 10−2
GO:0001501~skeletal system development103.33 × 10−2
Cellular Component GO:0005578~proteinaceous extracellular matrix195.11 × 10−6
GO:0005615~extracellular space194.03 × 10−3
GO:0044427~chromosomal part151.49 × 10−3
GO:0005694~chromosome157.06 × 10−3
GO:0044454~nuclear chromosome part82.25 × 10−3
GO:0000228~nuclear chromosome85.51 × 10−3
GO:0005657~replication fork65.21 × 10−5
GO:0000793~condensed chromosome64.05 × 10−2
GO:0042383~sarcolemma59.21 × 10−3
GO:0030018~Z disc44.69 × 10−2
Molecular FunctionGO:0019838~growth factor binding64.09 × 10−3
GO:0008094~DNA-dependent ATPase activity41.80 × 10−2
GO:0016875~ligase activity, forming carbon-oxygen bonds43.32 × 10−2
GO:0004812~aminoacyl-tRNA ligase activity43.32 × 10−2
GO:0005520~insulin-like growth factor binding34.08 × 10−2
KEGG Pathwaymmu03030:DNA replication81.35 × 10−6
mmu03430:Mismatch repair54.27 × 10−4
mmu04512:ECM-receptor interaction84.34 × 10−4
mmu03420:Nucleotide excision repair66.87 × 10−4
mmu03440:Homologous recombination49.92 × 10−3
mmu04510:Focal adhesion91.73 × 10−3
mmu03410:Base excision repair42.87 × 10−2
mmu00230:Purine metabolism74.71 × 10−2
Table 2. Correlation between physiological phenotypes and each module genes.
Table 2. Correlation between physiological phenotypes and each module genes.
ColorGene CountCorrelation Coefficient (R2)
blue 650.9218561
green 280.9203381
brown 600.8912573
turquoise810.8894315
yellow 510.889348
grey 580.8815949
Table 3. Enriched GOs and KEGG pathways for co-expression network genes.
Table 3. Enriched GOs and KEGG pathways for co-expression network genes.
ParameterTermCountp-Value
Biology ProcessGO:0007049~cell cycle90.012955
GO:0007155~cell adhesion80.02503
GO:0022610~biological adhesion80.025243
GO:0051301~cell division70.003211
GO:0000279~M phase70.003325
GO:0022403~cell cycle phase70.00678
GO:0022402~cell cycle process70.015593
GO:0006259~DNA metabolic process70.021132
GO:0006260~DNA replication60.001106
GO:0010817~regulation of hormone levels50.004264
GO:0009952~anterior/posterior pattern formation50.007972
GO:0000280~nuclear division50.016588
GO:0007067~mitosis50.016588
GO:0000087~M phase of mitotic cell cycle50.017769
GO:0048285~organelle fission50.018688
GO:0003002~regionalization50.024455
GO:0000278~mitotic cell cycle50.03702
GO:0042445~hormone metabolic process40.011174
GO:0030155~regulation of cell adhesion40.012565
GO:0001763~morphogenesis of a branching structure40.026619
GO:0035051~cardiac cell differentiation30.012461
KEGG pathwaymmu04512:ECM-receptor interaction40.010701
mmu04510:Focal adhesion40.009668
mmu03030:DNA replication30.015987
Table 4. Related miRNAs list.
Table 4. Related miRNAs list.
miRNAIDp-ValueFDR
mmu_TGCCTTA,MIR-124ADB_ID:5909.65 × 10−30.0014
mmu_GTGACTT,MIR-224DB_ID:5240.00020.0014
mmu_CTCTGGA,MIR-520ADB_ID:4840.00360.0126
mmu_ACCAAAG,MIR-9DB_ID:5880.00290.0126
mmu_ACTGAAA,MIR-30ADB_ID:4640.00650.0182
mmu_CTGAGCC,MIR-24DB_ID:5390.01070.0194
mmu_AACTGGA,MIR-145DB_ID:6140.01010.0194
mmu_AAGCACT,MIR-520FDB_ID:6150.01110.0194
Table 5. Related TFs list.
Table 5. Related TFs list.
TFIDp-ValueFDR
PAX4DB_ID:18308.59 × 10−62.58 × 10−5
MAZDB_ID:18157.72 × 10−62.58 × 10−5
MYCDB_ID:18194.99 × 10−62.58 × 10−5
NFATDB_ID:18221.40 × 10−53.15 × 10−5
FOXO4DB_ID:18013.72 × 10−56.70 × 10−5
SP1DB_ID:18372.00 × 10−43.00 × 10−4
LEF1DB_ID:18134.00 × 10−44.00 × 10−4
Table 6. Gene contributions to PC1-3.
Table 6. Gene contributions to PC1-3.
GeneContribution to PC1-3
Fndc50.97
Sertad40.97
Atp1b30.96
Fam57a0.95
P4htm0.89
Hoxd100.85
Psmc3ip0.85
Rab11a−0.85
Ypel5−0.9
Emx2−0.91
Jup−0.92
Gpc4−0.93
Slc25a33−0.94
Smarca1−0.94
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Talpur, H.S.; Rehman, Z.u.; Gouda, M.; Liang, A.; Bano, I.; Hussain, M.S.; FarmanUllah, F.; Yang, L. Molecular Genomic Study of Inhibin Molecule Production through Granulosa Cell Gene Expression in Inhibin-Deficient Mice. Molecules 2022, 27, 5595. https://doi.org/10.3390/molecules27175595

AMA Style

Talpur HS, Rehman Zu, Gouda M, Liang A, Bano I, Hussain MS, FarmanUllah F, Yang L. Molecular Genomic Study of Inhibin Molecule Production through Granulosa Cell Gene Expression in Inhibin-Deficient Mice. Molecules. 2022; 27(17):5595. https://doi.org/10.3390/molecules27175595

Chicago/Turabian Style

Talpur, Hira Sajjad, Zia ur Rehman, Mostafa Gouda, Aixing Liang, Iqra Bano, Mir Sajjad Hussain, FarmanUllah FarmanUllah, and Liguo Yang. 2022. "Molecular Genomic Study of Inhibin Molecule Production through Granulosa Cell Gene Expression in Inhibin-Deficient Mice" Molecules 27, no. 17: 5595. https://doi.org/10.3390/molecules27175595

APA Style

Talpur, H. S., Rehman, Z. u., Gouda, M., Liang, A., Bano, I., Hussain, M. S., FarmanUllah, F., & Yang, L. (2022). Molecular Genomic Study of Inhibin Molecule Production through Granulosa Cell Gene Expression in Inhibin-Deficient Mice. Molecules, 27(17), 5595. https://doi.org/10.3390/molecules27175595

Article Metrics

Back to TopTop