Transcriptomic analysis of granulosa cell populations proximal and distal to the germinal disc of chicken preovulatory follicles

DUSP8-attenuated ERK1/2 signaling mediates lipogenesis and steroidogenesis in chicken granulosa cells

The lipogenesis and steroidogenesis of granulosa cells are crucial during follicular development, yet it remains unclear whether dual-specificity phosphatase 8 (DUSP8) is involved. In this study, the specific role of DUSP8 in lipogenesis and steroidogenesis was investigated through culturing chicken granulosa cells in vitro. The results revealed that the expression levels of adipogenic genes were elevated after DUSP8 overexpression and reduced after knockdown. The same was observed for lipid deposition in granulosa cells. Meanwhile, the steroidogenic gene expression and progesterone synthesis were promoted after DUSP8 overexpression and inhibited after knockdown. In addition, we also found that DUSP8 blocked the phosphorylation of extracellular regulatory kinase 1/2 (ERK1/2). Based on the previous results that activated ERK1/2 signaling inhibited lipid deposition and progesterone synthesis in chicken granulosa cells, we demonstrated that DUSP8 promoted lipid deposition and progesterone synthesis through mediating the ERK1/2 signaling pathway. The results will improve our understanding of the molecular regulatory mechanisms regarding lipid metabolism and progesterone synthesis in chicken granulosa cells.

RNA-seq reveals role of cell-cycle regulating genes in the pathogenicity of a field very virulent infectious bursal disease virus

Infectious bursal disease virus (IBDV) infection causes highly contagious and immunosuppressive disease in poultry. The thymus, serving as the primary organ for T cell maturation and differentiation, plays an important role in the pathogenicity of IBDV in the infected chickens. However, there are no reports on the molecular pathogenesis of IBDV in the thymus currently. The aim of the study was to elucidate the molecular mechanisms underlying the pathogenicity of a field very virulent (vv) IBDV strain NN1172 in the thymus of SPF chickens using integrative transcriptomic and proteomic analyses. Our results showed that a total of 4,972 Differentially expressed genes (DEGs) in the thymus of NN1172-infected chickens by transcriptomic analysis, with 2,796 up-regulated and 2,176 down-regulated. Meanwhile, the proteomic analysis identified 726 differentially expressed proteins (DEPs) in the infected thymus, with 289 up-regulated and 437 down-regulated. Overall, a total of 359 genes exhibited differentially expression at both mRNA and protein levels, with 134 consistently up-regulated and 198 genes consistently down-regulated, as confirmed through a comparison of the RNA-seq and the proteomic datasets. The gene ontology (GO) analysis unveiled the involvement of both DEGs and DEPs in diverse categories encompassing cellular components, biological processes, and molecular functions in the pathological changes in IBDV-infected thymus. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the host mainly displayed severely disruption of cell survival/repair, proliferation and metabolism pathway, meanwhile, the infection triggers antiviral immune activation with a potential emphasis on the MDA5 pathway. Network inference analysis identified seven core hub genes, which include CDK1, TYMS, MCM5, KIF11, CCNB2, MAD2L1, and MCM4. These genes are all associated with cell-cycle regulating pathway and are likely key mediators in the pathogenesis induced by NN1172 infection in the thymus. This study discovered dominant pathways and genes which enhanced our understanding of the molecular mechanisms underlying IBDV pathogenesis in the thymus.

Machine learning identified MDK score has prognostic value for idiopathic pulmonary fibrosis based on integrated bulk and single cell expression data

Idiopathic pulmonary fibrosis (IPF) is a progressive and fatal lung disease that poses a significant challenge to medical professionals due to its increasing incidence and prevalence coupled with the limited understanding of its underlying molecular mechanisms. In this study, we employed a novel approach by integrating five expression datasets from bulk tissue with single-cell datasets; they underwent pseudotime trajectory analysis, switch gene selection, and cell communication analysis. Utilizing the prognostic information derived from the GSE47460 dataset, we identified 22 differentially expressed switch genes that were correlated with clinical indicators as important genes. Among these genes, we found that the midkine (MDK) gene has the potential to serve as a marker of Idiopathic pulmonary fibrosis because its cellular communicating genes are differentially expressed in the epithelial cells. We then utilized midkine and its cellular communication-related genes to calculate the midkine score. Machine learning models were further constructed through midkine and related genes to predict Idiopathic pulmonary fibrosis disease through the bulk gene expression datasets. The midkine score demonstrated a correlation with clinical indexes, and the machine learning model achieved an AUC of 0.94 and 0.86 in the Idiopathic pulmonary fibrosis classification task based on lung tissue samples and peripheral blood mononuclear cell samples, respectively. Our findings offer valuable insights into the pathogenesis of Idiopathic pulmonary fibrosis, providing new therapeutic directions and target genes for further investigation.

miR-128-3p regulates chicken granulosa cell function via 14-3-3β/FoxO and PPAR-γ/LPL signaling pathways

MicroRNAs (miRNAs) are class of 22 nt short RNA sequences which inhibit protein translation through binding to the 3'UTR of its target genes. The continuous ovulatory property of chicken follicle makes it a perfect model for studying granulosa cell (GC) functions. In this study, we found that large number of miRNAs including miR-128-3p, were differentially expressed in the GCs of F1 and F5 follicles of chicken. Subsequently, the results revealed that miR-128-3p inhibited proliferation, the formation of lipid droplets, and hormone secretion in chicken primary GCs through directly targeting YWHAB and PPAR-γ genes. To determine the effects of 14-3-3β (encoded by YWHAB) protein on GCs functions, we overexpressed or inhibited the expression of YWHAB, and the results showed that YWHAB inhibited the function of FoxO proteins. Collectively, we found that miR-128-3p was highly expressed in the chicken F1 follicles compared to the F5 follicles. In addition, the results indicated that miR-128-3p promoted GC apoptosis through 14-3-3β/FoxO pathway via repressing YWHAB, and inhibited lipid synthesis by impeding the PPAR-γ/LPL pathway, as well as reduced the secretion of progesterone and estrogen. Taken together, the results showed that miR-128-3p plays a regulatory role in chicken granulosa cell function via 14-3-3β/FoxO and PPAR-γ/LPL signaling pathways.

Quantitative Protein Analysis of ZPB2, ZPB1 and ZPC in the Germinal Disc and a Non-Germinal Disc Region of the Inner Perivitelline Layer in Two Genetic Lines of Turkey Hens That Differ in Fertility

The avian inner perivitelline layer (IPVL), containing the zona pellucida (ZP) family of proteins, surrounds the ovulated ovum. In mammalian species, ZP proteins serve as key component(s) in binding sperm and initiating the acrosome reaction. Sperm binding at the germinal disc (GD) region of the IPVL initiates fertilization in avian species, and the amount of sperm binding at the GD reflects female fertility. The current research determined whether reported differences in mRNA expression in two genetic lines of turkey hens (E, high fertility and F, low fertility) translated to the protein level. ZPB2 in the IPVL is greater in the GD region compared with the nongerminal disc (NGD) region, as indicated by both mRNA and protein expression. However, protein expressions of ZPB1 and ZPC in the IPVL of E- and F-line turkey hens was in contrast to previously reported mRNA expression. The results indicate that the mRNA expression of ZP proteins at their site of synthesis in E- and F-line hens often does not directly correlate with the IPVL abundance of these proteins. The greater protein concentration of ZPB2 in the GD region compared with the NGD regions suggests that this protein may be critical for sperm binding at the GD region.

Expression of Nerve Growth Factor and Its Receptor TrkA in the Reproductive System of Adult Zebrafish

Nerve growth factor (NGF), a member of the neurotrophin family, has emerged as an active mediator in different crucial events in the peripheral and central nervous system. At the same time, several studies showed that this neurotrophin can also play a role in non-neuronal tissues (e.g., among gonads). In spite of a large number of studies present in mammals, investigations devoted to NGF and its receptor TrkA in the reproductive system of other animal models, such as teleost fish, are scarce. To increase our knowledge of NGF and its receptor in a vertebrate gonads model, the present report describes the expression patterns of ngf and trka mRNA in the testis and ovary of adult zebrafish. By using chromogenic and fluorescence in situ hybridization, we demonstrate that in the testis of adult zebrafish, ngf and its receptor trka are mainly expressed in spermatogony B and spermatocytes. In the ovary of this fish, ngf and trka are expressed at different stages of oocyte development. Altogether, these results show that this neurotrophin and its receptor have an important role in the reproductive system that is conserved during vertebrate evolution.

Single-cell RNA-seq analysis and cell-cluster deconvolution of the human preovulatory follicular fluid cells provide insights into the pathophysiology of ovarian hyporesponse

Reduction in responsiveness to gonadotropins or hyporesponsiveness may lead to the failure of in vitro fertilization (IVF), due to a low number of retrieved oocytes. The ovarian sensitivity index (OSI) is used to reflect the ovarian responsiveness to gonadotropin stimulation before IVF. Although introduced to clinical practice already years ago, its usefulness to predict clinical outcomes requires further research. Nevertheless, pathophysiological mechanisms of ovarian hyporesponse, along with advanced maternal age and in younger women, have not been fully elucidated. Follicles consist of multiple cell types responsible for a repertoire of biological processes including responding to pituitary gonadotropins necessary for follicle growth and oocyte maturation as well as ovulation. Encouraging evidence suggests that hyporesponse could be influenced by many contributing factors, therefore, investigating the variability of ovarian follicular cell types and their gene expression in hyporesponders is highly informative for increasing their prognosis for IVF live birth. Due to advancements in single-cell analysis technologies, the role of somatic cell populations in the development of infertility of ovarian etiology can be clarified. Here, somatic cells were collected from the fluid of preovulatory ovarian follicles of patients undergoing IVF, and RNA-seq was performed to study the associations between OSI and gene expression. We identified 12 molecular pathways differentially regulated between hypo- and normoresponder patient groups (FDR<0.05) from which extracellular matrix organization, post-translational protein phosphorylation, and regulation of Insulin-like Growth Factor (IGF) transport and uptake by IGF Binding Proteins were regulated age-independently. We then generated single-cell RNA-seq data from matching follicles revealing 14 distinct cell clusters. Using cell cluster-specific deconvolution from the bulk RNA-seq data of 18 IVF patients we integrated the datasets as a novel approach and discovered that the abundance of three cell clusters significantly varied between hypo- and normoresponder groups suggesting their role in contributing to the deviations from normal ovarian response to gonadotropin stimulation. Our work uncovers new information regarding the differences in the follicular gene expression between hypo- and normoresponders. In addition, the current study fills the gap in understanding the inter-patient variability of cell types in human preovulatory follicles, as revealed by single-cell analysis of follicular fluid cells.

Characterization of Four Orphan Receptors (GPR3, GPR6, GPR12 and GPR12L) in Chickens and Ducks and Regulation of GPR12 Expression in Ovarian Granulosa Cells by Progesterone

The three structurally related orphan G protein-coupled receptors, GRP3, GPR6, and GPR12, are reported to be constitutively active and likely involved in the regulation of many physiological/pathological processes, such as neuronal outgrowth and oocyte meiotic arrest in mammals. However, the information regarding these orphan receptors in nonmammalian vertebrates is extremely limited. Here, we reported the structure, constitutive activity, and tissue expression of these receptors in two representative avian models: chickens and ducks. The cloned duck GPR3 and duck/chicken GPR6 and GPR12 are intron-less and encode receptors that show high amino acid (a.a.) sequence identities (66–88%) with their respective mammalian orthologs. Interestingly, a novel GPR12-like receptor (named GPR12L) sharing 66% a.a. identity to that in vertebrates was reported in the present study. Using dual-luciferase reporter assay and Western blot, we demonstrated that GPR3, GPR6, GPR12, and GPR12L are constitutively active and capable of stimulating the cAMP/PKA signaling pathway without ligand stimulation in birds (and zebrafish), indicating their conserved signaling property across vertebrates. RNA-seq data/qRT-PCR assays revealed that GPR6 and GPR12L expression is mainly restricted to the chicken brain, while GPR12 is highly expressed in chicken ovarian granulosa cells (GCs) and oocytes of 6 mm growing follicles and its expression in cultured GCs is upregulated by progesterone. Taken together, our data reveal the structure, function, and expression of GPR3, GPR6, GPR12, and GPR12L in birds, thus providing the first piece of evidence that GPR12 expression is upregulated by gonadal steroid (i.e., progesterone) in vertebrates.