An in silico approach to identify and prioritize miRNAs target sites polymorphisms in colorectal cancer and obesity

Genetic variants and haplotype structures of miRNA host genes in cancer and obesity

Cancer and obesity are two important public health problems. This study aimed to investigate the role of genetic variants and haplotypes of miRNA host genes in cancer and obesity. Data from the catalog of genome-wide association studies (GWAS) were used to find significant variants (index). Then, 1000-genome phase 3 data were used to find haplotypic variants (proxy) associated with these diseases. The candidate variants and haplotypes were identified from proxy and index variants. Finally, SNP function analysis was performed. All GWAS-significant cancer-associated miRNA host gene variants, including MIR4713HG, MIR663AHG, MIR99AHG and MIR4435-2HG, were also significantly associated with obesity. The rs703764 variant was common between cutaneous melanoma and obesity traits in the European population (P ≤ 5E-8). The rs2414098 variant was associated with endometrial cancer (P ≤ 5E-13), and the rs7173595 variant was associated with waist-hip ratio (P ≤ 5E-13) and new CGGCATCA haplotypic located at MIR4713HG was identified in the European population. In addition, the ATCTTGTT haplotype for rs17041868 in MIR4435-2HG was identified to be associated with obesity traits (waist-hip ratio and BMI) in the European population (P ≤ 5E-8). This study found that rs703764 is a common genetic marker between cancer and obesity. The CGGCATCA haplotype is common between endometrial cancer and waist-hip ratio. Also, ATCTTGTT haplotype is associated with obesity traits. These results indicate that the variants and haplotypes of miRNAs host genes play an important role between cancer and obesity in the European population. It is suggested to investigate the effect of these structures in other populations.Communicated by Ramaswamy H. Sarma.

In silico prioritisation of microRNA-associated common variants in multiple sclerosis

BACKGROUND

Genome-wide association studies (GWAS) have highlighted over 200 autosomal variants associated with multiple sclerosis (MS). However, variants in non-coding regions such as those encoding microRNAs have not been explored thoroughly, despite strong evidence of microRNA dysregulation in MS patients and model organisms. This study explores the effect of microRNA-associated variants in MS, through the largest publicly available GWAS, which involved 47,429 MS cases and 68,374 controls.

METHODS

We identified SNPs within the coordinates of microRNAs, ± 5-kb microRNA flanking regions and predicted 3'UTR target-binding sites using miRBase v22, TargetScan 7.0 RNA22 v2.0 and dbSNP v151. We established the subset of microRNA-associated SNPs which were tested in the summary statistics of the largest MS GWAS by intersecting these datasets. Next, we prioritised those microRNA-associated SNPs which are among known MS susceptibility SNPs, are in strong linkage disequilibrium with the former or meet a microRNA-specific Bonferroni-corrected threshold. Finally, we predicted the effects of those prioritised SNPs on their microRNAs and 3'UTR target-binding sites using TargetScan v7.0, miRVaS and ADmiRE.

RESULTS

We have identified 30 candidate microRNA-associated variants which meet at least one of our prioritisation criteria. Among these, we highlighted one microRNA variant rs1414273 (MIR548AC) and four 3'UTR microRNA-binding site variants within SLC2A4RG (rs6742), CD27 (rs1059501), MMEL1 (rs881640) and BCL2L13 (rs2587100). We determined changes to the predicted microRNA stability and binding site recognition of these microRNA and target sites.

CONCLUSIONS

We have systematically examined the functional, structural and regulatory effects of candidate MS variants among microRNAs and 3'UTR targets. This analysis allowed us to identify candidate microRNA-associated MS SNPs and highlights the value of prioritising non-coding RNA variation in GWAS. These candidate SNPs could influence microRNA regulation in MS patients. Our study is the first thorough investigation of both microRNA and 3'UTR target-binding site variation in multiple sclerosis using GWAS summary statistics.

Integrative analysis of microRNA-320a-related genes in osteoarthritis cartilage

Objectives

To investigate microRNA-320a-related differentially expressed genes (DEGs) and pathways in osteoarthritis (OA) by bioinformatic analysis.

Methods

The target genes of microRNA-320a were searched and collected from MiRTarBase microRNA Targets dataset, the TargetScan Predicted Nonconserved microRNA Targets dataset and the TargetScan Predicted Conserved microRNA Targets dataset. OA-related microRNAs and OA-related target genes were collected from GeneCards databases. The pathway enrichment analysis of miRNAs ware performed by Funrich analysis tool. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis was obtained from Database for Annotation, Visualization and Integrated Discovery (DAVID). GeneMANIA and STRING are used for protein-protein interaction (PPI) network analysis. Module analysis was performed by Cytoscape.

Results

A total of 176 OA related miRNAs were searched and collected for enrichment analysis, and microRNA-320a was one of OA related miRNAs. Enrichment pathway and analysis of 1721 miRNA-320a-related target genes from MiRTarBase and TargetScan were performed using the online tools Metascape. And results shown that the biological processes were remarkably enriched in chromatin organization, cellular response to DNA damage stimuli, mRNA metabolic process, protein ubiquitination, and regulation of cell adhesion. And then we analysed miRNA-320a-targeted OA genes via KEGG, GO enrichment and PPI Network. Our results showed that miRNA-320a played a role in OA through FoxO signaling pathway, PI3K-Akt signaling pathway, focal adhesion, MAPK signaling pathway, HIF-1 signaling pathway and cellular senescence. And we speculate that MAPK signaling pathway plays a key role in the effect of miRNA-320a on OA.

Conclusion

This study implied microRNA-320a-related DEGs and dysregulated pathways in OA. The aim is to screen miRNA-320a-related genes and pathways in OA and, eventually, to improve the understanding of underlying mechanisms of miRNA-320a in OA.

The role of inflammatory miRNA-mRNA interactions in PBMCs of colorectal cancer and obesity patients

Introduction

Inflammation is a critical hallmark in obesity and colorectal cancer (CRC). This study aimed to investigate effective microRNA (miRNA)–messenger RNA (mRNA) interactions on inflammatory networks involved in obesity and CRC.

Methods

The literature searches were applied to identify genes expression reported on peripheral blood mononuclear cells (PBMCs) and/or blood of CRC subjects and to find inflammatory miRNA  in blood samples. Furthermore, bioinformatics analysis was utilized to find inflammatory miRNA:mRNA interactions of the genes. Finally, a case‐control study was set to investigate the expression of LAMC1 and GNB3 genes besides miR‐10b, miR‐506‐3p, miR‐150‐5p, and miR‐124‐3p in CRC and control subjects.

Results

The expression of LAMC1 gene in healthy control groups was associated with body mass index (BMI) (p < .05). The level of miR‐10b (p < .001), miR‐506 (p < .001), and miR‐124 (p <. 001) were significantly increased in PBMCs of CRC patients, while they were not associated with BMI. The level of miR‐150 was associated with BMI in healthy subjects (p < .05).

Conclusions

The changes in the level of miR‐506 and miR‐124 in CRC patients may be associated with the regulatory role of these miRNAs on LAMC1 expression. The LAMC1 may be related to BMI, however, more observational studies on other populations are needed.

Association of miRNA targetome variants in LAMC1 and GNB3 genes with colorectal cancer and obesity

Background

Colorectal cancer (CRC) is one of the most common obesity‐associated cancers. Inflammation is also considered the most important factor between obesity and CRC. This study aimed to examine miRNAs binding sites variants on inflammatory genes identified using bioinformatics and systematic approach on clinical samples that were collected from CRC patients and controls.

Methods

The candidate variants related to CRC inflammatory genes were obtained from genome‐wide association studies and their population‐specific haplotypes. The variants were analyzed according to their genomic position on the miRNA targetome. Targetome variants in inflammation‐related genes were selected for genetic association study by TaqMan genotyping assay.

Results

The GG genotype of rs7473 decreased the risk of obesity (p < 0.05). Heterozygous genotype (GA) of rs1547715 decreased the risk of CRC (p < 0.05). In the rs7473/rs1547715 genotype and haplotype, the frequencies of AA/GA and GG/AA lessened in CRC and obesity, respectively (p < 0.05).

Conclusions

The variants of rs7473 and rs1547715 were associated with obesity and CRC, respectively. The above‐mentioned associations could be made based on the interactions of these variants with miRNAs.

miR-136 Suppresses the Aggressive Proliferation of Non-Small Cell Lung Cancer Through Restraining Histone Deacetylase 1 (HDAC1) and Phosphorylation of the Janus Kinase 2/Signal Transducer and Activator of Transcription 3 (Jak2/STAT3) Pathway

microRNA-136 can inhibit the proliferating activity of malignant cells and also participate in chemotherapy resistance of colorectal cancer via modulating HDAC1. This study assessed miR-136’s effect on NSCLC cell proliferation and underlying mechanisms. Tumor tissues and paracancerous tissues from NSCLC patients were collected to measure miR-136 and HDAC1 level. Cells were transfected with miR-136-mimics, miR-136-inhibitors or miR-136 mimics+HDAC1-OE followed by analysis of cell viability and apoptosis by CCK-8 method and flow cytometry, phosphorylation of Jak2/STAT3 by western blot. miR-136 was significantly downregulated in tumor tissues and NSCLC cells, accompanied by upregulated HDAC1. miR-136 overexpression suppressed HDAC1 expression, retarded phosphorylation and activation of Jak2/STAT3 signaling, reduced NSCLC cell viability and enhanced apoptosis. In addition, co-transfection of miR-136-mimics and HDAC1-OE reversed the inhibitory effects of miR-136 on NSCLC cells. In conclusion, miR-136 is reduced and HDAC1 is increased in NSCLC and miR-136 overexpression inhibited NSCLC cell proliferation and increased apoptosis possibly through regulating HDAC1/Jak2/STAT3 signal pathway, indicating that miR-136 might be a novel target for the treatment of NSCLC.

An in silico approach to identify and prioritize miRNAs target sites polymorphisms in colorectal cancer and obesity

Colorectal cancer (CRC) and obesity are linked clinical entities with a series of complex processes being engaged in their development. MicroRNAs (miRNAs) participate in these processes through regulating CRC and obesity‐related genes. This study aimed to develop an in silico approach to systematically identify and prioritize miRNAs target sites polymorphisms in obesity and CRC. Data from genome‐wide association studies (GWASs) were used to retrieve CRC and obesity‐associated variants. The polymorphisms that were resided in experimentally verified or computationally predicted miRNA target sites were retrieved and prioritized using a range of bioinformatics analyses. We found 6284 CRC and 38931 obesity unique variants. For CRC 33 haplotypes variants in 134 interactions were in miRNA targetome, while for obesity we found more than 935 unique interactions. Functionally prioritized SNPs revealed that, SNPs in 153 obesity and 50 CRC unique interactions were have disruptive effects on miRNA:mRNA integration by changing on target RNA secondary structure. Structural accessibility of target sites were decreased in 418 and 103 unique interactions and increased in 516 and 79 interactions, for obesity and CRC, respectively. The miRNA:mRNA hybrid stability was increased in 127 and 17 unique interactions and decreased in 33 and 24 interactions for the effect of obesity and CRC SNPs, respectively. In this study, seven SNPs with 15 interactions and three SNPs with four interactions were prioritized for obesity and CRC, respectively. These SNPs could be used for future studies for finding potential biomarkers for diagnoses, prognosis, or treatment of CRC and obesity.