Common miRNAs,candidate genes and their interaction network across four subtypes of epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is categorized into four major histological subtypes such as clear cell carcinoma (CCC), endometrioid carcinoma (EC), mucinous carcinoma (MC), and serous carcinoma (SC). Heterogeneity of the EOC leads to different clinical outcomes of the disease, although all the subtypes are originated from the same layer of tissue. Therefore, it is of interest to identify the common candidate genes, miRNA and their interaction network in four the subtypes of EOC. A comparative gene expression analysis identified 248 common differentially expressed genes (DEGs) in the four subtypes of EOC. Identified common DEGs were found to be enriched in cancer specific pathways. A protein-protein interaction (PPI) network of the common DEGs were constructed, and subsequent module and survival analyses identified seven key candidate genes (CCNB1, CENPM, CEP55, RACGAP1, TPX2, UBE2C, and ZWINT). We also documented 10 key candidate miRNAs (hsa-mir-16-5p, hsa-mir-23b-3p, hsa-mir-34a-5p, hsa-mir-103a-3p, hsa-mir-107, hsa-mir-124-3p, hsa-mir-129-2-3p, hsa-mir-147a, hsa-mir-205-5p, and hsa-mir-195-5p) linked to the candidate genes. These derived data find application in the understanding of EOC.     

microRNA-16-5p-containing exosomes derived from bone marrow-derived mesenchymal stem cells inhibit proliferation,migration, and invasion,while promoting apoptosis of colorectal cancer cells by downregulating ITGA2

Colorectal cancer (CRC) is a form of cancer developing from either the colon or rectum. Nowadays, research supports the functionality of exosome expressing microRNAs (miRNAs) as potential biomarker for various cancers including CRC. This study was performed with the intent of investigating the roles of both bone marrow-derived mesenchymal stem cells (BMSCs) and exosomal miR-16-5p in CRC by regulating integrin α2 (ITGA2). A microarray-based analysis was conducted to screen the CRC-associated differentially expressed genes (DEGs) as well as potential regulatory miRNAs. Next, the role of miR-16-5p in terms of its progression in association with CRC was determined. Subsequently, CRC cells were exposed to exosomes secreted by BMSCs transfected with miR-16-5p, isolated and cocultured with CRC cells in an attempt to identify the role of exosomes. Effects of BMSCs-derived exosomes overexpressing miR-16-5p on biological functions of CRC cells and tumorigenicity were all subsequently detected. Effects of miR-16-5p treated with CRC cells in regard to CRC in vivo were also measured. ITGA2 was overexpressed, while miR-16-5p was poorly expressed in CRC cells and miR-16-5p targeted ITGA2. The in vitro experiments revealed that the BMSCs-derived exosomes overexpressing miR-16-5p inhibited proliferation, migration, and invasion, while simultaneously stimulating the apoptosis of the CRC cells via downregulation of ITGA2. Furthermore, the results of in vivo experiments confirmed that the BMSCs-derived exosomes overexpressing miR-16-5p repressed the tumor growth of CRC. Collectively, BMSCs-derived exosomes overexpressing miR-16-5p restricted the progression of CRC by downregulating ITGA2.     

Profiling the circulating miRNAs in mice exposed to gram-positive and gram-negative bacteria by Illumina small RNA deep sequencing

Background

We profiled the expression of circulating microRNAs (miRNAs) in mice using Illumina small RNA deep sequencing in order to identify the miRNAs that may potentially be used as biomarkers to distinguish between gram-negative and gram-positive bacterial infections.

Results

Recombinant-specific gram-negative pathogen Escherichia coli (Xen14) and gram-positive pathogen Staphylococcus aureus (Xen29) were used to induce bacterial infection in mice at a concentration of 1 × 10⁸ bacteria/100 μL of phosphate buffered saline (PBS). Small RNA libraries generated from the serum of mice after exposure to PBS, Xen14, Xen29, and Xen14 + Xen29 via the routes of subcutaneous injection (I), cut wound (C), or under grafted skin (S) were analyzed using an Illumina HiSeq2000 Sequencer. Following exposure to gram-negative bacteria alone, no differentially expressed miRNA was found in the injection, cut, or skin graft models. Exposure to mixed bacteria induced a similar expression pattern of the circulating miRNAs to that induced by gram-positive bacterial infection. Upon gram-positive bacterial infection, 9 miRNAs (mir-193b-3p, mir-133a-1-3p, mir-133a-2-3p, mir-133a-1-5p, mir-133b-3p, mir-434-3p, mir-127-3p, mir-676-3p, mir-215-5p) showed upregulation greater than 4-fold with a p-value < 0.01. Among them, mir-193b-3p, mir-133a-1-3p, and mir-133a-2-3p presented the most common miRNA targets expressed in the mice exposed to gram-positive bacterial infection.

Conclusions

This study identified mir-193b-3p, mir-133a-1-3p, and mir-133a-2-3p as potential circulating miRNAs for gram-positive bacterial infections.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-014-0106-y) contains supplementary material, which is available to authorized users.     

Identification of novel biomarkers and key pathways of condyloma acuminata

Condyloma acuminata (CA) is a prevalent sexually transmitted disease, associated with human papilloma viruses (HPV) infections and host immune status. In this present study, we aimed to explore immune landscape and biomarkers for CA prevention and treatment. We obtained differentially expressed genes (DEGs) of CA vs normal tissues in GSE140662 and screened out hub genes from the protein–protein interaction (PPI) network. Hub genes were then subjected to microRNA (miRNA) analysis. Besides, CCK-8, transwell, flow cytometry assays were employed to assess the cell proliferation, migration and apoptosis in Hela cells. ImmuCellAI was firstly applied to identify immune cell infiltration levels of CA. We obtained 275 DEGs, 23 hub genes and key miRNAs. Subsequently, we verified four up-regulated hub genes IFIT1, IFI27, OASL, SAMD9L and down-regulated mir-146a-5p in CA tissues by RT-qPCR. Moreover, over-expression of miR-146a-5p reduced Hela cells proliferation, migration, blocked cell cycle and induced apoptosis. Up-regulated miR-146a-5p attenuated PI3K/AKT and activated p38/ MAPK signaling pathway. Proportions of Monocyte, NK cells, Gamma delta cells, Th17 cells were relatively low, while Th1 and CD8+ T cells were relatively high in CA skin. Our study revealed that mir-146a-5p contribute to CA progression through PI3K/AKT and p38/MAPK signaling pathway.     

Identification of Restricted Subsets of Mature microRNA Abnormally Expressed in Inactive Colonic Mucosa of Patients with Inflammatory Bowel Disease

Background

Ulcerative Colitis (UC) and Crohn''s Disease (CD) are two chronic Inflammatory Bowel Diseases (IBD) affecting the intestinal mucosa. Current understanding of IBD pathogenesis points out the interplay of genetic events and environmental cues in the dysregulated immune response. We hypothesized that dysregulated microRNA (miRNA) expression may contribute to IBD pathogenesis. miRNAs are small, non-coding RNAs which prevent protein synthesis through translational suppression or mRNAs degradation, and regulate several physiological processes.

Methodology/Findings

Expression of mature miRNAs was studied by Q-PCR in inactive colonic mucosa of patients with UC (8), CD (8) and expressed relative to that observed in healthy controls (10). Only miRNAs with highly altered expression (>5 or <0.2 -fold relative to control) were considered when Q-PCR data were analyzed. Two subsets of 14 (UC) and 23 (CD) miRNAs with highly altered expression (5.2->100 -fold and 0.05–0.19 -fold for over- and under- expression, respectively; 0.001<p≤0.05) were identified in quiescent colonic mucosa, 8 being commonly dysregulated in non-inflamed UC and CD (mir-26a,-29a,-29b,-30c,-126*,-127-3p,-196a,-324-3p). Several miRNA genes with dysregulated expression co-localize with acknowledged IBD-susceptibility loci while others, (eg. clustered on 14q32.31), map on chromosomal regions not previously recognized as IBD-susceptibility loci. In addition, in silico clustering analysis identified 5 miRNAs (mir-26a,-29b,-126*,-127-3p,-324-3p) that share coordinated dysregulation of expression both in quiescent and in inflamed colonic mucosa of IBD patients. Six miRNAs displayed significantly distinct alteration of expression in non-inflamed colonic biopsies of UC and CD patients (mir-196b,-199a-3p,-199b-5p,-320a,-150,-223).

Conclusions/Significance

Our study supports miRNAs as crucial players in the onset and/or relapse of inflammation from quiescent mucosal tissues in IBD patients. It allows speculating a role for miRNAs as contributors to IBD susceptibility and suggests that some of the miRNA with altered expression in the quiescent mucosa of IBD patients may define miRNA signatures for UC and CD and help develop new diagnostic biomarkers.     

Identifying potential biomarkers in hepatitis B virus infection and its response to the antiviral therapy by integrated bioinformatic analysis

The antiviral treatment efficacy varies among chronic hepatitis B (CHB) patients and the underlying mechanism is unclear. An integrated bioinformatics analysis was performed to investigate the host factors that affect the therapeutic responsiveness in CHB patients. Four GEO data sets (GSE54747, GSE27555, GSE66698 and GSE66699) were downloaded from the Gene Expression Omnibus (GEO) database and analysed to identify differentially expressed genes(DEGs). Enrichment analyses of the DEGs were conducted using the DAVID database. Immune cell infiltration characteristics were analysed by CIBERSORT. Upstream miRNAs and lncRNAs of hub DEGs were identified by miRWalk 3.0 and miRNet in combination with the MNDR platform. As a result, seventy-seven overlapping DEGs and 15 hub genes were identified including CCL5, CXCL9, MYH2, CXCR4, CD74, CCL4, HLA-DRB1, ACTA1, CD69, CXCL10, HLA-DRB5, HLA-DQB1, CXCL13, STAT1 and CKM. The enrichment analyses revealed that the DEGs were mainly enriched in immune response and chemokine signalling pathways. Investigation of immune cell infiltration in liver samples suggested significantly different infiltration between responders and non-responders, mainly characterized by higher proportions of CD8+ T cells and activated NK cells in non-responders. The prediction of upstream miRNAs and lncRNAs led to the identification of a potential mRNA-miRNA-lncRNA regulatory network composed of 2 lncRNAs (H19 and GAS5) and 5 miRNAs (hsa-mir-106b-5p, hsa-mir-17-5p, hsa-mir-20a-5p, hsa-mir-6720-5p and hsa-mir-93-5p) targeting CCL5 mRNA. In conclusion, our study suggested that host genetic factors could affect therapeutic responsiveness in CHB patients. The antiviral process might be associated with the chemokine-mediated immune response and immune cell infiltration in the liver microenvironment.     

Identification of a five-gene signature with prognostic value in colorectal cancer

Colorectal cancer (CRC) ranks as one of the most commonly diagnosed malignancies worldwide. Although mortality rates have been decreasing, the prognosis of CRC patients is still highly dependent on the individual. Therefore, identifying and understanding novel biomarkers for CRC prognosis remains crucial. The gene expression profiles of five-gene expression omnibus (GEO) data sets of CRC were first downloaded. A total of 352 consistent differentially expressed genes (DEGs) were identified for CRC and paired with normal tissues. Functional analysis including gene ontology and Kyoto encyclopedia of genes and genomes pathway enrichment revealed that these DEGs were related to metabolic pathways, tight junctions, and the cell cycle. Ten hub DEGs were identified based on the search tool for the retrieval of interacting genes database and protein–protein interaction networks. By using univariate Cox proportional hazard regression analysis, we found 11 survival-related genes among these DEGs. We finally established a five-gene signature (kinesin family member 15, N-acetyltransferase 2, glutathione peroxidase 3, secretogranin II, and chloride channel accessory 1) with prognostic value in CRC by step multivariate Cox regression analysis. Based on this risk scoring system, patients in the high-risk group had significantly poorer survival results compared with those in the low-risk group (log-rank test, p < 0.0001). Finally, we validated our gene signature scoring system in two independent GEO cohorts (GSE17536 and GSE33113). We found all five of the signature genes to be DEGs in The Cancer Genome Atlas database. In conclusion, our findings suggest that our five DEG-based signature can provide a novel biomarker with useful applications in CRC prognosis.     

Restoration of miR-29b exerts anti-cancer effects on glioblastoma

Background

Methylation plays an important role in the etiology and pathogenesis of colorectal cancer (CRC). This study aimed to identify aberrantly methylated-differentially expressed genes (DEGs) and pathways in CRC by comprehensive bioinformatics analysis.

Methods

Data of gene expression microarrays (GSE68468, GSE44076) and gene methylation microarrays (GSE29490, GSE17648) were downloaded from GEO database. Aberrantly methylated-DEGs were obtained by GEO2R. Functional and enrichment analyses of selected genes were performed using DAVID database. Protein–protein interaction (PPI) network was constructed by STRING and visualized in Cytoscape. MCODE was used for module analysis of the PPI network.

Results

Totally 411 hypomethylation-high expression genes were identified, which were enriched in biological processes of response to wounding or inflammation, cell proliferation and adhesion. Pathway enrichment showed cytokine–cytokine receptor interaction, p53 signaling and cell cycle. The top 5 hub genes of PPI network were CAD, CCND1, ATM, RB1 and MET. Additionally, 239 hypermethylation-low expression genes were identified, which demonstrated enrichment in biological processes including cell–cell signaling, nerve impulse transmission, etc. Pathway analysis indicated enrichment in calcium signaling, maturity onset diabetes of the young, cell adhesion molecules, etc. The top 5 hub genes of PPI network were EGFR, ACTA1, SST, ESR1 and DNM2. After validation in TCGA database, most hub genes still remained significant.

Conclusion

In summary, our study indicated possible aberrantly methylated-differentially expressed genes and pathways in CRC by bioinformatics analysis, which may provide novel insights for unraveling pathogenesis of CRC. Hub genes including CAD, CCND1, ATM, RB1, MET, EGFR, ACTA1, SST, ESR1 and DNM2 might serve as aberrantly methylation-based biomarkers for precise diagnosis and treatment of CRC in the future.

     

Identification of differentially expressed genes and biological characteristics of colorectal cancer by integrated bioinformatics analysis

Colorectal cancer (CRC) ranks as one of the most common malignant tumors worldwide. Its mortality rate has remained high in recent years. Therefore, the aim of this study was to identify significant differentially expressed genes (DEGs) involved in its pathogenesis, which may be used as novel biomarkers or potential therapeutic targets for CRC. The gene expression profiles of GSE21510, GSE32323, GSE89076, and GSE113513 were downloaded from the Gene Expression Omnibus (GEO) database. After screening DEGs in each GEO data set, we further used the robust rank aggregation method to identify 494 significant DEGs including 212 upregulated and 282 downregulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed by DAVID and the KOBAS online database, respectively. These DEGs were shown to be significantly enriched in different cancer-related functions and pathways. Then, the STRING database was used to construct the protein–protein interaction network. The module analysis was performed by the MCODE plug-in of Cytoscape based on the whole network. We finally filtered out seven hub genes by the cytoHubba plug-in, including PPBP, CCL28, CXCL12, INSL5, CXCL3, CXCL10, and CXCL11. The expression validation and survival analysis of these hub genes were analyzed based on The Cancer Genome Atlas database. In conclusion, the robust DEGs associated with the carcinogenesis of CRC were screened through the GEO database, and integrated bioinformatics analysis was conducted. Our study provides reliable molecular biomarkers for screening and diagnosis, prognosis as well as novel therapeutic targets for CRC.     

Inheritable changes in miRNAs expression in HeLa cells after X-ray and mitomycin C treatment

We identified 40 miRNAs with inherited aberrant expression by multiple parallel sequencing of human HeLa cells irradiated with X rays and mitomycin C. Twenty-two miRNAs were repressed and 15 miRNAs were induced after radiation and mytomycin C treatment. The expression of three miRNAs (miR-10b-5p, miR-148a-3p, and miR-340-5p) decreased after X-ray exposure and increased after mitomycin C treatment. The spectrum of aberrantly expressed miRNAs after X-ray and mitomycin C treatment is different, except for three miRNAs (mir-100-5p, miR-99b-5p, miR-501-3p), which showed the inherited decreased expression after both mutagens. It has been ascertained that for five miRNAs (miR-21-3p, miR-182-5p, miR-19b-3p, miR-30a-3p, and miR-30e-3p) with increased inherited expression, the targets are well-described tumor suppressor genes. For 9 miRNAs (miR-99b-5p, miR-148a-3p, miR-365a-3p, miR-193a-3p, miR-100-5p, miR-99a-5p, miR-29b-3p, miR-340-5p, and miR-23b-3p) with reduced inherited expression, the targets are oncogenes. The obtained results provide further support of the idea that induced epigenetic changes in the genome should be considered when assessing the long-term genetic effects of ionizing radiation and chemical compounds.     

Beta-adrenergic pathway activation enhances aggressiveness and inhibits stemness in head and neck cancer

Chronic stress leads to the activation of the beta-adrenergic pathway. Its activation has been implicated in the progression of different types of cancer but its role on head and neck squamous cell carcinomas (HNSCCs) remains undefined. The aim of this study was to investigate the influence of the beta-adrenergic pathway activation in the progression of HNSCCs and offer a panel of potential treatments for patients with the active beta-adrenergic pathway. Five hundred and twenty TCGA patients with primary HNSCCs were divided in two groups: ADRB2^low / SLC6A2^low and ADRB2^high / SLC6A2^high. Differentially expressed genes (DEGs) were identified through differential expression analysis. The association of clinicopathological and genomic features between the groups was analyzed using a bioinformatic approach. Potential drugs for treatment of HNSCC were identified based on the DEGs. There was association between ADRB2 and SLC6A2 expressions with age, race, tumor site, histologic grade, perineural invasion, and HPV p16 status. It was identified 898 DEGs between the groups. High ADRB2/SLC6A2 expression stimulated HNSCC proliferation, adhesion, invasion, and angiogenesis. On the other hand, genes related to cell stemness were downregulated in patients with activation of the beta- adrenergic pathway. Finally, 56 FDA-approved antineoplastic and immunotherapeutic drugs were identified as potential targets for the personalized treatment.     

Differential Expression of microRNAs and their Regulatory Networks in Skin Tissue of Liaoning Cashmere Goat during Hair Follicle Cycles

MicroRNAs (miRNAs) are endogenous small noncoding RNA molecules that negatively regulate gene expression. Herein, we investigated a selective number of miRNAs for their expression in skin tissue of Liaoning Cashmere goat during hair follicle cycles, and their intracellular regulatory networks were constructed based on bioinformatics analysis. The relative expression of six miRNAs (mir-103-3p, -15b-5p, 17-5p, -200b, -25-3p, and -30c-5p) at anagen phase is significantly higher than that at catagen and/or telogen phases. In comparison to anagen, the relative expression of seven miRNAs (mir-148a-3p, -199a-3p, -199a-5p, -24-3p, -30a-5p, -30e-5p, and -29a-3p) was revealed to be significantly up-regulated at catagen and/or telogen stages. The network analyses of miRNAs indicated those miRNAs investigated might be directly or indirectly involved in several signaling pathways through their target genes. These results provided a foundation for further insight into the roles of these miRNAs in skin tissue of Liaoning Cashmere goat during hair follicle cycles.     

Discovery of core gene families associated with liver metastasis in colorectal cancer and regulatory roles in tumor cell immune infiltration

In this study, we aimed to uncover genes that drive the pathogenesis of liver metastasis in colorectal cancer (CRC), and identify effective genes that could serve as potential therapeutic targets for treating with colorectal liver metastasis patients based on two GEO datasets. Several bioinformatics approaches were implemented. First, differential expression analysis screened out key differentially expressed genes (DEGs) across the two GEO datasets. Based on gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, we identified the enrichment functions and pathways of the DEGs that were associated with liver metastasis in CRC. Second, immune infiltration analysis identified key immune signature gene sets associated with CRC liver metastasis, among which two key immune gene families (CD and CCL) identified as key DEGs were filtered by protein-protein interaction (PPI) network. Some of the members in these gene families were associated with disease free survival (DFS) or overall survival (OS) in two subtypes of CRC, namely COAD and READ. Finally, functional enrichment analysis of the two gene families and their neighboring genes revealed that they were closely associated with cytokine, leukocyte proliferation and chemotaxis. These results are valuable in comprehending the pathogenesis of liver metastasis in CRC, and are of seminal importance in understanding the role of immune tumor infiltration in CRC. Our study also identified potentially effective therapeutic targets for liver metastasis in CRC including CCL20, CCL24 and CD70.     

miRNA expression profiling uncovers a role of miR-302b-3p in regulating skin fibroblasts senescence

Numbers of emerging evidence suggest that variable microRNA (miRNA) expression facilitates the aging process. In this study, we distinguished aberrant miRNA expression in aged skin and explored the biological functions and potential mechanism of upregulated miR-302b-3p. At first, miRNA microarray analysis was examined to explore miRNA expression profiling in the skin of aging mice model by D -galactose (d -gal) injection. We identified 29 aberrant miRNAs in aged mice skin. Next, KEGG enrichment analysis was conducted with DIANA-miPath v3.0, which was revealed that enrichment pathways involved in such processes as extracellular matrix-receptor interaction, MAPK signaling pathway, and mammalian target of rapamycin (mTOR) signaling pathway. The target genes of deregulated miRNAs were predicted from four bioinformatic algorithms (miRDB, Targetscan, miRwalk, and Tarbase). The interaction network of miRNAs and their targets were visualized using Cytoscape software. As a result, we found that some hub genes (including JNK2, AKT1/2/3, PAK7, TRPS1, BCL2L11, and IKZF2) were targeted by 12 potential miRNAs (including miR-302b-3p, miR-291a-5p, miR-139-3p, miR-467c-3p, miR-186-3p, etc.). Subsequently, we identified five upregulated miRNA via quantitative polymerase chain reaction and all of them were confirmed increased significantly in aged skin tissues compared with young control tissues. Among them, high expression of miR-302b-3p was verified in both aged skin tissues and senescence fibroblasts. Furthermore, miR-302b-3p mimic accelerated skin fibroblast senescence and suppressed the longevity-associated gene Sirtuin 1(Sirt1) expression, whereas miR-302b-3p inhibitor could delay skin fibroblast senescence and contribute Sirt1 expression. In addition, we demonstrated that c-Jun N-terminal kinase 2(JNK2) is a direct target of miR-302b-3p by a luciferase reporter assay. An inverse correlation was verified in fibroblasts between miR-302b-3p and JNK2. Most importantly, siRNA JNK2 confirmed that low expression of JNK2 could accelerate fibroblasts senescence. In conclusion, our results indicated that overexpressed miR-302b-3p plays an important biological role in accelerating skin aging process via directly targeting JNK2 gene.