scholarly journals Transcriptome Analyses Identify Potential Key microRNAs and Their Target Genes Contributing to Ovarian Reserve

2021 ◽  
Vol 22 (19) ◽  
pp. 10819
Author(s):  
Yoon-Young Kim ◽  
Kwang-Soo Kim ◽  
Yong-Jin Kim ◽  
Sung-Woo Kim ◽  
Hoon Kim ◽  
...  

Female endocrinological symptoms, such as premature ovarian inefficiency (POI) are caused by diminished ovarian reserve and chemotherapy. The etiology of POI remains unknown, but this can lead to infertility. This has accelerated the search for master regulator genes or other molecules that contribute as enhancers or silencers. The impact of regulatory microRNAs (miRNAs) on POI has gained attention; however, their regulatory function in this condition is not well known. RNA sequencing was performed at four stages, 2-(2 W), 6-(6 W), 15-(15 W), and 20-(20 W) weeks, on ovarian tissue samples and 5058 differentially expressed genes (DEGs) were identified. Gene expression and enrichment were analyzed based on the gene ontology and KEGG databases, and their association with other proteins was assessed using the STRING database. Gene set enrichment analysis was performed to identify the key target genes. The DEGs were most highly enriched in 6 W and 15 W groups. Figla, GDF9, Nobox, and Pou51 were significantly in-creased at 2 W compared with levels at 6 W and 20 W, whereas the expression of Foxo1, Inha, and Taf4b was significantly de-creased at 20 W. Ccnd2 and Igf1 expression was maintained at similar levels in each stage. In total, 27 genes were upregulated and 26 genes interacted with miRNAs; moreover, stage-specific upregulated and downregulated interactions were demonstrated. Increased and decreased miRNAs were identified at each stage in the ovaries. The constitutively expressed genes, Ccnd2 and Igf1, were identified as the major targets of many miRNAs (p < 0.05), and Fshr and Foxo3 interacted with miRNAs, namely mmu-miR-670-3p and mmu-miR-153-3p. miR-26a-5p interacted with Piwil2, and its target genes were downregulated in the 20 W mouse ovary. In this study, we aimed to identify key miRNAs and their target genes encompassing the reproductive span of mouse ovaries using mRNA and miRNA sequencing. These results indicated that gene sets are regulated in the reproductive stage-specific manner via interaction with miRNAs. Furthermore, consistent expression of Ccnd2 and Igf1 is considered crucial for the ovarian reserve and is regulated by many interactive miRNAs.

2021 ◽  
Vol 9 (1) ◽  
Author(s):  
E. Krasniqi ◽  
A. Sacconi ◽  
D. Marinelli ◽  
L. Pizzuti ◽  
M. Mazzotta ◽  
...  

Abstract Background In Western countries, ovarian cancer (OC) still represents the leading cause of gynecological cancer-related deaths, despite the remarkable gains in therapeutical options. Novel biomarkers of early diagnosis, prognosis definition and prediction of treatment outcomes are of pivotal importance. Prior studies have shown the potentials of micro-ribonucleic acids (miRNAs) as biomarkers for OC and other cancers. Methods We focused on the prognostic and/or predictive potential of miRNAs in OC by conducting a comprehensive array profiling of miRNA expression levels in ovarian tissue samples from 17 non-neoplastic controls, and 60 tumor samples from OC patients treated at the Regina Elena National Cancer Institute (IRE). A set of 54 miRNAs with differential expression in tumor versus normal samples (T/N-deregulated) was identified in the IRE cohort and validated against data from the Cancer Genoma Atlas (TCGA) related to 563 OC patients and 8 non-neoplastic controls. The prognostic/predictive role of the selected 54 biomarkers was tested in reference to survival endpoints and platinum resistance (P-res). Results In the IRE cohort, downregulation of the 2 miRNA-signature including miR-99a-5p and miR-320a held a negative prognostic relevance, while upregulation of miR-224-5p was predictive of less favorable event free survival (EFS) and P-res. Data from the TCGA showed that downregulation of 5 miRNAs, i.e., miR-150, miR-30d, miR-342, miR-424, and miR-502, was associated with more favorable EFS and overall survival outcomes, while miR-200a upregulation was predictive of P-res. The 9 miRNAs globally identified were all included into a single biologic signature, which was tested in enrichment analysis using predicted/validated miRNA target genes, followed by network representation of the miRNA-mRNA interactions. Conclusions Specific dysregulated microRNA sets in tumor tissue showed predictive/prognostic value in OC, and resulted in a promising biological signature for this disease.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Pusheng Quan ◽  
Kai Wang ◽  
Shi Yan ◽  
Shirong Wen ◽  
Chengqun Wei ◽  
...  

AbstractThis study aimed to identify potential novel drug candidates and targets for Parkinson’s disease. First, 970 genes that have been reported to be related to PD were collected from five databases, and functional enrichment analysis of these genes was conducted to investigate their potential mechanisms. Then, we collected drugs and related targets from DrugBank, narrowed the list by proximity scores and Inverted Gene Set Enrichment analysis of drug targets, and identified potential drug candidates for PD treatment. Finally, we compared the expression distribution of the candidate drug-target genes between the PD group and the control group in the public dataset with the largest sample size (GSE99039) in Gene Expression Omnibus. Ten drugs with an FDR < 0.1 and their corresponding targets were identified. Some target genes of the ten drugs significantly overlapped with PD-related genes or already known therapeutic targets for PD. Nine differentially expressed drug-target genes with p < 0.05 were screened. This work will facilitate further research into the possible efficacy of new drugs for PD and will provide valuable clues for drug design.


2021 ◽  
Vol 20 ◽  
pp. 153303382098682
Author(s):  
Zhipeng Zhu ◽  
Jiuhua Xu ◽  
Xiaofang Wu ◽  
Sihao Lin ◽  
Lulu Li ◽  
...  

Background: ADAMTS5 has different roles in multiple types of cancers and participates in various molecular mechanisms. However, the prognostic value of ADAMTS5 in patients with hepatocellular carcinoma (HCC) still remains unclear. We carried the study to evaluate the prognostic value and identified underlying molecular mechanisms in HCC. Methods: Firstly, the association of ADAMTS5 expression and clinicopathological parameters was evaluated by in GSE14520. Next, ADAMTS5 expression in HCC was performed using GSE14520, GSE36376, GSE76427 and The Cancer Genome Atlas (TCGA) profile. Furthermore, Kaplan-Meier analysis, Univariate and Multivariate Cox regression analysis, subgroup analysis was performed to evaluate the prognostic value of ADAMTS5 in HCC. Finally, GO enrichment analysis, gene set enrichment analysis (GSEA) and weighted gene co-expression network analysis (WGCNA) were performed to revealed underlying molecular mechanisms. Result: The expression of ADAMTS5 was positively correlated with the development of HCC. Next, high ADAMTS5 expression was significantly associated with poorer survival (all P < 0.05) and the impact of ADAMTS5 on all overall survival (OS), disease-free survival (DFS), relapse-free survival (RFS), disease specific survival (DSS) and progression free interval (PFI) was specific for HCC among other 29 cancer types. Subgroup analysis showed that ADAMTS5 overexpression was significantly associated with poorer OS in patients with HCC. Finally, ADAMTS5 might participate in the status conversion from metabolic-dominant to extracellular matrix-dominant, and the activation of ECM-related biological process might contribute to high higher mortality risk for patients with HCC. Conclusion: ADAMTS5 may play an important role in the progression of HCC, and may be considered as a novel and effective biomarker for predicting prognosis for patients with HCC.


2021 ◽  
Vol 12 ◽  
Author(s):  
Yani Dong ◽  
Likang Lyu ◽  
Daiqiang Zhang ◽  
Jing Li ◽  
Haishen Wen ◽  
...  

Long non-coding RNAs (lncRNAs) have been reported to be involved in multiple biological processes. However, the roles of lncRNAs in the reproduction of half-smooth tongue sole (Cynoglossus semilaevis) are unclear, especially in the molecular regulatory mechanism driving ovarian development and ovulation. Thus, to explore the mRNA and lncRNA mechanisms regulating reproduction, we collected tongue sole ovaries in three stages for RNA sequencing. In stage IV vs. V, we identified 312 differentially expressed (DE) mRNAs and 58 DE lncRNAs. In stage V vs. VI, we identified 1,059 DE mRNAs and 187 DE lncRNAs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses showed that DE mRNAs were enriched in ECM-receptor interaction, oocyte meiosis and steroid hormone biosynthesis pathways. Furthermore, we carried out gene set enrichment analysis (GSEA) to identify potential reproduction related-pathways additionally, such as fatty metabolism and retinol metabolism. Based on enrichment analysis, DE mRNAs with a potential role in reproduction were selected and classified into six categories, including signal transduction, cell growth and death, immune response, metabolism, transport and catabolism, and cell junction. The interactions of DE lncRNAs and mRNAs were predicted according to antisense, cis-, and trans-regulatory mechanisms. We constructed a competing endogenous RNA (ceRNA) network. Several lncRNAs were predicted to regulate genes related to reproduction including cyp17a1, cyp19a1, mmp14, pgr, and hsd17b1. The functional enrichment analysis of these target genes of lncRNAs revealed that they were involved in several signaling pathways, such as the TGF-beta, Wnt signaling, and MAPK signaling pathways and reproduction related-pathways such as the progesterone-mediated oocyte maturation, oocyte meiosis, and GnRH signaling pathway. RT-qPCR analysis showed that two lncRNAs (XR_522278.2 and XR_522171.2) were mainly expressed in the ovary. Dual-fluorescence in situ hybridization experiments showed that both XR_522278.2 and XR_522171.2 colocalized with their target genes cyp17a1 and cyp19a1, respectively, in the follicular cell layer. The results further demonstrated that lncRNAs might be involved in the biological processes by modulating gene expression. Taken together, this study provides lncRNA profiles in the ovary of tongue sole and further insight into the role of lncRNA involvement in regulating reproduction in tongue sole.


Blood ◽  
2018 ◽  
Vol 132 (23) ◽  
pp. 2470-2483 ◽  
Author(s):  
Shiro Tara ◽  
Yusuke Isshiki ◽  
Yaeko Nakajima-Takagi ◽  
Motohiko Oshima ◽  
Kazumasa Aoyama ◽  
...  

Abstract BCOR, encoding BCL-6 corepressor (BCOR), is X-linked and targeted by somatic mutations in various hematological malignancies including myelodysplastic syndrome (MDS). We previously reported that mice lacking Bcor exon 4 (BcorΔE4/y) in the hematopoietic compartment developed NOTCH-dependent acute T-cell lymphoblastic leukemia (T-ALL). Here, we analyzed mice lacking Bcor exons 9 and 10 (BcorΔE9-10/y), which express a carboxyl-terminal truncated BCOR that fails to interact with core effector components of polycomb repressive complex 1.1. BcorΔE9-10/y mice developed lethal T-ALL in a similar manner to BcorΔE4/y mice, whereas BcorΔE9-10/y hematopoietic cells showed a growth advantage in the myeloid compartment that was further enhanced by the concurrent deletion of Tet2. Tet2Δ/ΔBcorΔE9-10/y mice developed lethal MDS with progressive anemia and leukocytopenia, inefficient hematopoiesis, and the morphological dysplasia of blood cells. Tet2Δ/ΔBcorΔE9-10/y MDS cells reproduced MDS or evolved into lethal MDS/myeloproliferative neoplasms in secondary recipients. Transcriptional profiling revealed the derepression of myeloid regulator genes of the Cebp family and Hoxa cluster genes in BcorΔE9-10/y progenitor cells and the activation of p53 target genes specifically in MDS erythroblasts where massive apoptosis occurred. Our results reveal a tumor suppressor function of BCOR in myeloid malignancies and highlight the impact of Bcor insufficiency on the initiation and progression of MDS.


2013 ◽  
Vol 305 (1) ◽  
pp. G58-G65 ◽  
Author(s):  
Yu Fang ◽  
Hao Chen ◽  
Yuhui Hu ◽  
Zorka Djukic ◽  
Whitney Tevebaugh ◽  
...  

The barrier function of the esophageal epithelium is a major defense against gastroesophageal reflux disease. Previous studies have shown that reflux damage is reflected in a decrease in transepithelial electrical resistance associated with tight junction alterations in the esophageal epithelium. To develop novel therapies, it is critical to understand the molecular mechanisms whereby contact with a refluxate impairs esophageal barrier function. In this study, surgical models of duodenal and mixed reflux were developed in mice. Mouse esophageal epithelium was analyzed by gene microarray. Gene set enrichment analysis showed upregulation of inflammation-related gene sets and the NF-κB pathway due to reflux. Significance analysis of microarrays revealed upregulation of NF-κB target genes. Overexpression of NF-κB subunits (p50 and p65) and NF-κB target genes (matrix metalloproteinases-3 and -9, IL-1β, IL-6, and IL-8) confirmed activation of the NF-κB pathway in the esophageal epithelium. In addition, real-time PCR, Western blotting, and immunohistochemical staining also showed downregulation and mislocalization of claudins-1 and -4. In a second animal experiment, treatment with an NF-κB inhibitor, BAY 11-7085 (20 mg·kg−1·day−1 ip for 10 days), counteracted the effects of duodenal and mixed reflux on epithelial resistance and NF-κB-regulated cytokines. We conclude that gastroesophageal reflux activates the NF-κB pathway and impairs esophageal barrier function in mice and that targeting the NF-κB pathway may strengthen esophageal barrier function against reflux.


2021 ◽  
Vol 5 (Supplement_1) ◽  
pp. A752-A753
Author(s):  
Lei Su ◽  
Soravis Osataphan ◽  
Jessica Desmond ◽  
Rui Fang ◽  
Jeremy Chimene-Weiss ◽  
...  

Abstract We previously demonstrated that paternal obesity is associated with offspring metabolic risk during later life, and that paternal SGLT2i treatment improves offspring metabolic phenotypes. Since the placenta is a key determinant of prenatal growth and development, we hypothesized the placenta could mediate the impact of paternal obesity and paternal SGLT2i treatment. Male C57BL/6J mice were fed standard chow (Purina 9F) or 60% high-fat diet (HFD, D12492, Research Diet), or 60% HFD plus the SGLT2 inhibitor canagliflozin (CANA, 25 mg/kg/d) for 4 weeks before mating with chow-fed females. Placenta were collected on E16.5, and RNA-seq was performed on placenta from male offspring (paternal chow, pChow, n=4, pHFD, n=5, and pHFD+CANA, n=4), and differentially expressed genes were identified using Limma. Placenta weight was significantly lower in pHFD (0.089±0.004 g, 7 litters from 6 fathers) vs. both pChow (0.108±0.011 g, 4 litters, 4 fathers) and pHFD+CANA (0.107±0.013 g, 5 litters, 5 fathers)(p&lt;0.05). Litter size, fetal or liver weight, or fetal/placental weight ratio did not differ between groups. No genes were differentially expressed in pHFD vs. pChow (FDR&lt;0.1). Gene set enrichment analysis (GSEA) identified significance (FDR&lt;0.05, NES&gt;1.8) for gene sets in steroid metabolic, drug catabolic, and protein-containing complex remodeling processes. Genes responsible for enrichment included cholesterol biosynthesis (Hmgcs1), transport (Apob, Apoa1/2/4, Apom, Apoc1, Vldlr, Pcsk9) and steroid hormone biosynthesis genes (Hsd3b1, Cyp11b1), all upregulated in pHFD by 1.5-3-fold. These results suggest pHFD could potentially affect maternal and fetal cholesterol homeostasis. pHFD+CANA altered expression of 154 genes vs. pHFD (7 up-, 147 down, FDR &lt;0.1, FC &gt;|1.5|); 18 gene sets were downregulated by pHFD+CANA (GSEA NES&lt;-1.8 and FDR&lt;0.05), including the 3 sets upregulated by pHFD. ChEA3 enrichment analysis (ENCODE library) predicted regulation by transcription factors important for cholesterol and sterol biosynthesis (Srebf1/2), embryonic development (Foxa2), & glucose homeostasis (Hnf4g), suggesting these pathways could mediate the “rescue” effect of pHFD+CANA (FDR&lt;0.05). Expression of Foxa2 was significantly downregulated (4-fold) in pHFD+CANA vs. pHFD. We independently analyzed expression of the 78 detected imprinted genes. None were significantly different in pHFD, but both paternally expressed (Nnat) and maternally expressed genes (H19, Phlda2, Meg3, Meg8) were downregulated in pHFD+CANA vs. pHFD by 1.4 to 3.8 fold in pHFD+CANA (p&lt;0.001,FDR&lt;0.1). In summary, paternal SGLT2i reversed the impact of pHFD on placental weight. Robust impact of both pHFD and pSGLT2i on the transcriptome suggests that the placenta is a key mediator of paternal metabolic effects on offspring development and metabolic disease risk, potentially via modification of lipid transport.


2021 ◽  
Vol 12 ◽  
Author(s):  
Michal Marczyk ◽  
Agnieszka Macioszek ◽  
Joanna Tobiasz ◽  
Joanna Polanska ◽  
Joanna Zyla

A typical genome-wide association study (GWAS) analyzes millions of single-nucleotide polymorphisms (SNPs), several of which are in a region of the same gene. To conduct gene set analysis (GSA), information from SNPs needs to be unified at the gene level. A widely used practice is to use only the most relevant SNP per gene; however, there are other methods of integration that could be applied here. Also, the problem of nonrandom association of alleles at two or more loci is often neglected. Here, we tested the impact of incorporation of different integrations and linkage disequilibrium (LD) correction on the performance of several GSA methods. Matched normal and breast cancer samples from The Cancer Genome Atlas database were used to evaluate the performance of six GSA algorithms: Coincident Extreme Ranks in Numerical Observations (CERNO), Gene Set Enrichment Analysis (GSEA), GSEA-SNP, improved GSEA for GWAS (i-GSEA4GWAS), Meta-Analysis Gene-set Enrichment of variaNT Associations (MAGENTA), and Over-Representation Analysis (ORA). Association of SNPs to phenotype was calculated using modified McNemar’s test. Results for SNPs mapped to the same gene were integrated using Fisher and Stouffer methods and compared with the minimum p-value method. Four common measures were used to quantify the performance of all combinations of methods. Results of GSA analysis on GWAS were compared to the one performed on gene expression data. Comparing all evaluation metrics across different GSA algorithms, integrations, and LD correction, we highlighted CERNO, and MAGENTA with Stouffer as the most efficient. Applying LD correction increased prioritization and specificity of enrichment outcomes for all tested algorithms. When Fisher or Stouffer were used with LD, sensitivity and reproducibility were also better. Using any integration method was beneficial in comparison with a minimum p-value method in specific combinations. The correlation between GSA results from genomic and transcriptomic level was the highest when Stouffer integration was combined with LD correction. We thoroughly evaluated different approaches to GSA in GWAS in terms of performance to guide others to select the most effective combinations. We showed that LD correction and Stouffer integration could increase the performance of enrichment analysis and encourage the usage of these techniques.


Cells ◽  
2021 ◽  
Vol 10 (11) ◽  
pp. 3065
Author(s):  
Shreya Udawant ◽  
Carl Litif ◽  
Alma Lopez ◽  
Bonnie Gunn ◽  
Erin Schuenzel ◽  
...  

Glioblastoma (GBM) is the most lethal primary brain cancer that lacks effective molecular targeted therapies. The PI3K/AKT/mTOR pathway is activated in 90% of all Glioblastoma multiforme (GBM) tumors. To gain insight into the impact of the PI3K pathway on GBM metabolism, we treated U87MG GBM cells with NVP-BEZ235 (PI3K and mTOR a dual inhibitor) and identified differentially expressed genes with RNA-seq analysis. RNA-seq identified 7803 differentially regulated genes in response to NVP-BEZ235. Gene Set Enrichment Analysis (GSEA) identified two glycolysis-related gene sets that were significantly enriched (p < 0.05) in control samples compared to NVP-BEZ235-treated samples. We validated the inhibition of glycolytic genes by NVP-BEZ235 and examined the impact of the FOXO1 inhibitor (AS1842856) on these genes in a set of GBM cell lines. FOXO1 inhibition alone was associated with reduced LDHA expression, but not ENO1 or PKM2. Bioinformatics analyses revealed that PI3K-impacted glycolytic genes were over-expressed and co-expressed in GBM clinical samples. The elevated expression of PI3K-impacted glycolytic genes was associated with poor prognosis in GBM based on Kaplan–Meier survival analyses. Our results suggest novel insights into hallmark metabolic reprogramming associated with the PI3K-mTOR dual inhibition.


2009 ◽  
Vol 31 (5) ◽  
pp. 345-356
Author(s):  
Remond J. A. Fijneman ◽  
Lindsey K. Bade ◽  
Johannes R. Peham ◽  
Mark A. van de Wiel ◽  
Victor W. M. van Hinsbergh ◽  
...  

Background: The group IIA secretory phospholipase A2 gene, Pla2g2a, confers resistance to intestinal tumorigenesis in the ApcMin/+ mouse model. However, it is unclear how Pla2g2a exerts its tumor-suppressive effects and whether its mode of action depends on Apc-germline mutations.Methods: We tested whether expression of a Pla2g2a transgene provides protection against carcinogen-induced colon tumors, and examined whether the normal colon microenvironment is modulated by Pla2g2a expression.Results: Pla2g2a strongly inhibited colon tumorigenesis in mice following treatment with the DNA alkylating agent azoxymethane (AOM). Moreover, AOM-induced duodenal tumors were also attenuated by Pla2g2a expression. These tumors demonstrated upregulation of β-catenin, indicative of involvement of the Wnt signaling pathway. Comparison of genome-wide microarray expression profiles of healthy (non-pathologic) colon tissues from Pla2g2a-transgenic to non-transgenic mice revealed 382 genes that were differentially expressed, comprising clusters of genes involved in inflammation and microbial defense, cell signaling and cell cycle, transactivation, apoptosis and mitochondrial function, DNA repair, and lipid and energy metabolism. Pathway analysis using Gene Set Enrichment Analysis (GSEA) indicated that Pla2g2a suppresses the expression of interferon-induced genes.Conclusion: Our results demonstrate that Pla2g2a attenuates colon tumorigenesis independent of Apc-germline mutations, and reveal Pla2g2a target genes and pathways in non-pathologic colon microenvironment that influence conditions for colorectal cancer development.


Sign in / Sign up

Export Citation Format

Share Document