Integrative analysis in oral squamous cell carcinoma reveals DNA copy number-associated miRNAs dysregulating target genes

Identification of copy number alternation profiles in metastatic oral squamous carcinoma patients by using microarray-based comparative genomic hybridization: A study on Turkish patients

OBJECTIVE

Oral squamous cell carcinoma (OSCC) is a severe form of cancer affecting different anatomic sites of the oral cavity. OSCC ranks as the sixth most common cancer type with an increasing prevalence globally. However, the mechanisms of OSCC process at later stages are not well understood. In this study, we aimed to determine genetic alternations in metastatic OSCC patients to identify genomic changes occurred at metastatic phase of the disease.

MATERIAL AND METHODS

The Illumina CytoSNP-12 Array was used to determine copy number variations in OSCC cancer genome. Hybridization procedures were performed according to the manufacturer procedures (Illumina). Arrays were scanned on iScan System (Illumina). Data were analyzed using Illumina Genotyping module of Genome Studio software (version 1.2, Illumina). Multiple CNV algorithms and copy number alternations were accessed by Genome Studio. CNVs in whole genome were investigated by using a chromosomal heat map.

RESULTS

We reported that gains in 8q21.11-ter, 9p21.3, 13q14.11-ter, 13q13.3-ter and losses in 5q14.3-ter, 5q35 and 17p13.3-12 were associated with the development of OSCC. In addition, we also detected that deletion in 2q33.2-ter and 2q35-37.3 regions were also associated with OSCC metastasis process.

CONCLUSIONS

Our results were also showed that gains in 11q13.3-q13.4 and 2q13.2 chromosomal regions could promote the metastatic OSCC process. We believe that results of the study will help to find new biomarkers for diagnosis at later stage of OSCC.

miR-29a-3p promotes migration and invasion in ameloblastoma via Wnt/β-catenin signaling by targeting catenin beta interacting protein 1

BACKGROUND

Ameloblastoma (AB) is a common epithelial odontogenic tumor. The Wnt/β-catenin pathway has been found to be related to AB invasion.

METHODS

The alteration expression of microRNAs (miRNAs) and messenger RNAs (mRNAs) was performed by miRNA and mRNA microarray analysis and validated by quantitative real-time polymerase chain reaction (RT-PCR). The effects of miR-29a-3p on migration and invasion in AB cells were evaluated by a transwell assay. Bioinformatic prediction was conducted using the miRSystem and validated by quantitative RT-PCR, western blot, and a luciferase reporter assay.

RESULTS

miR-29a-3p was overexpressed in AB tissues, which promoted the migration and invasion of AB cells in vitro. Catenin beta interacting protein 1 (CTNNBIP1), a negative regulator of the Wnt/β-catenin pathway, was predicted to be a target of miR-29a-3p. miR-29a-3p inhibited the expression of CTNNBIP1 and promoted the expression of the downstream molecules of the Wnt/β-catenin pathway.

CONCLUSIONS

miR-29a-3p promoted migration and invasion in AB via Wnt/β-catenin signaling by targeting CTNNBIP1.

Monocyte Chemoattractant Protein 1 Promotes VEGF-A Expression in OSCC by Activating ILK and MEK1/2 Signaling and Downregulating miR-29c

Oral squamous cell carcinoma (OSCC) is an aggressive tumor that has a poor prognosis, with high levels of local invasion and lymph node metastasis. Vascular endothelial growth factor A (VEGF-A) plays essential roles in OSCC tumor angiogenesis and metastasis. Monocyte chemoattractant protein-1 (MCP-1, CCL2) is implicated in various inflammatory conditions and pathological processes, including oral cancer. The existing evidence has failed to confirm any correlation between MCP-1 or VEGF-A expression and OSCC angiogenesis. In this study, high expression levels of MCP-1 and VEGF-A were positively correlated with disease stage in patients with OSCC. In oral cancer cells, MCP-1 increased VEGF-A expression and subsequently promoted angiogenesis; miR-29c mimic reversed MCP-1 activity. We also found that MCP-1 modulated VEGF-A expression and angiogenesis through CCR2/ILK/MEK1/2 signaling. Ex vivo results of the chick embryo chorioallantoic membrane (CAM) assay revealed the angiogenic qualities of MCP-1, with increased numbers of visible blood vessel branches. Our data suggest that MCP-1 is a new molecular therapeutic target for the inhibition of angiogenesis and metastasis in OSCC.

Good or not good: Role of miR-18a in cancer biology

miR-18a is a member of primary transcript called miR-17-92a (C13orf25 or MIR17HG) which also contains five other miRNAs: miR-17, miR-19a, miR-20a, miR-19b and miR-92a. This cluster as a whole shows specific characteristics, where miR-18a seems to be unique. In contrast to the other members, the expression of miR-18a is additionally controlled and probably functions as its own internal controller of the cluster. miR-18a regulates many genes involved in proliferation, cell cycle, apoptosis, response to different kinds of stress, autophagy and differentiation. The disturbances of miR-18a expression are observed in cancer as well as in different diseases or pathological states. The miR-17-92a cluster is commonly described as oncogenic and it is known as 'oncomiR-1', but this statement is a simplification because miR-18a can act both as an oncogene and a suppressor. In this review we summarize the current knowledge about miR-18a focusing on its regulation, role in cancer biology and utility as a potential biomarker.

[MiR-140-3p Downregulation in Association with PDL-1 Overexpression in Many Cancers: A Review from the Literature Using Predictive Bioinformatics Tools]

Programmed death-ligand 1 (PD-L1) has been speculated to play a critical role in suppression of the immune system and it can be upregulated in cancer cells, which may allow cancers to evade the host immune system. MicroRNAs (miRNAs) are small non-coding RNA molecules (containing about 22 nucleotides), that function in RNA silencing and post-transcriptional regulation of gene expression. MiRNAs were found deregulated (upregulated or downregulated) and implicated in cancer development with various roles which depend on their gene target. Using targetscan web server prediction algorithm, we concluded that miR-140-3p is a targeting mirRNA with conserved consequential pairing of target region for PD-L1. Moreover, by reviewing all the available cancer studies in Pub/Medline about miR-140-3p, was found permanently down regulated. Furthermore, in recent immunotherapy related clinical trials in most cancers, evaluated PD-L1, it is found overexpressed. In the near future, in vitro or in vivo studies need to validate whether there is direct correlation between PD-L1 overexpression and miR-140-3p downregulation as targetscan performed algorithm predicted.

Differential expression of hsa-miR-221, hsa-miR-21, hsa-miR-135b, and hsa-miR-29c suggests a field effect in oral cancer

BACKGROUND

The theory of field effect suggests that the tumor-adjacent area, besides histopathologically normal, undergoes genetic and epigenetic changes that can eventually affect epithelial homeostasis, predisposing the patient to cancer development. One of the many molecular changes described in cancer are microRNAs (miRNAs), which regulates the expression of important genes during carcinogenesis. Thus, the aim of this study was to investigate the field effect in oral cancer.

METHODS

We investigated the differential expression profile of four miRNAs (hsa-miR-221, hsa-miR-21, hsa-miR-135b, and hsa-miR-29c) in cancerous oral tissue, in tumor-adjacent tissue and and in non-cancerous tissue samples from healthy volunteers.

RESULTS

Our results showed significant overexpression profiles of all four studied miRNAs in cancerous oral tissue compared to non-cancerous samples, as well as in tumor-adjacent tissue compared to cancer-free tissue. No significant difference was found when comparing the expression profile of cancerous and tissue-adjacent tissue groups. We found a negative correlation between the expression of hsa-miR-21 expression and STAT3 in oral squamous cell carcinoma.

CONCLUSION

These results suggest that the tissue adjacent to cancer cannot be considered a normal tissue because its molecular aspects are significantly altered. Our data corroborates the hypothesis of field cancerization.

Cigarette smoke and chewing tobacco alter expression of different sets of miRNAs in oral keratinocytes

Carcinogenic effect of tobacco in oral cancer is through chewing and/or smoking. Significant differences exist in development of oral cancer between tobacco users and non-users. However, molecular alterations induced by different forms of tobacco are yet to be fully elucidated. We developed cellular models of chronic exposure to chewing tobacco and cigarette smoke using immortalized oral keratinocytes. Chronic exposure to tobacco resulted in increased cell scattering and invasiveness in immortalized oral keratinocytes. miRNA sequencing using Illumina HiSeq 2500 resulted in the identification of 10 significantly dysregulated miRNAs (4 fold; p ≤ 0.05) in chewing tobacco treated cells and 6 in cigarette smoke exposed cells. We integrated this data with global proteomic data and identified 36 protein targets that showed inverse expression pattern in chewing tobacco treated cells and 16 protein targets that showed inverse expression in smoke exposed cells. In addition, we identified 6 novel miRNAs in chewing tobacco treated cells and 18 novel miRNAs in smoke exposed cells. Integrative analysis of dysregulated miRNAs and their targets indicates that signaling mechanisms leading to oncogenic transformation are distinct between both forms of tobacco. Our study demonstrates alterations in miRNA expression in oral cells in response to two frequently used forms of tobacco.

In-Silico Integration Approach to Identify a Key miRNA Regulating a Gene Network in Aggressive Prostate Cancer

Like other cancer diseases, prostate cancer (PC) is caused by the accumulation of genetic alterations in the cells that drives malignant growth. These alterations are revealed by gene profiling and copy number alteration (CNA) analysis. Moreover, recent evidence suggests that also microRNAs have an important role in PC development. Despite efforts to profile PC, the alterations (gene, CNA, and miRNA) and biological processes that correlate with disease development and progression remain partially elusive. Many gene signatures proposed as diagnostic or prognostic tools in cancer poorly overlap. The identification of co-expressed genes, that are functionally related, can identify a core network of genes associated with PC with a better reproducibility. By combining different approaches, including the integration of mRNA expression profiles, CNAs, and miRNA expression levels, we identified a gene signature of four genes overlapping with other published gene signatures and able to distinguish, in silico, high Gleason-scored PC from normal human tissue, which was further enriched to 19 genes by gene co-expression analysis. From the analysis of miRNAs possibly regulating this network, we found that hsa-miR-153 was highly connected to the genes in the network. Our results identify a four-gene signature with diagnostic and prognostic value in PC and suggest an interesting gene network that could play a key regulatory role in PC development and progression. Furthermore, hsa-miR-153, controlling this network, could be a potential biomarker for theranostics in high Gleason-scored PC.

Expression and Regulatory Network Analysis of miR-140-3p, a New Potential Serum Biomarker for Autism Spectrum Disorder

Given its prevalence and social impact, Autism Spectrum Disorder (ASD) is drawing much interest. Molecular basis of ASD is heterogeneous and only partially known. Many factors, including disorders comorbid with ASD, like TS (Tourette Syndrome), complicate ASD behavior-based diagnosis and make it vulnerable to bias. To further investigate ASD etiology and to identify potential biomarkers to support its precise diagnosis, we used TaqMan Low Density Array technology to profile serum miRNAs from ASD, TS, and TS+ASD patients, and unaffected controls (NCs). Through validation assays in 30 ASD, 24 TS, and 25 TS+ASD patients and 25 NCs, we demonstrated that miR-140-3p is upregulated in ASD vs.: NC, TS, and TS+ASD (Tukey's test, p-values = 0.03, = 0.01, < 0.0001, respectively). ΔCt values for miR-140-3p and YGTSS (Yale Global Tic Severity Scale) scores are positively correlated (Spearman r = 0.33; Benjamini-Hochberg p = 0.008) and show a linear relationship (p = 0.002). Network functional analysis showed that nodes controlled by miR-140-3p, especially CD38 and NRIP1 which are its validated targets, are involved in processes convergingly dysregulated in ASD, such as synaptic plasticity, immune response, and chromatin binding. Biomarker analysis proved that serum miR-140-3p can discriminate among: (1) ASD and NC (Area under the ROC curve, AUC: 0.70; sensitivity: 63.33%; specificity: 68%); (2) ASD and TS (AUC: 0.72; sensitivity: 66.66%; specificity: 70.83%); (3) ASD and TS+ASD (AUC: 0.78; sensitivity: 73.33%; specificity: 76%). Characterization of miR-140-3p network would contribute to further clarify ASD etiology. Serum miR-140-3p could represent a potential non-invasive biomarker for ASD, easy to test through liquid biopsy.

Correlation between the level of microRNA expression in peripheral blood mononuclear cells and symptomatology in patients with generalized anxiety disorder

This study investigated the correlation between the level of microRNA expression in peripheral blood mononuclear cells (PBMCs) and symptomatology in patients with generalized anxiety disorder (GAD). MicroRNA array was performed in peripheral blood mononuclear cells (PBMCs) obtained from GAD patients with gender, age, ethnicity-matched healthy controls. Then real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to verify the top 7 miRNAs with the highest fold-change values in 76 GAD patients and 39 healthy controls. It demonstrated that 5 miRNAs showed significantly differences in expression levels (P<0.01). These 5 GAD-associated miRNAs were finally selected into our study to analyze the association between the plasma level of miRNAs expression and symptomatology scores in Hamilton Anxiety Scale (HAMA). Results showed that the level of miR-4505 and miR-663 was negatively correlated with the total HAMA scores in GAD patients (r=0.2228, r=0.264 P<0.05). MiR-663 was selected into the regression equation of HAMA total scores and psychic anxiety symptomatology scores, and it could explain 5.3% of the HAMA total scores and 15.3% of the anxiety symptomatology scores. This study analyzed preliminarily possible circulating miRNAs expression changes in GAD patients, and the expression level of miR-663 highly correlated with psychic anxiety symptoms, further molecular mechanism of which needs to be explored.

Role of miR-140 in embryonic bone development and cancer

Bone is increasingly viewed as an endocrine organ with key biological functions. The skeleton produces hormones and cytokines, such as FGF23 and osteocalcin, which regulate an extensive list of homoeostatic functions. Some of these functions include glucose metabolism, male fertility, blood cell production and calcium/phosphate metabolism. Many of the genes regulating these functions are specific to bone cells. Some of these genes can be wrongly expressed by other malfunctioning cells, driving the generation of disease. The miRNAs are a class of non-coding RNA molecules that are powerful regulators of gene expression by suppressing and fine-tuning target mRNAs. Expression of one such miRNA, miR-140, is ubiquitous in chondrocyte cells during embryonic bone development. Activity in cells found in the adult breast, colon and lung tissue can silence genes required for tumour suppression. The realization that the same miRNA can be both normal and detrimental, depending on the cell, tissue and time point, provides a captivating twist to the study of whole-organism functional genomics. With the recent interest in miRNAs in bone biology and RNA-based therapeutics on the horizon, we present a review on the role of miR-140 in the molecular events that govern bone formation in the embryo. Cellular pathways involving miR-140 may be reactivated or inhibited when treating skeletal injury or disorder in adulthood. These pathways may also provide a novel model system when studying cancer biology of other cells and tissues.

The TF-miRNA Coregulation Network in Oral Lichen Planus

Oral lichen planus (OLP) is a chronic inflammatory disease that affects oral mucosa, some of which may finally develop into oral squamous cell carcinoma. Therefore, pinpointing the molecular mechanisms underlying the pathogenesis of OLP is important to develop efficient treatments for OLP. Recently, the accumulation of the large amount of omics data, especially transcriptome data, provides opportunities to investigate OLPs from a systematic perspective. In this paper, assuming that the OLP associated genes have functional relationships, we present a new approach to identify OLP related gene modules from gene regulatory networks. In particular, we find that the gene modules regulated by both transcription factors (TFs) and microRNAs (miRNAs) play important roles in the pathogenesis of OLP and many genes in the modules have been reported to be related to OLP in the literature.

Copy number variable microRNAs in schizophrenia and their neurodevelopmental gene targets

BACKGROUND

MicroRNAs (miRNAs) are key regulators of gene expression in the human genome and may contribute to risk for neuropsychiatric disorders. miRNAs play an acknowledged role in the strongest of genetic risk factors for schizophrenia, 22q11.2 deletions. We hypothesized that in schizophrenia there would be an enrichment of other rare copy number variants (CNVs) that overlap miRNAs.

METHODS

Using high-resolution genome-wide microarrays and rigorous methods, we compared the miRNA content of rare CNVs in well-characterized cohorts of schizophrenia cases (n = 420) and comparison subjects, excluding 22q11.2 CNVs. We also performed a gene-set enrichment analysis of the predicted miRNA target genes.

RESULTS

The schizophrenia group was enriched for the proportion of individuals with a rare CNV overlapping a miRNA (3.29-fold increase over comparison subjects, p < .0001). The presence of a rare CNV overlapping a miRNA remained a significant predictor of schizophrenia case status (p = .0072) in a multivariate logistic regression model correcting for total CNV size. In contrast, comparable analyses correcting for CNV size showed no enrichment of rare CNVs overlapping protein-coding genes. A gene-set enrichment analysis indicated that predicted target genes of recurrent CNV-overlapped miRNAs in schizophrenia may be functionally enriched for neurodevelopmental processes, including axonogenesis and neuron projection development. Predicted gene targets driving these results included CAPRIN1, NEDD4, NTRK2, PAK2, RHOA, and SYNGAP1.

CONCLUSIONS

These data are the first to demonstrate a genome-wide role for CNVs overlapping miRNAs in the genetic risk for schizophrenia. The results provide support for an expanded multihit model of causation, with potential implications for miRNA-based therapeutics.

MicroRNA gene polymorphisms and environmental factors increase patient susceptibility to hepatocellular carcinoma

BACKGROUND

Micro RNAs (miRNAs) are small RNA fragments that naturally exist in the human body. Through various physiological mechanisms, miRNAs can generate different functions for regulating RNA protein levels and balancing abnormalities. Abnormal miRNA expression has been reported to be highly related to several diseases and cancers. Single-nucleotide polymorphisms (SNPs) in miRNAs have been reported to increase patient susceptibility and affect patient prognosis and survival. We adopted a case-control research design to verify the relationship between miRNAs and hepatocellular carcinoma.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 525 subjects, including 377 controls and 188 hepatocellular carcinoma patients, were selected. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and real-time PCR were used to analyze miRNA146a (rs2910164), miRNA149 (rs2292832), miRNA196 (rs11614913), and miRNA499 (rs3746444) genetic polymorphisms between the control group and the case group. The results indicate that people who carry the rs3746444 CT or CC genotypes may have a significantly increased susceptibility to hepatocellular carcinoma (adjusted odds ratio [AOR] = 2.84, 95% confidence interval [CI] = 1.88-4.30). In addition, when combined with environmental risk factors, such as smoking and alcohol consumption, interaction effects were observed between gene polymorphisms and environmental factors (odds ratio [OR] = 4.69, 95% CI = 2.52-8.70; AOR = 3.38, 95% CI = 1.68-6.80).

CONCLUSIONS

These results suggest that a significant association exists between miRNA499 SNPs and hepatocellular carcinoma. Gene-environment interactions of miRNA499 polymorphisms, smoking, and alcohol consumption might alter hepatocellular carcinoma susceptibility.

Individual and combined effects of DNA methylation and copy number alterations on miRNA expression in breast tumors

Background

The global effect of copy number and epigenetic alterations on miRNA expression in cancer is poorly understood. In the present study, we integrate genome-wide DNA methylation, copy number and miRNA expression and identify genetic mechanisms underlying miRNA dysregulation in breast cancer.

Results

We identify 70 miRNAs whose expression was associated with alterations in copy number or methylation, or both. Among these, five miRNA families are represented. Interestingly, the members of these families are encoded on different chromosomes and are complementarily altered by gain or hypomethylation across the patients. In an independent breast cancer cohort of 123 patients, 41 of the 70 miRNAs were confirmed with respect to aberration pattern and association to expression. In vitro functional experiments were performed in breast cancer cell lines with miRNA mimics to evaluate the phenotype of the replicated miRNAs. let-7e-3p, which in tumors is found associated with hypermethylation, is shown to induce apoptosis and reduce cell viability, and low let-7e-3p expression is associated with poorer prognosis. The overexpression of three other miRNAs associated with copy number gain, miR-21-3p, miR-148b-3p and miR-151a-5p, increases proliferation of breast cancer cell lines. In addition, miR-151a-5p enhances the levels of phosphorylated AKT protein.

Conclusions

Our data provide novel evidence of the mechanisms behind miRNA dysregulation in breast cancer. The study contributes to the understanding of how methylation and copy number alterations influence miRNA expression, emphasizing miRNA functionality through redundant encoding, and suggests novel miRNAs important in breast cancer.

MicroRNA aberrances in head and neck cancer: pathogenetic and clinical significance

PURPOSE OF REVIEW

MicroRNAs (miRNAs) play crucial roles in modulating the neoplastic process of cancers including head and neck squamous cell carcinoma (HNSCC). miRNAs modulate pathogenesis by inhibiting target genes. Understanding how aberrant miRNAs are involved in HNSCC pathogenesis should help to validate potential clinical applications that target these entities.

RECENT FINDINGS

miR-21, miR-31, miR-504 and miR-10b are important oncogenic miRNAs that are involved in HNSCC and target tumour suppressor genes. The tumour suppressor roles of the let-7 family, the miR-99 family, miR-107, miR-133a, miR-137, miR-138 and miR-375 with respect to their targeting of oncogenes are unequivocal and have been confirmed by many studies. In addition, miR-21, let-7, miR-107, miR-138 and miR-200c seem to play complicated roles in regulating stemness or the epithelial-mesenchymal transition of tumour cells. The clinical implications of these tumour-associated miRNAs are generally in agreement with their functional roles.

SUMMARY

A number of pathways that become disregulated by aberrant miRNAs have been identified specifically for HNSCC. Analysis of these networks and their therapeutic interception might facilitate the prediction of disease status and help with the design of therapeutic trials.