Functional SNP in the microRNA-367 binding site in the 3'UTR of the calcium channel ryanodine receptor gene 3 (RYR3) affects breast cancer risk and calcification

Genetic polymorphism in untranslated regions of PRKCZ influences mRNA structure, stability and binding sites

BACKGROUND

Variations in untranslated regions (UTR) alter regulatory pathways impacting phenotype, disease onset, and course of disease. Protein kinase C Zeta (PRKCZ), a serine-threonine kinase, is implicated in cardiovascular, neurological and oncological disorders. Due to limited research on PRKCZ, this study aimed to investigate the impact of UTR genetic variants' on binding sites for transcription factors and miRNA. RNA secondary structure, eQTLs, and variation tolerance analysis were also part of the study.

METHODS

The data related to PRKCZ gene variants was downloaded from the Ensembl genome browser, COSMIC and gnomAD. The RegulomeDB database was used to assess the functional impact of 5' UTR and 3'UTR variants. The analysis of the transcription binding sites (TFBS) was done through the Alibaba tool, and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) was employed to identify pathways associated with PRKCZ. To predict the effect of variants on microRNA binding sites, PolymiRTS was utilized for 3' UTR variants, and the SNPinfo tool was used for 5' UTR variants.

RESULTS

The results obtained indicated that a total of 24 variants present in the 3' UTR and 25 variants present in the 5' UTR were most detrimental. TFBS analysis revealed that 5' UTR variants added YY1, repressor, and Oct1, whereas 3' UTR variants added AP-2alpha, AhR, Da, GR, and USF binding sites. The study predicted TFs that influenced PRKCZ expression. RNA secondary structure analysis showed that eight 5' UTR and six 3' UTR altered the RNA structure by either removal or addition of the stem-loop. The microRNA binding site analysis highlighted that seven 3' UTR and one 5' UTR variant altered the conserved site and also created new binding sites. eQTLs analysis showed that one variant was associated with PRKCZ expression in the lung and thyroid. The variation tolerance analysis revealed that PRKCZ was an intolerant gene.

CONCLUSION

This study laid the groundwork for future studies aimed at targeting PRKCZ as a therapeutic target.

Functional Analysis of BRCA1 3'UTR Variants Predisposing to Breast Cancer

Purpose

Breast Cancer (BC) is the main female cancer diagnosed worldwide, and it has been described that few genes, such as BRCA1, have a high penetrance for this type of cancer. In this manuscript, we were interested in evaluating the effect of 3'UTR variants on BRCA1 expression.

Patients and Methods

To accomplish this objective, Whole Exome Sequencing (WES) data of 400 patients with unselected BC was used to filter variants located in the region of interest of BRCA1 gene, finding two of them (c.*36C>G and c.*369_373del). miRGate and miRanda in silico tools were used to predict microRNA (miRNA) interaction.

Results

The two variants (c.*36C>G, c.*369_373del) were predicted to affect miRNA interaction. After cloning of BRCA1 3'UTR into pMIR-Report vector, the construct was transfected into two BC cell lines (MDA-MB-231 and MCF-7), and the variant c.*36C>G evidenced overexpression of reporter gene luciferase, showing that the transcript was not being degraded by the miRNA in MDA-MB-231 cells.

Conclusion

The variant seems to protect against Triple Negative BC probably due to the expression level of miRNA in this particular cell line (MDA-MB-231). This is consistent with the clinical history of the patients who harbor BC Hormone Receptors positive (HR+).

Development and validation of a novel prognostic lncRNA signature based on the APOBEC3 family genes in gastric cancer

Introduction

Gastric Cancer (GC) refers to a prevalent malignant cancer accompanied by a weak prognosis. The APOBEC3 family genes and lncRNAs are linked with cancer progression. Nevertheless, there is still a scarcity of data concerning the prognostic value of APOBEC3-related lncRNAs in GC.

Methods

We extracted the data from GC samples, including transcriptome as well as clinical data, obtained from the TCGA database. Then, we screened for lncRNAs that were correlated with the APOBEC3 family genes and constructed an APOBEC3-related lncRNA prognostic signature (LPS) by utilizing univariate Cox and lasso regression analysis. Furthermore, we validated our constructed signature and evaluated it thoroughly, including analysis of its function, immunity, mutations, and clinical applications via multiple methods, including Metascape, GSEA, and analyses including TIC and TME, immune checkpoints, CNV and SNPs, Kaplan-Meier survival curves, nomogram, decision tree and drug prediction analysis. Finally, we overexpressed LINC01094 to evaluate the impacts on the proliferation as well as migration with regards to KATO-2 cells.

Results

We selected eight lncRNAs for our APOBEC3-related LPS, which is demonstrated as a valuable tool in predicting the individual GC patients' prognosis. Subsequently, we segregated the samples into subgroups of high-as well as low-risk relying on the risk score with regards to APOBEC3-related LPS. By performing functional analysis, we have shown that immune-as well as tumor-related pathways were enriched in high- and low-risk GC patients. Furthermore, immune analysis revealed a robust correlation between the APOBEC3-related LPS and immunity. We found that immune checkpoints were significantly associated with the APOBEC3-related LPS and were greatly exhibited in GC tumor and high-risk samples. Mutational analysis suggested that the mutational rate was greater in low-risk samples. Furthermore, we predicted small molecular drugs displayed greater sensitivity in patients categorized as high-risk. Moreover, the immune response was also better in high-risk patients. Of these drugs, dasatinib was significant in both methods and might be considered a potential novel drug for treating high-risk GC patients. Finally, we found that LINC01094 has the potential to enhance the migration, proliferation as well as inhibit apoptosis of KATO-2 in GC cells. And Dasatinib has an inhibitory effect on the migration as well as proliferation in GC cells.

Conclusion

We created a novel APOBEC3-related LPS in predicting the prognosis with regards to individual GC patients. Importantly, this APOBEC3-related LPS was closely associated with immunity and might guide clinical treatment.

The Impact of microRNA SNPS on Breast Cancer: Potential Biomarkers for Disease Detection

Breast cancer is considered a significant health concern worldwide, with genetic predisposition playing a critical role in its etiology. Single nucleotide polymorphisms (SNPs), particularly those within the 3' untranslated regions (3'UTRs) of target genes, are emerging as key factors in breast cancer susceptibility. Specifically, miRNAs have been recognized as possible novel approach for biomarkers discovery for both prognosis and diagnosis due to their direct association with cancer progression. Regional disparities in breast cancer incidence underscore the need for precise interventions, considering socio-cultural and economic factors. This review explores into the differential effects of SNP-miRNA interactions on breast cancer risk, emphasizing both risk-enhancing and protective associations across diverse populations. Furthermore, it explores the clinical implications of these findings, highlighting the potential of personalized approaches in breast cancer management. Additionally, it reviews the evolving therapeutic prospect of microRNAs (miRNAs), extending beyond cancer therapeutics to encompass various diseases, indicative of their versatility as therapeutic agents.

Investigating the prognostic value of mTORC1 signaling in bladder cancer via bioinformatics evaluation

Bladder cancer, a prevalent and heterogeneous malignancy, necessitates the discovery of pertinent biomarkers to enable personalized treatment. The mammalian target of rapamycin complex 1 (mTORC1), a pivotal regulator of cellular growth, metabolism, and immune response, exhibits activation in a subset of bladder cancer tumors. In this study, we explore the prognostic significance of mTORC1 signaling in bladder cancer through the utilization of bioinformatics analysis. Our investigation incorporates transcriptomic, somatic mutation, and clinical data, examining the mTORC1 score of each sample, as well as the enrichment of differentially expressed genes (DEGs), differentiation characteristics, immunological infiltration, and metabolic activity. Our findings reveal that elevated mTORC1 levels serve as an adverse prognostic indicator for bladder cancer patients, exhibiting a significant association with Basal-type bladder cancer. Patients with heightened mTORC1 activation display heightened levels of pro-carcinogenic metabolism. Additionally, these individuals demonstrate enhanced response to immunotherapy. Finally, we develop an mTORC1-related signature capable of predicting the prognosis of bladder cancer patients.The signature offers novel mTORC1-related biomarkers and provides fresh insights into the involvement of mTORC1 in the pathogenesis of bladder cancer.

Identification of signature of tumor-infiltrating CD8 T lymphocytes in prognosis and immunotherapy of colon cancer by machine learning

BACKGROUND

To explore the specific marker of CD8+ T cell subsets which are closely related to the prognosis and immunotherapy of patients with colon cancer.

METHODS

18 kinds of immune cell expression profile data sets were obtained from GEO database. Compared with other immune cell types, the specific markers of CD8 (+) T cells (TI-CD8) in colorectal cancer were screened. Regression analyses were used to further screen prognostic related genes and construct a prognostic evaluation model. The patients were stratified and analyzed according to the risk scores, KRAS mutation status, stage, lymphatic infiltration and other indicators. The landscape of infiltration level, mutation and copy number variation of immune subsets in high and low TI-CD8Sig score groups were compared and analyzed. The difference of drug response between high and low TI-CD8Sig score groups was analyzed. Differential expression of the model genes was verified by the HPA database.

RESULTS

Six prognostic-related CD8T cell-specific gene targets were further screened, and the prognostic evaluation model was constructed. The AUC value of the model is >0.75. FAT3 and UNC13C showed a high mutation state in the low-risk group, while USH2A, MUC5B et al. specifically showed a high mutation state in the high-risk group. Compared with the low-risk group, the high-risk group had lower effective rate of drug response. The expression of PD-1 gene was positively correlated with the level of TI-CD8Sig score.

CONCLUSION

The risk assessment model based on CD8T cell-specific marker genes can effectively predict the prognosis and the drug response of patients with CRC.

Function and regulation of miR-186-5p, miR-125b-5p and miR-1260a in chordoma

BACKGROUND

The function and regulation of miRNAs in progression of chordoma were unclear.

METHODS

Five miRNAs were identified by the machine learning method from the miRNA expression array. CCk-8 assay, EDU assay, wound healing migration assay, and trans-well assay were used to reveal the effect of the miRNAs in chordoma cell lines. Moreover, bioinformation analysis and the mRNA expression array between the primary chordomas and recurrent chordomas were used to find the target protein genes of miRNAs. Furthermore, qRT-PCR and luciferase reporter assay were used to verify the result.

RESULTS

miR-186-5p, miR-30c-5p, miR-151b, and miR-125b-5p could inhibit proliferation, migration, and invasion of chordoma while miR-1260a enhances proliferation, migration, and invasion of chordoma. Recurrent chordoma has a worse disease-free outcome than the primary chordoma patients. AMOT, NPTX1, RYR3, and P2RX5 were the target protein mRNAs of miR-186-5p; NPTX1 was the target protein mRNAs of miR-125b-5p; and AMOT and TNFSF14 were the target protein mRNAs of miR-1260a.

CONCLUSIONS

miR-186-5p, miR-125b-5p, miR-1260a, and their target protein mRNAs including AMOT, NPTX1, RYR3, P2RX5, TNFSF14 may be the basement of chordoma research.

Calcium signals and potential therapy targets in ovarian cancer (Review)

Ovarian cancer (OC) is a deadly disease. The poor prognosis and high lethality of OC are attributed to its high degrees of aggressiveness, resistance to chemotherapy and recurrence rates. Calcium ion (Ca2+) signaling has received attention in recent years, as it appears to form an essential part of various aspects of cancer pathophysiology and is a potential therapeutic target for OC treatment. Disruption of normal Ca2+ signaling pathways can induce changes in cell cycle progression, apoptosis, proliferation and migration and invasion, leading to the development of the malignant phenotype of tumors. In the present review, the main roles of ion channel/receptor/pump‑triggered Ca2+ signaling pathways located at the plasma membrane and organelle Ca2+ transport in OC are summarized. In addition, the potential of Ca2+ signaling as a novel target for the development of effective treatment strategies for OC was discussed. Furthering the understanding into the role of Ca2+ signaling in OC is expected to facilitated the identification of novel therapeutic targets and improved clinical outcomes for patients.

Identification of Somatic Mutations in Plasma Cell-Free DNA from Patients with Metastatic Oral Squamous Cell Carcinoma

The accurate diagnosis and treatment of oral squamous cell carcinoma (OSCC) requires an understanding of its genomic alterations. Liquid biopsies, especially cell-free DNA (cfDNA) analysis, are a minimally invasive technique used for genomic profiling. We conducted comprehensive whole-exome sequencing (WES) of 50 paired OSCC cell-free plasma with whole blood samples using multiple mutation calling pipelines and filtering criteria. Integrative Genomics Viewer (IGV) was used to validate somatic mutations. Mutation burden and mutant genes were correlated to clinico-pathological parameters. The plasma mutation burden of cfDNA was significantly associated with clinical staging and distant metastasis status. The genes TTN, PLEC, SYNE1, and USH2A were most frequently mutated in OSCC, and known driver genes, including KMT2D, LRP1B, TRRAP, and FLNA, were also significantly and frequently mutated. Additionally, the novel mutated genes CCDC168, HMCN2, STARD9, and CRAMP1 were significantly and frequently present in patients with OSCC. The mutated genes most frequently found in patients with metastatic OSCC were RORC, SLC49A3, and NUMBL. Further analysis revealed that branched-chain amino acid (BCAA) catabolism, extracellular matrix-receptor interaction, and the hypoxia-related pathway were associated with OSCC prognosis. Choline metabolism in cancer, O-glycan biosynthesis, and protein processing in the endoplasmic reticulum pathway were associated with distant metastatic status. About 20% of tumors carried at least one aberrant event in BCAA catabolism signaling that could possibly be targeted by an approved therapeutic agent. We identified molecular-level OSCC that were correlated with etiology and prognosis while defining the landscape of major altered events of the OSCC plasma genome. These findings will be useful in the design of clinical trials for targeted therapies and the stratification of patients with OSCC according to therapeutic efficacy.

Clinical and genetic factors associated with radiographic damage in patients with ankylosing spondylitis

OBJECTIVES

To identify clinical and genetic factors associated with severe radiographic damage in patients with ankylosing spondylitis (AS).

METHODS

We newly generated genome-wide single nucleotide polymorphism data (833K) for 444 patients with AS. The severity of radiographic damage was assessed using the modified Stoke Ankylosing Spondylitis Spinal Score (mSASSS). To identify clinical and genetic factors associated with severe radiographic damage, multiple linear regression analyses were performed. Human AS-osteoprogenitor and control-osteoprogenitor cells were used for functional validation.

RESULTS

The significant clinical factors of final mSASSS were baseline mSASSS (β=0.796, p=3.22×10-75), peripheral joint arthritis (β=-0.246, p=6.85×10-6), uveitis (β=0.157, p=1.95×10-3), and smoking (β=0.130, p=2.72×10-2) after adjusting for sex, age and disease duration. After adjusting significant clinical factors, the Ryanodine receptor 3 (RYR3) gene was associated with severe radiographic damage (p=1.00×10-6). For pathway analysis, the PI3K-Akt signalling pathway was associated with severe radiographic damage in AS (p=2.21×10-4, false discovery rate=0.040). Treatment with rhodamine B, a ligand of RYR3, dose-dependently induced matrix mineralisation of AS osteoprogenitors. However, the rhodamine B-induced accelerated matrix mineralisation was not definitive in control osteoprogenitors. Knockdown of RYR3 inhibited matrix mineralisation in SaOS2 cell lines.

CONCLUSIONS

This study identified clinical and genetic factors that contributed to better understanding of the pathogenesis and biology associated with radiographic damage in AS.

Establishment and validation of an ultrasound-based nomogram with risk stratification for short disease-free survival in breast cancer

PURPOSE

This retrospective study aimed to develop and validate an Ultrasound (US)-based nomogram to predict short disease-free survival (short-DFS, less than 120 months DFS) in breast cancer (BC).

METHODS

Nomogram was established based on a training data of 311 BC patients by multivariable logistic regression, and were assessed by discrimination, calibration, and clinical usefulness. Risk stratification was performed by X-tile. An independent testing data of 200 patients with BC was used for external validation.

RESULTS

Nine predictors including three US features and six clinical parameters were screened into the nomogram by Lasso (log λ = -3.594) in training data. Better performance was obtained in the training data (C-index: 0.942) and testing data (C-index: 0.914). Calibration analysis indicated optimal agreement between nomogram predictions and actual observations (p = 0.67). Decision curve analysis showed a great clinical benefit (Youden index: 0.634). Three risk levels are low-risk (<184.0), moderate-risk (184.0-345.3) and high-risk (>345.3). Our nomograms had larger area under the receiver operating characteristic (ROC) curves compared with Magee Equation and Nottingham Prognostic models (0.942 vs. 0.824, 0.790).

CONCLUSION

The US-based nomogram and the practical score system facilitate individualized prediction of short-DFS to optimize clinical decisions and improve prognosis in patients with BC.

Recent Developments on the Roles of Calcium Signals and Potential Therapy Targets in Cervical Cancer

Intracellular calcium (Ca²⁺) concentration ([Ca²⁺]i) is implicated in proliferation, invasion, and metastasis in cancerous tissues. A variety of oncologic therapies and some candidate drugs induce their antitumor effects (in part or in whole) through the modulation of [Ca²⁺]i. Cervical cancer is one of most common cancers among women worldwide. Recently, major research advances relating to the Ca²⁺ signals in cervical cancer are emerging. In this review, we comprehensively describe the current progress concerning the roles of Ca²⁺ signals in the occurrence, development, and prognosis of cervical cancer. It will enhance our understanding of the causative mechanism of Ca²⁺ signals in cervical cancer and thus provide new sights for identifying potential therapeutic targets for drug discovery.

Ryanodine receptor 1-mediated Ca2+ signaling and mitochondrial reprogramming modulate uterine serous cancer malignant phenotypes

Background

Uterine serous cancer (USC) is the most common non-endometrioid subtype of uterine cancer, and is also the most aggressive. Most patients will die of progressively chemotherapy-resistant disease, and the development of new therapies that can target USC remains a major unmet clinical need. This study sought to determine the molecular mechanism by which a novel unfavorable prognostic biomarker ryanodine receptor 1 (RYR1) identified in advanced USC confers their malignant phenotypes, and demonstrated the efficacy of targeting RYR1 by repositioned FDA-approved compounds in USC treatment.

Methods

TCGA USC dataset was analyzed to identify top genes that are associated with patient survival or disease stage, and can be targeted by FDA-approved compounds. The top gene RYR1 was selected and the functional role of RYR1 in USC progression was determined by silencing and over-expressing RYR1 in USC cells in vitro and in vivo. The molecular mechanism and signaling networks associated with the functional role of RYR1 in USC progression were determined by reverse phase protein arrays (RPPA), Western blot, and transcriptomic profiling analyses. The efficacy of the repositioned compound dantrolene on USC progression was determined using both in vitro and in vivo models.

Results

High expression level of RYR1 in the tumors is associated with advanced stage of the disease. Inhibition of RYR1 suppressed proliferation, migration and enhanced apoptosis through Ca²⁺-dependent activation of AKT/CREB/PGC-1α and AKT/HK1/2 signaling pathways, which modulate mitochondrial bioenergetics properties, including oxidative phosphorylation, ATP production, mitochondrial membrane potential, ROS production and TCA metabolites, and glycolytic activities in USC cells. Repositioned compound dantrolene suppressed USC progression and survival in mouse models.

Conclusions

These findings provided insight into the mechanism by which RYR1 modulates the malignant phenotypes of USC and could aid in the development of dantrolene as a repurposed therapeutic agent for the treatment of USC to improve patient survival.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02419-w.

ALOX12: A Novel Insight in Bevacizumab Response, Immunotherapy Effect, and Prognosis of Colorectal Cancer

Colorectal cancer is a highly malignant cancer with poor prognosis and mortality rates. As the first biological agent approved for metastatic colorectal cancer (mCRC), bevacizumab was confirmed to exhibit good performance when combined with chemotherapy and immunotherapy. However, the efficacy of both bevacizumab and immunotherapy is highly heterogeneous across CRC patients with different stages. Thus, exploring a novel biomarker to comprehensively assess the prognosis and bevacizumab and immunotherapy response of CRC is of great significance. In our study, weighted gene co-expression network analysis (WGCNA) and the receiver operating characteristic (ROC) curves were employed to identify bevacizumab-related genes. After verification in four public cohorts and our internal cohort, ALOX12 was identified as a key gene related to bevacizumab response. Prognostic analysis and in vitro experiments further demonstrated that ALOX12 was closely associated with the prognosis, tumor proliferation, invasion, and metastasis. Multi-omics data analysis based on mutation and copy number variation (CNV) revealed that RYR3 drove the expression of ALOX12 and the deletion of 17p12 inhibited ALOX12 expression, respectively. Moreover, we interrogated the relationship between ALOX12 and immune cells and checkpoints. The results exhibited that high ALOX12 expression predicted a higher immune infiltration and better immunotherapy response, which was further validated in Tumor Immune Dysfunction and Exclusion (TIDE) and Subclass Mapping (SubMap) methods. Above all, our study provides a stable biomarker for clinical protocol optimization, prognostic assessment, precise treatment, and individualized treatment of CRC.

Neuronal microRNAs safeguard ER Ca2+ homeostasis and attenuate the unfolded protein response upon stress

Ca2+ is a critical mediator of neurotransmitter release, synaptic plasticity, and gene expression, but also excitotoxicity. Ca2+ signaling and homeostasis are coordinated by an intricate network of channels, pumps, and calcium-binding proteins, which must be rapidly regulated at all expression levels. Τhe role of neuronal miRNAs in regulating ryanodine receptors (RyRs) and inositol 1,4,5-triphosphate receptors (IP3Rs) was investigated to understand the underlying mechanisms that modulate ER Ca2+ release. RyRs and IP3Rs are critical in mounting and propagating cytosolic Ca2+ signals by functionally linking the ER Ca2+ content, while excessive ER Ca2+ release via these receptors is central to the pathophysiology of a wide range of neurological diseases. Herein, two brain-restricted microRNAs, miR-124-3p and miR-153-3p, were found to bind to RyR1-3 and IP3R3 3'UTRs, and suppress their expression at both the mRNA and protein level. Ca2+ imaging studies revealed that overexpression of these miRNAs reduced ER Ca2+ release upon RyR/IP3R activation, but had no effect on [Ca2+]i under resting conditions. Interestingly, treatments that cause excessive ER Ca2+ release decreased expression of these miRNAs and increased expression of their target ER Ca2+ channels, indicating interdependence of miRNAs, RyRs, and IP3Rs in Ca2+ homeostasis. Furthermore, by maintaining the ER Ca2+ content, miR-124 and miR-153 reduced cytosolic Ca2+ overload and preserved protein-folding capacity by attenuating PERK signaling. Overall, this study shows that miR-124-3p and miR-153-3p fine-tune ER Ca2+ homeostasis and alleviate ER stress responses.

Acute RyR1 Ca2+ leak enhances NADH-linked mitochondrial respiratory capacity

Sustained ryanodine receptor (RyR) Ca²⁺ leak is associated with pathological conditions such as heart failure or skeletal muscle weakness. We report that a single session of sprint interval training (SIT), but not of moderate intensity continuous training (MICT), triggers RyR1 protein oxidation and nitrosylation leading to calstabin1 dissociation in healthy human muscle and in in vitro SIT models (simulated SIT or S-SIT). This is accompanied by decreased sarcoplasmic reticulum Ca²⁺ content, increased levels of mitochondrial oxidative phosphorylation proteins, supercomplex formation and enhanced NADH-linked mitochondrial respiratory capacity. Mechanistically, (S-)SIT increases mitochondrial Ca²⁺ uptake in mouse myotubes and muscle fibres, and decreases pyruvate dehydrogenase phosphorylation in human muscle and mouse myotubes. Countering Ca²⁺ leak or preventing mitochondrial Ca²⁺ uptake blunts S-SIT-induced adaptations, a result supported by proteomic analyses. Here we show that triggering acute transient Ca²⁺ leak through RyR1 in healthy muscle may contribute to the multiple health promoting benefits of exercise.

Genetic variants of cell cycle pathway genes are associated with head and neck squamous cell carcinoma in the Chinese population

Genetic alterations in the cell cycle pathway are common in head and neck squamous cell carcinoma (HNSCC). We identified four novel HNSCC susceptibility loci (CDKN1C rs452338, CDK4 rs2072052, E2F2 rs3820028 and E2F2 rs2075993) through a two-stage matched case-control study. There was a combined effect among the four single nucleotide polymorphisms (SNPs), as the number of risk genotypes increased, the risk of HNSCC displayed an increasing trend (Ptrend < 0.001). And there were multiplicative interactions between rs452338 and rs2072052, rs2072052 and rs3820028, rs2072052 and rs2075993. Functional bioinformatics analysis and dual-luciferase reporter assay revealed that E2F2 rs2075993 T>C reduced the stability of E2F2 3'-UTR secondary structure and affected the binding of E2F2 to miR-940, which was up-regulated in HNSCC tumor tissues (P = 2.9e-8) and was correlated with poor overall survival of HNSCC (HR = 1.39, 95% CI = 1.02-1.90). In vitro assays, we discovered that the expression of miR-940 was regulated by METTL3, and miR-940 promoted the proliferation, migration and invasion, and inhibited the senescence and autophagy of tumor cells. In terms of mechanism, compared with rs2075993 allele T, we found that the protective variant rs2075993 allele C interfered with the translational inhibition of E2F2 by miR-940, resulting in increased expression of E2F2 protein, which further reduced the proliferation, migration, invasion, and increased the senescence of tumor cells.

Human microRNA similarity in breast cancer

MicroRNAs (miRNAs) play important roles in a variety of human diseases, including breast cancer. A number of miRNAs are up- and down-regulated in breast cancer. However, little is known about miRNA similarity and similarity network in breast cancer. Here, a collection of 272 breast cancer-associated miRNA precursors (pre-miRNAs) were utilized to calculate similarities of sequences, target genes, pathways and functions and construct a combined similarity network. Well-characterized miRNAs and their similarity network were highlighted. Interestingly, miRNA sequence-dependent similarity networks were not identified in spite of sequence–target gene association. Similarity networks with minimum and maximum number of miRNAs originate from pathway and mature sequence, respectively. The breast cancer-associated miRNAs were divided into seven functional classes (classes I–VII) followed by disease enrichment analysis and novel miRNA-based disease similarities were found. The finding would provide insight into miRNA similarity, similarity network and disease heterogeneity in breast cancer.

Ion Channels, Transporters, and Sensors Interact with the Acidic Tumor Microenvironment to Modify Cancer Progression

Solid tumors, including breast carcinomas, are heterogeneous but typically characterized by elevated cellular turnover and metabolism, diffusion limitations based on the complex tumor architecture, and abnormal intra- and extracellular ion compositions particularly as regards acid-base equivalents. Carcinogenesis-related alterations in expression and function of ion channels and transporters, cellular energy levels, and organellar H⁺ sequestration further modify the acid-base composition within tumors and influence cancer cell functions, including cell proliferation, migration, and survival. Cancer cells defend their cytosolic pH and HCO₃^- concentrations better than normal cells when challenged with the marked deviations in extracellular H⁺, HCO₃^-, and lactate concentrations typical of the tumor microenvironment. Ionic gradients determine the driving forces for ion transporters and channels and influence the membrane potential. Cancer and stromal cells also sense abnormal ion concentrations via intra- and extracellular receptors that modify cancer progression and prognosis. With emphasis on breast cancer, the current review first addresses the altered ion composition and the changes in expression and functional activity of ion channels and transporters in solid cancer tissue. It then discusses how ion channels, transporters, and cellular sensors under influence of the acidic tumor microenvironment shape cancer development and progression and affect the potential of cancer therapies.

Identifying Driver Genes for Individual Patients through Inductive Matrix Completion

MOTIVATION

The driver genes play a key role in the evolutionary process of cancer. Effectively identifying these driver genes is crucial to cancer diagnosis and treatment. However, due to the high heterogeneity of cancers, it remains challenging to identify the driver genes for individual patients. Although some computational methods have been proposed to tackle this problem, they seldom consider the fact that the genes functionally similar to the well-established driver genes may likely play similar roles in cancer process, which potentially promotes the driver gene identification. Thus, here we developed a novel approach of IMCDriver to promote the driver gene identification both for cohorts and individual patients.

RESULTS

IMCDriver first considers the well-established driver genes as prior information, and adopts the using multi-omics data (e.g., somatic mutation, gene expression and protein-protein interaction) to compute the similarity between patients/genes. Then, IMCDriver prioritizes the personalized mutated genes according to their functional similarity to the well-established driver genes via Inductive Matrix Completion. Finally, IMCDriver identifies the highly rank-ordered genes as the personalized driver genes. The results on five cancer datasets from TCGA show that our IMCDriver outperforms other existing state-of-the-art methods both in the cohort and patient-specific driver gene identification. IMCDriver also reveals some novel driver genes that potentially drive cancer development. In addition, even for the driver genes rarely mutated among a population, IMCDriver can still identify them and prioritize them with high priorities.

AVAILABILITY

Code available at https://github.com/NWPU-903PR/IMCDriver.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A ceRNA-associated risk model predicts the poor prognosis for head and neck squamous cell carcinoma patients

Head and neck squamous cell carcinoma (HNSCC) is one of the most malignant cancers with poor prognosis worldwide. Emerging evidence indicates that competing endogenous RNAs (ceRNAs) are involved in various diseases, however, the regulatory mechanisms of ceRNAs underlying HNSCC remain unclear. In this study, we retrieved differentially expressed long non-coding RNAs (DElncRNAs), messenger RNAs (DEmRNAs) and microRANs (DEmiRNAs) from The Cancer Genome Atlas database and constructed a ceRNA-based risk model in HNSCC by integrated bioinformatics approaches. Functional enrichment analyses showed that DEmRNAs might be involved in extracellular matrix related biological processes, and protein–protein interaction network further selected out prognostic genes, including MYL1 and ACTN2. Importantly, co-expressed RNAs identified by weighted co-expression gene network analysis constructed the ceRNA networks. Moreover, AC114730.3, AC136375.3, LAT and RYR3 were highly correlated to overall survival of HNSCC by Kaplan–Meier method and univariate Cox regression analysis, which were subsequently implemented multivariate Cox regression analysis to build the risk model. Our study provides a deeper understanding of ceRNAs on the regulatory mechanisms, which will facilitate the expansion of the roles on the ceRNAs in the tumorigenesis, development and treatment of HNSCC.

SNPs within microRNA binding sites and the prognosis of breast cancer

Single nucleotide polymorphisms (SNPs) within microRNA binding sites can affect the binding of microRNA to mRNA and regulate gene expression, thereby contributing to cancer prognosis. Here we performed a two-stage study of 2647 breast cancer patients to explore the association between SNPs within microRNA binding sites and breast cancer prognosis. In stage I, we genotyped 192 SNPs within microRNA binding sites using the Illumina Goldengate platform. In stage II, we validated SNPs associated with breast cancer prognosis in another dataset using the TaqMan platform. We identified 8 SNPs significantly associated with breast cancer prognosis in stage I (P<0.05), and only rs10878441 was statistically significant in stage II (AA vs CC, HR=2.21, 95% CI: 1.11-4.42, P=0.024). We combined the data from stage I and stage II, and found that, compared with rs10878441 AA genotype, CC genotype was associated with poor survival of breast cancer (HR=2.19, 95% CI: 1.30-3.70, P=0.003). Stratified analyses demonstrated that rs10878441 was related to breast cancer prognosis in grade II and lymph node-negative patients (P<0.05). The Leucine-rich repeat kinase 2 (LRRK2) rs10878441 CC genotype is associated with poor prognosis of breast cancer in a Chinese population and may be used as a potential prognostic biomarker for breast cancer.

• The LRRK2 rs10878441 CC genotype is associated with poor prognosis of breast cancer in a Chinese population.

• Stratified analyses demonstrated that rs10878441 was related to breast cancer prognosis in grade II patients and lymph node-negative patients.

miR-875-5p exerts tumor-promoting function via down-regulation of CAPZA1 in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC) is one of the leading causes of cancer deaths worldwide. Currently, efficient genetic markers for diagnosis and treatment of ESCC are lacking. MicroRNAs (miRNAs) are global genetic regulators that control cancer gene expression by binding to the 3'untranslated regions (3'UTRs) of targeting mRNAs. In addition, miRNAs function as oncogenes or tumor suppressors in the progression of tumors. In the current study, we found that hsa-miR-875-5p (miR-875-5p) exhibited amplification in ESCC according to the TCGA database. Then, xCELLigence Real-Time Cell Analyzer (RTCA)-MP system and colony formation assays were employed to detect cell proliferationand colony formationability. The results showed that miR-875-5p promoted the proliferation ESCC cells. Subsequently, transwell results indicated that miR-875-5p promoted the invasion and migration of ESCC cells. Furthermore, we showed that miR-875-5p was able to bind to CAPZA13'UTR, which contains the single nucleotide polymorphism (SNP), rs373245753, as reported in our previous study involving WGS and WES on ESCC. Subsequently, mRNA affinity pull-down assays verifiedthat the SNP disrupts miR-875-5p binding to CAPZA1. The current study is the first demonstration that miR-875-5p may function as an oncogene via down-regulation of CAPZA1 expression in ESCC.

Genetic Variability in the microRNA Binding Sites of BMPR1B, TGFBR1, IQGAP1, KRAS, SETD8 and RYR3 and Risk of Breast Cancer in Colombian Women

Background

Genetic variants in microRNA (miR) binding sites affect the regulation of miR-dependent gene expression and have been linked to the risk of a variety of cancers including breast cancer (BC). Most BC risk variants had been identified in women of European and Asian ancestry, but genetic data for Hispanic women are scarce. Here, we investigate the association between six variants in miR binding sites and BC risk in Colombian women.

Methods

We genotyped miR binding site variants in the BMPR1B, TGFBR1, IQGAP1, KRAS, SETD8 and RYR3 genes in 1022 BC cases and 1023 controls from the Colombian breast cancer case–control (Col-BCCC) study using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). Multiple logistic regression and permutation techniques were applied to assess the association between genetic variants and BC risk.

Results

We found no evidence of association between any of the six miR binding site variants and overall or estrogen receptor subtype-specific BC risk in Colombian women.

Conclusion

Our findings may point to ethnic differences in the association between genetic variability in miR binding sites and breast cancer risk.

Roles of NAD+ and Its Metabolites Regulated Calcium Channels in Cancer

Nicotinamide adenine dinucleotide (NAD⁺) is an essential cofactor for redox enzymes, but also moonlights as a regulator for ion channels, the same as its metabolites. Ca²⁺ homeostasis is dysregulated in cancer cells and affects processes such as tumorigenesis, angiogenesis, autophagy, progression, and metastasis. Herein, we summarize the regulation of the most common calcium channels (TRPM2, TPCs, RyRs, and TRPML1) by NAD⁺ and its metabolites, with a particular focus on their roles in cancers. Although the mechanisms of NAD⁺ metabolites in these pathological processes are yet to be clearly elucidated, these ion channels are emerging as potential candidates of alternative targets for anticancer therapy.

microRNA-related single-nucleotide polymorphisms and breast cancer

Breast cancer, as the most common cancer in women which affects patients both mentally and physically, requires great attention in many areas and many levels as this cancer is known to be multifactorial. Single-stranded molecules called microRNAs with near 22 nucleotides are seen to act in central dogma of molecular biology by inhibiting the translation process; it is demonstrated that any alteration in their sequence especially single-nucleotide polymorphisms (SNPs) may lead into increasing the breast cancer risk. miR-SNPs are considered to be the potential biomarkers for early detection of breast cancer. As a result, this review documents the well-known miR-SNPs that are known to be associated with breast cancer. In this regard, two principals were discussed: (a) SNPs in the target genes of microRNAs and the alteration in gene expression due to this phenomenon; (b) changes based on the SNPs in the microRNA coding region and the impact on their interaction with target messenger RNA.

Lessons from the Endoplasmic Reticulum Ca2+ Transporters-A Cancer Connection

Ca2+ is an integral mediator of intracellular signaling, impacting almost every aspect of cellular life. The Ca2+-conducting transporters located on the endoplasmic reticulum (ER) membrane shoulder the responsibility of constructing the global Ca2+ signaling landscape. These transporters gate the ER Ca2+ release and uptake, sculpt signaling duration and intensity, and compose the Ca2+ signaling rhythm to accommodate a plethora of biological activities. In this review, we explore the mechanisms of activation and functional regulation of ER Ca2+ transporters in the establishment of Ca2+ homeostasis. We also contextualize the aberrant alterations of these transporters in carcinogenesis, presenting Ca2+-based therapeutic interventions as a means to tackle malignancies.

Impact of NR5A2 and RYR2 3'UTR polymorphisms on the risk of breast cancer in a Chinese Han population

OBJECTIVES

The NR5A2 and RYR2 genes are important players in steroid metabolism and play an important role in cancer research. In this research, we want to evaluate the effect of NR5A2 and RYR2 polymorphisms on breast cancer (BC).

METHODS

Four single nucleotide polymorphisms on NR5A2 and RYR2 were selected to genotype by Agena MassARRAY in 379 BC patients and 407 healthy controls. Using the PLINK software to calculate the Odds ratio (OR) and 95% confidence intervals (CIs) via the logistic regression analysis to evaluate the risk for BC.

RESULTS

We found that NR5A2 rs2246209 significantly decreased the risk of BC with the AA genotype (OR 0.58, 95%CI 0.34-0.99, p = 0.049), and recessive model (OR 0.59, 95%CI 0.35-0.99, p = 0.046); rs12594 in the RYR2 gene significantly decreased the risk of BC in the GG genotype (OR 0.44, 95%CI 0.22-0.88, p = 0.020), and recessive model (OR 0.43, 95%CI 0.21-0.85, p = 0.016). Further stratification analysis showed that NR5A2 rs2246209 was related to a lower incidence of BC affected by age, lymph nodes metastasis, and tumor stage; RYR2 rs12594 was related to a decreased BC risk restricted by age, estrogen receptor (ER), progesterone receptor (PR), menopausal status, tumor size, and tumor stage. Rs12594 in the RyR2 gene remained significant on the genetic susceptibility of PR-positive BC after Bonferroni correction (p < 0.0125).

CONCLUSIONS

This study provides an evidence that NR5A2 rs2246209 and RYR2 rs12594 decreased the risk of breast cancer.

From Orai to E-Cadherin: Subversion of Calcium Trafficking in Cancer to Drive Proliferation, Anoikis-Resistance, and Metastasis

The common currency of epithelial differentiation and homeostasis is calcium, stored primarily in the endoplasmic reticulum, rationed according to need, and replenished from the extracellular milieu via store-operated calcium entry (SOCE). This currency is disbursed by the IP3 receptor in response to diverse extracellular signals. The rate of release is governed by regulators of proliferation, autophagy, survival, and programmed cell death, the strength of the signal leading to different outcomes. Intracellular calcium acts chiefly through intermediates such as calmodulin that regulates growth factor receptors such as epidermal growth factor receptor (EGFR), actin polymerization, and adherens junction assembly and maintenance. Here we review this machinery and its role in differentiation, then consider how cancer cells subvert it to license proliferation, resist anoikis, and enable metastasis, either by modulating the level of intracellular calcium or its downstream targets or effectors such as EGFR, E-cadherin, IQGAP1, TMEM16A, CLCA2, and TRPA1. Implications are considered for the roles of E-cadherin and growth factor receptors in circulating tumor cells and metastasis. The discovery of novel, cell type-specific modulators and effectors of calcium signaling offers new possibilities for cancer chemotherapy.

Genetic interactions among Pto-miR319 family members and their targets influence growth and wood properties in Populus tomentosa

MicroRNAs (miRNAs) play crucial roles in all aspects of plant growth and development, but the genetic interactions of miRNAs and their target genes in woody plants are largely unknown. Here, we integrated association genetics and expression profiling to decipher the allelic variations and interactions of the Pto-MIR319 family of miRNAs and 12 putative Pto-miR319 target genes related to wood formation in 435 unrelated individuals of Populus tomentosa Carrière (Chinese white poplar). Expression pattern analysis showed that among all pairings between expressions of pre-miRNA of Pto-MIR319 members and targets, 70.0% showed negative correlation of expression levels (r = - 0.944 to 0.674, P < 0.01) in eight tissues and organs of poplar, suggesting that Pto-miR319 may participate in the regulatory network of wood formation. Single SNP-based association studies identified 137 significant associations (P < 0.01, Q < 0.1), representing 126 unique SNPs from Pto-MIR319 members and their targets, with 10 tree growth traits, revealing that these genetic factors have common roles related to wood formation. Epistasis analysis uncovered 105 significant SNP-SNP associations (P < 0.01) influencing the 10 traits, demonstrating the close genetic interactions between Pto-MIR319 family members and the 12 Pto-miR319 target genes. Notably, one common SNP, in the precursor region of Pto-MIR319e, affected the stability of Pto-MIR319e's secondary structure by altering the stem-loop structure and minimum free energy, contributing to variations in the expression of Pto-MIR319e and Pto-miR319e target genes. This study enriches the understanding of the functions of miR319 family miRNAs in poplar and exemplifies a feasible approach to exploring the genetic effects underlying miRNA-mRNA interactions related to complex traits in trees.

Predictive role of single nucleotide polymorphism (rs11614913) in the development of breast cancer in Pakistani population

Aim: miRNAs play an important role in breast cancer (BC). Variations in miRNAs influence their maturation, expression and consequently regulation of their target genes. Materials & methods: In this study, single nucleotide polymorphism rs11614913 was genotyped in BC patients (n = 300) and 230 controls by employing tetra primer amplification refractory mutation system PCR and Sanger sequencing (Macrogen Korea). Results: A significant difference was observed in the genotypes through co-dominant (χ².#x00A0;= 42.03; p < 0.0001), additive (odds ratio [OR] = 0.6441 [0.4887-0.8490, 95% confidence interval]; p < 0.0019), dominant (OR = 0.3996 [0.2809-0.5686], p < 0.0001) and recessive (OR = 0.2993 [0.1220-0.7347], p < 0.009) statistical models showed decreased risk association of C allele with BC. Conclusion: Females having CT genotype are at higher risk of BC as compared with those having CC genotype.

A STAT3 of Addiction: Adipose Tissue, Adipocytokine Signalling and STAT3 as Mediators of Metabolic Remodelling in the Tumour Microenvironment

Metabolic remodelling of the tumour microenvironment is a major mechanism by which cancer cells survive and resist treatment. The pro-oncogenic inflammatory cascade released by adipose tissue promotes oncogenic transformation, proliferation, angiogenesis, metastasis and evasion of apoptosis. STAT3 has emerged as an important mediator of metabolic remodelling. As a downstream effector of adipocytokines and cytokines, its canonical and non-canonical activities affect mitochondrial functioning and cancer metabolism. In this review, we examine the central role played by the crosstalk between the transcriptional and mitochondrial roles of STAT3 to promote survival and further oncogenesis within the tumour microenvironment with a particular focus on adipose-breast cancer interactions.

The Role of Vitamin D and Sunlight Incidence in Cancer

BACKGROUND

Vitamin D (VD) deficiency affects individuals of different ages in many countries. VD deficiency may be related to several diseases, including cancer.

OBJECTIVE

This study aimed to review the relationship between VD deficiency and cancer.

METHODS

We describe the proteins involved in cancer pathogenesis and how those proteins can be influenced by VD deficiency. We also investigated a relationship between cancer death rate and solar radiation.

RESULTS

We found an increased bladder cancer, breast cancer, colon-rectum cancer, lung cancer, oesophagus cancer, oral cancer, ovary cancer, pancreas cancer, skin cancer and stomach cancer death rate in countries with low sunlight. It was also observed that amyloid precursor protein, ryanodine receptor, mammalian target of rapamycin complex 1, and receptor for advanced glycation end products are associated with a worse prognosis in cancer. While the Klotho protein and VD receptor are associated with a better prognosis in the disease. Nfr2 is associated with both worse and better prognosis in cancer.

CONCLUSION

The literature suggests that VD deficiency might be involved in cancer progression. According to sunlight data, we can conclude that countries with low average sunlight have high cancers death rate. New studies involving transcriptional and genomic data in combination with VD measurement in long-term experiments are required to establish new relationships between VD and cancer.

A robust two-gene signature for glioblastoma survival prediction

Glioblastoma multiforme (GBM) is a highly malignant brain tumor. We explored the prognostic gene signature in 443 GBM samples by systematic bioinformatics analysis, using GSE16011 with microarray expression and corresponding clinical data from Gene Expression Omnibus as the training set. Meanwhile, patients from The Chinese Glioma Genome Atlas database (CGGA) were used as the test set and The Cancer Genome Atlas database (TCGA) as the validation set. Through Cox regression analysis, Kaplan-Meier analysis, t-distributed Stochastic Neighbor Embedding algorithm, clustering, and receiver operating characteristic analysis, a two-gene signature (GRIA2 and RYR3) associated with survival was selected in the GSE16011 dataset. The GRIA2-RYR3 signature divided patients into two risk groups with significantly different survival in the GSE16011 dataset (median: 0.72, 95% confidence interval [CI]: 0.64-0.98, vs median: 0.98, 95% CI: 0.65-1.61 years, logrank test P < .001), the CGGA dataset (median: 0.84, 95% CI: 0.70-1.18, vs median: 1.21, 95% CI: 0.95-2.94 years, logrank test P = .0017), and the TCGA dataset (median: 1.03, 95% CI: 0.86-1.24, vs median: 1.23, 95% CI: 1.04-1.85 years, logrank test P = .0064), validating the predictive value of the signature. And the survival predictive potency of the signature was independent from clinicopathological prognostic features in multivariable Cox analysis. We found that after transfection of U87 cells with small interfering RNA, GRIA2 and RYR3 influenced the biological behaviors of proliferation, migration, and invasion of glioblastoma cells. In conclusion, the two-gene signature was a robust prognostic model to predict GBM survival.

Calcium transport and signalling in breast cancer: Functional and prognostic significance

Comprised of a complex network of numerous intertwining pathways, the Ca²⁺ signalling nexus is an essential mediator of many normal cellular activities. Like many other such functions, the normal physiological activity of Ca²⁺ signalling is frequently co-opted and reshaped in cases of breast cancer, creating a potent oncogenic drive within the affected cell population. Such modifications can occur within pathways mediating either Ca²⁺ import (e.g. TRP channels, ORAI-STIM1) or Ca²⁺ export (e.g. PMCA), indicating that both increases and decreases within cellular Ca²⁺ levels have the potential to increase the malignant potential of a cell. Increased understanding of these pathways may offer clinical benefit in terms of both prognosis and treatment; patient survival has been linked to expression levels of certain Ca²⁺ transport proteins, whilst selective targeting of these factors with novel anti-cancer agents has demonstrated a variety of anti-tumour effects in in vitro studies. In addition, the activity of several Ca²⁺ signalling pathways has been shown to influence chemotherapy response, suggesting that a synergistic approach coupling traditional chemotherapy with Ca²⁺ targeting agents may also improve patient outcome. As such, targeted modulation of these pathways represents a novel approach in precision medicine and breast cancer therapy.

MicroRNA-binding site polymorphisms and risk of colorectal cancer: A systematic review and meta-analysis

Genetic variations in miRNAs binding site might participate in cancer risk. This study aimed to systematically review the association between miRNA-binding site polymorphisms and colorectal cancer (CRC). Electronic literature search was carried out on PubMed, Web of Science (WOS), Scopus, and Embase. All types of observational studies till 30 November 2018 were included. Overall 85 studies (21 SNPs) from two systematic searches were included analysis. The results showed that in the Middle East population, the minor allele of rs731236 was associated with decreased risk of CRC (heterozygote model: 0.76 [0.61-0.95]). The minor allele of rs3025039 was related to increased risk of CRC in East Asian population (allelic model: 1.25 [1.01-1.54]). Results for rs3212986 were significant in overall and subgroup analysis (P < .05). For rs1801157 in subgroup analysis the association was significant in Asian populations (including allelic model: 2.28 [1.11-4.69]). For rs712, subgroup analysis revealed a significant (allelic model: 1.41 [1.23-1.61]) and borderline (allelic model: 0.92 [0.84-1.00]) association in Chinese and Czech populations, respectively. The minor allele of rs17281995 increased risk of CRC in different genetic models (P < .05). Finally, rs5275, rs4648298, and rs61764370 did not show significant associations. In conclusion, minor allele of rs3025039, rs3212986, and rs712 polymorphisms increases the risk of CRC in the East Asian population, and heterozygote model of rs731236 polymorphism shows protective effect in the Middle East population. In Europeans, the minor allele of rs17281995 may increase the risk of CRC, while rs712 may have a protective effect. Further analysis based on population stratifications should be considered in future studies.

MiR-124-3p inhibits the migration and invasion of Gastric cancer by targeting ITGB3

BACKGROUND

Gastric cancer is one of the major malignant tumors in the world. Integrins expressed in cancer cells can promote tumor progression and migration. MiRNAs can inhibit the expression of target genes by directly binding to their mRNAs and can affect various important biological processes. The aim of this study was to investigate the role of miR-124- 3p and ITGB3 in gastric cancer.

METHODS

RT-PCR and western blot are used to detect the expression of miR-124-3p, ITGB3 and integrin β3 in gastric cancer tissues and cells. The wound healing, CCK-8 assay, transwell migration and invasion assay were performed to determine the cell proliferation, migration and invasion. What's more, bioinformatics prediction and luciferase assay was conducted to demonstrated the binding efficiency between miR-124-3p and ITGB3.

RESULTS

We verified that ITGB3 and miR-124-3p changes the migration and invasion of gastric cancer cells in vitro. The overexpression or silencing of miR-124-3p inhibited or promoted the proliferation, migration and invasion of both selected gastric cancer cells, and ITGB3 is just the reverse. Meanwhile, we validated that ITGB3 is the target of miR-124-3p by bioinformatics prediction and luciferase assay. Lastly, the expression of ITGB3 in 40 pairs of gastric cancer tissues were significantly higher than that in the adjacent normal tissues, while the expression level of miR-124-3p was significantly decreased in cancer tissues.

CONCLUSIONS

miR-124-3p inhibits the migration and invasion of Gastric cancer by targeting ITGB3 in gastric cancer cells. Our results suggested that miR-124-3p and ITGB3 may reasonably serve as a promising therapeutic target.

Therapeutic Potential of the miRNA-ATM Axis in the Management of Tumor Radioresistance

The ataxia-telangiectasia mutated (ATM) protein kinase is widely known for its function as a chief mobilizer of the DNA damage response (DDR) upon DNA double-strand breaks. ATM orchestrates the DDR by modulating the expression of various miRNAs through several mechanisms. On the other hand, a set of miRNAs contribute to tight regulation of ATM by directly targeting the 3'-untranslated region of ATM mRNA. This review addresses the therapeutic application and molecular mechanisms that underlie the intricate interactions between miRNAs and ATM. It also describes therapeutic delivery of miRNAs in different environments such as hypoxic tumor microenvironments.

GOLGA7 rs11337, a Polymorphism at the MicroRNA Binding Site, Is Associated with Glioma Prognosis

MicroRNAs bind to the 3' untranslated regions of mRNAs, affecting translation, tumorigenesis, and apoptosis. This study evaluated the role of TYMS (rs1059394, C > T, and rs2847153, G > A), RYR3 (rs1044129, G > A), KIAA0423 (rs1053667, T > C), and GOLGA7 (rs11337, G > T) polymorphisms for assessment of glioma risk and prognosis among the Chinese Han population. Five single-nucleotide polymorphisms were assessed in 605 glioma patients and 1,300 controls. We found a significant correlation between rs1059394 and glioma susceptibility in the homozygote and dominant genetic models (TT versus CC, odds ratio [OR] = 0.71, 95% confidence interval [CI] = 0.52-0.97, p = 0.03; CT+TT versus CC, OR = 0.74, 95% CI = 0.55-0.99, p = 0.04). The results of the Kaplan-Meier and log rank tests revealed that the rs11337 GG genotype correlated with better overall survival of glioma patients (p = 0.017) than the GT genotype. Multivariate Cox regression analysis results also showed that the rs11337 GT genotype correlated with worse overall survival (p = 0.017, hazard ratio [HR] = 1.25, 95% CI = 1.04-1.5) than the GG genotype. These results suggest that GOLGA7 (rs11337) polymorphism may play a role in the prognosis of glioma patients and that TYMS (rs1059394) is associated with glioma risk.

lnc-3215 Suppression Leads to Calcium Overload in Selenium Deficiency-Induced Chicken Heart Lesion via the lnc-3215-miR-1594-TNN2 Pathway

Selenium deficiency has been proven to induce calcium disorders in the chicken heart. However, detailed regulatory mechanisms, e.g., the long noncoding RNA (lncRNA)-microRNA (miRNA)-mRNA regulatory axis, have not yet been described. Here, we point out lnc-2315, miR-1594, and Troponin T (TNNT2) based on the results of lncRNA and miRNA comparative genomics group analysis of Se-deficient chicken hearts compared with control hearts. We employed lnc-3215 and TNNT2 knockdown, miR-1594 knockdown, and overexpression models in the chicken embryos in vivo, and lnc-3215, miR-1594, and TNNT2 knockdown and overexpression models in cardiomyocytes in vitro. The dual-luciferase reporter assay and quantitative real-time PCR were used to confirm the relationships between miR-1594 and TNNT2, lnc-3215, and miR-1594 in cardiomyocytes. Our results revealed that TNNT2 suppression induced cardiac calcium overload in vivo and in vitro. miR-1594 activates cardiac calcium overload by targeting TNNT2. Moreover, we found that lnc-3215 regulates miR-1594, and thus influences the TNNT2 expression in vivo and in vitro; these conclusions were verified by gene knockdown in chicken embryos. Our present study revealed a novel regulatory model of a calcium program, which comprises lnc-3215, miR-1594, and TNNT2 in the chicken heart. Our conclusions may provide a feasible diagnostic tool for Se-deficient cardiomyocytes injury.

Combined MicroRNA In Situ Hybridization and Immunohistochemical Detection of Protein Markers

MicroRNAs are short (18-23 nucleotides) noncoding RNAs involved in posttranscriptional regulation of gene expression through their specific binding to the 3'UTR of mRNAs. MicroRNAs can be detected in tissues using specific locked nucleic acid (LNA)-enhanced probes. The characterization of microRNA expression in tissues by in situ detection is often crucial following a microRNA biomarker discovery phase in order to validate the candidate microRNA biomarker and allow better interpretation of its molecular functions and derived cellular interactions. The in situ hybridization data provides information about contextual distribution and cellular origin of the microRNA. By combining microRNA in situ hybridization with immunohistochemical staining of protein markers, it is possible to precisely characterize the microRNA-expressing cells and to identify the potential microRNA targets. This combined technology can also help to monitor changes in the level of potential microRNA targets in a therapeutic setting. In this chapter, we present a fluorescence-based detection method that allows the combination of microRNA in situ hybridization with immunohistochemical staining of one and, in this updated version of the paper, two protein markers detected with primary antibodies raised in the same host species.

SNP rs4937333 in the miRNA-5003-Binding Site of the ETS1 3'-UTR Decreases ETS1 Expression

Mutations in and reduced expression of the ETS1 gene may be associated with systemic lupus erythematosus (SLE). Here, we report a replication study to investigate associations of eight ETS1 single-nucleotide polymorphisms in the 3′-untranslated region (3′-UTR) with SLE and their regulation of ETS1 expression in a study population. We found that the rs4937333 T allele was associated with a significantly increased risk of SLE (odds ratio: 1.800, 95% confidence interval: 1.02–3.157, P = 0.040) and with dramatically reduced levels of ETS1 in B cells from SLE subjects. Functionally, the rs4937333 T allele alters the binding affinity between miR-5003 and its ETS1 3′-UTR target, thus enhancing suppression of ETS1 expression. Furthermore, immunoglobulin M-secreting plasmacytes were significantly reduced among B cells with the rs4937333 C allele versus the T allele according to FACS and ELISA. Additionally, miR-5003 expression was higher in B cells than in T cells from SLE patients, and a negative correlation between miR-5003 and ETS1 was found, especially in B cells with the T allele. These findings suggest that the rs4937333 T allele is a risk factor for susceptibility to SLE in the studied population. The rs4937333 T allele may enhance the binding of miR-5003 to ETS1, which probably promotes the involvement of ETS1 in the differentiation of B cells into plasmacytes.

Evidence of functional ryanodine receptors in rat mesenteric collecting lymphatic vessels

In the current study, the potential contributions of ryanodine receptors (RyRs) to intrinsic pumping and responsiveness to substance P (SP) were investigated in isolated rat mesenteric collecting lymphatic vessels. Responses to SP were characterized in lymphatic vessels in the absence or presence of pretreatment with nifedipine to block L-type Ca²⁺ channels, caffeine to block normal release and uptake of Ca²⁺ from the sarcoplasmic reticulum, ryanodine to block all RyR isoforms, or dantrolene to more selectively block RyR1 and RyR3. RyR expression and localization in lymphatics was also assessed by quantitative PCR and immunofluorescence confocal microscopy. The results show that SP normally elicits a significant increase in contraction frequency and a decrease in end-diastolic diameter. In the presence of nifedipine, phasic contractions stop, yet subsequent SP treatment still elicits a strong tonic contraction. Caffeine treatment gradually relaxes lymphatics, causing a loss of phasic contractions, and prevents subsequent SP-induced tonic contraction. Ryanodine also gradually diminishes phasic contractions but without causing vessel relaxation and significantly inhibits the SP-induced tonic contraction. Dantrolene treatment did not significantly impair lymphatic contractions nor the response to SP. The mRNA for all RyR isoforms is detectable in isolated lymphatics. RyR2 and RyR3 proteins are found predominantly in the collecting lymphatic smooth muscle layer. Collectively, the data suggest that SP-induced tonic contraction requires both extracellular Ca²⁺ plus Ca²⁺ release from internal stores and that RyRs play a role in the normal contractions and responsiveness to SP of rat mesenteric collecting lymphatics.NEW & NOTEWORTHY The mechanisms that govern contractions of lymphatic vessels remain unclear. Tonic contraction of lymphatic vessels caused by substance P was blocked by caffeine, which prevents normal uptake and release of Ca²⁺ from internal stores, but not nifedipine, which blocks L-type channel-mediated Ca²⁺ entry. Ryanodine, which also disrupts normal sarcoplasmic reticulum Ca²⁺ release and reuptake, significantly inhibited substance P-induced tonic contraction. Ryanodine receptors 2 and 3 were detected within the smooth muscle layer of collecting lymphatic vessels.

Genetic variants in the calcium signaling pathway genes are associated with cutaneous melanoma-specific survival

Remodeling or deregulation of the calcium signaling pathway is a relevant hallmark of cancer including cutaneous melanoma (CM). In this study, using data from a published genome-wide association study (GWAS) from The University of Texas M.D. Anderson Cancer Center, we assessed the role of 41,377 common single-nucleotide polymorphisms (SNPs) of 167 calcium signaling pathway genes in CM survival. We used another GWAS from Harvard University as the validation dataset. In the single-locus analysis, 1830 SNPs were found to be significantly associated with CM-specific survival (CMSS; P ≤ 0.050 and false-positive report probability ≤ 0.2), of which 9 SNPs were validated in the Harvard study (P ≤ 0.050). Among these, three independent SNPs (i.e. PDE1A rs6750552 T>C, ITPR1 rs6785564 A>G and RYR3 rs2596191 C>A) had a predictive role in CMSS, with a meta-analysis-derived hazards ratio of 1.52 (95% confidence interval = 1.19-1.94, P = 7.21 × 10-4), 0.49 (0.33-0.73, 3.94 × 10-4) and 0.67 (0.53-0.86, 0.0017), respectively. Patients with an increasing number of protective genotypes had remarkably improved CMSS. Additional expression quantitative trait loci analysis showed that these genotypes were also significantly associated with mRNA expression levels of the genes. Taken together, these results may help us to identify prospective biomarkers in the calcium signaling pathway for CM prognosis.

Association of MIR146A rs2910164 variation with a predisposition to sporadic breast cancer in a Pakistani cohort

Single-nucleotide polymorphisms (SNPs) in genes coding for microRNAs (miRNAs) play a pivotal role in the progression of breast cancer (BC). We investigated the association of miR-146a rs2910164 GC polymorphism with the risk of BC in the Pakistani population. The miR-146a rs2910164 polymorphism was genotyped in 300 BC cases and 300 age- and gender-matched healthy controls using T-ARMS-PCR. Genotype and allele frequencies were calculated and the association between genotypes and the risk of BC was calculated by odds ratio (OR) and confidence interval (95%). A significant difference in genotypic frequencies (χ²  = 63.10; P = <0.0001) and allelic frequencies (OR = 0.3955 (0.3132-0.4993); P = < 0.0001) was observed between cases and controls. Furthermore, we also found that miR-146 rs2910164 CC homozygote increased the risk of BC in the dominant (OR = 0.2397 (0.1629-0.3526); P = 0.0001; GG vs. GC + CC) and recessive (OR = 2.803 (1.865-4.213); P = <0.0001; CC vs. GC + GG) inheritance models. In summary, miR-146a rs2910164 GC is significantly associated with BC in the Pakistani population. To our knowledge, this is the first study that assessed MIR146a rs2910164 G > C SNP in Pakistani population. By analyzing the secondary structure of MIR146A variant, a significant structural modification was noted. Study with a larger sample size is needed to further confirm of these findings.

A DAAM1 3'-UTR SNP mutation regulates breast cancer metastasis through affecting miR-208a-5p-DAAM1-RhoA axis

Background

Dishevelled-associated activator of morphogenesis 1 (DAAM1) is a member of microfilament-related formins and mediates cell motility in breast cancer (BrCa). However, the genetic mutation status of DAAM1 mRNA and its correlation with pathological characteristics are still unclearly. Methods: A patient cohort and BrCa cells were recruited to demonstrate the role of functional SNP in microRNA-208a-5p binding site of DAAM1 3′-UTR and underlying mechanism in BrCa metastasis.

Methods

A patient cohort and BrCa cells were recruited to demonstrate the role of functional SNP in microRNA-208a-5p binding site of DAAM1 3′-UTR and underlying mechanism in BrCa metastasis.

Results

The expression and activation of DAAM1 increased markedly in lymphnode metastatic tissues. A genetic variant (rs79036859 A/G) was validated in the miR-208a-5p binding site of DAAM1 3′-UTR. The G genotype (AG/GG) was a risk genotype for the metastasis of BrCa by reducing binding affinity of miR-208a-5p for the DAAM1 3′-UTR. Furthermore, the miR-208a-5p expression level was significantly suppressed in lymphnode metastatic tissues compared with that in non-lymphnode metastatic tissues. Overexpression of miR-208a-5p inhibited DAAM1/RhoA signaling pathway, thereby leading to the decrease of the migratory ability.

Conclusion

Overall, the rs79036859 G variant of DAAM1 3′-UTR was identified as a relevant role in BrCa metastasis via the diversity of miR-208a-5p binding affinity.

Electronic supplementary material

The online version of this article (10.1186/s12935-019-0747-8) contains supplementary material, which is available to authorized users.

Comparison of breast cancer risk factors among molecular subtypes: A case-only study

Epidemiological studies have a clear definition of the risk factors for breast cancer. However, it is unknown whether the distribution of these factors differs among breast cancer subtypes. We conducted a hospital‐based case‐only study consisting of 8067 breast cancer patients basing on the Tianjin Cohort of Breast Cancer Cases. Major breast cancer subtypes including luminal A, luminal B, human epidermal growth factor receptor 2 (HER2)‐enriched and basal‐like were defined by estrogen receptor, progesterone receptor, HER2, and Ki‐67 status. Variables including demographic characteristics, reproductive factors, lifestyle habits, imaging examination, and clinicopathologic data were collected for patients. Chi‐square test and one‐way analysis of variance were used to compare the distributions of variables among the four breast cancer subtypes. Multivariate logistic regression was used to estimate the odds ratios and associated 95% confidence intervals where luminal A patients served as the reference group. Overall, more commonality rather than heterogeneity on the distributions of factors was found between the four molecular subtypes of breast cancer. The proportion of overweight and obesity were lower in HER2‐enriched subtype. Women with age at menarche ≤13 years were more likely to be found in basal‐like subtype. Postmenopausal women were more frequent in HER2‐enriched and basal‐like subtypes. Women with benign breast disease and higher breast density were more common in HER2‐enriched subtype. Risk factor scoring showed that total risk scores were similar among the four subtypes. HER2‐enriched and basal‐like subtypes were more frequently diagnosed with large tumors. Calcification was more likely to be found in luminal B and HER2‐enriched subtypes, whereas less distributed in basal‐like subtype. Most of the breast cancer risk factors were similarly distributed among the four major breast cancer subtypes; commonality is predominant.

Potential dysfunctional effects of synonymous variants: Insights from an exhaustive in silico analysis of the ABCB4 gene

The multiple drug resistance 3 (MDR3) protein is a canalicular phospholipid translocator involved in the bile secretion and encoded by the ABCB4 gene. Its deficiency is related to a large spectrum of liver diseases. Taking into account the increased evidence about the involvement of synonymous variants in inherited diseases, this study aims to explore the putative effects of silent genetic variants on the ABCB4 expression. We performed an exhaustive computational approach using ESE finder, RegRNA 2.0, MFOLD, SNPfold, and %MinMax software added to the measurement of the Relative Synonymous Codon Usage. This analysis included 216 synonymous variants distributed throughout the ABCB4 gene. Results have shown that 11 synonymous coding SNPs decrease the ESE activity, while 8 of them change the codon frequency. Besides, the c.24C>T variation, located 21 nucleotides downstream the start A (Adenine) U (Uracil) G (Glutamine) AUG causes an increase in the local stability. Moreover, the computational analysis of the 3'UTR region showed that six of the eight variants located in this region affected the Wild Type (WT) pattern of the miRNA targets sites and/or their proper display. The 26 sSNPs retained as putatively functional possessed a very low allele frequency, supporting their pathogenicity. In conclusion, the obtained results suggest that some synonymous SNPs in the ABCB4 gene, considered up to now as neutral, may be involved in the MDR3 deficiency.