Tumor suppression by miR-31 in esophageal carcinoma is p21-dependent

A bioinformatics analysis for diagnostic roles of the E2F family in esophageal cancer

BACKGROUND

Esophageal cancer (EC) is the eighth most commonly occurring cancer worldwide and the sixth leading cause of cancer-related deaths. The therapeutic effect of EC patients is not ideal, and new biomarkers are needed to guide diagnosis and prognosis of EC patients. E2F family transcription factors are among the most important links in the cell cycle regulatory network. E2Fs dysregulation not only promotes the early stages of tumor development but also the progression of benign tumors to malignant tumors. E2F is expected to be a new biomarker. The prognostic significance of the E2F family in EC requires further research.

METHODS

We analyzed The Cancer Genome Atlas (TCGA), Gene Expression Profiling Interactive Analysis (GEPIA), and GeneMANIA databases to obtain RNA-sequencing data and clinical data. The clinical data included age, gender, race, stage, type, status, etc. The prognosis outcome included overall survival (OS) and progression-free interval (PFI). Subsequently, we conducted further research on gene expressions, enrichment analysis, interaction network, and prognostic values by R software, containing ggplot2, ComplexHeatmap, DESeq2, pROC R package, based on t-test, Wilcoxon rank sum test, Spearman rank correlation analysis, log-rank test and COX model.

RESULTS

We found that mRNA transcription levels of E2F1, E2F3-8 were more highly expressed in esophageal carcinoma (ESCA) tissues than in normal tissues. E2F8 expression was correlated with tumor stage [Pr(>F)=0.00856]. E2F-related genes played a role in development and differentiation, and were prevalent in the endoplasmic reticulum lumen, Golgi lumen, and lipoprotein particle, catalyzing translation activities and lipid metabolism. Each gene was found to be related to each other to some degree. The GeneMANIA network analysis revealed links between E2Fs and other genes. We compared the correlations between 24 kinds of tumor-infiltrating immune cells and E2Fs. E2F1 (AUC =0.945, CI: 0.890-1.000) and E2F7 (AUC =0.958, CI: 0.920-0.996) exhibited higher predictive power accuracy. However, only E2F7 was closely related to OS [HR =1.91 (1.16-3.16), P=0.011].

CONCLUSIONS

We discover that E2F7 is a prognostic biomarker. E2F family may take part in the development of EC through lipid metabolism pathways, which is helpful to predict the prognosis of EC patients and guide accurate diagnosis and treatment.

HiFreSP: A novel high-frequency sub-pathway mining approach to identify robust prognostic gene signatures

With the increasing awareness of heterogeneity in cancers, better prediction of cancer prognosis is much needed for more personalized treatment. Recently, extensive efforts have been made to explore the variations in gene expression for better prognosis. However, the prognostic gene signatures predicted by most existing methods have little robustness among different datasets of the same cancer. To improve the robustness of the gene signatures, we propose a novel high-frequency sub-pathways mining approach (HiFreSP), integrating a randomization strategy with gene interaction pathways. We identified a six-gene signature (CCND1, CSF3R, E2F2, JUP, RARA and TCF7) in esophageal squamous cell carcinoma (ESCC) by HiFreSP. This signature displayed a strong ability to predict the clinical outcome of ESCC patients in two independent datasets (log-rank test, P = 0.0045 and 0.0087). To further show the predictive performance of HiFreSP, we applied it to two other cancers: pancreatic adenocarcinoma and breast cancer. The identified signatures show high predictive power in all testing datasets of the two cancers. Furthermore, compared with the two popular prognosis signature predicting methods, the least absolute shrinkage and selection operator penalized Cox proportional hazards model and the random survival forest, HiFreSP showed better predictive accuracy and generalization across all testing datasets of the above three cancers. Lastly, we applied HiFreSP to 8137 patients involving 20 cancer types in the TCGA database and found high-frequency prognosis-associated pathways in many cancers. Taken together, HiFreSP shows higher prognostic capability and greater robustness, and the identified signatures provide clinical guidance for cancer prognosis. HiFreSP is freely available via GitHub: https://github.com/chunquanlipathway/HiFreSP.

m6 A demethylase ALKBH5 promotes proliferation of esophageal squamous cell carcinoma associated with poor prognosis

Esophageal squamous cell carcinoma (ESCC) is one of the most fatal types of malignant tumors worldwide. Epitranscriptome, such as N⁶ -methyladenosine (m⁶ A) of mRNA, is an abundant post-transcriptional mRNA modification and has been recently implicated to play roles in several cancers, whereas the significance of m⁶ A modifications is virtually unknown in ESCC. Analysis of tissue microarray of the tumors in 177 ESCC patients showed that higher expression of m⁶ A demethylase ALKBH5 correlated with poor prognosis and that ALKBH5 was an independent prognostic factor of the survival of patients. There was no correlation between the other demethylase FTO and prognosis. siRNA knockdown of ALKBH5 but not FTO significantly suppressed proliferation and migration of human ESCC cells. ALKBH5 knockdown delayed progression of cell cycle and accumulated the cells to G0/G1 phase. Mechanistically, expression of CDKN1A (p21) was significantly up-regulated in ALKBH5-depleted cells, and m⁶ A modification and stability of CDKN1A mRNA were increased by ALKBH5 knockdown. Furthermore, depletion of ALKBH5 substantially suppressed tumor growth of ESCC cells subcutaneously transplanted in BALB/c nude mice. Collectively, we identify ALKBH5 as the first m⁶ A demethylase that accelerates cell cycle progression and promotes cell proliferation of ESCC cells, which is associated with poor prognosis of ESCC patients.

Role of MicroRNAs in Treatment Response in Prostate Cancer

Prostate cancer (PCa) is the most common non-skin cancer in men worldwide, resulting in significant mortality and morbidity. Depending on the grade and stage of the cancer, patients may be given radiation therapy, hormonal therapy, or chemotherapy. However, more than half of these patients develop resistance to treatment, leading to disease progression and metastases, often with lethal consequences. MicroRNAs (miRNAs) are short, non-coding RNAs, which regulate numerous physiological as well as pathological processes, including cancer. miRNAs mediate their regulatory effect predominately by binding to the 3'-untranslated region (UTR) of their target mRNAs. In this review, we will describe the mechanisms by which miRNAs mediate resistance to radiation and drug therapy (i.e. hormone therapy and chemotherapy) in PCa, including control of apoptosis, cell growth and proliferation, autophagy, epithelial-to-mesenchymal transition (EMT), invasion and metastasis, and cancer stem cells (CSCs). Furthermore, we will discuss the utility of circulating miRNAs isolated from different body fluids of prostate cancer patients as non-invasive biomarkers of cancer detection, disease progression, and therapy response. Finally, we will shortlist the candidate miRNAs, which may have a role in drug and radioresistance, that could potentially be used as predictive biomarkers of treatment response.

LRP5 regulates the expression of STK40, a new potential target in triple-negative breast cancers

Triple-negative breast cancers (TNBCs) account for a large proportion of breast cancer deaths, due to the high rate of recurrence from residual, resistant tumor cells. New treatments are needed, to bypass chemoresistance and improve survival. The WNT pathway, which is activated in TNBCs, has been identified as an attractive pathway for treatment targeting. We analyzed expression of the WNT coreceptors LRP5 and LRP6 in human breast cancer samples. As previously described, LRP6 was overexpressed in TNBCs. However, we also showed, for the first time, that LRP5 was overexpressed in TNBCs too. The knockdown of LRP5 or LRP6 decreased tumorigenesis in vitro and in vivo, identifying both receptors as potential treatment targets in TNBC. The apoptotic effect of LRP5 knockdown was more robust than that of LRP6 depletion. We analyzed and compared the transcriptomes of cells depleted of LRP5 or LRP6, to identify genes specifically deregulated by LRP5 potentially implicated in cell death. We identified serine/threonine kinase 40 (STK40) as one of two genes specifically downregulated soon after LRP5 depletion. STK40 was found to be overexpressed in TNBCs, relative to other breast cancer subtypes, and in various other tumor types. STK40 depletion decreased cell viability and colony formation, and induced the apoptosis of TNBC cells. In addition, STK40 knockdown impaired growth in an anchorage-independent manner in vitro and slowed tumor growth in vivo. These findings identify the largely uncharacterized putative protein kinase STK40 as a novel candidate treatment target for TNBC.

LncRNA ATB promotes proliferation and metastasis in A549 cells by down-regulation of microRNA-494

Lung cancer is a commonly diagnosed disease with poor prognosis. Novel therapeutic targets and deep understanding of the regulatory mechanisms in lung cancer are of great importance. We aimed to figure out the functional roles of lncRNA-activated by transforming growth factor-β (ATB) in A549 cells as well as the underlying molecular mechanisms. ATB was non-physiologically expressed in A549 cells after cell transfection. Then, cell proliferation, expressions of proteins related to proliferation and epithelial-mesenchymal transition (EMT), migration, and invasion were measured by BrdU incorporation assay, Western blot analysis, and Transwell assay, respectively. Afterwards, miR-494 expression in transfected A549 cells was determined by quantitative reverse transcription PCR. Meanwhile, effects of miR-494 overexpression on ATB-overexpressed cells were assessed. Finally, the phosphorylation levels of AKT and key kinases in the Janus-activated kinase (JAK)/signal transducer and activator of transcription-3 (STAT3) pathway were detected by Western blot analysis. ATB overexpression promoted proliferation, migration, and invasion of A549 cells. Meanwhile, EMT of A549 cells was also enhanced. ATB silence showed the opposite influence. Expression of miR-494 was negatively regulated by ATB. Following experiments showed ATB-induced alterations of proliferation, migration, invasion, and EMT were all reversed by miR-494 overexpression. Finally, we proved that ATB increased phosphorylated levels of AKT, JAK1, and STAT3, and those increases were all reversed by miR-494 overexpression. We interestingly figured out that ATB promoted proliferation, migration, invasion, and EMT through down-regulating miR-494 in A549 cells. Moreover, ATB might activate AKT and the JAK/STAT3 pathway via down-regulating miR-494.

Deregulation of microRNA‑31a‑5p is involved in the development of primary hypertension by suppressing apoptosis of pulmonary artery smooth muscle cells via targeting TP53

The present study aimed to identify the association between microRNA (miRNA/miR)-31a-5p and the development of hypertension, and its potential molecular mechanism. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) and western blot analyses were performed to validate the candidate miRNA and genes involved in hypertension, following which an online miRNA database search, luciferase assay, and RT-qPCR and western blot analyses were performed to confirm the interaction between miR-31a-5p and TP53. A MTT assay and flow cytometric analysis were utilized to determine the effect of miR-31a-5p on cell growth and apoptosis. The results revealed that miR-31a-5p and TP53 were the candidate miRNA and gene regulating hypertension, and that TP53 was the virtual target gene of miR-31a-5p with a binding site located in the TP53 3′ untranslated region (3′UTR). It was confirmed by luciferase activity that miR-31a-5p markedly reduced the luciferase activity of the Luc-wild-type-TP53-3′UTR, whereas the mutated putative miR-31a-5p binding located on the TP53-3′UTR was found to eliminate such an inhibitory effect. miR-31a-5p had no effect on specificity protein 1, E2F transcription factor 2 or forkhead box P3 luciferase activity. Smooth muscle cells collected from spontaneously hypertensive rats treated with gold nano-particles containing anti-rno-miR-31a-5p exhibited a lower growth rate and a higher apoptotic rate. The results of the RT-qPCR and western blot analyses showed that miR-31a-5p negatively regulated the expression of TP53, and transfection with the hsa-miR-31a-5p mimic significantly promoted cell growth and inhibited cell apoptosis, whereas transfection with the anti-hsa-miR-31a-5p mimic significantly suppressed cell growth and induced cell apoptosis. Taken together, these findings indicated that miR-31a-5p is involved in hypertension via the accelerated proliferation of arterial smooth muscle cells and inhibition of apoptosis through targeting TP53.

Circulating miR-31 as an effective biomarker for detection and prognosis of human cancer: a meta-analysis

PURPOSE

Circulating miR-31 was found to be associated with cancers detection and prognosis. The present meta-analysis aimed to explore the effect of circulating miR-31 on cancer detection and prognosis.

METHOD

The studies were accessed using multiple databases. RevMan5.3, Meta-DiSc 1.4, and STATA14.0 were used to estimate the pooled effects, heterogeneity among studies, and publication bias.

RESULTS

A total of 14 studies with 1397 cancer patients and 1039 controls were included. For the 12 prognostic tests, the adjusted pooled-AUC was 0.79 (95% CI: 0.73-0.86) as the pooled sensitivity, specificity, positive likelihood ratio (PLR), negative likelihood ratio (NLR), diagnostic odd ratio (DOR) from 10 tests was 0.79 (95% CI: 0.76-0.82), 0.79 (95% CI: 0.76-0.82), 3.81 (95% CI: 2.90-5.01), 0.26 (95% CI: 0.20-0.35), and 16.81 (95% CI: 9.67-29.25), respectively. For the 5 prognosis analyses, the pooled HR (hazard ratio) of overall survival (OS) was 1.55 (95% CI 1.30-1.86) for high versus low circulating miR-31 expression. However, high expression of circulating miR-31 did not significantly increase the risk of poor differentiation (pooled OR=1.39, 95% CI: 0.56-3.47) and LNM (pooled OR=3.46, 95% CI: 0.96-12.42) in lung cancer.

CONCLUSION

Circulating miR-31 is an effective biomarker and could be used as a component of miRs signature for cancer detection and prognosis surveillance.

Function of miR-25 in the invasion and metastasis of esophageal squamous carcinoma cells and bioinformatical analysis of the miR-106b-25 cluster

MicroRNAs (miRNAs/miRs) are a class of small, non-coding RNA molecules that serve a key function in carcinogenesis and tumor progression. Recent evidence indicates that miRNAs may act as powerful regulators of migration and invasion. The present study aimed to investigate the effect of miR-25 on the invasion and metastasis of KYSE-150 and EC109 esophageal squamous cell carcinoma (ESCC) cells, and predict the mechanism of this effect by bioinformatically analyzing the miR-106b-25 cluster. In order to alter the expression of miR-25 in the two cell lines, a miR-25 inhibitor or mimic were transfected into the cells, which were then studied via Transwell migration and invasion assays. Subsequently, the target genes of the miR-106b-25 cluster were predicted using miRanda, PicTar, TargetScan and miRTarbase, and the functions of the target genes were predicted via Gene Ontology term and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. Then, a protein-protein interaction (PPI) network was produced using the Search Tool for the Retrieval of Interacting Genes. The results revealed that overexpressing miR-25 led to significantly increased cell migration and invasion in KYSE150 and EC109 cells. Suppressing miR-25 resulted in significantly decreased cell migration and invasion in KYSE150 cells, while the result was not significant in EC109 cells. Target genes of the miR-106b-25 cluster were significantly enriched in the biological process regulation of cellular metabolic process and several cancer-associated pathways, such as those for glioma and melanoma. The PPI network revealed that PTEN, TP53, MDM2, E2F1, PRMT5, MCM2, RB1, CDKN1A, SHAD7 and EZH2 may serve core roles within the network and associate with one another during the pathogenesis of ESCC. These results indicate that a high expression of miR-25 promotes the invasion and metastasis of ESCC cells, while the influence of low expression of miR-25 differs with cells with different degrees of differentiation. Invasion and metastasis are not effected in cells with poor differentiation, while they were decreased in well differentiated cells. Furthermore, PTEN, TP53, MDM2, E2F1, PRMT5, MCM2, RB1, CDKN1A, SHAD7 and EZH2 may be targeted by the miR-106b-25 cluster, and act together to regulate the development of ESCC.

Molecular interplay of pro-inflammatory transcription factors and non-coding RNAs in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma is the sixth most common cancer in the developing world. The aggressive nature of esophageal squamous cell carcinoma, its tendency for relapse, and the poor survival prospects of patients diagnosed at advanced stages, represent a pressing need for the development of new therapies for this disease. Chronic inflammation is known to have a causal link to cancer pre-disposition. Nuclear factor kappa B and signal transducer and activator of transcription 3 are transcription factors which regulate immunity and inflammation and are emerging as key regulators of tumor initiation, progression, and metastasis. Although these pro-inflammatory factors in esophageal squamous cell carcinoma have been well-characterized with reference to protein-coding targets, their functional interactions with non-coding RNAs have only recently been gaining attention. Non-coding RNAs, especially microRNAs and long non-coding RNAs demonstrate potential as biomarkers and alternative therapeutic targets. In this review, we summarize the recent literature and concepts on non-coding RNAs that are regulated by/regulate nuclear factor kappa B and signal transducer and activator of transcription 3 in esophageal cancer progression. We also discuss how these recent discoveries can pave way for future therapeutic options to treat esophageal squamous cell carcinoma.

Senescence-Associated MicroRNAs

Senescent cells arise as a consequence of cellular damage and can have either a detrimental or advantageous impact on tissues and organs depending on the specific cell type and metabolic state. As senescent cells accumulate in tissues with advancing age, they have been implicated in many age-related declines and diseases. The major facets of senescence include two pathways responsible for establishing and maintaining a senescence program, p53/CDKN1A(p21) and CDKN2A(p16)/RB, as well as the senescence-associated secretory phenotype. Numerous MicroRNAs influence senescence by modulating the abundance of key senescence regulatory proteins, generally by lowering the stability and/or translation of mRNAs that encode such factors. Accordingly, understanding the molecular mechanisms by which MicroRNAs influence senescence will enable diagnostic and therapeutic opportunities directed at senescent cells. Here, we review senescence-associated (SA)-MicroRNAs and discuss their implications in senescence-relevant pathologies.

Low p21 level is necessary for the suppressive effects of micoRNA-31 on glioma cell migration and invasion

MicroRNAs (miRNAs), a kind of endogenous non-coding RNAs, regulate gene expression through binding to the 3'-untranslational region (UTR) of target messenger RNAs (mRNAs) and act as endogenous agents of RNA interference, resulting in either mRNA degradation or translational repression. MiR-31 has been demonstrated to be associated with the development and progression of glioma. However, the underlying molecular mechanism remains largely unclear. In the present study, we demonstrated that miR-31 only inhibited the cell migration and invasion, as well as the expression of a known miR-31 target oncogene radixin, in U251 glioma cells that expressed low level of p21; however, miR-31 showed no above effects on glioma SHG44 cells that highly expressed p21. Moreover, upregulation of p21 in U251 cells reversed the suppressive effects of miR-31 on the cell migration and invasion, suggesting that low p21 level is necessary for the miR-31-mediated inhibitory effects on glioma. Furthermore, analysis for 35 glioma specimens showed that the expression of radixin was negatively correlated with the miR-31 level in glioma tissues with low p21 expression; however, no such correlation was found in glioma tissues with high p21 level, further supporting that the low p21 level is necessary for the suppressive effect of miR-31 on the expression of its target oncogenes. In summary, our study demonstrates that the suppressive effect of miR-31 on glioma cell migration and invasion is p21-dependent, and suggests that miR-31 may be used for the treatment of patients with p21-deficent glioma.

Identification of differentially expressed microRNAs between Bacillus thuringiensis Cry1Ab-resistant and -susceptible strains of Ostrinia furnacalis

The Asian corn borer (ACB), Ostrinia furnacalis (Guenée), can develop strong resistance to Cry1Ab, the most widely commercialized Cry toxin for Bt maize worldwide. It is essential to understand the mechanism of resistance for management of this species, but information on the post-transcriptional regulation of Bt resistance in this target insect is limited. In the present study, RNA was extracted from the ACB in various larval stages (1–5 instar) from Cry1Ab-sensitive (ACB-BtS) and -resistant (ACB-AbR) strains, each of which included two biological replicates. Using Illumina sequencing, a total of 23,809,890 high-quality reads were collected from the four ACB libraries. The numbers of known microRNAs (miRNAs) were 302 and 395 for ACB-BtS and 268 and 287 for ACB-AbR. Using Mireap software, we identified 32 and 16 potential novel miRNAs for ACB-BtS and 18 and 22 for ACB-AbR. Among them, 21 known and 1 novel miRNAs had significantly different expression between ACB-BtS and ACB-AbR. Several miRNAs were observed to target potential Bt receptor genes, such as aminopeptidase N and cadherin-like protein. The glycosylphosphatidylinositol-anchor biosynthetic process and ABC transporters pathway were identified through Gene Ontology and KEGG pathway analysis of target genes of the differentially expressed miRNAs.

Prognostic value of cancer stem cell marker CD133 expression in esophageal carcinoma: A meta-analysis

CD133 has been identified as a putative neoplastic stem cell marker in esophageal carcinoma. However, the prognostic value of CD133 overexpression in patients with esophageal carcinoma remains controversial. A meta-analysis of previous studies was performed, in order to assess the association of CD133 overexpression with the clinicopathological characteristics of esophageal carcinoma patients. A total of 7 studies, including 538 patients, were subjected to the final analysis. Our results indicated that a positive CD133 expression was significantly associated with lymph node metastasis [odds ratio (OR)=3.09, 95% confidence interval (CI): 1.93-4.95; P<0.00001], clinical stage (OR=4.26, 95% CI: 1.55-11.73; P=0.005) and histopathological grade (OR=2.40, 95% CI: 1.16-4.94; P=0.02). There was no statistically significant association of CD133 with depth of invasion (OR=1.89, 95% CI: 0.42-8.43; P=0.41). Based on the results of this study, we concluded that CD133 is an efficient prognostic factor in esophageal carcinoma. Higher CD133 expression is significantly associated with lymph node metastasis, clinical stage and histopathological grade.

Expression of TP53, BCL-2, and VEGFA Genes in Esophagus Carcinoma and its Biological Significance

BACKGROUND

The pathogenesis of esophagus carcinoma involves a cascade process consisting of multiple factors and accumulation of gene mutations. It is known that vascular endothelial growth factor (VEGF) mainly regulates de novo vascular formation while B-cell lymphoma-2 (BCL-2) gene exerts a tumor-suppressing effect. The prominent expression of VEGFA and BCL-2 genes, along with the most famous tumor-suppressor gene, TP53, raise the possibly of gene interaction. This study therefore investigated the effect and correlation of TP53, BCL-2, and VEGFA genes on cell proliferation and apoptosis of esophagus carcinoma.

MATERIAL AND METHODS

A total of 30 male rats were prepared by subcutaneous injection of methyl-benzyl-nitrosamine (MBNA) to induce esophagus cancer, along with 30 controlled rats which received saline instead. After 4, 10, 20, or 30 weeks, rats were sacrificed to observe the morphological changes of esophageal mucosa. Cell apoptosis was quantified by terminal deoxynucleotidyl transferase-mediated dUTP nick-end-labeling (TUNEL) assay. Immunohistochemical (IHC) staining was employed to examine the expression of TP53, BCL-2 and VEGFA genes.

RESULTS

With the progression of cancer, pathological damages of esophageal tissue aggravated while the cancer cell apoptosis gradually decreased compared to controlled animals. Protein levels of p53, Bcl-2, and VEGF in the model group were significantly elevated at each time point. Positive correlations existed between p53 and Bcl-2 or VEGF.

CONCLUSIONS

Abnormally elevated expression of TP53, BCL-2, and VEGFA genes may participate in the proliferation of esophagus cancer cells in a synergistic manner.

Differential Expression of miR-34c and Its Predicted Target Genes in Testicular Tissue at Different Development Stages of Swine

To verified the target genes of miR-34c, bioinformatics software was used to predict the targets of miR-34c. Three possible target genes of miR-34c related to spermatogenesis and male reproductive development: zinc finger protein 148 (ZNF148), kruppel-like factor 4 (KLF4), and platelet-derived growth factor receptor alpha (PDGFRA) were predicted. Then, the expression of miR-34c and its target genes were detected in swine testicular tissue at different developmental stages by quantitative polymerase chain reaction. The results suggested that the expression of PDGFRA has the highest negative correlation with miR-34c. Then immunohistochemical staining was done to observe the morphology of swine testicular tissue at 2-days and 3, 4, 5-months of age, which indicated that PDGFRA was mainly expressed in the support cells near the basement membrane during the early development stages of testicular tissue, but that the expression of PDGFRA was gradually reduced in later stages. Therefore, western blot analyzed that the highest expression of PDGFRA was generated in 2-days old testicular tissues and the expression levels reduced at 3 and 4-months old, which correlated with the results of immunohistochemical staining. In conclusion, PDGFRA is a target gene of miR-34c.

Differential Expression of MicroRNAs in Papillary Thyroid Carcinoma and Their Role in Racial Disparity

Objective

MicroRNAs (miRNAs) are known to play important roles in the diagnosis and prognosis of papillary thyroid cancer (PTC), and they are useful in developing targeted therapies. However, there have been no studies on the existence of racial differences in miRNAs expression that could explain differential overall survival of PTC patients. Expression analysis of miRNAs in major racial groups would be important for optimizing personalized treatment strategies. In the current study, we assessed the differential expression of 8 miRNAs between normal and tumor tissues, and also assessed racial differences between African American (AA) and Caucasian American (CA).

Methods

First, the miRNA expression profiling was performed using formalin-fixed paraffin embedded (FFPE) tissue sections of tumor containing over 70% tumor cells. Normal and tumor sections of thyroid tissues were studied from AA and CA patients. The miRNA microarray profiling was done using miRBase version 18 (LC Sciences, Houston, TX, USA). Quantitative real-time PCR (qRT-PCR) was used to validate expression of 8 selected miRNAs.

Results

Ingenuity pathway analysis showed involvement of target genes, such as Ras and NF-κB. Deregulated miRNAs such as miR-221 and miR-31 were found to be statistically significant between the two races. Using qRT-PCR, we found that miR-21, miR-146b, miR-221, miR-222, miR-31, and miR-3613 were up-regulated while miR-138 and miR-98 were down-regulated in tumors compared to normal tissues.

Conclusion

Though sample size was small, we found several deregulated miRNAs having racial differences. The differential expression of miRNAs suggest that these miRNAs and their target genes could be useful to gain further mechanistic insight of PTC and their clinical implications, including miRNA replacement therapy or their knockdown strategies.